|
（PARI）polisprimitive（poli）=np=2^poldegree（poly）-1；if（类型（（x^np-1）/poli）！=“t_POL”，返回（0））；对于步骤（k=np-1，1，-1，如果（类型（（x^k-1）/poli）==“t_POL”，返回（0））；）；返回（1）；
列表（nn）={对于（n=0，nn，poli=Mod（1，2）*（x^（2*n+2）+x+1）；如果（polisirreducible（poli）&&polisprimitive（poli\\米歇尔·马库斯，2013年5月27日
（鼠尾草）
定义is_primitive（p）：
d=2^（p.degree（））-1
如果不是p.divides（x^d-1）：返回False
对于k in（d//q对于q in d.prime_factors（））：
如果p.divides（x^k-1）：返回False
return True
P.<x>=GF（2）[]
对于范围（11000）内的n：
p=x^（2*n+2）+x+1
如果p.为_irreducible（）且为_primitive（p）：
打印（n）
| 