【BZOJ 2004】: [Hnoi2010]Bus 公交线路

题目链接：

　　TP

题解：

　　  所以说，超显眼的数据范围啊。

　　很显然我们对于每个P的区间都是要有k个站被bus停留，然后考虑转移的话应该是把这k个站里的某个bus往前走，那么转移也很显然了，n的范围很大，所以直接上矩阵。

代码：

#define Troy 

#include <bits/stdc++.h>

using namespace std;

inline int read(){
    int s=,k=;char ch=getchar();
    while(ch<''|ch>'')    ch=='-'?k=-:,ch=getchar();
    while(ch>&ch<='')    s=s*+(ch^),ch=getchar();
    return s*k;
}

const int N=,mod=;

int n,k,p,BGM,stk[N],m;

struct Matrix{
    int a[N][N];
    Matrix(){memset(a,,sizeof(a));}
    inline void evoid(){
        for(int i=;i^m;++i)
            a[i][i]=;
    }
    inline friend Matrix operator *(const Matrix &x,const Matrix &y){
        Matrix c;
        for(int i=;i^m;++i)
            for(int j=;j^m;++j)
                for(int k=;k^m;++k)
                    c.a[i][j]=(c.a[i][j]+x.a[i][k]*y.a[k][j]%mod)%mod;
        return c;
    }
    inline friend Matrix operator ^(Matrix x,int b){
        Matrix ret;ret.evoid();
        while(b){
            if(b&)ret=ret*x;
            b>>=;x=x*x;
        }return ret;
    }
};

inline bool move_to(int to,int from){
    from^=<<p-;
    from<<=;
    int deta=from^to;
    return (deta&(-deta))==deta;
}

int main(){
    n=read(),k=read(),p=read();
    register int i,j;
    for(i=(<<p-);i^(<<p);++i){
        int x=i;
        for(j=;x;x^=(x&(-x)),++j);
        if(j==k){
            if(i==(<<p)-(<<p-k))
                BGM=m;
            stk[m++]=i;
        }
    }
    Matrix ans,t;
    ans.a[BGM][]=;
    for(i=;i^m;++i)
        for(j=;j^m;++j){
                t.a[i][j]=move_to(stk[j],stk[i]);
        }
    ans=(t^(n-k))*ans;
    printf("%d\n",ans.a[BGM][]);
}