BZOJ1875 HH去散步

这题 我承认 我过于NAIVE了 我把不能走回边 搞成了 走到一个点必须休息一下。。。

　　发现它的边数非常少 所以 按边建立矩阵 就可以十分科学的完成本题 

　　矩乘我一定要封装进p,q

BZOJ 1875

#include<bits/stdc++.h>
#define mod 45989
#define ll long long
#define breaks printf("!
");
using namespace std;
 
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct edge
{
    int u,v,next;
}vs[300];
int st[80],ee=1;
int n,m,t,A,B;
struct matrix
{
    int mat[130][130];
    void clr()
    {
        for(int i=1;i<=ee;i++)
            for(int j=1;j<=ee;j++)
        mat[i][j]=0;
    }
}M;
matrix operator*(matrix a,matrix b)
{
    matrix tmp; tmp.clr();
    for(int i=1;i<=ee;i++)
        for(int j=1;j<=ee;j++)
            for(int k=1;k<=ee;k++)
    (tmp.mat[i][j]+=1ll*a.mat[i][k]*b.mat[k][j]%mod)%=mod;
    return tmp;
}
matrix operator^(matrix x,int p)
{
    matrix tmp;tmp.clr();
    for(int i=1;i<=ee;i++)
        tmp.mat[i][i]=1;
    for(;p;p>>=1,x=x*x)
        if(p&1)tmp=tmp*x;
    return tmp;
}
inline void addedge(int u,int v)
{
    vs[++ee].v=v;vs[ee].u=u;
    vs[ee].next=st[u];st[u]=ee;
}
int main()
{
//    freopen("a.out","w",stdout);
    n=read();m=read();t=read();
    A=read();B=read();
    for(int i=1;i<=m;i++)
    {
        int a=read(),b=read();
        addedge(a,b);addedge(b,a);
    }
    matrix xx; xx.clr();
    for(int i=st[A];i;i=vs[i].next)
        xx.mat[1][i]++;
    for(int i=2;i<=ee;i++)
        for(int j=2;j<=ee;j++)
    if(vs[i].v==vs[j].u)
    if(i!=(j^1))M.mat[i][j]++;
   // xx.print(); M.print();
    xx=xx*(M^(t-1)); int ans=0;
    for(int i=st[B];i;i=vs[i].next)
        (ans+=xx.mat[1][i^1])%=mod;
    cout<<ans<<endl;
    return 0;
}