a89[n_]：=a89[n]=乘积[{p，e}=pe；其中[p<3&e==1，1，p==2&e>1，0，Mod[p，4]==1、2，Mod[p，4]==3，0，True，a89[p^e]]，{pe，FactorInteger[n]}]；
a86[n_]：=a86[n]=乘积[{p，e}=pe；其中[p==1|p==3&e==1，1，p==3&r&e>1，0，Mod[p，3]==1、2，Mod[p,3]==2，0，True，a86[p^e]]，{pe，FactorInteger[n]}]；
a1615[n_]：=n和[MoebiusMu[d]^2/d，{d，除数[n]}]；
a1616[n_]：=和[EulerPhi[GCD[d，n/d]]，{d，除数[n]}]；
a1617[n]：=1+a1615[n]/12-a89[n]/4-a86[n]/3-a1616[n]/2；
