u[n_,k_,r]:=(r*二项式[k*n+r,n]/(k*n+r));
e[n,k_]:=和[u[j,k,1+(n-2*j)*k/2],{j,0,n/2}]
c[n_,k_]:=如果[n==0,1,(除数和[n,EulerPhi[n/#]*二项式[k*#,#]&]+除数和[CCD[n-1,k],Euler Phi[#]*二项式[n*k/#,(n-1)/#]&])/(k*n)-二项式[k*n,n]/(n*(k-1)+1)];
T[n_,k_]:=(1/2)*(c[n,k]+如果[n==0,1,如果[OddQ[k],如果[奇Q[n],2*u[商[n,2],k,(k+1)/2],u[n/2,k,1]+u[n/2-1,k,k]],e[n,k]+如果[奇Q[n],u[商[n,2],k,k/2]]]/2])/。空->0;
a[n_]:=T[n,4];
表[a[n],{n,0,30}](*Jean-François Alcover公司2018年6月14日,继Andrew Howroyd和A303929型*)
|