题目描述
给出A,B,C三个01矩阵,问A*B是否等于C.
原题:
输入
本题有多组数据. 第一行的ca表示一共有ca组数据. 对于每组数据,第一行有一个整数N(1<=N<=1000) 然后就是A,B,C三个矩阵.
原:
输出
对于每组数据,输出一行. 如果A*B=C,就输出“YES”,否则输出“NO”。
原:
样例输入
2
2
10
01
11
10
11
01
3
111
111
000
100
010
001
111
111
000
样例输出
NO
YES
这题,看到一串串00000,11111你也许会想到异或,但看完题目描述后,就会发现事实并不是这样。读完题后,首先我们应注意一个式子——
Ci,j∑k=1n=Ai,k∗Bk,j
因为这个式子,我们不能用暴力了,只好另辟行径。
那么我们可以这么做,random出一个由01串D,组成的一维数组,然后把它和A矩阵处理成一个维的数组,然后把这个一维数组再一次与B矩阵处理,而C矩阵也与D处理一下,完成所有处理后,就变成了一个由A,B矩阵处理出的一维数组a和一个C矩阵处理出的一维数组b,这三次处理都只用了
n2
,所以是可行的,最后再把a,b比较一下即可。注意要进行多次random,以免出现特殊情况。
总体来说,方法就是:
D:=random(2);
A*B*D是否=C*D
这题其实算是一道概率题,因为random,最多可能只需4到5次就可以求出答案。
欲知详情,请点击这里!!!
另赠代码一份
var
bz:boolean;
a,b,c:array[1..1000,1..1000] of longint;
d,s,x,ans:array[1..1000] of longint;
o:char;
i,ii,j,k,t,n:longint;
procedure sample;
var
i:longint;
begin
for i:=1 to n do d[i]:=random(2);
end;
begin
readln(t);
randomize;
for ii:=1 to t do
begin
readln(n);
bz:=true;
for i:=1 to n do
begin
for j:=1 to n do
begin
read(o);
val(o,a[i,j]);
end;
readln;
end;
for i:=1 to n do
begin
for j:=1 to n do
begin
read(o);
val(o,b[i,j]);
end;
readln;
end;
for i:=1 to n do
begin
for j:=1 to n do
begin
read(o);
val(o,c[i,j]);
end;
readln;
end;
k:=0;
while (bz=true) and (k<4) do
begin
inc(k);
sample;
fillchar(s,sizeof(s),0);
fillchar(x,sizeof(x),0);
fillchar(ans,sizeof(ans),0);
for i:=1 to n do
for j:=1 to n do
s[i]:=(s[i]+a[j,i]*d[j]) mod 2;
for i:=1 to n do
for j:=1 to n do
x[i]:=(x[i]+b[j,i]*s[j]) mod 2;
for i:=1 to n do
for j:=1 to n do
ans[i]:=(ans[i]+c[j,i]*d[j]) mod 2;
for i:=1 to n do if ans[i]<>x[i] then bz:=false;
end;
if bz=true then writeln('YES') else writeln('NO');
end;
end.