bzoj1875 [SDOI2009]HH去散步

Description

---

HH有个一成不变的习惯，喜欢饭后百步走。所谓百步走，就是散步，就是在一定的时间 内，走过一定的距离。 但
是同时HH又是个喜欢变化的人，所以他不会立刻沿着刚刚走来的路走回。 又因为HH是个喜欢变化的人，所以他每
天走过的路径都不完全一样，他想知道他究竟有多 少种散步的方法。 现在给你学校的地图（假设每条路的长度都
是一样的都是1），问长度为t，从给定地 点A走到给定地点B共有多少条符合条件的路径

数据保证Ai ！= Bi,但不保证任意两个路口之间至多只有一条路相连接。 
路口编号从0到N -1。 
同一行内所有数据均由一个空格隔开，行首行尾没有多余空格。没有多余空行。 
答案模45989。
N ≤ 20，M ≤ 60，t ≤ 2^30，0 ≤ A,B

来自 https://www.lydsy.com/JudgeOnline/problem.php?id=1875

Solution

---

仔细审题发现并不能来回走，因此直接上矩阵乘法是会wa的

考虑dp，设f[I,j]表示走了i步到达了第j条边的终点的方案数，那么f[I+1,j]=Σf[I,k] (满足e[k].y=e[j].x)
第一维可以去掉，那么转移矩阵就很显然了。注意边是双向的

Code

---

#include <stdio.h>
#include <string.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))

const int MOD=45989;
const int N=125;

int rc[N][N],f[N];
int n,m,t;
int ls[N],edCnt=1;

struct edge {int x,y,next;} e[N];

struct mat {
    int a[N][N];
    mat operator *(mat b) {
        mat c; fill(c.a,0);
        rep(i,1,edCnt) rep(j,1,edCnt) {
            rep(k,1,edCnt) c.a[i][j]=(c.a[i][j]+a[i][k]*b.a[k][j]%MOD)%MOD;
        }
        return c;
    }
} A,ret;

void add_edge(int x,int y) {
    e[++edCnt]=(edge) {x,y,ls[x]}; ls[x]=edCnt;
    e[++edCnt]=(edge) {y,x,ls[y]}; ls[y]=edCnt;
}

mat ksm(mat x,int dep) {
    ret=x; dep-=1;
    while (dep) {
        if (dep&1) ret=ret*x;
        x=x*x; dep/=2;
    }
    return ret;
}

int main(void) {
    int st,ed;
    scanf("%d%d%d%d%d",&n,&m,&t,&st,&ed);
    st++; ed++;
    rep(i,1,m) {
        int x,y; scanf("%d%d",&x,&y);
        x++; y++;
        add_edge(x,y);
    }
    rep(i,2,edCnt) rep(j,2,edCnt) {
        if (e[i].x==e[j].y&&i!=(j^1)) A.a[j][i]=1;
    }
    A=ksm(A,t-1);
    for (int i=ls[st];i;i=e[i].next) f[i]=1;
    int ans=0;
    rep(i,1,edCnt) if (e[i].y==ed) {
        rep(j,1,edCnt) ans=(ans+f[j]*A.a[j][i]%MOD)%MOD;
    }
    printf("%d\n", ans);
    return 0;
}