%I#7 2014年12月12日05:04:40

%S 3,5,7,11,15,19,23,25,27,31,35,39,43,45,47,51,55,59,63,65,67,71,75,79，

%电话83,85,87,91,95,99103105107111115119123125127131135139，

%U 143145147151155159163165167171175179183185187191195199

%N数N，25除以和[Prime[k]^N，{k，1，N}]。

%C a（n）小于等于a（7）=23与A007665[n+1]=河内塔（带5个桩）重合。a（n）似乎包括所有A007665[n]={1、3、5、7、11、15、19、23、27、31、39、47、55、63、71、79、87、95、103、111、127、143、159、175、191、207、223、239、255、271、287、303、319、335、351、383、415、447、479、511、543、575、607、639、671、703、735、767、799…}，A007665[1]=1。

%C该序列中的素数包括5个素数以及形式为4k+3、A002145[n]的所有素数。术语包括10k+5（非负k），A017329[n]形式的所有数字。

%e有25个素数p<100，p（n）={2，3，5，7，11，13，17，19，23，29，31，37，41，43，47，53，59，61，67，71，73，79，83，89，97}。

%e a（1）=因为25除以总和[p（n）^3，{n，1,25}]=2^3+3^3+…+89^3+97^3=A098999[25]，不除和[p（n）^1，{n，1,25}]=A007504[25]和和[p。

%e下一个a（2）=5，因为25除以总和[p（n）^5，{n，1,25}]=A122103[25]，而不除以总和[p（n）*4，{n、1,25}]=A122102[25]。

%t选择[Range[300]，IntegerQ[Sum[Prime[k]^#1，{k，1,25}]/25]&]

%o（PARI）对于（n=1100，如果（总和（k=1,25，素数（k）^n）%25==0，打印1（n，“，”））；

%o打印；print（“不使用质数的替代方法：”）；

%o表示（n=1100，m=（n-1）%6；打印1（（n-m）*3+（n-m+如果（m>1，（m-1）*12-1，m*6-1））/3，“，”））\\_K.Spage_，2009年10月23日

%Y参考A007504、A024450、A098999、A122102、A122103、A007665。

%Y参考A002144、A002145、A017329。

%K非n

%O 1,1号机组

%A_Alexander Adamchuk，2006年8月21日，2006年9月18日，2006日9月21日