[SDOI2009] HH去散步 (矩阵乘法)

link

$solution:$

将边化为点后重新建矩阵，跑$T-1$幂即可(因为跑的是新边)。

最后直接找与$x,y$所相连的边即可。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define int long long
#define mod 45989
using namespace std;
inline int read(){
    int f=,ans=;char c=getchar();
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}
    return f*ans;
}
const int MAXN=;
struct Matrix{
    int a[MAXN<<][MAXN<<];
    void init(){memset(a,,sizeof(a));return;}
}F,G;
int cnt,n,m,A,B,T,sta[MAXN][MAXN<<],End[MAXN<<];
Matrix operator*(Matrix x1,Matrix x2){
    Matrix x3;x3.init();
    for(int i=;i<=(m<<);i++)
        for(int k=;k<=(m<<);k++)
            for(int j=;j<=(m<<);j++) x3.a[i][j]+=x1.a[i][k]*x2.a[k][j],x3.a[i][j]%=mod;
    return x3;
}
Matrix ksm(Matrix a,int b){
    Matrix ans;ans.init();
    for(int i=;i<=(m<<);i++) ans.a[i][i]=;
    while(b){
        if(b&) ans=ans*a;
        a=a*a,b>>=;
    }return ans;
}int ans;
int uu[MAXN<<],vv[MAXN<<];
bool check(int idi,int idj){
    if(vv[idi]==uu[idj]) return ;
    return ;
}
signed main(){
    n=read(),m=read(),T=read(),A=read()+,B=read()+;
    for(int i=;i<=m;i++){
        int u=read()+,v=read()+;
        if(v==B) End[++End[]]=i;
        if(u==B) End[++End[]]=i+m;
        sta[u][++sta[u][]]=i;
        uu[i]=u,vv[i]=v;
        sta[v][++sta[v][]]=i+m;
        uu[i+m]=v,vv[i+m]=u;
    }F.init(),G.init();
    for(int i=;i<=(m<<);i++){
        for(int j=;j<=(m<<);j++) if(check(i,j)&&(i+m!=j&&j+m!=i)) G.a[i][j]=;
    }
    //print(G);
    //printf("=========================\n");
    F=ksm(G,T-);
    //print(F);
    //for(int i=1;i<=End[0];i++) printf("End(%d):%d\n",i,End[i]);
    //for(int i=1;i<=sta[A][0];i++) printf("sta(%d):%d\n",i,sta[A][i]);
    for(int i=;i<=sta[A][];i++)
        for(int j=;j<=End[];j++) ans+=F.a[sta[A][i]][End[j]],ans%=mod;
    printf("%lld\n",ans);
}