BZOJ 2440

最新推荐文章于 2019-04-27 17:59:28 发布

NeighThorn

最新推荐文章于 2019-04-27 17:59:28 发布

阅读量375

点赞数 1

分类专栏： BZOJ 模板文章标签： vijos 莫比乌斯反演

本文链接：https://blog.csdn.net/neighthorn/article/details/51868079

版权

模板同时被 2 个专栏收录

54 篇文章 0 订阅

订阅专栏

BZOJ

45 篇文章 1 订阅

订阅专栏

传送门

求第k个不含完全平方数的数

分析：
同vijos1889
但是需要注意的是：Notice！！！
vijos求得是第k个含有完全平方数的数，那么就要从2开始循环

int check(int x){
    int sum=0,lala=sqrt(x);
    for(int i=2;i<=lala;i++)
        sum-=mu[i]*x/(i*i);
    return sum; 
}

很明显如果1也算完全平方数的话说，所有数都含有了
此题需要从1开始循环

int check(int x){
    int sum=0,m=sqrt(x);
    for(int i=1;i<=m;i++)
        sum+=mu[i]*x/(i*i);
    return sum;
}

代码如下：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#define int long long
using namespace std;
const int maxn=100000+5;
int cas,k,m,prime[maxn],vis[maxn],mu[maxn],cnt;
inline int read(){
    char ch=getchar();
    int f=1,x=0;
    while(!(ch>='0'&&ch<='9')){
        if(ch=='-')
            f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
        x=x*10+ch-'0',ch=getchar();
    return f*x;
}
int check(int x){
    int sum=0,m=sqrt(x);
    for(int i=1;i<=m;i++)
        sum+=mu[i]*x/(i*i);
    return sum;
}
signed main(void){
    cas=read();
    cnt=0,memset(vis,0,sizeof(vis)),mu[1]=1;
    for(int i=2;i<=50000;i++){
        if(!vis[i])
            prime[++cnt]=i,mu[i]=-1;
        for(int j=1;j<=cnt&&i*prime[j]<=50000;j++){
            vis[i*prime[j]]=1;
            if(i%prime[j]==0){
                mu[i*prime[j]]=0;
                break;
            }
            else
                mu[i*prime[j]]=-mu[i];
        }
    }
    while(cas--){
        k=read();
        int l=k,r=1644934081LL,ans;
        while(l<=r){
            int mid=(l+r)>>1;
            if(check(mid)>=k)
                ans=mid,r=mid-1;
            else
                l=mid+1;
        }
        cout<<ans<<endl;
    }
    return 0;
}