跟着大哥混系列。。
模板题: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个变元就可以得到全部解。