HDU 5778 abs（数论）

Description
给定一个数x，求正整数y≥2，使得满足以下条件：
1.y-x的绝对值最小
2.y的质因数分解式中每个质因数均恰好出现2次。
Input
第一行输入一个整数T（1≤T≤50)
每组数据有一行，一个整数x（1≤x≤10^​18​​ )
Output
对于每组数据，输出一行y-x的最小绝对值
Sample Input
5
1112
4290
8716
9957
9095
Sample Output
23
65
67
244
70
Solution
由于10^9内相邻素数的差值不会超过300，y作为某个素数的平方是满足第二个条件的，而必然存在相邻两素数p1,p2，使得p1*p1<=x<=p2*p2，故直接枚举sqrt(x)附近的数判断其是否无平方因子即可
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
#define maxn 111111
int prime[maxn],res,is_prime[maxn];
void get_prime(int n)
{
    memset(is_prime,0,sizeof(is_prime));
    res=0;
    for(int i=2;i<n;i++)
        if(!is_prime[i])
        {
            prime[res++]=i;
            for(int j=2*i;j<n;j+=i)is_prime[j]=1;
        }
}
int T;
ll x;
int check(ll n)
{
    for(int i=0;i<res&&prime[i]*prime[i]<=n;i++)
    {
        if(n%prime[i]==0)
        {
            n/=prime[i];
            if(n%prime[i]==0)return 0;
        }
    }
    return 1;
}
int main()
{
    get_prime(50000);
    scanf("%d",&T);
    while(T--)
    {
        scanf("%I64d",&x);
        if(x<=4)printf("%d\n",4-x);
        else
        {
            ll ans=x,y=(ll)sqrt(1.0*x);
            for(int i=0;;i++)
            {
                ll temp=y+i;
                if(temp*temp>=x&&check(temp))
                {
                    ans=min(ans,temp*temp-x);
                    break;
                }
            }
            for(int i=0;;i++)
            {
                ll temp=y-i;
                if(temp*temp<=x&&check(temp))
                {
                    ans=min(ans,x-temp*temp);
                    break;
                }
            }
            printf("%d\n",ans);
        }
    }
    return 0;
}