2021年5月18日火曜日

210518

PARI


Sum_{n = 1..24} n^2 = 70^2(6)

n-gonal number であり、n-gonal pyramidal number である数について
比較的小さなもの(ただし0, 1 を除く)を出力してみた。

a(n) = for(k=2, n, if(ispolygonal(m=k*(k+1)*((n-2)*k-n+5)/6, n), print(n, ": ", [k, m])));
for(n=3, 100000, a(n))

出力結果
3: [3, 10]
10: [5, 175]
14: [6, 441]
30: [17, 23001]
43: [33, 245905]
50: [34, 314755]
88: [15, 48280]
276: [26, 801801]
322: [28, 1169686]
823: [113, 197427385]
2378: [103, 432684460]
2386: [420, 29437553530]
9525: [2169, 16195753597485]
16420: [6936, 913053565546276]
19605: [1191, 5519583702676]
31265: [259, 90525801730]
31368: [14858, 17147031694579605]
83135: [1310, 31148407558500]

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。