矩阵乘法
// C = A * B
void mul(int c[][N],int a[][N],int b[][N])
{
static int t[N][N]; //作为一个临时的数组
memset(t,0,sizeof t); //先清空
for(int i = 0;i < N;i++) //三次循环
for(int j = 0;j < N;j++)
for(int k = 0;k < N;k++)
t[i][j] = (t[i][j] + (LL)a[i][k] * b[k][j]);
memcpy(c,t,sizeof t);
}