Time Limit: 20 Sec Memory Limit: 64 MB
Submit: 2333 Solved: 1204
[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

不能沿著剛剛走來的路走回
然後就很煩啊

解決方案：把雙向邊拆成兩條單向邊後按邊建矩陣就不會走回去啦啦啦~(≧▽≦)/~

#include<iostream>
#include<cstdio>
#define M 45989
#define LL long long 
using namespace std;

int i,m,n,j,k,a,b,t, x,y,s[1000][2],ans;

struct vv
{
    int g[200][200];
} ch,f,cs,l;

vv cheng(vv a,vv b)
{
    vv d=ch;
    for(int i=1;i<=m+m;i++)
        for(int j=1;j<=m+m;j++) 
            for(int k=1;k<=m+m;k++) 
                d.g[i][j]=(int)(d.g[i][j]+(LL)a.g[i][k]*b.g[k][j])%M; 
    return d;
}


vv ksm(vv a,int x)
{
    vv c=cs;
    for(x;x>1;x>>=1)
    {
        if(x&1) c=cheng(c,a);
        vv v=cheng(a,a);
        a=v;
    }
    return cheng(a,c);
}

int main()
{
    scanf("%d%d%d%d%d",&n,&m,&t,&a,&b);
    for(i=1;i<=m+m;i++) cs.g[i][i]=1;
    for(i=1;i<=m;i++) scanf("%d%d",&s[i][0],&s[i][1]), s[i+m][0]=s[i][1], s[i+m][1]=s[i][0];
    for(i=1;i<=m+m;i++)
        for(j=1;j<=m+m;j++)
            if(s[i][1]==s[j][0] && i-j!=m && j-i!=m) f.g[i][j]=1;
    vv d=ksm(f,t-1);
    for(i=1;i<=m+m;i++)
        if(s[i][0]==a) l.g[1][i]=1;
    l=cheng(l,d);
    for(i=1;i<=m+m;i++) if(s[i][1]==b) ans=(ans+l.g[1][i])%M; 
    printf("%d",ans);
}