高斯消元法

跟着大哥混系列。。

模板题：luogu3389

//普通高斯消元法模板
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
using namespace std;
#define eps 1e-6
int n;
double a[110][110];
double ans[110];
int r;
void swap(int x,int y)
{
	double t;
	for (int i=1;i<=n+1;i++)
	{
		t = a[x][i];a[x][i] = a[y][i];a[y][i] = t;
	}
}
void Debug()
{
	for (int i=1;i<=n;i++)
	{
		for(int j=1;j<=n+1;j++)
			printf("%.2f ", a[i][j]);
		printf("\n");
	}
	cout<<endl;
}
int Gauss()//返回0表示无解，否则返回增广矩阵的秩
{
	//消元
	r = n;//rank of A'
	double t;
	for (int i = 1;i<=n;i++)
	{
		int pos = i;
		for (int j=i+1;j<=n;j++)
			if (fabs(a[pos][i]) < fabs(a[j][i])) pos = j;//找出系数最大的，一，便于判断是否全部为0，二，让系数尽量小，防止爆double
		if (fabs(a[pos][i]) < eps) 
		{
			for (int j=i;j<=n;j++)
				if (fabs(a[j][n+1]) >= eps) return 0;//无解
			r--;continue;//rank减小
		}
		if (pos != i) swap(pos,i);
		for (int j=i+1;j<=n;j++)
			if (fabs(a[j][i]) >= eps)
			{
				t = a[j][i] / a[i][i];
				for (int k=i;k<=n+1;k++)
					a[j][k] -= t * a[i][k];
			}

		//Debug();

	}
	//回代
	for (int i=n;i>=1;i--)
	{
		ans[i] = a[i][n+1] / a[i][i];
		for (int j = i-1;j>=1;j--)
			a[j][n+1] -= ans[i] * a[j][i];
	}
	return r;//返回增广矩阵的秩

}
int main()
{
	scanf("%d", &n);
	for (int i=1;i<=n;i++)
		for (int j=1;j<=n+1;j++)
			scanf("%lf", &a[i][j]);
	if (Gauss() != n)
		printf("No Solution\n");
	else 
		for (int i=1;i<=n;i++)
			printf("%.2f\n", ans[i]);
	return 0; 
}

poj1222

棋盘，与每个点有关的只有上下左右

//高斯消元解异或方程组
//poj1222
//注意建立方程组的时候，对于越过边界的点一定要特殊处理，这个和普通的矩阵不一样，会多算出其他的变元，自己真的没看出来
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
using namespace std;
int T,n,res,T1;
int ans[35];
int a[35][35];
void swap(int x,int y)
{
	int t;
	for (int i=1;i<=n+1;i++)
	{
		t = a[x][i];a[x][i] = a[y][i];a[y][i] = t;
	}
}
void Gauss()
{
	//n=30
		//消元
		for (int i=1;i<=n;i++)
		{
			int pos = i;
			while (pos <= n && a[pos][i] != 1) pos++;
			if (pos != i) swap(i,pos);
			for (int j = i+1;j <= n;j++)
				if (a[j][i] == 1)
					for (int k = i;k<=n+1;k++)
						a[j][k] ^= a[i][k];
		}
		//回代
		for (int i = n;i>=1;i--)
		{
			ans[i] = a[i][n+1];
			for (int j = i-1;j>=1;j--)
				if (a[j][i]) 
				{
					a[j][i] ^= a[i][i];
					a[j][n+1] ^= a[i][n+1];
				}
		}

}
int main()
{
	scanf("%d", &T);
	T1 = T;
	while (T--)
	{
		cout<<"PUZZLE #"<<T1-T<<endl;
		memset(a,0,sizeof(a));
		//构造方程
		n = 0;
		for (int i=1;i<=6;i++)
		{
			n++;
			a[n][i] = a[n][i-1] = a[n][i+1] = a[n][i+6] = 1;
			if (i == 1) a[n][i-1] = 0;
			if (i == 6) a[n][i+1] = 0;
		}

		for (int i=2;i<=5;i++)
			for (int j=1;j<=6;j++)
			{
				n++;
				a[n][(i-1)*6+j] = a[n][(i-2)*6+j] = a[n][(i-1)*6+j-1] = a[n][i*6+j] = a[n][(i-1)*6+j+1] = 1;
				if (j == 1) a[n][(i-1) * 6 + j - 1] = 0;
				if (j == 6) a[n][(i-1) * 6 + j + 1] = 0;
			}
		
		for (int i=1;i<=5;i++)
			for (int j=1;j<=6;j++)
				scanf("%d", &a[(i-1)*6 + j][n+1]);//初始状态


		Gauss();
		//输出答案
		for (int i=1;i<=5;i++)
		{
			for (int j=1;j<=6;j++)
				cout<<ans[(i-1) * 6 + j]<<' ';
			cout<<endl;
		}
	}

}

poj3185

对于高斯消元法的解法，题目要求最小的翻碗的次数，就要在所有的解中找到次数最小的。-->枚举自由变元

任意n-r个变元都可以是自由变元，每个自由变元都有0/1两种取值，因此只要随便枚举其中n-r个变元就可以得到全部解。