1. 程式人生 > >bzoj 4903: [Ctsc2017]吉夫特【lucas+狀壓dp】

bzoj 4903: [Ctsc2017]吉夫特【lucas+狀壓dp】

二進制 clu 根據 main namespace ios ctsc ace 狀壓dp

首先根據lucas,
\[ C_n^m\%2=C_{n\%2}^{m\%2}*C_{n/2}^{m/2} \]
讓這個式子的結果為計數的情況只有n&m==m,因為m的每一個為1的二進制位都需要n中這一位為1,否則結果就是0
所以枚舉子集,設f[i]為以i開頭的合法子序列個數,dp的時候枚舉子集從後往前dp即可

#include<iostream>
#include<cstdio>
using namespace std;
const int N=300005,mod=1e9+7;
int n,a[N],p[N],f[N],ans;
int read()
{
    int r=0,f=1;
    char p=getchar();
    while(p>‘9‘||p<‘0‘)
    {
        if(p==‘-‘)
            f=-1;
        p=getchar();
    }
    while(p>=‘0‘&&p<=‘9‘)
    {
        r=r*10+p-48;
        p=getchar();
    }
    return r*f;
}
void jia(int &x,int y)
{
    x+=y;
    x>=mod?x-=mod:0;
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
    {
        a[i]=read();
        p[a[i]]=i;
    }
    for(int i=1;i<=233333;i++)
        if(p[i])
        {
            for(int j=i;j;j=(j-1)&i)
                if(p[j]>p[i])
                    jia(f[i],f[j]);
            jia(ans,f[i]++);
        }
    printf("%d\n",ans);
    return 0;
}

bzoj 4903: [Ctsc2017]吉夫特【lucas+狀壓dp】