矩阵（二分法）

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <string.h>
#include <queue>
#include <vector>    //poj3233 矩阵（二分法）
using namespace std;
int n, m;
struct matrix
{
	int a[35][35];
	void init()
	{
		for(int j=0; j<35; ++j)
		{
			for(int t=0; t<35; ++t)
				a[j][t]=0;
		}
	}
}p;

matrix add(matrix a1, matrix b1)
{
	int t, j;
	matrix q;
	q.init();
	for(j=0; j<n; ++j)
	{
		for(t=0; t<n; ++t)
			q.a[j][t]=(a1.a[j][t]+b1.a[j][t])%m;
	}
	return q;
}

matrix mul(matrix a1, matrix b1)   //a1取行，b1取列
{
	int t, j, g;
	matrix q;
	q.init();
	for(j=0; j<n; ++j)   //行
	{
		for(t=0; t<n; ++t)  //列
		{
			for(g=0; g<n; ++g)
			{
				q.a[j][t]+=a1.a[j][g]*b1.a[g][t];
				q.a[j][t]=q.a[j][t]%m;
			}
		}
	}
	return q;
}

matri