sicily 1211 商人的宣传

Description

 Bruce是K国的商人，他在A州成立了本身的公司，此次他的公司生产出了一批性能很好的产品，准备宣传活动开始后的第L天到达B州进行新品拍卖，期间Bruce打算将产品拿到各个州去作推销宣传,以增长其影响力。

K国有不少个州,每一个州都与其余一些州相邻,可是K国对商人做宣传却有一些很奇怪的规定：
一、 商人只能从某些州到达另一些州，即连通路线是单向的，并且有些州多是到达不了的。
二、 商人不容许在同一个州连续宣传两天或以上，天天宣传完必须离开该州。
三、 商人能够屡次来到同一个州进行宣传。

"我必须找出一条影响力最大的路线才行"，Bruce想，"但我首先必须知道到底有多少这种符合规定的宣传路线可供我选择。"如今Bruce把任务交给了你。而且出于考虑之后的须要，你还要帮他算出给出的两州之间的路线的总数。

Input

输入文件第一行包含三个整数n,m,L（1≤n,L≤100），分别表示K国的州数、连通路线的数量，以及多少天后必须到达B州。接下来有m行，每行一队整数x,y（1≤x,y≤n），表示商人能从x州到达y州。
第m+2行为一个整数q（1≤q≤100），表示Bruce有q个询问。下面q行每行两个整数A,B（1≤A,B≤n），即A、B州的位置。

Output

输出文件包含q行，每行一个整数t，为所求的从A州到B州知足上述规定的路线总数。
输入数据中的询问将保证答案t在长整数范围内，即t<2 ³¹。

Sample Input

4 5 6
1 2
2 3
3 4
4 1
2 4
2
1 4
4 2

Sample Output

2
1

分析：

比较偷懒的作法是直接使用矩阵乘法。具体作法是，能够初始化地图矩阵m1全部节点为0，而后一旦联通则设为1，那么2天内由A到B的方法为M1*M1后获得的矩阵的第a行第b列，递推能够知道第l天方法数为矩阵吗m1连乘l-1次所得矩阵。这种方法写法简便，可是速度可能稍慢。

代码：

// Problem#: 1211
// Submission#: 1899573
// The source code is licensed under Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License
// URI: http://creativecommons.org/licenses/by-nc-sa/3.0/
// All Copyright reserved by Informatic Lab of Sun Yat-sen University
#include <iostream>
#include <cstring>
using namespace std;

#define MAX 100

int main() {
    int a, b;
    int m1[MAX][MAX] = {0};
    int m2[MAX][MAX] = {0};
    int m3[MAX][MAX] = {0};
    int n, m, l, q;
    cin >> n >> m >> l;
    while (m--) {
        cin >> a >> b;
        ++m1[a - 1][b - 1];
        ++m2[a - 1][b - 1];
    }
    for (int c = 0; c < l - 1; ++c) {
        memset(m3, 0, sizeof(m3));
        for (int i = 0; i < n; ++i) {
            for (int j = 0; j < n; ++j)
                for(int k = 0; k < n; ++k)
                    m3[i][j] += m1[i][k] * m2[k][j]; 
        }
        memcpy(m1, m3, sizeof(m1));
    }
    cin >> q;
    while (q--) {
        cin >> a >> b;
        cout << m3[a - 1][b - 1] << endl;
    }
    return 0;
}