码迷,mamicode.com
首页 > 其他好文 > 详细

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

时间:2016-07-31 14:22:08      阅读:151      评论:0      收藏:0      [点我收藏+]

标签:

分析: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+5;
bool v[N];
int prime[N>>1],cnt;
void getprime(){
  for(int i=2;i*i<=N-5;++i)
    if(!v[i])
      for(int j=i*i;j<=N-5;j+=i)
        v[j]=1;
  for(int i=2;i<=N-5;++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=-1,r=-1;
    for(int k=x;k>=2;--k){
     bool flag=false;
     for(int i=1;i<=cnt&&1ll*prime[i]*prime[i]<=k;++i){
      if(k%(1ll*prime[i]*prime[i])==0){flag=true;break;}
     }
     if(!flag){l=k;break;}
    }
    for(int k=x+1;;++k){
       bool flag=false;
     for(int i=1;i<=cnt&&1ll*prime[i]*prime[i]<=k;++i){
      if(k%(1ll*prime[i]*prime[i])==0){flag=true;break;}
      
    }if(!flag){r=k;break;}
   }
   LL ret;
   if(l==-1)ret=fun(r*r,n);
   else ret=min(fun(n,l*l),fun(r*r,n));
    printf("%I64d\n",ret);
  }
  return 0;
}
View Code

 

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

标签:

原文地址:http://www.cnblogs.com/shuguangzw/p/5722916.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!