BZOJ1875: [SDOI2009]HH去散步

矩乘，因为不能走刚走过的边，所以矩阵里面的位置代表边，一个条边的两端开两个点放进矩阵里，这个点的值代表刚走过这条边走到这个点的方案数

---

code

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 130;
const int Mod = 45989;
struct mat
{
    int a[maxn][maxn];
}st,one,zero,zz;
int n,m,t,sq,eq;
int d[maxn][maxn],num[maxn];
int edge[maxn][2];
int N;

mat operator *(mat x,mat y)
{
    mat r;
    for(int i=0;i<N;i++)
    {
        for(int j=0;j<N;j++)
        {
            r.a[i][j]=0;
            for(int k=0;k<N;k++)
            {
                r.a[i][j]+=x.a[i][k]*y.a[k][j]%Mod;
            }
            r.a[i][j]%=Mod;
        }
    }
    return r;
}
mat pw(int k)
{
    mat r=zz,nw=st;
    for(;k;k>>=1,nw=nw*nw)
        if(k&1)r=r*nw;
    return r;
}

int main()
{
    scanf("%d%d%d%d%d",&n,&m,&t,&sq,&eq); sq++; eq++;
    if(!t){printf("%d\n",sq==eq);return 0;}
    N=m*2+1;

    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y); x++; y++;
        edge[i][0]=x; edge[i][1]=y;
        num[x]++; d[x][num[x]]=i;
        num[y]++; d[y][num[y]]=i;
    }
    for(int i=0;i<N;i++) 
        for(int j=0;j<N;j++)
            st.a[i][j]=zero.a[i][j]=0,zz.a[i][i]=1;

    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=num[i];j++)
        {
            int k=d[i][j];
            int x=(edge[k][1]==i);
            int idx=k*2-1+x;
            for(int l=1;l<=num[i];l++)
                if(j!=l)
                {
                    int ky=d[i][l];
                    int y=(edge[ky][0]==i);
                    int idy=ky*2-1+y;
                    st.a[idx][idy]++;
                }
        }
    }

    t--;

    one=zero;
    for(int i=1;i<=num[sq];i++)
    {
        int k=d[sq][i];
        int x=(edge[k][0]==sq);
        one.a[1][k*2-1+x]++;
    }
    mat r=one*pw(t);
    int ans=0;
    for(int i=1;i<=num[eq];i++)
    {
        int k=d[eq][i];
        int x=(edge[k][1]==eq);
        ans+=r.a[1][k*2-1+x];
    }
    printf("%d\n",ans%Mod);

    return 0;
}