洛谷 P4318 完全平方数【莫比乌斯】【数论分块】【二分】

---

题目：

传送门

---

题意：

定义$n$为无平方因子数，当且仅当不存在$d>1$满足$d^2|n$
求第$k$个无平方因子数

---

分析：

因为范围可能会很大，所以我们不能一个个去枚举，我们采用二分的方法找到无平方因子数恰好等于$k$的数
现在问题来到怎么求$n$的无平方因子数
这个数值显然可以表示成$n-p_i^2|n+(p_i\ast p_j{)}^2|n\dots \dots$，这个就是经典的容斥模型啦
再看到容斥的系数，我们发现恰好是等于$\mu (n)$的
所以我们可以用线性筛预处理出所有$\mu (n)$
在本菜的代码中我在计算时还用了数论分块，不用其实影响也不大，大家就看着切就好了

---

代码：

#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<cmath>
#define LL long long 
#define LZX Mu
using namespace std;
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
LL m[40050],tf[40050],p[40050],cnt=0,n=read();
bool check(LL x,LL k)
{
	LL ans=0,p;
	for(LL i=1;i*i<=x;i=p+1)
	{
		p=min((LL)sqrt(x/(x/(i*i))),(LL)sqrt(x));
		ans+=x/(i*i)*(m[p]-m[i-1]);
	}
	return ans>=k;
}
void ef(LL x)
{
	LL l=x,r=x<<1;
	while(l<=r)
	{
		LL mid=(l+r)>>1;
		if(check(mid,x)) r=mid-1;
		else l=mid+1;
	}
	if(check(l,x)) printf("%lld\n",l);
	else printf("%lld\n",r);
	return;
}
int main()
{
	m[1]=1;
	for(LL i=2;i<=40005;i++)
	{
		if(!tf[i]) p[++cnt]=i,m[i]=-1;
		for(LL j=1;j<=cnt&&p[j]*i<=40005;j++)
		{
			tf[p[j]*i]=1;
			if(!(i%p[j])) break;
			m[i*p[j]]=-m[i];
		}
 	}
 	for(LL i=1;i<=40005;i++) 
	   m[i]+=m[i-1];
	 while(n--)
	 {
	 	LL x=read();
	 	ef(x);
	 }
	return 0;
}