# BZOJ 2440: [中山市选2011]完全平方数

/*
BZOJ 2440: [中山市选2011]完全平方数

二分答案+莫比乌斯函数
*/
#include <cstdio>
#include <iostream>

#define rg register
{
rg char c = getchar ();
for (n = 0; !isdigit (c); c = getchar ());
for (; isdigit (c); n = n * 10 + c - '0', c = getchar ());
}

#define Max 1000006
int mu[Max], p[Max]; bool is[Max];

void Euler (int N)
{
rg int i, j, k; int C = 0; mu[1] = 1;
for (i = 2; i <= N; ++ i)
{
if (!is[i]) p[++ C] = i, mu[i] = -1;
for (j = 1; j <= C; ++ j)
{
if ((k = i * p[j]) > N) break;
is[k] = true;
if (i % p[j] == 0) { mu[k] = 0; break; }
else mu[k] = -mu[i];
}
}
}

#define INF 2e9
typedef long long LL;
LL Calc (LL N)
{
LL res = 0;
for (rg int i = 1; i * i <= N; ++ i)
res += N / (i * i) * mu[i];
return res;
}
int main (int argc, char *argv[])
{
int T, K; read (T); rg LL l, r, m, res;
Euler (Max - 1);
for (; T; -- T)
{
l = K, r = INF;
for (; l <= r; )
{
m = l + r >> 1;
if (Calc (m) >= K) r = m - 1, res = m;
else l = m + 1;
}
printf ("%lld\n", res);
}
return 0;
}

myj 吊打我Orz，xxy 捆起来打我Orz，myl 文化课上天Orz， lrh 姿势水平敲高Orz， hkd 特别胖Orz%%%，cys 智商感人Orz，syl zz专业Orz，我没有学上， 我们未来一片光明
posted @ 2017-12-20 13:54  ZlycerQan  阅读(167)  评论(0编辑  收藏  举报