hdu5778 枚举

题目链接：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;
}