2023年9月8日金曜日

230908

PARI


A357782, A357783 and A357784

出力してみた。

a(n) = if(n==0, 1, 2 * sum(k=0, n-1, binomial(n-1,k)*c(k)));
b(n) = if(n==0, 0,     sum(k=0, n-1, binomial(n-1,k)*a(k)));
c(n) = if(n==0, 0,     sum(k=0, n-1, binomial(n-1,k)*b(k)));

for(n=0, 20, print1(a(n),", "))
for(n=0, 20, print1(b(n),", "))
for(n=0, 20, print1(c(n),", "))

出力結果
1, 0, 0, 2, 12, 50, 184, 686, 2996, 16642, 110328, 784190, 5645876, 40685762, 296458344, 2226254766, 17564381332, 147289101090, 1312394060536, 12305546886398, 119906479624084, 
0, 1, 1, 1, 3, 21, 131, 705, 3515, 17389, 91739, 547889, 3746227, 28241373, 224124083, 1821051233, 15023818091, 126366334125, 1094358852075, 9858890038513, 92983173940419, 
0, 0, 1, 3, 7, 17, 61, 343, 2231, 14301, 88561, 542011, 3397483, 22638993, 164336085, 1299899087, 10991061663, 97070035205, 881323166809, 8173386231395, 77489746906355, 

0 件のコメント:

コメントを投稿

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