%I#7 2021年1月25日14:41:26

%第7,8,24,25,29,31,36,42,49,59,65,66,69,73,74,76,80,83,93,94,99102页，

%电话：104110117118121122123124126127138140144149150152161，

%电话：163168169174175178181185188190195203209210212213

%N五边形数（A000326）的指数，是其他两个正五边形数字的总和。

%F A000326（a（n））=A000326。

%e a（1）=7，因为P（7）=70是最小的五边形数，可以写为其他两个五边形数字的和，P（7。

%o（PARI）P（n）=n*（3*n-1）>>1/*也称为A000326*/

%o是Pent（t）=P（平方（t<<1\3）+1）==t

%o表示（i=1999，表示（j=1，（i+1）\sqrt（2），isPent（P（i）-P（j））&print1（i“，”）||next（2））

%o/*以下速度要快得多，代价是实现sum2sqr（），参见A133388*/

%o A136116next（i）=i=6*i-1；直到（0，对于（j=2，#t=sum2sqr（（i+=6）^2+1），t[j]%6==[5,5]&&break（2）））；i\6+1

%o A136116vect（n，i=0）=矢量（n，j，i=A136116next（i））

%o A136116（n，i=0）=直到（！n--，i=A136116next（i））；2007年12月25日，i.F.Hasler

%Y参考A000326，A136112-A136118。

%K nonn公司

%O 1,1号机组

%A _M.F.Hasler_，2007年12月15日；2007年12月25日更正