【中山市选2008】矩阵

题目描述

给出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\sum_{k=1}^n=Ai,k*Bk,j$$

因为这个式子，我们不能用暴力了，只好另辟行径。
那么我们可以这么做，random出一个由01串D,组成的一维数组，然后把它和A矩阵处理成一个维的数组，然后把这个一维数组再一次与B矩阵处理，而C矩阵也与D处理一下，完成所有处理后，就变成了一个由A,B矩阵处理出的一维数组a和一个C矩阵处理出的一维数组b，这三次处理都只用了$n^2$,所以是可行的，最后再把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.