HDOJ 5778 abs

最新推荐文章于 2020-04-20 22:53:03 发布

RaAlGhul

最新推荐文章于 2020-04-20 22:53:03 发布

阅读量467

点赞数 1

分类专栏： ACM_数论初步

本文链接：https://blog.csdn.net/RaAlGhul/article/details/52078154

版权

ACM_数论初步专栏收录该内容

26 篇文章 0 订阅

订阅专栏

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5778

题意：给我们一个数，判断离它最近的那个能够分拆成质因子集每一个质因子的个数都是2的形式的数。

比赛的时候觉得好难的一道题，比赛完了再看真是觉得懊恼……这么简单的一道题居然没有想出来。

首先我们要求的数y = p1*p1 * p2*p2 * p3*p3……pk*pk = (p1*p2*p3……*pk) * (p1*p2*p3……*pk) = z*z。这样我们就可以把求y转换成求z，接下来我们就可以对x开根号之后向两边进行枚举即可求得答案。

#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 100000+5;
bool vis[maxn];
typedef long long LL;
LL prime[maxn], num;
void getPrime()
{
	for(LL i=2; i<maxn; i++)
	{
		if(!vis[i]) prime[num++] = i;
		for(LL j=1; i*j<maxn; j++)
			vis[i*j] = 1;
	}
}
bool judge(LL x)
{
	if(x < 2) return false;
	for(LL i=0; i<num; i++)
	{
		LL tmp = prime[i]*prime[i];
		if(tmp > x) return true;
		if(x % tmp == 0) return false;
	}
	return true;
}
LL solve(LL x)
{
	LL L = sqrt(x);
	LL R = L+1;
	while(1)
	{
		if(abs(L*L - x) < abs(R*R - x))
			if(judge(L))
				return abs(L*L-x);
			else L--;
		else
			if(judge(R))
				return abs(R*R-x);
			else R++;
	}
}
int main()
{
	getPrime();
//	printf("%d\n", num);
//	for(int i=0; i<num; i++) printf("%d ", prime[i]);
	int T;
	scanf("%d", &T);
	while(T--)
	{
		LL x;
		scanf("%I64d", &x);
		printf("%I64d\n", solve(x));
	}
	return 0;
}