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]；
