洛谷4318 完全平方数

$题意:$
$求第k个不含平方因子的数$


$题解:$

$我们可以题目转为求给定n，求n以内的无平方因子数的个数，然后二分得到结果$
$求n以内的无平方因子数的个数，即求$
$$\sum _{k=1}^n{\mu }^2(k),这里\mu 是莫比乌斯函数$$

$$\begin{gathered} 我们考虑一个素数p，那么p^{的}倍数都含有平方因子，个数是\lfloor \frac{n}{p^2}\rfloor ,需要从答 \\ 案中去掉。但是这样会多去掉一些数，比如对于不同的素数p_1,p_2,p_1^2\cdot p_2^2的 \\ 倍数就被去掉了两次，个数是\lfloor \frac{n}{p_1^2\cdot p_2^2}\rfloor ,需要加回来,这里就是使用了容斥原理. \end{gathered}$$

$$\begin{gathered} 如果d是s个不同素数的乘积,那么对答案的贡献就是(-1{)}^s\lfloor \frac{n}{d^2}\rfloor ,如果d含有平方因 \\ 子，d对答案没有贡献,也就是说容斥的系数就是Mobius函数 \end{gathered}$$

$因此可以有下面的结论，答案就是$
$$\sum _{k=1}^n{\mu }^2(k)=\sum _{d=1}^{\sqrt{n}}\mu (d)\lfloor \frac{n}{d^2}\rfloor$$

$枚举d，预处理莫比乌斯函数即可$


$代码:$

#include <bits/stdc++.h>
using namespace std;


#define Debug(x) cout<<#x<<"="<<x<<endl;
typedef long long LL;

const int N = 100100;
int n, mu[N+50], p[N+50];
int tot;
bool used[N+50];

//预处理出素数 莫比乌斯函数
void init() {
	mu[1] = 1;
	for (int i = 2; i <= N; i++) {
		if (!used[i]) {
			mu[i] = -1;
			p[++tot] = i;
		}
		for (int j = 1; j <= tot && p[j] * i <= N; j++) {
			used[p[j] * i] = true;
			if (i % p[j] == 0) {
				mu[i * p[j]] = 0;
				break;
			}
			mu[i * p[j]] = mu[i] * (-1);
		}
	}
	
}


LL check(LL n) {
	LL res = 0;
	for (int i = 1; i * i <= n; i++) {
		res += mu[i] * n / (i * i);
	}
	return res;
}


LL cal(LL k) {
	LL l = 1, r = k * 2;
	LL ans;
	while (l <= r) {
		int mid = (l + r) / 2;
		int number = check(mid);
		if (number >= k) r = mid - 1, ans = mid;
		else l = mid + 1;
	}
	return ans;
}


int main() {
	int T;
	cin >> T;
	init();
	while (T--) {
		int x;
		cin >> x;
		cout << cal(x) << endl;
	}
	return 0;
}