【数论】poj1777 Vivian's Problem（梅森素数）

题意：给定k个数，然后为每个数添加一个幂ei（0=<ei<=10），最后k项累乘的结果为M，如果M的所有因子的和可以写成2^x,求x的最大值，如果没有条件满足，输出NO。

刚看到题没什么思路，解题需要一些结论。

1：梅森素数

所谓梅森数，是指形如2p－1的一类数，其中指数p是素数，常记为Mp 。如果梅森数是素数，就称为梅森素数。

（下文中将梅森素数的指数P称为梅森指数）

（推荐看一下百度百科，多了解一些）

2：

关于梅森素数，有一个重要的定理：“一个数能够写成几个不重复的梅森素数的乘积” 等价于 “这个数的约数和是2的幂次”，注意不能重复，

比如说3是梅森素数，9就不满足约数和为2的幂，

并且这个幂次正好等于作为因子的梅森素数的梅森指数的和。

比如 3(2的2次幂-1)X7(2的3次幂-1) =21;

21的因数和1+3+7+21=32=2^5;

解题思路：

知道上面两个结论，这题就很明确了。只要找到题中所给的梅森素数，或梅森素数的累乘，则一定满足题干要求。

不是梅森数则幂次为0，梅森数不能重复作为因子出现，幂次为1。如果没有则输出NO。

代码如下：

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <cmath>typedef long long ll;using namespace std;#define MAXN 105int mp[8] = {2,3,5,7,13,17,19,31} ;int mason[10];int ans, cnt, status[MAXN];int k, a;void get(int a) {    //判断是否梅森素数，同时将状态存入status数组    status[cnt] = 0;    for (int i = 0; i < 8 && a >= mason[i]; i++)        if (a % mason[i] == 0) {            status[cnt] |= 1<<i;  //状压，位数代表对应的梅森指数            a /= mason[i];            if (a % mason[i] == 0)                return;        }    if (a == 1)        cnt++;}int  calc(int st) {    int res = 0;    for (int i = 0; i < 8; i++)        if (st & (1<<i))  //因数和S 与 梅森素数对应的指数之和x 关系： S=2^x；求出X即可。            res += mp[i];    return res;}void dfs(int cur, int st) {    ans = max(ans, calc(st));    for (int i = cur; i < cnt; i++)        if ((st & status[i]) == 0)  //如果作为因数的梅森素数没有重复，则相乘            dfs(i+1, st | status[i]);}int main() {    for (int i = 0; i < 8; i++)        mason[i] = (1 << mp[i]) - 1;    while (scanf("%d", &k) != EOF) {        cnt = 0;        for (int i = 0; i < k; i++) {            scanf("%d", &a);            get(a);        }        if (cnt == 0) {            printf("NO\n");        } else {            ans = 0;            dfs(0, 0);            printf("%d\n", ans);        }    }    return 0;}