完全数

基本理论：

完全数，又称为完美数或完备数，是一些特殊的自然数，它所有的真因子(即除了自身以外的约数)的和，恰好等于它本身。完美数是非常少见的。迄今为止，也不过发现29个。其他的数就当然不是完美数。不完美数可以分为两类——盈数和亏数。如果一个数的所有真因子和大于该数本身，那么就称它为盈数；如果一个数所有真因子和小于该数本身，那么就称为亏数。

基本理论：

(1)如果n是一个正整数且 σ(n)=2n,那么n就称为完全数。(σ是希腊字母，英文表达sigma，汉语译音为“西格玛”。)

(2)正整数n是一个偶完全数，当且仅当n=(2^m-1)*2^(m-1),其中m>=2是使得2^m-1是一个素数的整数。

梅森素数中的一个定理：如果m是一个正整数且2^m-1是一个素数，那么m必是素数。

(3)完全数都能写成连续自然数的和。

(4)每个完全数的全部因数倒数之和都是2，因此每个完全数都是调和数。

(5)完全数都是以6或者8结尾。

下面是前十二个完全数：

1……6

2……28

3……496

4……8,128

5……33,550,336

6……8,589,869,056

7……137,438,691,328

8……2,305,843,008,139,952,128

9……2,658,455,991,569,831,744,654,692,615,953,842,176

10……191,561,942,608,236,107,294,793,378,084,303,638,130,997,321,548,169,216

11……13,164,036,458,569,648,337,239,753,460,458,722,910,223,472,318,386,943,117,783,728,128

12……14,474,011,154,664,524,427,946,373,126,085,988,481,573,677,491,474,835,889,066,354,349,131,199,152,128

例题一：完美数(sdut 1220)

Problem Description

任何一个自然数的约数中都有1和它本身，我们把小于它本身的因数叫做这个自然数的真约数。
如6的所有真约数是1、2、3，而且6=1＋2＋3。像这样，一个数所有真约数的和正好等于这个数，通常把这个数叫做完美数。
古希腊人非常重视完美数。毕达哥拉斯发现它之后，人们就开始了对完美数的研究。
现在要求输出所有在m和n范围内的完美数。

Input

输入数据有多组，每组占一行，包括两个整数m和n（1≤m≤n≤99999999）。
输入以0 0结束。

Output

对于每个测试实例，要求输出所有在给定范围内的完美数，就是说，输出的完美数必须大于等于m，并且小于等于n，如果有多个，则要求从小到大排列在一行内输出，之间用一个空格隔开；
如果给定的范围内不存在完美数，则输出No；
每个测试实例的输出占一行。

Example Input

1 100
0 0

Example Output

6 28

分析：

由于1<=m<=n<=99 999 999,所以数据比较小，打表到第五个完全数即可；

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int main()
{
    int p[8]={6,28,496,8128,33550336};
    int m,n;
    while(scanf("%d%d",&m,&n)!=EOF)
    {
        if(m==0&&n==0)
            break;
        bool flag=0;
        int i;
        int c;
        for(i=0;i<5;i++)
            if(p[i]>=m&&p[i]<=n)
        {
            c=i;
            cout<<p[i];
            flag=1;
            break;
        }
        if(flag==0)
        {
            cout<<"No"<<endl;
            continue;
        }
        for(i=c+1;i<5;i++)
            if(p[i]>=m&&p[i]<=n)
        {
            cout<<" "<<p[i];
        }
        cout<<endl;
    }
    return 0;
}

例题二：完美数(ecnu 1197)

题意：给定一个数，判断是盈数，亏数还是完全数。

输入：输入数据有多组。

           每组数据输入仅一行，一个整数N(1<=N<=1 000 000 000)

           当N=0时输入结束

输出：

对于每组测试数据，若N是完美数，输出0；若N是盈数，输出1；若N是亏数，输出-1；

输入样例：

4

6

36

0

输出样例：

-1

0

1

分析:

对于输入的数，计算其因子和，当N 比较大时可以不必计算全部，计算一部分的因子和，在循环中判断当前因子和大于2N时就可以直接判断是盈数，不必计算全部。

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int main()
{
    long long n,m,s;///s记录的是因数和
    while(scanf("%d",&n)!=EOF)
    {
        if(n==0)
            break;
            m=n,s=1;
        for(int i=2;i*i<=n;i++)
            if(n%i==0)
        {
            long long a=1;///a记录某个素因子的幂
            while(n%i==0)
            {
                a*=i;
                n/=i;
            }
            s*=(a*i-1)/(i-1);
            if(s>2*m)///如果因子和大于2N
                break;
        }
        if(s>2*m)
        {
            cout<<1<<endl;
            continue;
        }
        if(n>1)
            s*=(n+1);
        s-=m
        if(s>m)
            cout<<1<<endl;
        else if(s==m)
        cout<<0<<endl;
        else
            cout<<-1<<endl;
    }
    return 0;
}