%I#30 2022年8月2日17:17:25
%S 0,1,0,4,12,0,95212,0239510657,0126068375941,0494335720513691,
%电话:0271911476
%数字和为N的N位素数的个数。
%Ca(3n)=0.-_Sascha Kurz,2002年8月23日
%e a(2)=1,因为只有两位数的素数和2是11。a(5)=12:素数为10103、10211、10301、11003、12011、12101、13001、20021、20201、21011、21101和30011。
%t表[Length[Select[Prime[Range[PrimePi[10^(n-1)]+1,PrimePi[10 ^n]],Total[Integer Digits[#]]==n&]],{n,8}](*Harvey P.Dale_2011年8月9日*)
%o(Python)
%o来自sympy import isprime
%o定义nextsod(n,基础):
%o c,b,w=0,基数,0
%o为True时:
%o d=n%b
%o如果d+1<b和c:
%o收益率(n+1)*b**w+((c-1)%(b-1)+1)*b**
%o c+=d;n//=b;w+=1
%o定义a(n):
%o如果n%3==0:返回0
%o c,t=0,10**(n-1)+((n-1)%9+1)*10**((n-1)//9)-1
%o当t<10**n时:
%o如果是质数(t):c+=1
%o t=下一个sod(t,10)
%o返回c
%o打印([a(n)代表范围(1,13)中的n)]#_Michael S.Branicky_,2022年8月2日
%Y参考A069710。
%K基,nonn
%O 1,4型
%A _Amarnath Murthy,2002年8月18日
%E更多来自_Sascha Kurz的条款,2002年8月23日
%修正了E a(4)和a(5),并修正了示例,作者:Harvey P.Dale_,2011年8月9日
%2011年8月10日,Donovan Johnson_修正了E a(11)和a(13)-a(15)
%E a(16)-a(18)摘自_David Radcliffe_,2015年5月5日
%E a(19),来自Michael S.Branicky,2022年8月2日
|