HDU 4920 水

矩阵乘法

因为答案要MOD3，所以矩阵中会有很多值为0，对这些不乘就行了，，，，，，，这样也能水过。。。

BUT ：

这样写会超时：
            for (int i=1; i<=n; i++)
                for (int j=1; j<=n; j++)
                    for (int k=1; k<=n; k++)
                        c[i][j]+=a[i][k]*b[k][j];

这样写就能过:

        for (int k=1; k<=n; k++)
            for (int i=1; i<=n; i++)
                for (int j=1; j<=n; j++)
                    c[i][j]+=a[i][k]*b[k][j];

BECAUSE：

我们知道内存中二维数组是以行为单位连续存储的，逐列访问将会每次跳1000*4(bytes)。根据cpu cache的替换策略，将会有大量的cache失效。
时间居然会相差很多。 可见利用好cpu cache优化我们的程序，是非常有必要掌握的技能。
平时写程序时，也应当尽量使cpu对内存的访问，是尽可能连续的。

#include "stdio.h"
#include "string.h"
int a[810][810],b[810][810],c[810][810];
int main()
{
    int n,i,j,k;

    while (scanf("%d",&n)!=EOF)
    {
        for (i=1;i<=n;i++)
            for (j=1;j<=n;j++)
            {
                scanf("%d",&a[i][j]);
                a[i][j]%=3;
            }

        for (i=1;i<=n;i++)
            for (j=1;j<=n;j++)
            {
                scanf("%d",&b[i][j]);
                b[i][j]%=3;
            }

        memset(c,0,sizeof(c));

        for (k=1;k<=n;k++)
            for (i=1;i<=n;i++)
            if (a[i][k]!=0)
                for (j=1;j<=n;j++)
                {
                    c[i][j]+=a[i][k]*b[k][j]%3;

                }
        for (i=1;i<=n;i++)
        {
            for (j=1;j<n;j++)
                printf("%d ",c[i][j]%3);
            printf("%d\n",c[i][j]%3);
        }


    }
    return 0;
}