BJ模拟:and(容斥)

传送门

题解:
这里写图片描述

我反过来用类似记忆化搜索的方式实现的,具体看代码。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair <int,int> pii;

const int RLEN=1<<18|1;
inline char nc() {
    static char ibuf[RLEN],*ib,*ob;
    (ib==ob) && (ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
    return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
    char ch=nc(); int i=0,f=1;
    while(!isdigit(ch)) {if(ch=='-')f=-1; ch=nc();}
    while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=nc();}
    return i*f;
}
inline void W(LL x) {
    static int buf[50];
    if(!x) {putchar('0'); return;}
    if(x<0) {putchar('-'); x=-x;}
    while(x) {buf[++buf[0]]=x%10; x/=10;}
    while(buf[0]) putchar(buf[buf[0]--]+'0');
}
const int N=1<<13, mod=1e9+7;
inline int add(int x,int y) {return (x+y>=mod) ? (x+y-mod) : (x+y);}
inline int dec(int x,int y) {return (x-y<0) ? (x-y+mod) : (x-y);}
int n, a[N]; LL f[2][N+50][N+50];
int main() {
    n=rd();
    for(int i=1; i<=n; i++) a[i]=rd();
    int now=1, lim=(1<<13)-1;
    for(int i=1; i<=n; i++) a[i]^=lim;
    f[now][lim][0]++;
    for(int p=n;p>=1;p--) {
        now^=1;
        for(int i=0;i<=lim;i++)
            for(int j=i;;j=(j-1)&i) {
                f[now][i][j]=f[now^1][i|a[p]][j-(j&a[p])];
                f[now][i][j]-=f[now^1][i|a[p]][(j-(j&a[p]))|(a[p]-(a[p]&i))];
                f[now][i][j]+=f[now^1][i][j];
                if(!j) break;
            }
    }
    printf("%d\n",f[now][0][0]%mod);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值