2016.5.28【初中部 NOIP提高组】模拟赛C

Space_chicken

于 2016-05-28 15:59:41 发布

阅读量395

点赞数 2

分类专栏：成长总结

本文链接：https://blog.csdn.net/Space_chicken/article/details/51525309

版权

总结同时被 2 个专栏收录

8 篇文章 0 订阅

订阅专栏

成长

2 篇文章 0 订阅

订阅专栏

今天做题时，本来已经做得比较认真了，但还是做的不好。
T1:在做这题的时候，因为之前做过的一道类似的题，我甚至还想用异或可不可以做，但纠结了很久，然后又想了很多其他奇奇怪怪的算法，但没有一种是可行的……然后只好打了一个暴力。
T2:想完T1这题更让我一脸懵逼，花了半天终于读懂题后，我有不断在纸上模拟各种方法，但无济于事。到最后，比赛快结束是，我才发现原来打个搜索就可以拿分了，而最可气的是，比赛结束后，有同学的搜索还过了。
T3:在被T1,T2狂虐之后，看到T3的我，仿佛看到了曙光一般——终于可以AC一题了？一阵惊喜后，各种方法与程序实现方面的问题向我袭来，又把我推入了万丈深渊，比如说：如果有两个以上三角形重叠会怎样？……但最后我才发现这些其实都不用担心的。
无论怎样，吸取这次教训，下星期要更好发挥。

【中山市选2008】矩阵

题目描述

给出A,B,C三个01矩阵,问A*B是否等于C.

输入

本题有多组数据. 第一行的ca表示一共有ca组数据. 对于每组数据,第一行有一个整数N(1<=N<=1000) 然后就是A,B,C三个矩阵.

输出

对于每组数据,输出一行. 如果A*B=C,就输出“YES”,否则输出“NO”。

样例输入

111

000

100

010

001

111

000

样例输出

YES

这题，看到一串串00000，11111你也许会想到异或，但看完题目描述后，就会发现事实并不是这样。读完题后，首先我们应注意一个式子——

Ci,j∑k=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.