2440: [中山市选2011]完全平方数

题目链接

题目大意：求第k个不含有平方数因子(这里不包括1)的数

题解：明显没法直接求，二分一下，需要求1–mid中不含有平方数因子的数的个数

$无平方因子->\mu(i)!=0~~莫比乌斯函数大力容斥一波$

$枚举\sqrt{mid}中的所有质数，Num=0个质数平方的倍数的数量(1的倍数)-1个质数平方的倍数的数量(9,25...的倍数)+2个质数平方的倍数的数量(36,100...的倍数)...$

$对于一个数x，x的倍数对答案的贡献符号为 μ(x)(偶+奇-) ，出现次数自然是[mid/x]$

于是有

$$cnt=\sum\limits^{\lfloor{\sqrt{mid}}\rfloor}_{i=1}(\mu(i)\lfloor{\frac{mid}{i^2}}\rfloor)$$

我的收获：莫比乌斯函数

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
 
#define ll long long
const int N=44723;
 
int T,k;
int p[N+5],u[N+5];
bool vis[N+5];
ll ans;
 
void euler()
{
    u[1]=1;
    for(int i=2;i<=N;i++){
        if(!vis[i]) p[++p[0]]=i,u[i]=-1;
        for(int j=1;j<=p[0]&&p[j]*i<=N;j++){
            vis[p[j]*i]=1;
            if(i%p[j]==0){u[p[j]*i]=0;break;}
            u[p[j]*i]=-u[i];
        }
    }
}
 
bool check(ll x)
{
    int ret=0;
    for(ll i=1;i*i<=x;i++)
        ret+=u[i]*(x/(i*i));
    return ret>=k;
}
 
void work()
{
    cin>>k;
    for(ll l=1,r=k<<1;l<=r;){
        ll mid=l+r>>1;
        if(check(mid)) r=mid-1,ans=mid;
        else l=mid+1;
    }
    printf("%lld\n",ans);
}
 
int main()
{
    euler();
    cin>>T;
    while(T--) work();
    return 0;
}