完数

Write a program to output the first 555 perfect numbers. A perfect number is defined to be a positive integer where the sum of its positive integer divisors excluding the number itself equals the number.

For example: 1+2+3=61+ 2 + 3 = 61+2+3=6, and 666 is the first perfect number.

There is no input for this problem.

求前五个完数；

其实前4个都很好求，但是第五个很容易超时；

欧拉曾发现如此一个性质：若p是素数，且pow(2,p)-1也是素数，那么(pow(2,p)-1)*pow(2,p-1)就是一个完数。如果pow(2,p)-1是素数，称它为梅森素数，而梅森素数至今只发现了50个，所以完数的个数也不会大于50。

解题思路已经很清晰了，先做一次素数筛，再把满足p是素数，且pow(2,p)-1也是素数的所有p筛出来，但是这个素数筛并不需要筛到1e12。我们让pow(2,p)-1最大到1e7即可，那么这个式子(pow(2,p)-1)*pow(2,p-1)一定就不会比1e12小。

pow(2,p)-1最大到1e7，那么p不会比24大，所以先做一次1e7的素数筛，然后最多再24次就能把所有在1e12以内的完数找出来，每次询问的时候在表内查询即可。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e7+5;
int  vis[maxn+1];//0就代表i是素数或者1
void sieve(int n)
{
    int m=(int)sqrt(n+0.5);
    memset(vis,0,sizeof(vis));
    for(int i=2; i<=m; i++) if(!vis[i])
            for(int j=i*i; j<=n; j+=i) vis[j]=1; 
}
vector<ll> ans;
int main()
{
    sieve(maxn);
    ans.clear();
    for(int i=2; i<=24; i++)
    {
        int tem=(1<<i)-1;
        if(tem>maxn)break;
        if(!vis[i]&&!vis[tem])
        {
            ll tt=(ll)tem*(ll)(1<<(i-1));
            ans.push_back(tt);
        }
    }
    ll n;
    int flag=0;
    for(n=2;n<2147483647;n++){
        bool ok=false;
        for(int i=0; i<ans.size(); i++)
            if(ans[i]==n)
            {
                ok=true;
                break;
            }		

        if(ok){
        	cout<<n<<endl;
        	++flag;
		}
		if(flag==5)
			break;
    }
    return 0;
}