mobius函数的计算

计算n的mobius函数u(n)

定义：将n分解成素数幂的乘积：n=p1^i1*p2^i2*...*pr^ir

u(n) = { 1 (n==1)

           0  (乘积式中至少有一个ir>1)

          (-1)^r (i1==i2==...==ir==1)

算法实现

const int n = 1 << 20;
int mu[n];
int getMu(int n)
{
 for (int i = 1;i <= n;i++)
 {
  int target = i == 1 ? 1 : 0;
  int delta = target - mu[i];
  mu[i] = delta;
  for(int j = i+i;j <= n;j += i)mu[j] += delta;
 }
}