【BZOJ 2440】[中山市选2011]完全平方数 莫比乌斯反演+容斥

考虑二分答案，所以需要的问题就是1->n以内有多少个数不是完全平方数的倍数即u[i]!=0

ans=n-是完全平方数的倍数的个数

即

ans=n-只有一个质因数的个数+有两个质因数的个数-有三个质因数的个数。。。。。。

考虑枚举完全平方数的根，然后x/(i*i)就得到在x内它的倍数有多少个，如果这个数已经是完全平方数的倍数u[i]==0对答案并没有影响，否则恰好是莫比乌斯函数

#include<cstdio>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
LL T,k,pri[100021],vis[100021],u[100021],cnt;
void pre(){
	u[1]=1;
	for(LL i=2;i<=100000;i++){
		if(!vis[i])pri[++cnt]=i,u[i]=-1;
		for(LL j=1;j<=cnt&&pri[j]*i<=100000;j++){
			vis[pri[j]*i]=1;
			u[pri[j]*i]=-u[i];
			if(i%pri[j]==0){u[pri[j]*i]=0;break;}
		}
	}
}
bool cheak(LL x){
	LL ans=0;
	for(LL i=1;i*i<=x;i++){
		ans+=u[i]*(x/(i*i));
	}return ans>=k;
}
int main(){
	pre();
	scanf("%lld",&T);
	while(T--){
		scanf("%lld",&k);
		LL l=1,r=1ll<<31;
		while(l<r){
			LL mid=l+r>>1ll;
			if(cheak(mid))r=mid;
			else l=mid+1;
		}
		printf("%lld\n",l);
	}
	return 0;
}