码迷,mamicode.com
首页 > 其他好文 > 详细

HDU 4920 水

时间:2014-08-06 14:48:51      阅读:222      评论:0      收藏:0      [点我收藏+]

标签:style   blog   color   io   for   2014   ar   时间   

矩阵乘法

因为答案要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;
}




HDU 4920 水,布布扣,bubuko.com

HDU 4920 水

标签:style   blog   color   io   for   2014   ar   时间   

原文地址:http://blog.csdn.net/u011932355/article/details/38400471

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!