【JZOJ 4937】与运算

最新推荐文章于 2020-03-29 20:50:18 发布

HOWARLI

最新推荐文章于 2020-03-29 20:50:18 发布

阅读量795

点赞数

分类专栏： DP 文章标签：位运算状态压缩

本文链接：https://blog.csdn.net/HOWARLI/article/details/54428454

版权

DP 专栏收录该内容

50 篇文章 0 订阅

订阅专栏

Description

这里写图片描述

Solution

肯定先把所有的数转成二进制，
设一个 $g_i$ 表示i这个二进制被多少个数包含，也就是这个数前面可以放多少个数（满足 $a_i$ and i=i），
这样就可以状态压缩DP了，
当前做到了i，枚举在i的基础上减掉一位转移到j，那么转移式为： $f_j=f_i+j*(g_j-g_i)$

复杂度： $O(m\log(m))$

Code

#include <iostream>
#include <cstdio>
#include <cstdlib>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long LL;
const int N=8148576+1;
int read(int &n)
{
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
int m,n;
int er[22];
LL g[N];
LL f[N],ans;
void divide(int l,int r)
{
    if(l==r)return;
    int t=(l+r)/2;
    divide(l,t);divide(t+1,r);
    fo(i,l,t)g[i]+=g[i+t-l+1];
}
int main()
{
    freopen("and.in","r",stdin);
    freopen("and.out","w",stdout);
    int mx=0,q;
    er[1]=1;fo(i,2,21)er[i]=er[i-1]<<1;
    read(n);
    fo(i,1,n)g[read(q)]++,mx=max(mx,q);
    for(m=1;er[m]<=mx;m++);
    divide(0,er[m]-1);
    fod(i,er[m]-1,2)
    {
        fo(j,1,m-1)if(!(er[j]&i))
        {
            f[i]=max(f[i],f[i+er[j]]+i*(g[i]-g[i+er[j]]));
        }
        ans=max(ans,f[i]);
    }
    printf("%lld\n",ans);
    return 0;
}