题意
传送门 P4318 完全平方数
题解
考虑二分答案,问题转化为区间 [ 1 , n ] [1,n] [1,n] 内非完全平方数的数量。设全集 U U U 为区间 [ 1 , n ] [1,n] [1,n] 的元素, A i A_i Ai 代表区间内第 i i i 个素数的平方数的整数倍代表的子集,答案为 ∣ U ∣ − ∣ ∪ A i ∣ |U|-|\cup A_i| ∣U∣−∣∪Ai∣,根据容斥原理,将其展开后观察每一项,发现 ∣ A i ∩ A j ⋯ ∩ A k ∣ |A_i\cap A_j\dots \cap A_k| ∣Ai∩Aj⋯∩Ak∣ 前面的系数对应的和数 i × j × k i\times j \times k i×j×k 的莫比乌斯函数,全集可以看做零个素数的平方的整数倍的集合。那么答案为 ∑ i = 1 ⌊ n ⌋ μ ( i ) ⌊ n i 2 ⌋ \sum\limits_{i=1}^{\lfloor \sqrt{n}\rfloor} \mu (i) \lfloor \frac{n}{i^2}\rfloor i=1∑⌊n⌋μ(i)⌊i2n⌋.
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <map>
#include <vector>
using namespace std;
#define maxn 1000005
typedef long long ll;
int T, K, prime[maxn >> 1], used[maxn], mu[maxn];
void sieve(int n)
{
mu[1] = 1;
int p = 0;
for (int i = 2; i < maxn; ++i)
{
if (!used[i])
{
prime[p++] = i;
mu[i] = -1;
}
for (int j = 0; j < p && i * prime[j] < maxn; ++j)
{
int k = i * prime[j];
used[k] = 1;
if (i % prime[j] == 0)
{
mu[k] = 0;
break;
}
else
{
mu[k] = -mu[i];
}
}
}
}
ll count(ll n)
{
ll res = 0;
for (ll i = 1LL; i * i <= n; ++i)
res += mu[i] * n / (i * i);
return res;
}
int main()
{
sieve(maxn);
scanf("%d", &T);
while (T--)
{
scanf("%d", &K);
ll lb = 0, ub = 1e10;
while (ub - lb > 1)
{
ll m = (lb + ub) >> 1;
if (count(m) < K)
lb = m;
else
ub = m;
}
printf("%lld\n", ub);
}
return 0;
}