|
b1[m_，n_]：=总和[EulerPhi[c]*EulerPhi[d]*2^（m*n/LCM[c，d]），{c，除数[m]}，{d，除数[n]}]/（4*m*n）；
b2a[m_，n_]：=如果[OddQ[m]，2^（（m+1）*n/2）/（4*n），（2^；
b2b[m_，n_]：=除数和[n，如果[#>=2，EulerPhi[#]*2^（（m*n）/#），0]&]/（4*n）；
b2c[m_，n_]：=如果[OddQ[m]，总和[If[OddQ[n/GCD[j，n]]，2^（（m+1）*GCD[j，n]/2）-2^（m*GCD[j，n]），0]，{j，1，n-1}]/（4*n）2^（m*GCD[j，n]+1），0]，{j，1，n-1}]/（8*n）]；
b2[m，n]：=b2a[m，n]+b2b[m，n]+b2c[m，n]；
b3[m，n]：=b2[n，m]；b4oo[m，n]：=2^（（m*n-3）/2）；
b4eo[m，n]：=3*2^（m*n/2-3）；b4ee[m，n]：=7*2^（m*n/2-4）；
a[m_，n_]：=模[{b}，If[OddQ[m]，If[CodQ[n]，b=b4oo[m，n]，b=b4eo[m，n]]，If[OddQ[n]、b=b4meo[m、n]，b=b4ee[m，n]]；b+=b1[m，n]+b2[m，n]+b3[m，n]；返回[b]]；
a[0]=1；a[n]：=a[n，n]；
| 