[莫比乌斯函数] BZOJ 2986 Non-Squarefree Numbers & BZOJ 2440 [中山市选2011]完全平方数

不知不觉就发现了2440的双倍经验

考虑二分 然后用容斥求出数量 

x以内的无平方因子数=无需是任何质数的倍数的数的数量(即x)-是至少一个质数平方倍数的数的数量+是至少两个质数平方倍数的数的数量-是至少三个质数平方倍数的数的数量...

发现每个数前的符号与莫比乌斯函数一致

然后么 复杂度分析 Orz

2986的代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;

inline char nc(){
	static char buf[100000],*p1=buf,*p2=buf;
	if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
	return *p1++;
}

inline void read(ll &x){
	char c=nc(),b=1;
	for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
	for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int maxn=500000;
int mobius[maxn+5],sum[maxn+5];
int vst[maxn+5],prime[maxn+5],num=0;

inline void Pre()
{
	mobius[1]=1;
	for (int i=2;i<=maxn;i++)
	{
		if (!vst[i])
			mobius[i]=-1,prime[++num]=i;
		for (int j=1;j<=num && prime[j]*i<=maxn;j++)
		{
			vst[i*prime[j]]=1;
			if (i%prime[j]==0)
			{
				mobius[i*prime[j]]=0; break;
			}
			else
				mobius[i*prime[j]]=-mobius[i];
		}
	}
}

ll K;

inline ll Cal(ll x)
{
	ll ret=0;
	ll t=sqrt(x);
	for (ll i=1;i<=t;i++)
		ret+=x/(i*i)*mobius[i];
	return x-ret;
}

int main()
{
	ll Q;
	ll L,R,MID;
	freopen("t.in","r",stdin);
	freopen("t.out","w",stdout);
	Pre();
	Q=1;
	while (Q--)
	{
		read(K);
		L=0; R=210000000000LL;
		while (L+1<R)
			if (Cal(MID=(L+R)>>1)>=K)
				R=MID;
			else
				L=MID;
		printf("%lld\n",R);
	}
	return 0;
}