bzoj2396: 神奇的矩阵（矩阵乘法+随机化）

　　这题n三方显然会GG...

　　运用矩阵乘法的性质A*B*R=A*(B*R)=C*R，于是随机化出一个一列的R，就可以把复杂度降低成n方...大概率是不会错的

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=,inf=1e9;
typedef int mtx[maxn][maxn];
int n,m,x,y,z,tot;
mtx f,a,b,c;
void read(int &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
void mul(mtx &a,mtx b)
{
    mtx c;memset(c,,sizeof(c));
    for(int i=;i<=n;i++)
    for(int j=;j<=;j++)
    for(int k=;k<=n;k++)
    c[i][j]+=a[i][k]*b[k][j];
    memcpy(a,c,sizeof(c));
}
int main()
{
    srand();
    for(int i=;i<=;i++)f[i][]=+rand()%;
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
        read(a[i][j]);
        for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
        read(b[i][j]);
        for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
        read(c[i][j]);
        mul(b,f);mul(c,f);mul(a,b);
        int flag=;
        for(int i=;i<=n;i++)
        if(a[i][]!=c[i][])
        {
            flag=;
            break;
        }
        printf("%s\n",flag?"No":"Yes");
    }
    return ;
}