bzoj1875 HH去散步

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。

Output

一行，表示答案。

Sample Input

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

Sample Output

4

HINT

對於30%的數據，N ≤ 4，M ≤ 10，t ≤ 10。 對於100%的數據，N ≤ 20，M ≤ 60，t ≤ 2^30，0 ≤ A,B

Source

Day1

矩陣乘法，跟hdu的那道題一樣，模數變一下而已。

http://blog.csdn.net/lzr010506/article/details/51136389

換了一種寫法。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define mod 45989
#define N 22
using namespace std;
int n,m,t,a,b;
struct Matrix
{
    int map[N * 6][N * 6];
    void clear() {memset(map,0,sizeof(map));}
    friend Matrix operator * (Matrix &a,Matrix &b)
    {
        Matrix ret;
        for(int i = 0; i <= 2 * m; i ++)
            for(int j = 0; j <= 2 * m; j ++)
            {
                ret.map[i][j] = 0;
                for(int k = 0; k <= 2 * m; k ++)      
                    ret.map[i][j] = ret.map[i][j] + a.map[i][k] * b.map[k][j] % mod;
                ret.map[i][j] %= mod;
            }
        return ret;
    }
    friend Matrix operator ^ (Matrix &a,int b)
    {
        Matrix ret;
        ret.clear();
        for(int i = 0; i <= m * 2; i ++) ret.map[i][i] = 1;
        while(b)
        {
            if(b & 1) ret = ret * a;
            a = a * a;
            b >>= 1;
        }
        return ret;
    }
}ori,base;
int head[N],cnt;
struct node
{
    int from,to,next;
}edge[N * 6];
void init()
{
    memset(head,-1,sizeof(head));
    cnt = 1;
}
void edgeadd(int from,int to)
{
    edge[cnt].from = from;
    edge[cnt].to = to;
    edge[cnt].next = head[from];
    head[from] = cnt ++;
}
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;
}
int main()
{
    init();
    n = read();
    m = read();
    t = read();
    a = read();
    b = read();
    a ++;
    b ++;
    for(int i = 1; i <= m; i ++)
    {
        int x, y;
        x = read();
        y = read();
        x ++;
        y ++;
        edgeadd(x, y);
        edgeadd(y, x);
    }
    for(int i = head[a]; i != -1; i = edge[i].next)
        base.map[0][i] ++;
    for(int i = 1; i < cnt; i ++)
    {
        int x = edge[i].from, y = edge[i].to;
        for(int j = head[y]; j != -1; j = edge[j].next)
            if(j != (i + ((i & 1)? 1: -1)))
                ori.map[i][j] ++;
    }
    ori = ori ^ (t - 1);
    base = base * ori;
    int sum = 0;
    for(int i = 1; i < cnt; i ++)
        if(edge[i].to == b)
            sum = (sum + base.map[0][i])%mod;
    printf("%d\n",sum);
    return 0;
}