【bzoj4887】：[Tjoi2017]可乐 矩阵乘法,快速幂

【bzoj4887】：[Tjoi2017]可乐

题目大意：一张无相连通图(n<=30)，从1号点开始走，每秒可以走到相邻的点也可以自爆，求第t秒(t<=1e6)后所有的方案数是多少对2017取模

恩。。就是一个矩阵快速幂。。矩阵就是原图的邻接矩阵。。然后f[i][i]也是1。。

但是这是不会自爆的情况下的矩阵，算上自爆的话要把每次转移的结果求和。。蒟蒻想了半天。。

然后发现其实只要再加一行一列，然后f[n+1][i]=1,就可以了。。

意会一下好了。。矩阵什么的感觉讲不清楚啊。。

/* http://www.cnblogs.com/karl07/ */
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;

struct MAT{
    int m[35][35],x,y,sum;
    MAT () {for (int i=0;i<35;i++) for (int j=0;j<35;j++) m[i][j]=0; sum=x=y=0;}
}a,b;
int n,m,t;
int mp[35][35];

MAT operator *(const MAT &a,const MAT &b){
    MAT c;
    c.x=b.x;c.y=a.y;
    for (int i=1;i<=c.x;i++){
        for (int j=1;j<=c.y;j++){
            c.m[i][j]=0;
            for (int k=1;k<=a.x;k++){
                c.m[i][j]+=a.m[k][j]*b.m[i][k]%2017;
                c.m[i][j]%=2017;
            }
            c.sum+=c.m[i][j];
            c.sum%=2017;
        }
    }
    return c;
} 

int main(){
    scanf("%d%d",&n,&m);
    a.x=n+1,a.y=1,a.m[1][1]=1;
    b.x=b.y=n+1;
    for (int i=1;i<=m;i++){
        int p,q;
        scanf("%d%d",&p,&q);
        b.m[p][q]++;
        b.m[q][p]++;
    }
    for (int i=1;i<=n+1;i++){
        b.m[i][i]=1;
        b.m[n+1][i]=1;
    }
    scanf("%d",&t);
    while (t){
        if (t&1) a=a*b;
        b=b*b;
        t=(t>>1);
    }
    printf("%d
",a.sum);
    return 0;
}

啊。。终于不是游记了。。