莫欺少年穷,修仙之旅在这开始—>算法基础课题解
可参考: 高斯消元解线性方程组
#include<bits/stdc++.h>
using namespace std;
const int N = 110;
int n;
int a[N][N];
int gauss()
{
//枚举所有列
int r,c;
for(r=c=0;c<n;c++)
{
//找到当前列最大的行
int t=r;
for(int i=r;i<n;i++)
if(a[i][c])
t=i;
if(!a[t][c]) continue;
//与当前行交换
for(int i=c;i<=n;i++) swap(a[t][i],a[r][i]);
//让当前行后面的所有行的当前列变为 0
for(int i=r+1;i<n;i++)
if(a[i][c])
for(int j=c;j<=n;j++)
a[i][j]^=a[r][j];
r++;
}
if(r<n)
{
for(int i=r;i<n;i++)
if(a[i][n])
return 2;
return 1;
}
for(int i=n-1;i>=0;i--)
for(int j=i+1;j<n;j++)
a[i][n]^=a[i][j]*a[j][n];
return 0;
}
int main()
{
cin>>n;
for(int i=0;i<n;i++)
for(int j=0;j<=n;j++)
cin>>a[i][j];
int t=gauss();
//唯一解
if(!t) for(int i=0;i<n;i++) cout<<a[i][n]<<endl;
//无数解
else if(t==1) puts("Multiple sets of solutions");
//无解
else puts("No solution");
return 0;
}