bzoj1297 [SCOI2009]迷路——拆点+矩阵快速幂

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1297

一看感觉是矩阵快速幂之类的，但边权不好处理啊；

普通的矩阵快速幂只能处理边权为1的，所以想办法把边权处理成1；

仔细一看还有一个条件是边权小于10；

所以拆点！把一个点拆成10个点表示到它不同的距离，那么和它相连的那些点就可以跟某个距离的点连边权为1的边；

虽然没有自己想出来，不过1A还是极好的！(因为太简单了)

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,T,w[][],tot,id[][],mod=;
struct Matrix{
    int a[][];
    Matrix(){memset(a,,sizeof a);}
    Matrix operator * (const Matrix &y) const
    {
        Matrix x;
        for(int i=;i<=tot;i++)
            for(int k=;k<=tot;k++)
                for(int j=;j<=tot;j++)
                    (x.a[i][j]+=a[i][k]*y.a[k][j])%=mod;
        return x;
    }
    void init(){for(int i=;i<=tot;i++)a[i][i]=;}
}f,ans;
Matrix pw(Matrix x,int k)
{
    Matrix ret; ret.init();
    for(;k;k>>=,x=x*x)
        if(k&)ret=ret*x;
    return ret;
}
int main()
{
    scanf("%d%d",&n,&T);
    char ch[];
    for(int i=;i<n;i++)
    {
        scanf("%s",&ch);
        for(int j=;j<n;j++)
            w[i][j]=ch[j]-'';
    }
    for(int i=;i<n;i++)
        for(int j=;j<=;j++)
        {
            id[i][j]=++tot;
            if(j)f.a[tot-][tot]++;
        }
    for(int i=;i<n;i++)
        for(int j=;j<n;j++)
            if(w[i][j])f.a[id[i][w[i][j]-]][id[j][]]++;
    ans.a[][id[][]]=;
    Matrix fn=pw(f,T);
    ans=ans*fn;
    printf("%d",ans.a[][id[n-][]]);
    return ;
}