EulerInvTransform[{}]={};EulerInvTransform[seq_]:=模块[{final={}},对于[i=1,i<=长度[seq],i++,AppendTo[final,i*seq[[i]]-总和[final[[d]]*seq[[i-d]],{d,i-1}]];
表[Sum[MoebiusMu[i/d]*final[[d]],{d,Divisors[i]}]/i,{i,Length[seq]}]];
EulerInvTransform[表[DivisorSigma[1,n],{n,30}]]
|