【bzoj1875】【JZYZOJ1354】[SDOI2009]HH去散步 矩阵快速幂 点边转换

http://172.20.6.3/Problem_Show.asp?id=1354

http://www.lydsy.com/JudgeOnline/problem.php?id=1875 

题意：求从起始点到终点走过的路程为x有多少种走法，不保证两点间只有一条路，刚走完一条路时不能原路返回，每条路长度为1。

 

把每条路更换为正反两条边，再造一条出到起点的边s造一条入为终点的边t，将每一条边作为邻接矩阵中的点，矩乘即可

注意题意告诉我们一条路中的正反两条边不能相连。

 

代码

#include<cstdio>
#include<cstring>
#include<iostream>
const int maxn=10010;
const double eps=1e-8;
const int modn=45989;
int n,m,d,s,t;
struct mat{
    long long e[130][130];
    mat(){ memset(e,0,sizeof(e)); }
};mat a;
struct nod{
    int y;
    int f;
    int rev;
    int next;
}e[200];
int tot=0;
int head[50]={};
mat Mul(mat x,mat y){
    mat z;
    for(int i=0;i<=m*2+1;i++){
        for(int j=0;j<=m*2+1;j++){
            for(int k=0;k<=m*2+1;k++){
                z.e[i][j]+=x.e[i][k]*y.e[k][j];
                z.e[i][j]%=modn;
            }
        }
    }
    return z;
}
mat Pow(mat x,int k){
    mat z;
    for(int i=0;i<=m*2+1;i++){
        z.e[i][i]=1;
    }
    while(k>0){
        if(k&1){
            z=Mul(z,x);
        }
        x=Mul(x,x);
        k/=2;
    }
    return z;
}
void init(int x,int y,int rev){
    e[++tot].y=y;
    e[tot].next=head[x];
    e[tot].rev=rev;
    head[x]=tot;
}
void dfs(int x,int v){
    int y;
    if(x==t){
        a.e[v][2*m+1]=1;
    }
    for(int i=head[x];i;i=e[i].next){
        y=e[i].y;
        if(v!=e[i].rev){
            a.e[v][i]=1;
        }
        if(!e[i].f){
            e[i].f=1;
            dfs(y,i);
        }
    }
}
int main(){
    scanf("%d%d%d%d%d",&n,&m,&d,&s,&t);
    s+=1,t+=1;
    int x,y;
    for(int i=1;i<=m;i++){
        scanf("%d%d",&x,&y);
        init(x+1,y+1,tot+2);
        init(y+1,x+1,tot);
    }dfs(s,0);
    mat z=Pow(a,d+1);
    printf("%lld
",z.e[0][2*m+1]);
    return 0;
}