2021年2月23日火曜日

210223

PARI


Number of perfect matchings in the graph C_{2n} X C_k

kが3以上のとき、計算してみた。

a341533(n, k) = round(sqrt(prod(a=1, n,   prod(b=1, k,   4*sin((2*a-1)*Pi/(2*n))^2+4*sin((2*b-1)*Pi/k)^2))));
a341738(n, k) = round(sqrt(prod(a=1, n,   prod(b=1, k-1, 4*sin((2*a-1)*Pi/(2*n))^2+4*sin(2*b*Pi/k)^2))));
a341739(n, k) = round(     prod(a=1, n-1, prod(b=1, k,   4*sin(a*Pi/n)^2          +4*sin((2*b-1)*Pi/(2*k))^2)));
T(n, k) = a341533(n, k)/2+a341738(n, k)+2*((k+1)%2)*a341739(n, ceil(k/2));

K=7;
for(n=2, 10, for(k=3, K, print1(T(n, k), ", ")); print)

出力結果
50, 272, 722, 3108, 10082, 
224, 3108, 9922, 90176, 401998, 
1058, 39952, 155682, 3113860, 19681538, 
5054, 537636, 2540032, 114557000, 1034315998, 
24200, 7379216, 41934482, 4357599552, 55820091938, 
115934, 102144036, 694861522, 169094614280, 3044533460992, 
555458, 1418981392, 11527389122, 6640849944580, 166779871224962, 
2661344, 19742287908, 191304901282, 262671237617216, 9152970837103102, 
12751250, 274848860432, 3175220160032, 10433179552323108, 502711247500143362, 

0 件のコメント:

コメントを投稿

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