(蟒蛇)
导入sympy
从sympy import isprime
def poly4(x):
..如果为isprime(x):
....f=x**4-x**3-x**2-x-1
....如果是主(f**4-f**3-f**2-f-1):
......返回True
..返回False
x=1
当x<10**5时:
..如果是poly4(x):
....打印(x**4-x**3-x**2-x-1)
..x+=1
(平价)s=[];对于素数(p=2700,n=p^4-p^3-p^2-p-1;如果(isprime(n^4-n^3-n^2-n-1),s=concat(s,n));s\\科林·巴克2014年2月11日
|