%I#20 2024年5月2日17:33:24

%电话：974871455732261275393220861155412342959105184724639193，

%电话：4037309370697317318519014728913345633017910000453154075723，

%电话：130702609304154189143584109183298237794511671058720273528452035676204134592184511589077963965154024740819880414772484029

%N最小奇素数基数q，使得p^5除以q^（p-1）-1，其中p=素数（N）。

%H Robert Israel，n的表，n=1..500的a（n）</a>

%H W.Keller和J.Richstein，<a href=“http://www1.uni-hamburg.de/RRZ/W.Keller/FermatQuotient.html“>可被p整除的费马商。

%p f：=proc（n）局部p，k，j，q，R；

%p p：=i素数（n）；

%p R:=排序（映射(rhs@操作，[msolve（q^（p-1）-1，p^5）]）；

%p代表0 do中的k

%R do中j的p

%pq:=k*p^5+j；

%p如果isprime（q），则返回q-fi；

%日期

%日期

%p端程序：

%p映射（f，[1..100]美元）；#_罗伯特·伊斯雷尔，2019年4月11日

%t Do[p=素数[n]；q=2；而[PowerMod[q，p-1，p^5]！=1，q=下一素数[q]]；打印[q]，{n，100}]（*_Ryan Propper_，2007年3月31日*）

%o（PARI）{a（n）=局部（p，x，y）；如果（n==1，返回（97））；p=prime（n）；x=znprimroot（p^5）^（p^4）；vecsort（向量（p-1，i，y=lift（x^i）；while（！isprime（y），y+=p^5

%o（Python）

%o从itertools导入计数

%o从sympy导入nthroot_mod，isprime，prime

%o定义A125646（n）：

%o m=（p:=素数（n））**5

%o r=排序（nthroot_mod（1，p-1，m，all_roots=真））

%o对于计数（0，m）中的i：

%o对于in r：

%o如果isprime（i+a）：返回i+a#_Chai Wah Wu_，2024年5月2日

%Y参考A125609、A125610、A125611、A125612、A125632、A125633、A125634、A125635、A125636、A125637、A125645、A125647、A125648、A125649。

%K nonn公司

%O 1,1号机组

%A_Alexander Adamchuk，2006年11月29日

%E 2007年3月31日_Ryan Propper_的更多条款

%E更多条款来自2007年5月30日的Max Alekseyev