（矩阵快速幂）hdu2157 How many ways??

题意大概是给出一个图的邻接矩阵，求出从a点走到b点长度恰好为k的路的条数。很基础的离散数学知识，这种路的条数即是邻接矩阵乘k次幂后对应位置a[i][j]的值，只不过数据范围较大，需要用矩阵快速幂来实现。其实是套模板稍加修改就可以过的一道题，但因为刚接触，天生愚钝，WA了无数次调了一个多小时才勉强ac，伤不起......

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

const int mod = 1000;
const int maxn = 20;
const int maxm = 100;

int n, m, T, A, B, k;

struct mat {
    ll a[maxn][maxn];
    mat() {
        memset(a, 0, sizeof(a));
    }
};

void print(mat ans) {
    for(int i = 0;i < n;i++) {
        for(int j = 0;j < n;j++)
            printf("%d ", ans.a[i][j]);
        printf("\n");
    }
    return;
}

mat init0() {
    mat tt;
    for(int i = 0;i < n;i++)
        for(int j = 0;j < n;j++) tt.a[i][j] = 0;
    return tt;
}

mat init1() {
    mat tt;
    for(int i = 0;i < n;i++)
        for(int j = 0;j < n;j++) tt.a[i][j] = (i==j);
    return tt;
}

mat mul(mat x, mat y) {
    mat tt;
    for(int i = 0;i < n;i++) {
        for(int j = 0;j < n;j++) {
            tt.a[i][j] = 0;
            for(int k = 0;k < n;k++) {
                tt.a[i][j] += x.a[i][k] * y.a[k][j];
            }
            tt.a[i][j] %= mod;
        }
    }
    return tt;
}

mat power(mat tt, int b) {
    mat res = init1();
    while(b) {
        if(b & 1) {
            res = mul(res, tt);
        }
        b >>= 1;
        tt = mul(tt, tt);
    }
    return res;
}

int main() {
    while(scanf("%d %d", &n, &m) != EOF && (n||m)) {
        int l, r;
        mat gra = init0();
        for(int i = 0;i < m;i++) {
            scanf("%d %d", &l, &r);
            gra.a[l][r] = 1;
        }
        scanf("%d", &T);
        while(T--) {
            scanf("%d %d %d", &A, &B, &k);
            mat b = power(gra, k);
            printf("%d\n", b.a[A][B]);
        }
    }
    return 0;
}