求第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;
}
by >o< neighthorn