紫书习题 10-8 UVa 10622（gcd）

最新推荐文章于 2020-01-27 10:07:07 发布

weixin_34336292

最新推荐文章于 2020-01-27 10:07:07 发布

阅读量83

点赞数

原文链接：http://www.cnblogs.com/sugewud/p/9819483.html

版权

本文介绍了一种利用质因数分解求解整数因子个数的最大公约数(GCD)的方法。通过预处理得到一定范围内的所有质数，并以此为基础进行质因数分解，最终求得每个质因数的次数的GCD。文章提供了详细的C++实现代码。

摘要由CSDN通过智能技术生成

把这个数质因数分解然后求因子个数的gcd就ok了。

一些细节

（1）这道题的质因数不需要存下来，每一次做完取一次gcd就ok了

（2）判断奇偶用ans & 1的时候要加括号，位运算要注意括号

（3）这道题在32位带符号整数范围内，也就是说用int可以了

（3）注意ans最后的处理，看代码

#include<cstdio>
#include<vector>
#include<cmath>
#include<cstring>
#define REP(i, a, b) for(int i = (a); i < (b); i++)
using namespace std;

const int MAXN = 1123456;
bool is_prime[MAXN];
vector<int> prime;

void init()
{
	memset(is_prime, true, sizeof(is_prime));
	is_prime[0] = is_prime[1] = false;
	REP(i, 2, MAXN)
	{
		if(is_prime[i]) prime.push_back(i);
		REP(j, 0, prime.size())
		{
			if(i * prime[j] >= MAXN) break;
			is_prime[i * prime[j]] = false;
			if(i % prime[j] == 0) break;
		}
	}
}

int gcd(int a, int b) { return !b ? a : gcd(b, a % b); }

int solve(int n)
{
	int t = n;
	if(n < 0) n = -n;
	int ans = 0;
	REP(i, 0, prime.size())
		if(n % prime[i] == 0)
		{
			int cnt = 0;
			while(n % prime[i] == 0)
			{
				n /= prime[i];
				cnt++;
			}
			ans = gcd(ans, cnt);
			if(n == 1) break;
		}
	
	if (ans == 0) ans = 1; //可能是一个非常大的质数，所以要加这个判断 
	if (t < 0) 
		while ((ans & 1) == 0) //加括号 
			ans >>= 1;        //是负数的话答案必须是奇数，因为质因数里面只有2的偶数 
	return ans;               //所以把2除完剩下一定是奇数 
}

int main()
{
	init();
	int n;
	while(~scanf("%d", &n) && n)
		printf("%d\n", solve(n));
	return 0;
}

转载于:https://www.cnblogs.com/sugewud/p/9819483.html