莫比烏斯反演

標籤：target   分享   int   logs   bsp   https   函數   class   技術   

學習連結：https://www.cnblogs.com/peng-ym/p/8647856.html

莫比烏斯函數線性篩代碼：

 1 bool not_prime[maxn]; 2 int prime[maxn]; 3 int Mob[maxn]; 4 void Mobius_sieve(int n) 5 { 6     int tot = 0; 7     not_prime[1] = 1; 8     Mob[1] = 1; 9     for(int i = 2; i <= n; i++)10     {11         if(!not_prime[i])prime[tot++] = i, Mob[i] = - 1;12         for(int j = 0; j < tot && 1LL * prime[j] * i <= n; j++)13         {14             not_prime[prime[j] * i] = 1;//每個合數x由它最小素因子prime[j]篩掉15             Mob[i * prime[j]] = (i % prime[j] ? -Mob[i]: 0);16             if(i % prime[j] == 0)break;//如果i % prime[j] == 0,不停止迴圈17             //那麼接下來將用prime[j+1]篩去i*prime[j+1]，但實際上應該用prime[i]篩去,因為i%prime[j]==018         }19     }20 }

 

莫比烏斯反演