给定一个有向图，问从A点恰好走k步（允许重复经过边）到达B点的方案数---矩阵乘法

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <string.h>
#include <vector>
#include <queue>    //hdu2157
#define mod 1000 
using namespace std;
/*
题目大意：给定一个有向图，问从A点恰好走k步（允许重复经过边）到达B点的方案数mod p的值
把给定的图转为邻接矩阵，即A(i,j)=1当且仅当存在一条边i->j。令C=A*A，那么C(i,j)=ΣA(i,k)*A(k,j)，
实际上就等于从点i到点j恰好经过2条边的路径数（枚举k为中转点。
类似地，C*A(这里的C已经是A*A)的第i行第j列就表示从i到j经过3条边的路径数。同理，如果要求经过k步的 路径数，我们只需要二分求出A^k即可。
注意：矩阵多次幂可以防止爆栈（相比二分法求矩阵多次幂）
*/
int n;
struct matrix
{
	int a[25][25];
	void init()
	{
		for(int j=0; j<25; ++j)
		{
			for(int t=0; t<25; ++t)
				a[j][t]=0;
		}
	}
}p;

matrix mul(matrix a1, matrix b1)
{
	matrix q;
	q.init();
	int t, j, k;
	for(j=0; j<n; ++j)
	{
		for(t=0; t<n; ++t)
		{
			for(k=0; k<n; &#