%I#17 2021年6月8日02:23:42

%S 1,3,21,57219399903153326674161707717157185432912732547，

%电话：506736591974109961411137991114681134463194691227829304311，

%电话：3523173839073891934630715166216729877934986351789877993279912524611353639

%N将k编号为k^3除以4^（k^2）-1。

%C自2020年8月13日起（开始）

%C除1外，所有项都可以被3整除，但不能被5或9整除。

%C所有>3的项都可以被7、19和73中的至少一项整除。

%C所有条款都是自由的吗？（结束）

%H Robert Israel，n表，n=1..500的a（n）（前49个术语来自Robert Price）

%p过滤器：=n->4&^（n^2）-1模（n^3）=0:

%p选择（过滤器，[1，序列（i，i=3..10^6,6）]）；#_罗伯特·伊斯雷尔，2020年8月13日

%t k=2；Do[p=质数[k]；如果[InterQ[（PowerMod[p+1，n^2，n^3]-1）/n^3]，打印[｛k，p，n｝]]，｛n，1200000｝]

%t k=2；Do[p=质数[k]；如果[IntegerQ[（PowerMod[p+1，n^2，n^3]-1）/n^3]，打印[{k，p，n}]]，{n，1000000}]（*_Robert G.Wilson v_，2007年4月6日*）

%t加入[{1}，选择[Range[3000000]，PowerMod[4，#^2，#^3]==1&]]（*_Robert Price_，2020年3月31日*）

%Y参考A014945（数字k除以4^k-1）。

%Y参考A127104（数字k，k^2除以4^k-1）。

%Y参考A128678（数字k，k^3除以4^（k^2）+1）。

%K nonn公司

%O 1,2号机组

%A_Alexander Adamchuk，2007年4月3日

%E更多条款摘自_Robert G.Wilson v_，2007年4月6日