题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5778
过了很久才来写题解,以前不会做,看了题解,有一个地方不是很懂,一直不是很明白为什么只需要找出10^5内的素数就可以了,现在有遇到这到题,突然想通了,来说一下我的理解吧。
题目要求给定x 要找到y-x的绝对值最小,y满足以下条件,y的质因数分解中所有的质因数都恰好只出现了两次,那么很容易就知道y是一个平方数,就可以转换为求sqrt(y)-sqrt(x)的最小值,那么x就缩小到了10^9级别,然后枚举x左右的满足条件的y 左边和右边第一个满足条件的y的选最小的一个就可以了。
我的问题就出在如何判断y满足条件,因为y已经开了根号,那么其实对于其质因数分解,每个质因数只出现一次,之前一直看不懂为什么是10^5以内的质数,因为我们判断的时候,如果当前数能被这个素数整除,那么把它除以这个素数,如果还能被这个数整除,那么肯定不满足条件,之所以打表10^5是因为如果y能被10^5以后的数整除一次,如果还能整除,那么y就已经大于10^10了否则是不可能被10^5以后的素数整除两次的。ok 基本就这么多 看代码。
#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<string>
#include<stack>
#include<queue>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
#define LL long long
LL ABS(LL x)
{
return x<0?-x:x;
}
const int N = 100005;
int prime[N];
int cnt = 0;
void Init()
{
bool isNotPrime[N] = {0};
cnt = 0;
for(int i = 2;i<N;i++)
{
if(!isNotPrime[i])
{
prime[cnt++] = i;
//cout << i << endl;
}
for(int j = 0;j<cnt && i*prime[j]<N;j++)
{
isNotPrime[i*prime[j]] = true;
if(!(i%prime[j]))
break;
}
}
}
bool isAns(LL x)
{
for(int i = 0;i<cnt;i++)
{
if(x%prime[i]==0)
x = x/prime[i];
if(x%prime[i]==0)
return false;
if(x==1)
return true;
}
return true;
}
LL solve(LL x,int ch)
{
LL temp = sqrt(x);
if(ch>0)
temp++;
//printf("%d\n",temp);
while(!isAns(temp))
{
temp += ch;
}
return ABS(temp*temp-x);
}
int main()
{
int t;
scanf("%d",&t);
Init();
while(t--)
{
LL x;
scanf("%lld",&x);
if(x<4)
{
printf("%lld\n",ABS(4-x));
continue;
}
LL ans = solve(x,1);
LL tans = solve(x,-1);
//cout << ans << " " << tans << endl;
printf("%lld\n",min(ans,tans));
}
return 0;
}