%I#29 2024年3月22日17:21:33
%S 2,5,31,43,53,79,83,89,97109131139151199229233239283313317,
%电话:359367389433443479487569571577601617641643659677769,
%电话:79782385385792993794197110191049106910691087109311117116311711181123112491283
%N素数p,使得多项式x^4-x^3-x^2-x-1是不可约模p。
%H Robert Israel,n的表,a(n)表示n=1..10000</a>
%p p:=x^4-x^3-x^2-x-1:
%p选择(p->Irreduc(p)mod p,[seq(ithprime(i),i=1..1000)]);
%o(Python)
%o从itertools导入岛
%o从sympy导入Poly,nextprime
%o从sympy.abc导入x
%o def A370830_gen():#术语生成器
%o p=2
%o当为True时:
%o如果Poly(x*(x*,x*(x-1)-1)-1,x,模量=p)。is_irreducible:
%o产量p
%o p=下一素数(p)
%o A370830_list=list(岛屿(A370830-gen(),20))#_Chai Wah Wu_,2024年3月14日
%Y A106283的后续序列。参见106309。
%K nonn公司
%O 1,1
%以色列罗贝尔特,2024年3月13日
|