bzoj2052: Pku1777 Vivian

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2052

2052: Pku1777 Vivian

Time Limit: 10 Sec  Memory Limit: 256 MB
Submit: 28  Solved: 14
[Submit][Status][Discuss]

Description

Input

 

Output

 

Sample Input

 

Sample Output

 

HINT

 

Source

梅森素数。。。

没写这道题之前根本不知道梅森素数。。。好像如果f(n)=sigma(d|n,d)=2^k，n一定是不重复的梅森素数之积，且k=梅森指数之和。。。

具体看这篇吧：http://blog.csdn.net/acm_cxlove/article/details/7860735

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int k,cnt,pp[8]={2,3,5,7,13,17,19,31},sum[256],ans;
long long p[8]={3,7,31,127,8191,131071,524287,2147483647};
bool v[256];
int calc(int x){
    int y=0;
    for(int i=0;i<8;i++)
    if(x%p[i]==0){
        x/=p[i];y=y+(1<<i);
    }
    if(x!=1)return -1;
    return y;
}
int main(){
    //freopen("data.in","r",stdin);
    //freopen("data.out","w",stdout);
    for(int i=0;i<256;i++)
    for(int j=0;j<8;j++)
    if(i&(1<<j))sum[i]+=pp[j];
    while(~scanf("%d",&k)){
        memset(v,0,sizeof(v));v[0]=1;
        for(int i=1;i<=k;i++){
            int x;scanf("%d",&x);
            x=calc(x);
            if(x==-1)continue;
            for(int j=255;j;j--)
            if((j&x)==x&&v[j^x])v[j]=1;
        }
        ans=0;
        for(int i=0;i<256;i++)
        if(v[i]&&sum[i]>ans)ans=sum[i];
        if(ans)printf("%d
",ans);
        else printf("NO
");
    }
    return 0;
}