poj 3318 Matrix Multiplication

http://poj.org/problem?id=3318

矩阵A＊矩阵B是否等于矩阵C

#include <cstdio>
#include <cstring>
#include <time.h>
#include <algorithm>
#define maxn 1010
using namespace std;

int a[maxn][maxn],b[maxn][maxn],c[maxn][maxn],d[maxn];
int n;
int b1[maxn],c1[maxn],a1[maxn];

int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
            scanf("%d",&a[i][j]);
        }
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
            {
                scanf("%d",&b[i][j]);
            }
        }
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
            {
                scanf("%d",&c[i][j]);
            }
        }
        srand((unsigned int)time(0));
        for(int i=0; i<n; i++)
        {
            d[i]=rand()%100;
        }
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
            {
                b1[i]+=b[i][j]*d[j];
                c1[i]+=c[i][j]*d[j];
            }
        }
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
            {
                a1[i]+=a[i][j]*b1[j];
            }
        }
        bool flag=false;
        for(int i=0; i<n; i++)
        {
            if(a1[i]!=c1[i])
            {
                flag=true;
                break;
            }
        }
        if(!flag) printf("YES
");
        else printf("NO
");
    }
    return 0;
}