|
功能[l]：=
如果[Range[Length[l]].l>长度[l]，{}，len=Length[1]；
选择[排列[l]，#。范围[len]==len&]]
函数b[s_，m_]：=模块[{l=0}，
如果[m+s==0，1，
如果[m s==0，0，
如果[m>=s，
如果[m>s，0，1]，
如果[m==1，CatalanNumber[s]^2，
如果[s-m==1,4m，
我=
扁平[Map[functiony，Integer Partitions[s+m，{s}]-1]，1]；
地图[Times@@#&，
地图[Map[r，Range[1，s]]^#&，
l] ]。（映射[Times@@#&，映射[Factorial，l]]）^（-1）*m！]
]
]
]
]
]
a[n_，s_]：=和[二项式[n，m]b[s，m]，{m，1，s}]
b[s，m]：=
如果[s+m>0，表1[[s+1，m+1]]，如果[s+m==0，1，0]]
f[n，k_]：=
k（2n+2k-
4)! （2 n+k-1）/（（n+k-1）！（n+k-2）！不！（n+k）！）-
求和[table2[[s+1]]a[k+s，n-s]，{s，0，n-1}]
r[n]：=（二项式[2 n，n]）^2/（n^2+2 n+1）
答案[n]：=f[n，3]
指数=24；
table1=表[functionb[s，m]，{s，0，index}，{m，0，index}]；
表2=范围[索引]；
对于[i=2，i<=索引，i++，表2[[i]]=f[i-1，3]]；
f[索引-1，3]
| 