|
（PARI）
卡（n）=我的（c=1，s=0）；s=σ（n）-n；while（isprime（s），c++；s=西格玛（n，c）-n^c）；c（c）--
a（n）=我的（x=0）；对于（k=1，+oo，x=card（k））；如果（x>=n，返回（k））
（Python）
从itertools导入计数
从数学导入prod
来自sympy import isprime，factor
对于计数（2）中的m：
f=因子（m）.items（）
如果全部（map（isprime，（prod（（p**（（e+1）*i）-1）//（p**i-1）for p，e in f）-m**i for i in range（1，n+1）））：
| 