51nod 1204 莫比乌斯函数之和

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/AC_jie/article/details/82655631

今天来做莫比乌斯，这个题就是裸的莫比乌斯，但是由于如果 想打表做的话，并不是网上写的那样就做不了，虽然说不能完全打出素数表，但是仔细想一个的话，我们可以先做一个1e6的素数表，然后进行素数分解，如果不断进行分解后，N的值还是不为1的话，这就说明N是一个很大的素数，因为在给定的数据范围内每个数都被它的最小素因子筛去，如果一个合数的最小素因子是1e6这个级别的，那样这个数就超出了数据范围，综上，肯定存在一个我们用1e6的素数表筛不到的大素数。这样单独判断一个就好

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
typedef long long LL;
const int MAXN = 1e7 + 10;
LL prime[MAXN];
void getPrime()
{
    memset(prime,0,sizeof(prime));
    for(int i = 2;i < MAXN;++i)
    {
        if(!prime[i]) prime[++prime[0]] = i;
        for(int j = 1;j <= prime[0] && prime[j] < MAXN / i; ++j)
        {
            prime[prime[j] * i] = 1;
            if(i % prime[j] == 0) break;
        }
    }
}
int miu(int n)
{
    int ans = 1;
    for(int i = 1;i <= prime[0] && n != 1;++i)
    {
        int k = 0;
        while(n % prime[i] == 0)
        {
            if(k == 1) return 0;
            n /= prime[i];
            k++;
            ans *= -1;
        }
    }
    if(n > 1) ans *= -1; // 判断是否是一个大素数
    return ans;

}
int main()
{
    ios::sync_with_stdio(false);
    getPrime();
   // cout << prime[prime[0]] << endl;
    LL n;
    while(cin >> n)
    {
        cout << miu(n) << endl;
    }
    return 0;
}

反思

1. 筛法写的不熟。
2. 主观想象没有推理。记忆空间有误。