[TJOI2017]可乐

题目大意：
　　一个$n(nle30)$个点$m(mle100)$条边的无向图，每次可以选择移动到相邻点、原地不动或原地自爆。一开始在编号为$1$的点，问$t(tle10^6)$秒后有多少种可能的行为方案？

思路：
　　新建结点编号$0$表示原地自爆，则在原图基础上对于每个点加上自环和通向$0$的有向边，用邻接矩阵$T$表示。用矩阵$F$记录方案数，其中$F_{0,i}$表示最后一步在$i$处的方案数，初始状态$F_{0,1}=1$。用矩阵快速幂算出$F'=F imes T^t$，答案即为$sum F'_{0,i}$。

#include<cstdio>
#include<cctype>
#include<cstring>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int N=31,mod=2017;
int tmp[N][N];
struct Matrix {
    int n,m,val[N][N];
    void resize(const int &n,const int &m) {
        this->n=n;
        this->m=m;
    }
    Matrix &operator *= (const Matrix another) {
        for(register int i=0;i<n;i++) {
            for(register int j=0;j<another.m;j++) {
                tmp[i][j]=0;
                for(register int k=0;k<m;k++) {
                    (tmp[i][j]+=val[i][k]*another.val[k][j]%mod)%=mod;
                }
            }
        }
        memcpy(val,tmp,sizeof val);
        return *this;
    }
};
Matrix f,t;
int main() {
    const int n=getint(),m=getint();
    f.resize(1,n+1);
    t.resize(n+1,n+1);
    for(register int i=0;i<=n;i++) {
        t.val[i][i]=t.val[i][0]=1;
    }
    for(register int i=0;i<m;i++) {
        const int u=getint(),v=getint();
        t.val[u][v]=t.val[v][u]=1;
    }
    f.val[0][1]=1;
    for(register int k=getint();k;k>>=1,t*=t) {
        if(k&1) f*=t;
    }
    int ans=0;
    for(register int i=0;i<=n;i++) {
        (ans+=f.val[0][i])%=mod;
    }
    printf("%d
",ans);
    return 0;
}