VIJOS P1889 天真的因数分解

描述
小岛: 什么叫做因数分解呢?
doc : 就是将给定的正整数n, 分解为若干个素数连乘的形式.
小岛: 那比如说 n=12 呢?
doc : 那么就是 12 = 2 X 2 X 3 呀.
小岛: 呜呜, 好难, 居然素数会重复出现, 如果分解后每一个素数都只出现一次, 我就会.
wish: 这样来说, 小岛可以正确分解的数字不多呀.
doc : 是呀是呀.
wish: 现在问题来了, 对于给定的k, 第 k 个小岛无法正确分解的数字是多少?

---

【分析】
容斥原理+莫比乌斯反演
其实就是把莫比乌斯函数中1和-1的值颠倒
代码里有详解
我自己对于最后的二分阶段还有些问题

---

【代码】

//VIJOS 1889 天真的因数分解 
//莫比乌斯反演！（即1和-1交换） 
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define fo(i,j,k) for(i=j;i<=k;i++)
#define LL long long
using namespace std;
const int n=160000;
long long k;
int miu[250050],b[250050],p[250050];
inline void gm()  //线性筛法求素数，并求出莫比乌斯函数的值（反的） 
{
    int i,j;
    for(i=2;i<=n;i++)
    {
        if(!b[i])
        {
            p[++p[0]]=i;
            miu[i]=1;
        }
        for(j=1;j<=p[0] && i*p[j]<=n;j++)
        {
            b[i*p[j]]=1;
            if(i%p[j]==0)
            {
                miu[i*p[j]]=0;
                break; 
            } //有平方因子的特判 
            else miu[i*p[j]]=-miu[i];  //miu[i*p[j]]比miu[i]多一个质因子p[j] 
        }
    }
}
LL chai(LL mid)
{
    LL t=sqrt(mid),sum=0,i;
    fo(i,1,t) sum+=mid/(i*i)*miu[i];
    return sum;
}
int main()
{
    scanf("%lld",&k);
    gm();
    LL l=k,r=25505460948LL,i,j,ans;
    while(l<=r)
    {
        LL mid=(l+r)>>1;
        if (chai(mid)>=k) ans=mid,r=mid-1;
        else l=mid+1;
    }
    printf("%lld\n",ans);
    return 0;
}