1,1

因为g^3=g+1，我们有g^4=g^2+g，g^5=g^3+g^2，g^6=g^4+g^3。。。，g^（k+3）=g^（k+1）+g^k。因此，使用g和g^2，我们可以计算本原根的所有幂，类似于A003147号，其中g^（k+2）=g^，（k+1）+g^k（参见Shanks参考）。

罗伯特·伊斯雷尔，n=1..10000时的n，a（n）表

D.柄，斐波那契原根,文章结尾，光纤。夸脱。，10 (1972), 163-168, 181.

过滤器：=proc（p）局部x，r；

如果不是isprime（p），则返回false fi；

对于地图中的r（t->rhs（op（t）），[msolve（x^3-x-1，p）]）do

如果numtheory:-order（r，p）=p-1，则返回真fi

od；

假

结束进程：

选择（过滤器，[seq（i，i=3..2000，2）]）#罗伯特·伊斯雷尔2017年10月2日

selQ[p_]：=AnyTrue[PrimitiveRootList[p]，Mod[#^3-#-1，p]==0&]；

选择[Prime[Range[200]]，selQ](*Jean-François Alcover公司2020年7月29日*）

（PARI）

Z（r，p）=znorder（Mod（r，p））==p-1；\\r是否为本原根mod p

Y（p）=对于（r=2，p-2，如果（Z（r，p）&&Mod（r^3-r-1，p）==0，返回（1）））；0; \\ 测试p

对于素数（p=2，10^3，如果（Y（p），打印1（p，“，”））；

囊性纤维变性。A003147号（本原根g，使得g^2=g+1 mod p）。

囊性纤维变性。A293201型（本原根g使得g^3=g^2+g+1 mod p）。

囊性纤维变性。A104217号.

上下文中的序列：A354357型 A156104号 A191080号*A067830号 A191019号 A106862号

相邻序列：A293197型 A293198型 A293199型*A293201型 A293202型 A293203型

非n

乔格·阿恩特2017年10月2日

经核准的