luogu4159 迷路 (矩阵加速)

考虑如果只有距离为1的边，那我用在时间i到达某个点的状态数矩阵 乘上转移矩阵（就是边的邻接矩阵），就能得到i+1时间的

然后又考虑到边权只有1~9，那可以把边拆成只有距离为1的

具体做法是一个点拆成9个然后串联

 #include<bits/stdc++.h>
 #define pa pair<ll,ll>
 #define CLR(a,x) memset(a,x,sizeof(a))
 #define MP make_pair
 using namespace std;
 typedef long long ll;
 const int maxn=,maxm=,P=;

 inline char gc(){
     return getchar();
     static const int maxs=<<;static char buf[maxs],*p1=buf,*p2=buf;
     return p1==p2&&(p2=(p1=buf)+fread(buf,,maxs,stdin),p1==p2)?EOF:*p1++;
 }
 inline ll rd(){
     ll x=;char c=gc();bool neg=;
     while(c<''||c>''){if(c=='-') neg=;c=gc();}
     while(c>=''&&c<='') x=(x<<)+(x<<)+c-'',c=gc();
     return neg?(~x+):x;
 }

 int id[maxn][maxn],N,M,T;
 int trans[maxm][maxm],tmp[maxm][maxm],ans[maxm][maxm];

 int main(){
     //freopen("","r",stdin);
     int i,j,k=;
     N=rd();T=rd();
     for(i=;i<=N;i++){
         for(j=;j<=;j++){
             id[i][j]=++M;
             if(j!=) trans[id[i][j]][id[i][j-]]=;
         }
     }
     for(i=;i<=N;i++){
         char s[];scanf("%s",s+);
         for(j=;j<=N;j++){
             if(s[j]!='') trans[id[i][]][id[j][s[j]-'']]=;
         }
     }
     T--;
     memcpy(ans,trans,sizeof(ans));
     while(T){
         if(T&){
             CLR(tmp,);
             for(i=;i<=M;i++){
                 for(j=;j<=M;j++){
                     for(k=;k<=M;k++){
                         tmp[i][j]=(tmp[i][j]+ans[i][k]*trans[k][j])%P;
                     }
                 }
             }
             memcpy(ans,tmp,sizeof(tmp));
         }
         CLR(tmp,);
         for(i=;i<=M;i++){
             for(j=;j<=M;j++){
                 for(k=;k<=M;k++){
                     tmp[i][j]=(tmp[i][j]+trans[i][k]*trans[k][j])%P;
                 }
             }
         }
         memcpy(trans,tmp,sizeof(tmp));
         T>>=;
     }
     printf("%d\n",ans[id[][]][id[N][]]);
     return ;
 }