bzoj 2004

这几天就顾着看以前的题，不写的码手还是会生疏啊

状态压缩+矩阵乘法

发现在长度为P的区间中符合题意的很有限，所以先压缩完建出图，用矩阵加速转移

为什么有状态压缩大家都要弄个DP上去呢，感觉只有压缩的部分没有DP呀。。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<iostream>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define mem(a) memset(a,0,sizeof(a))
#define inf 30031
#define ll long long
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define NM 200
using namespace std;
int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int c[NM],n,m,p,cnt,ans;
struct mat{
    int a[NM][NM];
}s,t;
mat operator*(const mat&x,const mat&y){
    mat s;mem(s.a);
    inc(i,1,cnt)inc(j,1,cnt)inc(k,1,cnt)
    (s.a[i][j]+=x.a[i][k]*y.a[k][j])%=inf;
    return s;
}
void dfs(int x,int t,int num){
    if(x==p){
        if(num==m)c[++cnt]=t;return;
    }
    dfs(x+1,t<<1,num);dfs(x+1,t<<1|1,num+1);
}
void out(mat t){
    inc(i,1,cnt){
        inc(j,1,cnt)printf("%d ",t.a[i][j]);putchar('
');
    }putchar('
');
}
int main(){
    freopen("data.in","r",stdin);
    n=read();m=read();p=read();
    dfs(1,1,1);
    inc(i,1,cnt)inc(j,1,cnt){
        int x=(c[i]<<1)^succ(p)^c[j];
        if(x==lowbit(x))t.a[i][j]=1;
    }
//    out(t);
    inc(i,1,cnt)s.a[i][i]=1;
    p=n-m;
    for(mat V=t;p;p>>=1,V=V*V)
    if(p&1)s=s*V;
//    out(s);
    printf("%d
",s.a[cnt][cnt]);
    return 0;
}