(最大值)a(n):=3*n*求和(二项式(n+1,j)*求和(二项式(j,i-j)*2^(4*j-2*n-i-2)*6^(n+i+1)*二项(n-j+1,3*j-n-i-2),i,j,n+j)/6^(3*j),j,0,(n+1))/*弗拉基米尔·克鲁奇宁,2011年11月21日*/
(Python)
#打印序列中的第一个num_entries条目
导入数学,sympy;x=符号('x')
k=3;num_entries=64
P=范围(k+1);eP=总和(P]中d的[x**d/math.阶乘(d));r=[0,1];curr_pow=eP
对于范围(1,num_entries-1)中的术语:
curr_pow=(curr_pow*eP).expand()
r.append(curr_pow.coeff(x**term)*math.factorial(term))