【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[][],x,y,sum;
     MAT () {for (int i=;i<;i++) for (int j=;j<;j++) m[i][j]=; sum=x=y=;}
 }a,b;
 int n,m,t;
 int mp[][];

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

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

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