BZOj-1875: [SDOI2009]HH去散步 (矩阵快速幂)

1875: [SDOI2009]HH去散步

Time Limit: 20 Sec  Memory Limit: 64 MB
Submit: 1999  Solved: 980
[Submit][Status][Discuss]

Description

HH有个一成不变的习惯，喜欢饭后百步走。所谓百步走，就是散步，就是在一定的时间 内，走过一定的距离。 但

是同时HH又是个喜欢变化的人，所以他不会立刻沿着刚刚走来的路走回。 又因为HH是个喜欢变化的人，所以他每

天走过的路径都不完全一样，他想知道他究竟有多 少种散步的方法。 现在给你学校的地图（假设每条路的长度都

是一样的都是1），问长度为t，从给定地 点A走到给定地点B共有多少条符合条件的路径

Input

第一行：五个整数N，M，t，A，B。

N表示学校里的路口的个数

M表示学校里的 路的条数

t表示HH想要散步的距离

A表示散步的出发点

B则表示散步的终点。

接下来M行

每行一组Ai，Bi，表示从路口Ai到路口Bi有一条路。

数据保证Ai ！= Bi,但不保证任意两个路口之间至多只有一条路相连接。 

路口编号从0到N -1。 

同一行内所有数据均由一个空格隔开，行首行尾没有多余空格。没有多余空行。 

答案模45989。

N ≤ 20，M ≤ 60，t ≤ 2^30，0 ≤ A,B

Output

一行，表示答案。

Sample Input

4 5 3 0 0
0 1
0 2
0 3
2 1
3 2

Sample Output

4

HINT

 

Source

Day1

laj一开始把题目看错了= = 不能走回头路哇……laj就说这题怎么可能这么傻逼呢……

然后就是很神的做法咯……把边存到矩阵里收尾相接的两条边在矩阵里连边，然后求t-1次幂，在开另一个矩阵，剩下的开代码qwq

laj又忘记膜辣……wa了一发 _(:зゝ∠)_

#include "bits/stdc++.h"
using namespace std;
typedef long long LL;
const int MAX=205,mod=45989;
int n,m,t,st,ed;
struct Mat{
    int x,y;
    int mat[MAX][MAX];
    Mat(){x=y=0;memset(mat,0,sizeof(mat));}
    friend Mat operator * (Mat cc,Mat tt) {
        Mat an;int i,j,k;
        an.x=cc.x,an.y=tt.y;
        for (i=1;i<=cc.x;i++)
            for (k=1;k<=cc.y;k++)
                if (cc.mat[i][k])
                    for (j=1;j<=an.y;j++)
                        an.mat[i][j]=(an.mat[i][j]+cc.mat[i][k]*tt.mat[k][j])%mod;
        return an;
    }
    friend Mat operator ^(Mat cc,int tt) {
        Mat an;int i,j;
        an.x=cc.x,an.y=cc.y;
        for (i=1;i<=cc.x;i++) an.mat[i][i]=1;
        while (tt){
            if (tt&1) an=an*cc;
            cc=cc*cc,tt>>=1;
        }
        return an;
    }
    void out(){
        int i,j;
        for (i=1;i<=x;i++){
            for (j=1;j<=y;j++)
                printf("%d ",mat[i][j]);
            printf("\n");
        }
        printf("\n");
    }
}a,b;
struct Edge{int x,y;}edge[1005];
int tot,head[MAX],adj[1005],next[1005];
void addedge(int u,int v){tot++;adj[tot]=v,next[tot]=head[u],head[u]=tot,edge[tot].x=u,edge[tot].y=v;}
inline int read(){
    int an=0,x=1;char c=getchar();
    while (c<'0' || c>'9') {if (c=='-') x=-1;c=getchar();}
    while (c>='0' && c<='9') {an=(an<<3)+(an<<1)+c-'0';c=getchar();}
    return an*x;
}
int main(){
    freopen ("walk.in","r",stdin);freopen ("walk.out","w",stdout);
    int i,j,x,y;tot=1;
    n=read(),m=read(),t=read(),st=read(),ed=read();st++,ed++;
    for (i=1;i<=m;i++){
        x=read(),y=read();x++,y++;
        addedge(x,y),addedge(y,x);
    }
    a.x=a.y=b.x=b.y=tot;
    for (i=head[st];i;i=next[i]) a.mat[1][i]++;
    for (i=2;i<=tot;i++)
        for (j=2;j<=tot;j++)
            if (edge[i].y==edge[j].x && i!=(j^1))
                b.mat[i][j]++;
    a=a*(b^(t-1));
    int ans=0;
    for (i=head[ed];i;i=next[i])
        ans+=a.mat[1][i^1];
    printf("%d",ans%mod);
    return 0;
}