UVa 11466 - Largest Prime Divisor

最新推荐文章于 2016-07-28 08:38:00 发布

原创最新推荐文章于 2016-07-28 08:38:00 发布 · 1.1k 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

解题报告同时被 2 个专栏收录

1435 篇文章

订阅专栏

数论

191 篇文章

订阅专栏

本文介绍了一种利用素数表和短除法快速找出整数n的最大素数因子的方法，适用于负数和正数，代码实现简洁高效。

題目：給你一個整數n（不超過14位），求出他的最大的素數因子，如果只有一個素數因子輸出-1。

分析：數論。直接打表計算10^7內的所有素數因子，然後用短除法除n，記錄最大的因子即可。

如果最後下的數字不是1，則它就是最大的素數因子。

說明：注意n可能為負數。

#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>

using namespace std;

int visit[10000001];
int prime[10000001];

int main()
{
    memset(visit, 0, sizeof(visit));  
    int count = 0;  
    for (int i = 2; i < 10000001; ++ i) {  
        if (!visit[i]) prime[count ++] = i;  
        for (int j = 0; j < count && i*prime[j] < 10000001; ++ j) 
            visit[i*prime[j]] = 1;  
    }
	
	long long n,m;
	while (cin >> n && n) {
		if (n < 0) n = -n;
		int size = 0;
		for (int i = 0; i < count; ++ i) {
			if (n < prime[i]) break;
			if (n%prime[i] == 0) {
				m = prime[i];
				++ size;
				while (n%m == 0) n /= m;
			}
		}
		if (n != 1) {
			m = n;
			++ size;
		}
		
		if (size > 1) 
			cout << m << endl;
		else cout << "-1" << endl;
	}
    return 0;
}