10.13NOIP模拟赛

T1 魔界大门

题意：

小蛮，龙幽，姜云凡拜入蜀山门下，多日之后闲来无事，忽然想去魔界查访神农鼎的下落，当他们到魔界大门时，却被守门卫士执宿，疚业拦住，无奈之中只能开打，胜利之后卫士消失，魔界大门却迟迟没有动静，这是眼尖的小蛮发现门上有许多的数字，门的上方写着0/1，懂得魔族语言的龙幽反复观察门上的文字，发现门上给出的是三个N*N 的矩阵A,B,C，当A*B=C 他们只要把门上的1全部按下去就能打开大门，否则需要按下所有的0。四个人的数学都不好，所以请你帮助他们判断A*B 是否与C 相等。

提示：矩阵乘法+随机化

考试的时候memset写错导致丢掉暴力40分QAQ

算法1：暴力矩阵乘法 期望得分40

但是O(n^3)算法铁定超时，所以我们想一想，为什么要全算完呢？

算法2：随机化坐标 我们可以rand()矩阵中的坐标，只计算这么一个坐标，之后与目标矩阵中的对应点比较，相同继续rand，不同直接跳出，但是也不能一直算，所以期望得分70~80

既然可以随机化坐标，那为什么不能随机化矩阵呢？

算法3：随机化矩阵 我们一次rand出一个n*1的矩阵，之后乘，每次O(n^2)效率，随机2~3次即可解决，期望得分100

 

#include<cstdio>
#include<cstring>
#include<cstdlib>

using namespace std;

const int MAXN=1020;

unsigned int T,n,a[MAXN][MAXN],b[MAXN][MAXN],c[MAXN][MAXN],d[MAXN][MAXN];
unsigned int A[MAXN][MAXN],B[MAXN][MAXN],C[MAXN][MAXN],t[MAXN];

bool init() {
    memset(A,0,sizeof(A));
	memset(B,0,sizeof(B));
	memset(d,0,sizeof(d));
	int i,j,k;
	scanf("%d",&n);
	for(i=1;i<=n;++i)
	    for(j=1;j<=n;++j)
	        scanf("%d",&a[i][j]);
	for(i=1;i<=n;++i)
	    for(j=1;j<=n;++j)
	        scanf("%d",&b[i][j]);
	for(i=1;i<=n;++i)
	    for(j=1;j<=n;++j)
	        sca