【BZOJ】2440 [中山市选2011]完全平方数 容斥

题目传送门

杜教筛求$\mu^2$的前缀和什么的真的不会啊，还好有数学拯救我于水火之中。

考虑容斥：x以内的无平方因子数=0个质数乘积的平方的倍数的数的数量(1的倍数)-每个质数的平方的倍数的数的数量(9的倍数,25的倍数,…)+每2个质数乘积的平方的倍数的数的数量(36的倍数,100的倍数,…)-…

容易发现每一项$a$的正负恰好是$\mu(a)$（其实是根据莫比乌斯函数的性质），那么$[1,x]$间无平方因子的数的个数就是$\sum_{i=1}^{\lfloor \sqrt x \rfloor}\mu(i)\lfloor \frac{x}{i^2} \rfloor$。于是预处理后二分即可。

p.s.RE什么鬼……l+r爆int了……

附上AC代码：

#include <cstdio>
using namespace std;

const int N=5e4+10;
int t,n,p[N],miu[N],num,l,r,mid,ans;
bool b[N];

inline void init(int n){
    miu[1]=1;
    for (int i=2; i<=n; ++i){
        if (!b[i]) p[++num]=i,miu[i]=-1;
        for (int j=1; j<=num&&p[j]*i<=N; ++j){
            b[p[j]*i]=1;
            if (i%p[j]==0) {miu[p[j]*i]=0;break;}
            else miu[p[j]*i]=-miu[i];
        }
    }
    return;
}

inline int check(int mid){
    int ret=0;
    for (int i=1; i*i<=mid; ++i) ret+=miu[i]*(mid/i/i);
    return ret;
}

int main(void){
    init(5e4);
    for (scanf("%d",&t); t; --t){
        scanf("%d",&n),l=1,r=2e9,ans=0;
        while (l<=r)
            if (check(mid=l+(r-l>>1))<n) l=mid+1;
            else r=mid-1,ans=mid;
        printf("%d\n",ans);
    }
    return 0;
}