POJ 1777 Vivian's Problem(梅森素数)

本文探讨了使用梅森素数解决特定数学问题的算法思路,通过状态压缩动态规划方法,实现对给定数值进行有效处理,以达到题目要求的最优解。文章详细介绍了梅森素数的概念,以及如何利用这些素数的特性进行问题求解。

题目链接:http://poj.org/problem?id=1777

题意:给出n个数Ai,给每个数字加上一个指数ei(0<=e1<=10),使得:

思路:(1)梅森素数:2^i-1的素数叫做梅森素数;

(2)“一个数能够写成几个不重复的梅森素数的乘积” 等价于 “这个数的约数和是2的幂次”。

由于题目给的数据范围内的梅森素数只有8个,可以状态压缩。

 

int mason[8]={3,7,31,127,8191,131071,524287,2147483647};
int mp[8]={2,3,5,7,13,17,19,31};
vector<int> V;
int n,f[1<<10];

int ok(int t)
{
    int ans=0,i;
    FOR0(i,8) if(t%mason[i]==0)
    {
        ans|=(1<<i);
        t/=mason[i];
        if(t%mason[i]==0) return 0;
    }
    if(t!=1) return 0;
    return ans;
}


int cal(int t)
{
    int cnt=0,i;
    FOR0(i,8) if(t&(1<<i)) cnt+=mp[i];
    return cnt;
}

int main()
{
    Rush(n)
    {
        V.clear();
        int i,x;
        FOR0(i,n)
        {
            RD(x);
            if(x=ok(x)) V.pb(x);
        }
        n=SZ(V);
        if(!n)
        {
            puts("NO");
            continue;
        }
        clr(f,0);
        f[0]=1;
        int j;
        FOR0(i,(1<<8)) if(f[i]) FOR0(j,n) if(!(i&V[j]))
        {
            f[i|V[j]]=1;
        }
        int ans=0;
        FOR0(i,(1<<8)) if(f[i]) upMax(ans,cal(i));
        PR(ans);
    }
}

  

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值