hdu 4920 快速矩阵相乘 以后得换一种写法了

【题意】：求两个矩阵相乘的结果 

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define mod 3

int a[808][802],b[808][802];
int c[808][808],n;

void mul()
{
    for(int i=0; i<n; i++)
        for(int j=0; j<n; j++)
            if(a[i][j])
            {
                for(int k=0; k<n; k++)
                    //if(b[j][k])    加了这一句居然就超时了 实在不晓得是为什么
                    c[i][k]=(c[i][k]+a[i][j]*b[j][k]);
                //取模运算是最耗时的运算 经计算这里不会超过数据范围 所以不要一步一取模，最后一步来取余就可以了
            }
}

int main()
{
    int i,j,m,t;
    while(~scanf("%d",&n))
    {
        for(i=0; i<n; i++)
            for(j=0; j<n; j++)
               {
                    scanf("%d",&a[i][j]);
                    a[i][j]%=3;
               }
        for(i=0; i<n; i++)
            for(j=0; j<n; j++)
            {
                scanf("%d",&b[i][j]);
                b[i][j]%=3;
            }
        memset(c,0,sizeof(c));
        mul();
        for(i=0; i<n; i++)
        {
            printf("%d",c[i][0]%3);
            for(j=1; j<n; j++)
                printf(" %d",c[i][j]%3);
            printf("
");
        }
    }
    return 0;
}