高斯消元总结

#include<bits/stdc++.h>
using namespace std;
int n;
double a[105][105];
bool f;
int main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n+1;j++){
			cin>>a[i][j];
		}
	}
	for(int i=1;i<=n;i++){
		int maxx=i;
		for(int j=i+1;j<=n;j++){
			if(fabs(a[maxx][i])<fabs(a[j][i]))maxx=j;
		}
		swap(a[i],a[maxx]);
		if(fabs(a[i][i])<1e-9)continue;
		double chu=a[i][i];
		for(int j=i;j<=n+1;j++){
			a[i][j]/=chu;
		}
		for(int j=1;j<=n;j++){
			if(i==j)continue;
			double chu=a[j][i]/a[i][i];
			for(int k=i;k<=n+1;k++){
				a[j][k]-=a[i][k]*chu;
			}
		}
//		for(int j=1;j<=n;j++){
//			for(int k=1;k<=n+1;k++){
//				cout<<a[j][k]<<" ";
//			}
//			cout<<endl;
//		}
//		cout<<endl;
	}
	for(int i=1;i<=n;i++){
		int cnt=0;
		for(int j=1;j<=n;j++){
			if(fabs(a[i][j])<1e-9)cnt++;
		}
		if(cnt==n){
			if(fabs(a[i][n+1])<1e-9){
				f=true;
			}
			else{
				cout<<-1;
				return 0;
			}
		}
	}
	if(f==true){
		cout<<0;
		return 0;
	}
	for(int i=1;i<=n;i++){
		printf("x%d=%.2lf\n",i,a[i][n+1]);
	}
	return 0;
} 

	#include<bits/stdc++.h>
using namespace std;
const int maxn=3005;
int n,m,ans;
char c[1005];
bitset<1005>a[maxn];
int main(){
	cin>>n>>m;
	for(int i=1;i<=m;i++){
		scanf("%s",c+1);
		for(int j=1;j<=n;j++){
			if(c[j]=='1')a[i][j]=1;
		}
		int x;
		cin>>x;
		a[i][n+1]=x;
	}
	for(int i=1;i<=n;i++){
		int pos=i;
		while(pos<=m&&(!a[pos][i]))pos++;
		if(pos>m){
			puts("Cannot Determine");
			return 0;
		}
		swap(a[i],a[pos]);
		ans=max(ans,pos);
		for(int j=1;j<=m;j++){
			if(i==j||(!a[j][i]))continue;
			a[j]^=a[i];
		}
//		if(i==n)cout<<pos<<endl;
	}
	cout<<ans<<endl;
	for(int i=1;i<=n;i++){
		if(a[i][n+1])puts("?y7M#");
		else puts("Earth");
	}
	return 0;
}