|
mob[m_，n_]：=如果[Mod[m，n]==0，MoebiusMu[m/n]，0]；
EULERi[b_]：=模[{a，c，i，d}，c={}；对于[i=1，i<=长度[b]，i++，c=Append[c，i*b[i]]-求和[c[[d]]*b[[i-d]]，{d，1，i-1}]]；a={}；对于[i=1，i<=长度[b]，i++，a=Append[a，（1/i）*Sum[mob[i，d]*c[[d]]，{d，1，i}]]；返回[a]]；
b[n_，i_]：=b[n，i]=如果[n==0，{0}，如果[i<1，{}，展平@表[Map[Function[{p}，p+j*x^i]，b[n-i*j，i-1]]，{j，0，n/i}]]；
g[n_，k_]：=g[n，k]=总和[Sum[2^Sum[Sum[GCD[i，j]*系数[s，x，i]*系数[t，x，j]，{j，1，指数[t，x]}]，{i，1，指标[s，x]{]/乘积[i^系数[s、x、i]*系数值[s，x，i]！，{i，1，指数[s，x]}]/乘积[i^系数[t，x，i]*系数[t、x、i]！，{i，1，指数[t，x]}]，{t，b[n+k，n+k]}，{s，b[n，n]}]；
A[n_，k_]：=g[最小值[n，k]，绝对值[n-k]]；
b[d_]：=总和[A[n，d-n]，{n，0，d}]；
| 