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

最新推荐文章于 2020-09-19 00:25:36 发布

原创最新推荐文章于 2020-09-19 00:25:36 发布 · 1.3k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#省选 #bzoj

OI 专栏收录该内容

324 篇文章

订阅专栏

本文介绍了一种基于Mu函数的高效算法实现，该算法通过预处理得到Mu函数值，并利用这些值来解决特定数学问题。文章提供了完整的C++代码实现，可用于计算给定数值下Mu函数对应的解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门：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;
}