isPent(t)=P(平方(t<<1\3)+1)==t
对于(i=1299,对于(j=1,(i+1)\sqrt(2),isPent(P(i)-P(j))&&print1(P(i)“,”)||next(2))
/*下面的速度要快得多,代价是实现sum2sqr(),cf。A133388号*/
A136117next(i)=i=平方(i\3*2)*6+5;直到(0,对于(j=2,#t=sum2sqr((i+=6)^2+1),t[j]%6==[5,5]&&break(2)));i^2\24
A136117vect(n,i)=矢量(n,j,i=A136117next(i))/*第二个参数默认为0,但允许从其他位置开始*/
|