描述
小岛: 什么叫做因数分解呢?
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;
}