HDU 4920 居然会超时

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

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

 1 #include <cstdio>
 2 #include <iostream>
 3 #include <sstream>
 4 #include <cmath>
 5 #include <cstring>
 6 #include <cstdlib>
 7 #include <string>
 8 #include <vector>
 9 #include <map>
10 #include <set>
11 #include <queue>
12 #include <stack>
13 #include <algorithm>
14 using namespace std;
15 #define ll long long
16 #define _cle(m, a) memset(m, a, sizeof(m))
17 #define repu(i, a, b) for(int i = a; i < b; i++)
18 #define repd(i, a, b) for(int i = b; i >= a; i--)
19 #define sfi(n) scanf("%d", &n)
20 #define pfi(n) printf("%d\n", n)
21 #define MAXN 100010
22 const int N = 810;
23 int n;
24 int a[N][N], b[N][N], c[N][N];
25 int scan()
26 {
27     int res = 0 , ch;
28     while( !( ( ch = getchar() ) >= ‘0‘ && ch <= ‘9‘ ) )
29     {
30         if( ch == EOF )  return 1 << 30 ;
31     }
32     res = ch - ‘0‘ ;
33     while( ( ch = getchar() ) >= ‘0‘ && ch <= ‘9‘ )
34         res = res * 10 + ( ch - ‘0‘ ) ;
35     return res ;
36 }
37
38 int main()
39 {
40     while(~scanf("%d", &n))
41     {
42         repu(i,0,n)
43         repu(j,0,n)
44         {
45             a[i][j] = scan() % 3;
46         }
47         repu(i,0,n)
48         repu(j,0,n)
49         {
50             b[i][j] = scan() % 3;
51         }
52         memset(c, 0, sizeof(c));
53         repu(i,0,n)
54         {
55             repu(j,0,n)
56             {
57                 repu(k,0,n)
58                 {
59                     //c[i][j] += a[i][k] * b[k][j];///会超时
60                     //c[i][k] += a[i][j] * b[j][k];///不会超时
61                     c[j][k] += a[j][i] * b[i][k];///不会超时
62                 }
63             }
64         }
65         for(int i = 0; i < n; i++)
66         {
67             for(int j = 0; j < n; j++)
68             {
69                 c[i][j] %= 3;
70                 if(j) printf(" %d", c[i][j]);
71                 else printf("%d", c[i][j]);
72             }
73             puts("");
74         }
75     }
76     return 0;
77 }

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

不过我还有一点疑惑就是为什么同样是三层循环（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];///不会超时