%I#31 2021年6月17日09:13:46

%S 1,2,3,4,5,0,6,7,8,0,9,0,10,0,11,12,0,13,0,0,0，14,0,15,0,16,0,17,0，

%T 18,19,0,0,0，0,0,10,00,02,0,21,0,22,0,00,23,0,24,0,,0,025,0,0-0,0.0,0,26,0，

%U 27,0,2,28,0,0,29,0,0，0,30,0,31,0,0.0,0,032,0,33,0,34,0,0-0,0,10,35,0

%如果N=A000961（k）（素数的幂），则N a（N）=k；如果N不在A000961。

%C该名称已被编辑，以澄清指数k指的是A000961（“素数幂”={1}U A246655），而不是正确素数幂的列表A246655_M.F.Hasler，2021年6月16日

%H Reinhard Zumkeller，<a href=“/A095874/b095874.txt”>n的表，a（n）表示n=1.10000</a>

%F a（n）=和{1<=k<=n}A010055（k）；[由M.F.Hasler_于2021年6月15日更正]

%F a（n）=A065515（n）*（A065516（n）-A065515（n-1））。

%F a（n）=A065515（n）*A069513（n）.-_M.F.Hasler，2021年6月16日

%t连接[{1}，模块[{k=2}，表[If[PrimePowerQ[n]，k；k++，0]，{n，2100}]]（*哈维·P·戴尔，2020年8月15日*）

%o（哈斯克尔）

%o a095874 n | y==n=长度xs+1

%o |否则=0

%o其中（xs，y:ys）=跨度（<n）a000961_list

%o--_Reinhard Zumkeller_，2012年2月16日，2011年6月26日

%o（PARI）a（n）=if（isprimepower（n），sum（i=1，logint（n，2），primepi（sqrtnint（n，i））+1，n==1）\\_Charles R Greathouse IV_，2015年4月29日

%o（PARI）{M95874=Map（）；A095874（n，k）=if（mapisdefined（M95874，n，&k），k，isprimepower（n），mapput（M9587，n，k=sum（i=1，指数（n）），primepi（sqrtnint（n，i）））+1）；k，n==1）}\\带记忆的变量，可能对计算A097621，A344826和相关变量有用。如果已知n是素数的幂，则可以省略“isprimepower（n）”（可能需要因式分解）和“，n==1”，即为A000961求左逆_M.F.Hasler，2021年6月15日

%Y参见A000961（右反转）、A049084、A097621。

%K nonn公司

%O 1,2号机组

%A _Reinhard Zumkeller，2004年6月10日

%E编辑：M.F.Hasler，2021年6月15日