uvaoj 10791 Minimum Sum LCM 类似欧拉定理的写法

uvaoj 10791 Minimum Sum LCM 类似欧拉定理的写法

给定一个数n，求一个集合，这个集合的最小公倍数是n，并且这个集合中数的和最小。

处理好集中情况：

1）n为1时，应该输出2。

2）n为素数时，应该输出n+1。

3）n为某个素数的次方，即只有一个质因子，应该输出n+1。

4）n为2^31-1时，这个数为质数，应该输出2^31，所以要考虑用long long。

5）其他情况下，答案为质因子的相应次方加起来。

代码如下：

/*************************************************************************
	> File Name: 10791.cpp
	> Author: gwq
	> Mail: gwq5210@qq.com 
	> Created Time: 2015年01月06日 星期二 20时48分01秒
 ************************************************************************/

#include <cmath>
#include <ctime>
#include <cctype>
#include <climits>
#include <cstdio>
#include <cstdlib>
#include <cstring>

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <string>
#include <vector>
#include <sstream>
#include <iostream>
#include <algorithm>

#define INF (INT_MAX / 10)
#define clr(arr, val) memset(arr, val, sizeof(arr))
#define pb push_back
#define sz(a) ((int)(a).size())

using namespace std;
typedef set<int> si;
typedef vector<int> vi;
typedef map<int, int> mii;
typedef long long ll;

const double esp = 1e-5;

int main(int argc, char *argv[])
{
	ll n;
	ll c = 0;
	while (scanf("%lld", &n) != EOF) {
		if (n == 0) {
			break;
		}
		//ll m = n;
		ll m = (int)sqrt(n + 0.5);
		ll a = n;
		ll res = 0;
		ll flag = 0;
		// 这里i使用int爆了
		//for (ll i = 2; i * i <= m; ++i) {
		for (int i = 2; i <= m; ++i) {
			if (n % i == 0) {
				flag++;
				ll tmp = 1;
				while (n % i == 0) {
					n /= i;
					tmp *= i;
				}
				res += tmp;
			}
		}
		if (n == a) {
			res = n + 1;
		} else if (flag == 1 || n > 1) {
			res += n;
		}
		printf("Case %lld: %lld\n", ++c, res);
	}
	return 0;
}

参考：http://blog.csdn.net/goomaple/article/details/8550381