和与或(dp)

最新推荐文章于 2024-05-30 16:17:47 发布

FAUX123455

最新推荐文章于 2024-05-30 16:17:47 发布

阅读量60

点赞数

分类专栏： DP 文章标签：算法动态规划

本文链接：https://blog.csdn.net/weixin_50935546/article/details/125886812

版权

DP 专栏收录该内容

6 篇文章 0 订阅

订阅专栏

题目链接.
假设 $s u m = a [1] + a [2] + ... + a [n] = a [1] ∣ a [2] ∣...∣ a [n]$ ;

若要满足上述条件， $s u m$ 二进制的每一位1只能来自于一个a[i],否则就会导致 $a [1] + a [2] + ... + a [n] > a [1] ∣ a [2] ∣...∣ a [n]$

$d p [i] [j]$ i指 $s u m$ 二进制的第 $i$ 位， $j$ 为 $n$ 个数的第 $i$ 位是否提供 $1$ 的二进制序列的状态压缩。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll mod=1e9+9;

int n;
ll r[12];
int c[12][62];
ll dp[62][1<<10+5];

ll dfs(int x,int pos)
{
    if(dp[x][pos])return dp[x][pos];
    if(x==0)return 1;
    int res=pos;
    for(int i=1;i<=n;i++)
        res=res|(c[i][x]<<(i-1));
    dp[x][pos]=dfs(x-1,res);
    for(int i=1;i<=n;i++)
    {
        int p=(pos>>(i-1))&1;
        int q=(p==1?1:c[i][x]);
        if(q)
            dp[x][pos]=(dp[x][pos]+dfs(x-1,res^((p==0)<<(i-1))))%mod;
    }
    return dp[x][pos]%mod;
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);cout.tie(0);

    cin>>n;
    for(int i=1;i<=n;i++)cin>>r[i];
    
    //预处理，c[i][j]记录r[i]的二进制第j位的值
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=60;j++)
        {
            c[i][j]=(r[i]>>(j-1))&1;
        }
    }
    cout<<dfs(60,0);

    return 0;
}