【莫比乌斯函数】BZOJ2440(中山市选2011)[完全平方数]题解

题目概述

求第 $k$ 个不是完全平方数（此题中 $1$ 不是完全平方数）倍数的数。

解题报告

显然是求 $\sum_{i=1}^{n}\mu^2(i)\ge k$ 的最小的 $n$ ，好像有个性质：

$$\mu^2(n)=\sum_{d^2|n}\mu(d)$$
然而我并不会证明QAQ，就当做是结论记一下吧。那么：

$$\sum_{i=1}^{n}\mu^2(i)=\sum_{i=1}^{n}\sum_{d^2|i}\mu(d)=\sum_{d=1}^{\lfloor\sqrt x\rfloor}\mu(d)\lfloor{x\over d^2}\rfloor$$
所以先用线性筛筛出莫比乌斯函数，二分之后直接暴枚验证就行了。

示例程序

#include<cstdio>
#include<cmath>
using namespace std;
typedef long long LL;
const int maxn=100000;

int te,n,mu[maxn+5],p[maxn+5];bool pri[maxn+5];

void Make()
{
    pri[1]=true;mu[1]=1;
    for (int i=2;i<=maxn;i++)
    {
        if (!pri[i]) p[++p[0]]=i,mu[i]=-1;
        for (int j=1,t;j<=p[0]&&(t=i*p[j])<=maxn;j++)
        {
            pri[t]=true;mu[t]=-mu[i];
            if (i%p[j]==0) {mu[t]=0;break;}
        }
    }
}
inline bool check(int mid)
{
    int S=sqrt(mid);LL ans=0;
    for (int i=1;i<=S;i++) ans+=mu[i]*mid/(i*i);
    return ans>=n;
}
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    for (Make(),scanf("%d",&te);te;te--)
    {
        scanf("%d",&n);int L=1,R=2e9;
        for (int mid=L+(R-L>>1);L<=R;mid=L+(R-L>>1))
            if (check(mid)) R=mid-1; else L=mid+1;
        printf("%d\n",L);
    }
    return 0;
}