%I#4 2015年5月29日06:29:36
%S 5,4,7,12,19,17,25,20,10,28,45,42,39,17,37,21,36,35,13,33,65,28,67,32,
%电话:52,40,74,31,70,85,35,16,60,70,77,68,42,30105,76,59,26,74,49115,19,
%U 125115102110,92,56103,29145100114,77,92,47,63108152,95,22116号
%N最小指数m>0,使A136117(N)-A000326(m)再次为五边形数。
%e a(1)=5是最小整数m>0,因此A136117(1)-P(m)是五边形数,即P(7)-P。
%e a(2)=4是最小整数m>0,因此A136117(2)-P(m)是五边形数,即P(8)-P。
%o(PARI)A136118vect(n,i=-1)=向量(n,k,直到(0,对于(j=2,#n=sum2sqr((i+=6)^2+1),n[j]%6==[5,5]||接下来;n=n[j];断裂(2));n[1]\6+1)/*这使用sum2sqr(),参见A133388。下面是一些简单但速度慢得多的代码*/
%o我的(P=A000326(n)=n*(3*n-1)/2,是Pent(t)=P(平方(t*2\3)+1)==t);对于(i=1299,对于(j=1,(i+1)\sqrt(2),是Pent(P(i)-P(j))&print1(j“,”)||next(2))
%Y参考A000326,A136112-A136117。
%K nonn公司
%O 1,1
%A _M.F.Hasler_,2007年12月25日
|