[BZOJ 2440]完全平方数

题目描述：

求第K个无平方因子的数

题目分析：

无平方因子数，即分解质因数后所有质因数的次数都为1的数.
首先转化为判定性问题，即1-n中有多少个数为无平方因子数
二分答案.
根据容斥原理可得
n以内的无平方因子数=
0个质数乘积的平方的倍数的数的个数(即N)
- 1个质数乘积的平方的倍数的数的个数(如4的倍数,9的倍数)
+ 2个质数乘积的平方的倍数的数的个数(如36的倍数,100的倍数)
- 3个质数乘积的平方的倍数的数的个数…
每个枚举的乘积的贡献恰为μ(a) 例如μ(3)=-1,9对答案的贡献即为负数
μ(6)=1,36对答案的贡献即为负数
F(x)=Σ(μ(i)*(x/(i*i)))(i=[1,sqrt(x)])

题目链接：

BZOJ 2440

Ac 代码:

#include <cstdio>
#include <iostream>
#include <cmath>
#define sqr(x) x*x
#define ll long long
const int maxm=100011;
int mu[maxm],prime[maxm],cnt,n;
bool vis[maxm];
void getmu()
{
    mu[1]=1,vis[1]=1;
    for(int i=2;i<=maxm;i++)
    {
        if(!vis[i])
        {
            prime[++cnt]=i;
            //printf("%d\n",i);
            mu[i]=-1;
        }
        for(int j=1;j<=cnt&&prime[j]*i<=maxm;j++)
        {
            vis[prime[j]*i]=1;
            if(i%prime[j]==0)
            {
                mu[prime[j]*i]=0;
                break;
            }
            mu[prime[j]*i]=-mu[i];
        }
    }
}
int sum(ll mid)
{
    ll ans=0;
    ll div=std::sqrt(mid);
    for(int i=1;i<=div;i++) 
     ans+=mu[i]*(mid/(i*i));
    return ans;
}
void work()
{
    scanf("%d",&n);
    ll l=1,r=(1LL<<31)-1;
    ll ans=1;
    while(l<=r)
    {
        ll mid=(l+r)>>1;
        if(sum(mid)>=n) ans=mid,r=mid-1;
        else l=mid+1;
    }
    printf("%lld\n",ans);
}
int main()
{
    getmu();
    int Q;
    scanf("%d",&Q);
    while(Q--)
     work();
    return 0;
}