安德鲁·霍罗伊德,<a href=“/A002830号/b002830号_1.txt“>n、a(n)表,n=0..200</a>
检验过的
经核准的
OEIS服务器:已将新的b文件安装为b002830.txt。旧的b文件现在是b002830_1.txt。
提出
编辑
seq(n)={Vec(substpol(prod(k=1,2*n,g(2*n、k)),x^2,x))}\\安德鲁·霍罗伊德2017年12月14日;于2002年5月更新, 2023
1,1,5, 16, 86, 448, 3580, 34981, 448628, 6854130, 121173330, 2403140605, 52655943500, 1260724587515, 32726520985365, 915263580719998, 27432853858637678, 877211481667946811, 29807483816421710806, 1072542780403547030073, 40739888428757581326987
1,2
0,3
安德鲁·霍罗伊德,<a href=“/A002830号/b002830号_1.txt“>表格n,a(n)代表n=10..30200</a>
通用条款:exp(总和(总和_{k个 >= 1} F(x^k)/k, k个 >= 1)其中F(x)是指A002831号. -肖恩·欧文2014年9月9日
permcount(v)={my(m=1,s=0,k=0,t);对于(i=1,#v,t=v[i];k=if(i>1&&t==v[i-1],k+1,1);m*=t*k;s+=t);s!/m}
b(k,q个第页)={如果(q个k个%2,如果(k个第页%2,0,我的(j=k个第页/2);q个k个^j*(2*j)/(j!*2^j)),总和(j=0,k个第页\2,二项式(k个, 第页, 2*j)*q个k个^j*(2*j)/(j!*2^j))}
pm(v)={my(p=总和(i=1,#v,x^v[i]));prod(i=1,极度(p),b(polceoff(p,i),i))}
g(n,k)={和(r=0,n\k,x^(k*r)*b(k,r)^3/(r!*k^r))+O(x*x^n)}
一序列(n) = {我的Vec公司(substpol公司(秒=0); 部分触头(第页k个=1, 2*n,秒+=permcount公司(第页)*下午(第页)^三); 秒/克(2个*个,k个)), x^2, x))!} \\安德鲁·霍罗伊德2017年12月14日; 已更新 五月 02, 2023
a(0)=1由安德鲁·霍罗伊德2023年5月2日
安德鲁·霍罗伊德:更好的程序和更大的b文件。