UVa 10490 - Mr. Azad and his Son!!!!!

标签：style   blog   io   os   sp   for   on   2014   log   

题目：给你一个数k判断2^(k-1)*(2^k-1)是不是完全数（真因数之和和自身相等），不是判断k是不是素数。

分析：数论。欧拉证明了所有偶完全数都满足式子2^(k-1)*(2^k-1)；其中2^k-1为素数时，上式为完全数。

            满足2^k-1形式的素数叫梅森素数，这里打表计算50000内的素数判断2^k-1是不是素数即可。

           （如果存在，超过50000的素数因子只能有一个，而且一定是素数）

说明：又见数论(⊙o⊙)。

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

using namespace std;

int visit[50000];
int prime[50000];
int perfect[32];

int main()
{
	memset(visit, 0, sizeof(visit));
	visit[0] = visit[1] = 1;
	int count = 0;
	for (int i = 2 ; i < 50000 ; ++ i)
		if (!visit[i]) {
			prime[count ++] = i;
			for (int j = i+i ; j < 50000 ; j += i)
				visit[j] = 1;
		}
	memset(perfect, 0, sizeof(perfect));
	for (int i = 2 ; i < 32 ; ++ i) {
		int V = (1<<i)-1,flag = 0;
		for (int j = 0 ; j < count ; ++ j)
			while (V%prime[j] == 0) {
				V /= prime[j];
				flag ++;
			}
		if (V != 1) flag ++;
		if (flag == 1) 
			perfect[i] = 1;
	}
	
	int n;
	while (scanf("%d",&n) && n)
		if (perfect[n])
			printf("Perfect: %lld!\n",(1LL<<(n-1LL))*((1LL<<n)-1LL));
		else if (!visit[n])
			printf("Given number is prime. But, NO perfect number is available.\n");
		else printf("Given number is NOT prime! NO perfect number is available.\n");
    return 0;
}

UVa 10490 - Mr. Azad and his Son!!!!!

标签：style   blog   io   os   sp   for   on   2014   log   

原文地址：http://blog.csdn.net/mobius_strip/article/details/41559815