VIJOS P1889 天真的因数分解

描述

小岛: 什么叫做因数分解呢?
doc : 就是将给定的正整数n, 分解为若干个素数连乘的形式.
小岛: 那比如说 n=12 呢?
doc : 那么就是 12 = 2 X 2 X 3 呀.
小岛: 呜呜, 好难, 居然素数会重复出现, 如果分解后每一个素数都只出现一次, 我就会.

wish: 这样来说, 小岛可以正确分解的数字不多呀.
doc : 是呀是呀.
wish: 现在问题来了, 对于给定的k, 第 k 个小岛无法正确分解的数字是多少?

格式

输入格式

输入只有一行, 只有一个整数 k.

输出格式

输出只有一行, 只有一个整数, 表示小岛无法正确分解出来的第k个数字.

样例1

样例输入1

10

样例输出1

27

限制

对于30%的数据, k <= 2,000,000
对于100%的数据, 1 <= k <= 10,000,000,000

提示

前 10 个小岛无法正确分解出来的数字依次是: 4 8 9 12 16 18 20 24 25 27

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

莫比乌斯函数+容斥原理～

循环时ij要设成long long型～

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
#define ll long long
#define kk 250050

const int n=160000;

ll k;
int u[kk],q[kk];
bool b[kk];

inline void ya()
{
    for(ll i=2;i<=n;i++)
    {
        if(!b[i])
        {
            q[++q[0]]=i;u[i]=1;
        }
        for(ll j=1;j<=q[0] && i*q[j]<=n;j++)
		{
            b[i*q[j]]=1;
            if(!(i%q[j]))
            {
                u[i*q[j]]=0;break;
            }
            u[i*q[j]]=-u[i];
        }
    }
}

ll qiu(ll mid)
{
    ll z=sqrt(mid),sum=0;
    for(ll i=1;i<=z;i++) sum+=u[i]*mid/(i*i);
    return sum;
}

int main()
{
    scanf("%lld",&k);
    ya();
    ll mid,ans,l,r;
    l=k;r=25505460948LL;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(qiu(mid)>=k) ans=mid,r=mid-1;
        else l=mid+1;
    }
    printf("%I64d\n",ans);
    return 0;
}