登录
a(n)是最小整数,其适当除数的i次幂之和是1≤i≤n的素数,如果不存在这样的数,则为-1。
1

%I#30 2022年11月16日06:58:01

%S 4,498188298829191220291052282761881号

%N a(N)是最小整数,其适当除数的第i次幂之和是1≤i≤N的素数,如果不存在这样的数,则为-1。

%e4是一个项,因为4的严格除数是{1,2},1^1+2^1=3和1^2+2^2=5是素数,并且没有数字<4具有这个性质。

%o(PARI)

%o卡(n)=我的(c=1,s=0);s=σ(n)-n;while(isprime(s),c++;s=西格玛(n,c)-n^c);c(c)--

%o a(n)=我的(x=0);对于(k=1,+oo,x=card(k));如果(x>=n,返回(k))

%o(Python)

%o从itertools导入计数

%o来自math导入prod

%o来自sympy import isprime,factor

%o定义A358199(n):

%o对于计数(2)中的m:

%o f=因子(m).items()

%o如果全部(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))):

%o返回m#_Chai Wah Wu_,2022年11月15日

%Y参考A001065(s1)、A180852、A357324、A057709、A000203、A001157、A001158。

%Y A037020的子序列。

%K nonn,更多

%O 1,1号机组

%A _Jean-Marc Rebert,2022年11月2日