bzoj1875 [SDOI2009]HH去散步——矩阵快速幂

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

有个限制是不能走回头路，比较麻烦；

所以把矩阵中的元素设成边的经过次数，单向边之间就好转移了；

最后从单向边的经过次数得到点的路径方案数。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int const mod=;
int n,m,T,A,B,head[],ct=,tot;
struct N{
    int hd,to,next;
    N(int h=,int t=,int n=):hd(h),to(t),next(n) {}
}edge[];
struct Matrix{
    int a[][];
    Matrix(){memset(a,,sizeof a);}//
    Matrix operator * (const Matrix &y) const
    {
        Matrix x;
        for(int i=;i<=ct;i++)
            for(int k=;k<=ct;k++)
                for(int j=;j<=ct;j++)
                    (x.a[i][j]+=a[i][k]*y.a[k][j])%=mod;
        return x;
    }
    void init()
    {
        for(int i=;i<=ct;i++)a[i][i]=;
    }
}tr,ans;
Matrix operator ^ (Matrix x,int y)
{
    Matrix ret; ret.init();
    for(int i=y;i;i>>=,x=x*x)
        if(i&)ret=ret*x;//i 而不是 y!!!
    return ret;
}
void add(int x,int y){edge[++ct]=N(x,y,head[x]); head[x]=ct;}
int main()
{
    scanf("%d%d%d%d%d",&n,&m,&T,&A,&B);
    for(int i=,x,y;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        add(x,y); add(y,x);
    }
    for(int i=head[A];i;i=edge[i].next) ans.a[][i]++;
    for(int i=;i<=ct;i++)//
        for(int j=;j<=ct;j++)
            if(edge[i].to==edge[j].hd && i!=(j^)) tr.a[i][j]++;//
    ans=ans*(tr^(T-));
    for(int i=head[B];i;i=edge[i].next)
        (tot+=ans.a[][i^])%=mod;
    printf("%d",tot);
    return ;
}