bzoj2396 神奇的矩阵（随机化）

Time Limit: 5 Sec  Memory Limit: 512 MB

给出三个行数和列数均为N的矩阵A、B、C，判断A*B=C是否成立。

题目可能包含若干组数据。
    对于每组数据，第一行一个数N，接下来给出三个N*N的矩阵，依次为A、B、C三个矩阵。

对于每组数据，若A*B=C成立，则输出Yes，否则No。每个答案占一行。

Sample Input

1
2
2
100

Sample Output

No

HINT

对于90%的数据，N不超过100；

对于100%的数据，N不超过1000，矩阵中的数字大于等于0小于1000，数据组数不超过5组。

---

poj原题......就是改了下数据范围

详见：blog：poj3318 Matrix Multiplication

根据这题的数据范围$O(n^2),n<=1000$，每组大概判断17次，错误率$1/2^{17}$

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#define rint register int
using namespace std;
#define N 1005
int i,j,k,n,a[N][N],b[N][N],c[N][N],d[N],e[N],f[N];
int main(){
    srand();
    while(scanf("%d",&n)!=EOF){
        for(i=;i<=n;++i)
            for(j=;j<=n;++j)
                scanf("%d",&a[i][j]);
        for(i=;i<=n;++i)
            for(j=;j<=n;++j)
                scanf("%d",&b[i][j]);
        for(i=;i<=n;++i)
            for(j=;j<=n;++j)
                scanf("%d",&c[i][j]);
        for(k=;k<=;++k){
            for(i=;i<=n;++i) d[i]=rand()&;
            for(e[i=]=;i<=n;e[++i]=)
                for(j=;j<=n;++j)
                    e[i]+=c[i][j]*d[j];
            for(f[i=]=;i<=n;f[++i]=)
                for(j=;j<=n;++j)
                    f[i]+=b[i][j]*d[j];
            for(i=;i<=n;++i) d[i]=f[i];
            for(f[i=]=;i<=n;f[++i]=)
                for(j=;j<=n;++j)
                    f[i]+=a[i][j]*d[j];
            for(i=;e[i]==f[i]&&i<=n;++i);
            if(i<=n) break;
        }puts(k>?"Yes":"No");
    }return ;
}