[DP 分块子集枚举] UOJ#300.【CTSC2017】吉夫特

最新推荐文章于 2019-03-04 22:32:00 发布

LowestJN

最新推荐文章于 2019-03-04 22:32:00 发布

阅读量432

点赞数 1

分类专栏： DP

本文链接：https://blog.csdn.net/Coldef/article/details/79090203

版权

DP 专栏收录该内容

84 篇文章 0 订阅

订阅专栏

DP转移很好写

暴力是 $O(3^{18})$ 的

把每个二进制数分成前9位和后9位，分块统计复杂度就是 $O(6^9)$ 了

#include <cstdio>
#include <iostream>
#include <algorithm>

using namespace std;

const int N=333333,P=1e9+7;

int n,ans,a[N],f[N];

inline int Query(int x){
  int a=x&~511,b=(~x)&511,ret=0;
  for(int i=b;i;i=(i-1)&b)
    ret=(ret+f[a|((~i)&511)])%P;
  return (ret+f[a|511])%P;
}

inline void Add(int x,int y){
  int a=x>>9,b=x&511;
  for(int i=a;i;i=(i-1)&a)
    f[i<<9|b]=(f[i<<9|b]+y)%P;
  f[b]=(f[b]+y)%P;
}

int main(){
  freopen("1.in","r",stdin);
  freopen("1.out","w",stdout);
  scanf("%d",&n);
  for(int i=1;i<=n;i++) scanf("%d",&a[i]);
  for(int i=1;i<=n;i++){
    int cur=Query(a[i])+1;
    ans=(ans+cur)%P;
    Add(a[i],cur);
  }
  printf("%d\n",(ans+P-n)%P);
  return 0;
}