Description
Solution
肯定先把所有的数转成二进制,
设一个
gi
表示i这个二进制被多少个数包含,也就是这个数前面可以放多少个数(满足
ai
and i=i),
这样就可以状态压缩DP了,
当前做到了i,枚举在i的基础上减掉一位转移到j,那么转移式为:
fj=fi+j∗(gj−gi)
复杂度: O(mlog(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;
}