HDU 5778 abs (BestCoder Round #85 C)素数筛+暴力

分析:y是一个无平方因子数的平方,所以可以从sqrt(x)向上向下枚举找到第一个无平方因子比较大小

大家可能觉得这样找过去暴力,但实际上无平方因子的分布式非常密集的,相关题目,可以参考

CDOJ:无平方因子数 http://acm.uestc.edu.cn/#/problem/show/618

这个题和CDOJ的题虽然不一样,但是可以从CDOJ发现这种数是很多的

官方题解:官方题解说这个无平方因子的枚举量在logn级别,可见非常小

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <map>
#include <queue>
#include <vector>
using namespace std;
typedef long long LL;
const int N = 1e5+;
bool v[N];
int prime[N>>],cnt;
void getprime(){
for(int i=;i*i<=N-;++i)
if(!v[i])
for(int j=i*i;j<=N-;j+=i)
v[j]=;
for(int i=;i<=N-;++i)
if(!v[i])prime[++cnt]=i;
}
LL fun(LL x,LL y){
if(x>y)swap(x,y);
return y-x;
}
int main(){
int T;
getprime();
scanf("%d",&T);
while(T--){
LL n;scanf("%I64d",&n);
LL x=(LL)(sqrt(1.0*n)+0.5),l=-,r=-;
for(int k=x;k>=;--k){
bool flag=false;
for(int i=;i<=cnt&&1ll*prime[i]*prime[i]<=k;++i){
if(k%(1ll*prime[i]*prime[i])==){flag=true;break;}
}
if(!flag){l=k;break;}
}
for(int k=x+;;++k){
bool flag=false;
for(int i=;i<=cnt&&1ll*prime[i]*prime[i]<=k;++i){
if(k%(1ll*prime[i]*prime[i])==){flag=true;break;} }if(!flag){r=k;break;}
}
LL ret;
if(l==-)ret=fun(r*r,n);
else ret=min(fun(n,l*l),fun(r*r,n));
printf("%I64d\n",ret);
}
return ;
}
上一篇:HDU 4859 海岸线(最大流最小割)


下一篇:hdu 5101 Select(Bestcoder Round #17)