提出
经核准的
编辑
罗伯特·C·莱昂斯谢谢你,米歇尔。
罗伯特·C·莱昂斯:我还需要做其他事情来批准这个序列吗?
(PARI)i=0;a=2;打印1(a,“,”);while(i<100,对于素数(p=2,10^6,if(Mod(a,p^2)^(p-1)==1&p%2=0&&((a-1)%p^2)&&(a+1)%p*2),打印1(p,“,”);i++;a=p;中断({n=1}))\\米歇尔·马库斯2023年1月21日
米歇尔·马库斯:是的,这与A244550不同;所以我只是把pari带到了那里,并把它改成在这里工作(但用同样的外观,以便看到不同之处)
(Python)
从sympy导入nextprime
从gmpy2导入powmod
最大n=45
a=2
seq=[a]
对于范围(2,max_n+1)中的i:
p=2
为True时:
p_squared=p*p
如果powmod(a,p-1,p_squared)==1和(a-1)%p_squaled!=0和(a+1)%p_squared!=0:
顺序追加(p)
a=p
打破
其他:
p=下一素数(p)
打印(seq)
囊性纤维变性。A178871号,17.94万澳元,A179678号,A244550型,A297846型.
罗伯特·C·莱昂斯:我添加了一个Python程序。它使用维基百科文章中指定的定义/算法计算术语(https://en.wikipedia.org/wiki/Wiefrich_prime#Wieferich_sequence). 程序的输出与维基百科文章中的顺序相匹配。
米歇尔·马库斯:你看到A244550了吗?我想知道为什么你会得到不同的东西
罗伯特·C·莱昂斯:明天,我将尝试编写一个Python程序。
囊性纤维变性。A178871号,A179400个,A179678号,A297846型.
罗伯特·C·莱昂斯:我添加了外部参照。我没有节目。谢谢。
囊性纤维变性。A178871号,A179400个,179678英镑.
米歇尔·马库斯:我有些不明白;我们看到a(4)=20771和a(5)=18043;所以20771应该是A179400的术语,18043应该是A1179678的术语??
罗伯特·C·莱昂斯我认为你是对的,米歇尔。我无法验证(甚至无法理解:-)A179400和A179678的PARI程序。
米歇尔·马库斯:实际上我不知道
米歇尔·马库斯:是否可以添加外部参照A178871、A297846?
米歇尔·马库斯:A179400和A179678:都表示“素数序列递增”;但在这里没有增加
米歇尔·马库斯你有节目吗?
2, 1093,5,20771,18043,5,20771,18043,5,20771,18043,5,20771,18043,5,20771,18043,5,20771,18043,5,20771,18043,5,20771,18043,5,20771,18043, 5, 20771, 18043, 5, 20771, 18043, 5, 20771, 18043, 5, 20771, 18043, 5, 20771, 18043
这个启动 具有 一(三),这个该序列是周期性的,具有以下周期,即威弗里奇三联体:52077118043。