HDU 4920 居然会超时

题意：求两个n*n的矩阵相乘的结果，得出的每个元素%3；

分析：2000ms然后n的范围是800，我们自己估算的时间复杂度并不会超时，但是结果就是超时了。

 #include <cstdio>
 #include <iostream>
 #include <sstream>
 #include <cmath>
 #include <cstring>
 #include <cstdlib>
 #include <string>
 #include <vector>
 #include <map>
 #include <set>
 #include <queue>
 #include <stack>
 #include <algorithm>
 using namespace std;
 #define ll long long
 #define _cle(m, a) memset(m, a, sizeof(m))
 #define repu(i, a, b) for(int i = a; i < b; i++)
 #define repd(i, a, b) for(int i = b; i >= a; i--)
 #define sfi(n) scanf("%d", &n)
 #define pfi(n) printf("%d\n", n)
 #define MAXN 100010
 const int N = ;
 int n;
 int a[N][N], b[N][N], c[N][N];
 int scan()
 {
     int res =  , ch;
     while( !( ( ch = getchar() ) >= '' && ch <= '' ) )
     {
         if( ch == EOF )  return  <<  ;
     }
     res = ch - '' ;
     while( ( ch = getchar() ) >= '' && ch <= '' )
         res = res *  + ( ch - '' ) ;
     return res ;
 }

 int main()
 {
     while(~scanf("%d", &n))
     {
         repu(i,,n)
         repu(j,,n)
         {
             a[i][j] = scan() % ;
         }
         repu(i,,n)
         repu(j,,n)
         {
             b[i][j] = scan() % ;
         }
         memset(c, , sizeof(c));
         repu(i,,n)
         {
             repu(j,,n)
             {
                 repu(k,,n)
                 {
                     //c[i][j] += a[i][k] * b[k][j];///会超时
                     //c[i][k] += a[i][j] * b[j][k];///不会超时
                     c[j][k] += a[j][i] * b[i][k];///不会超时
                 }
             }
         }
         for(int i = ; i < n; i++)
         {
             for(int j = ; j < n; j++)
             {
                 c[i][j] %= ;
                 if(j) printf(" %d", c[i][j]);
                 else printf("%d", c[i][j]);
             }
             puts("");
         }
     }
     return ;
 }

我当时找稀疏矩阵的时候就找到的是直接矩阵相乘的模板，但是还是超时了，后来才知道我的输入挂需要重新更新了。

不过我还有一点疑惑就是为什么同样是三层循环（i，j，k），但是跑出来的时间不一样的呢。。。。。。。

c[i][j] += a[i][k] * b[k][j];///会超时
c[i][k] += a[i][j] * b[j][k];///不会超时
c[j][k] += a[j][i] * b[i][k];///不会超时