VIJOS 1889 天真的因数分解

最新推荐文章于 2018-05-10 11:08:00 发布

Orion_Rigel

最新推荐文章于 2018-05-10 11:08:00 发布

阅读量369

点赞数

分类专栏：莫比乌斯文章标签： vijos

本文链接：https://blog.csdn.net/Orion_Rigel/article/details/51868434

版权

莫比乌斯专栏收录该内容

2 篇文章 0 订阅

订阅专栏

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

【题目分析】
感人的莫比乌斯反演。听起来难，但是写起来还是比较好写的。（代码借鉴自神犇hzw）莫比乌斯函数其实就是统计了质因子的个数。（1 -奇数 2 - 偶数）然后这样子再用容斥定理。然后二分就可以得到答案。

【代码】

#include <cstdio>
#include <cmath> 
long long ans,k,T,n=160000;
int tot,mu[200005],p[200005],ma[200005];
inline void getmu()
{
    for (int i=2;i<=n;++i)
    {
        if (!ma[i]) p[++tot]=i,mu[i]=1;
        for (int j=1;j<=tot&&p[j]*i<=n;++j)
        {
            ma[i*p[j]]=1;
            if (i%p[j]==0){mu[i*p[j]]=0;break;}
            else mu[i*p[j]]=-mu[i];
        }
    }
}
long long te(long long x)
{
    long long sum=0,t=sqrt(x);
    for (long long i=1;i<=t;++i) sum+=x/(i*i)*mu[i];
    return sum;
}
int main()
{
    getmu();
    scanf("%lld",&k);
    long long l=k,r=25505460948LL;
    while (l<=r)
    {
        long long mid=(l+r)>>1;
        if (te(mid)>=k) ans=mid,r=mid-1;
        else l=mid+1;
    }
    printf("%lld\n",ans);
}