（矩阵快速幂）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;
}