A-Matrix Equation_第 45 届国际大学生程序设计竞赛（ICPC）亚洲区域赛（济南）（重现赛）

         自己写一个2维矩阵或者3维矩阵就可以发现对于每一列来说都是独立的，每一列的n个Cij都是都关系的，这就构成了一个n元一次方程组，其实这就是解一下这个方程组，但是他是问的有多少个矩阵，对于这个方程组构造出的矩阵来说，有多少自由元就说明有多少个Cij是可以任意取值的，也就是有1，2两种选择，不是自由元的就只能有1种；自由元是指系数都是0的未知数，自由元的个数就是高斯消元后0行的个数，所以用高斯消元就可以求出来；

        又因为有运算是模2的，所以也就可以看成是异或型矩阵，那么运算方式就要稍微的发生以下变化，比如

上述图片来自2020icpc-济南站_琅歌的博客-CSDN博客 

#include 
using namespace std;
#define endl '\n'
#define int long long
const int mod=998244353;
const int inf=1e18;
const int N = 2e5+100;
const double eps=1e-8;
int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int getinv(int a){return qpow(a,mod-2LL);}
int n;
int a[220][220],b[220][220],d[220][220],fac[N];
int gauss()
{int newline=1;for(int k=1;k<=n;k++){int p=newline;for(int i=newline+1;i<=n;i++)if(d[i][k]>d[p][k]) p=i;if(d[p][k]==0) continue;swap(d[newline],d[p]);for(int i=newline+1;i<=n;i++){if(d[i][k]==0) continue;for(int j=1;j<=n;j++)d[i][j]^=d[newline][j];}newline++;}newline--;return n-newline;
}
void print()
{cout<<"--------------------\n";for(int i=1;i<=n;i++){for(int k=1;k<=n;k++)cout<>n;fac[0]=1;for(int i=1;i<=100000;i++) fac[i]=fac[i-1]*2LL%mod;for(int i=1;i<=n;i++)for(int j=1;j<=n;j++) cin>>a[i][j],d[i][j]=a[i][j];for(int i=1;i<=n;i++)for(int j=1;j<=n;j++) cin>>b[i][j];int ans=1;for(int j=1;j<=n;j++){for(int i=1;i<=n;i++) for(int k=1;k<=n;k++) d[i][k]=a[i][k];for(int i=1;i<=n;i++)d[i][i]=a[i][i]^b[i][j];//print();int res=gauss();//print();//cout<

#include 
using namespace std;
#define endl '\n'
#define int long long
const int mod=998244353;
const int inf=1e18;
const int N = 2e5+100;
const double eps=1e-8;
int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int sgn(double x)
{if(fabs(x)>n;for(int i=1;i<=n;i++)for(int j=1;j<=n+1;j++) cin>>a[i][j];gauss();system("pause");return 0;
}

#include 
using namespace std;
#define endl '\n'
#define int long long
const int mod=998244353;
const int inf=1e18;
const int N = 4e5+100;
const double eps=1e-8;
int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int sgn(double x)
{if(fabs(x)>n>>m;for(int i=1;i<=n;i++) a[i][i]=a[i][n+1]=1;for(int i=1;i<=m;i++){int u,v;cin>>u>>v;a[u][v]=a[v][u]=1;}//print();gauss();ans=inf;//print();dfs(n,0);cout<

#include 
using namespace std;
#define endl '\n'
#define int long long
const int mod=998244353;
const int inf=1e18;
const int N = 4e5+100;
const double eps=1e-10;
int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int sgn(double x)
{if(fabs(x)v[105];
int gauss()
{int nl=1;for(int k=1;k<=n;k++){int p=nl;for(int i=nl+1;i<=n;i++) if(fabs(a[p][k])1) return 0;return ma;
}
void print()
{cout<<"-----------------------------\n";for(int i=1;i<=n;i++){for(int j=1;j<=n+1;j++)cout<>n;for(int i=1;i<=n+1;i++){int m;cin>>m;for(int j=1;j<=m;j++){int x;cin>>x;v[i].push_back(x);}cin>>b[i];}int ans=-1,mark=0;for(int i=1;i<=n+1;i++){for(int j=1;j<=n;j++)for(int k=1;k<=n+1;k++) a[j][k]=0;int cnt=0;for(int j=1;j<=n+1;j++){if(i==j)continue;cnt++;for(int k=0;k1){ans=-1;break;}}}if(ans==-1) cout<<"illegal\n";else cout<

#include 
using namespace std;
#define endl '\n'
#define int long long
const int mod=998244353;
const int inf=1e18;
const int N = 4e5+100;
const double eps=1e-10;
int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int sgn(double x)
{if(fabs(x)a[2005];
bool gauss()
{ans=0;int nl=1;for(int k=1;k<=n;k++){int p=nl;for(int i=nl;i<=m;i++) if(a[i][k]==1){p=i;break;}if(a[p][k]==0) return 0;ans=max(ans,p);swap(a[p],a[nl]);for(int i=1;i<=m;i++){if(i==nl||a[i][k]==0) continue;// for(int j=k;j<=n+1;j++)// a[i][j]^=a[nl][j];a[i]^=a[nl];}nl++;}nl--;return nl>=n;
}
void print()
{cout<<"-----------------------------\n";for(int i=1;i<=n;i++){for(int j=1;j<=n+1;j++)cout<>n>>m;for(int i=1;i<=m;i++){cin>>(s[i]+1);for(int j=1;j<=n;j++) a[i][j]=(s[i][j]-'0');int x;cin>>x;a[i][n+1]=x;}if(gauss()){cout<

#include 
using namespace std;
#define endl '\n'
#define int long long
const int mod=998244353;
const int inf=1e18;
const int N = 4e5+100;
const double eps=1e-10;
int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int sgn(double x)
{if(fabs(x)>n>>m>>p>>q;P=p/q;for(int i=1;i<=m;i++){int u,v;cin>>u>>v;in[u]++;in[v]++;e[u][v]=e[v][u]=1;}for(int i=1;i<=n;i++) a[i][i]=1;for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){if(e[i][j]) a[i][j]-=a[j][j]*(1.0-P)/in[j];}}a[1][n+1]=1;gauss();system("pause");return 0;
}