【BZOJ】【P2440】【中山市选2011】【完全平方数】【题解】【莫比乌斯反演】

最新推荐文章于 2018-10-18 00:00:00 发布

iamzky

最新推荐文章于 2018-10-18 00:00:00 发布

阅读量1.2k

点赞数

分类专栏： OI 文章标签：省选 bzoj

本文链接：https://blog.csdn.net/iamzky/article/details/27706483

版权

OI 专栏收录该内容

324 篇文章 0 订阅

订阅专栏

传送门：www.lydsy.com/JudgeOnline/problem.php?id=2440

还不是很会，先放着代码

Code：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=1e5+10;
typedef long long lld;
lld u[maxn],prime[maxn];
lld T;
bool ok[maxn];
void getmu(){
	u[1]=1;
	for(int i=2;i<maxn;i++){
		if(!ok[i]){
			prime[++prime[0]]=i;
			u[i]=-1;
		}
		for(int j=1;j<=prime[0];j++){
			if(i*prime[j]>=maxn)break;
			ok[i*prime[j]]=1;
			if(i%prime[j])
				u[i*prime[j]]=-u[i];
			else{
				u[i*prime[j]]=0;
				break;
			}
		}
	}
}
lld calc(lld x){
	lld sum=0,sq=sqrt(x);
	for(int i=1;i<=sq;i++)
		sum+=u[i]*(x/(i*i));
	return sum;
}
lld work(lld x){
	lld l=0,r=x*2;
	while(l<r){
		lld mid=(l+r)>>1;
		if(calc(mid)>=x)
			r=mid;
		else
			l=mid+1;
	}return l;
}
int main(){
	scanf("%lld",&T);
	getmu();
	while(T--){
		lld x;scanf("%lld",&x);
		printf("%lld\n",work(x));
	}
	return 0;
}