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
#include
#include
#include
#include
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;iif(!is_prime[i])
        {
            prime[res++]=i;
            for(int j=2*i;j1;
        }
}
int T;
ll x;
int check(ll n)
{
    for(int i=0;iif(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;
}

你可能感兴趣的:(HDU,数论)