记录编号 |
358227 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[vijos1889]天真的因数分解 |
最终得分 |
100 |
用户昵称 |
sxysxy |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
0.561 s |
提交时间 |
2016-12-15 07:32:35 |
内存使用 |
8.89 MiB |
显示代码纯文本
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <list>
#include <queue>
#include <vector>
#include <cmath>
#include <iostream>
using namespace std;
#define MAXN 1000001
int miu[MAXN];
bool vis[MAXN];
int prime[MAXN];
int tot;
void pre()
{
miu[1] = 1;
for(int i = 2; i <= MAXN; i++)
{
if(!vis[i])
{
miu[i] = 1;
prime[tot++] = i;
}
for(int j = 0; j < tot; j++)
{
int k = prime[j]*i;
if(k > MAXN)break;
vis[k] = true;
if(i%prime[j])miu[k] = -miu[i];
else break;
}
}
//for(int i = 1; i <= MAXN; i++)miu[i] += miu[i-1];
}
typedef long long LL;
LL calc(LL x)
{
LL ans = 0;
LL t = sqrt(x+1.0);
for(LL i = 2; i <= t; i++)
ans += x/(i*i)*miu[i];
return ans;
}
int main()
{
freopen("naive.in", "r", stdin);
freopen("naive.out", "w", stdout);
LL k;
scanf("%lld", &k);
pre();
LL l = k, r = 1e11;
while(l < r)
{
LL m = (l+r)>>1;
if(calc(m) >= k)r = m;
else l = m+1;
}
printf("%lld\n", r);
return 0;
}