12
2
2015
0

【vijos1052】【高斯消元】贾老二算算术

描述

贾老二是个品学兼优的好学生,但由于智商问题,算术学得不是很好,尤其是在解方程这个方面。虽然他解决 2x=2 这样的方程游刃有余,但是对于 {x+y=3 x-y=1} 这样的方程组就束手无策了。于是他要你来帮忙。前提是一次方程组且保证在integer的范围内可以处理所有问题。

格式

输入格式

第一行一个数字N(1≤N≤100)表示要求的未知数的个数,同时也是所给的方程个数。

第2到N+1行,每行N+1个数。前N个表示第1到N个未知数的系数。第N+1个数表示N个未知数乘以各自系数后的加和。(保证有唯一整数解)

输出格式

一行N个数,表示第1到N个未知数的值。

样例1

样例输入1[复制]

 
2
1 1 3
1 –1 1

样例输出1[复制]

 
2 1

来源

ASYZOI贾老二专题


高斯消元模板题~

用一节课重新理解高斯消元,重新写高斯消元模板~

果然就算品学兼优也无法弥补智商的缺陷么T_T?NO NO NO

这道题 看题+打代码+提交+修改+AC=25min~

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#define eps 1e-3
#define for1(i,n) for(int i=1;i<=(n);i++)
using namespace std;
const int M=100+5;
int n;
double a[M][M];
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
    return x*f;
}
void Gauss_elimination()
{
	int r;
	for(int i=1;i<=n;i++)
	{
		r=i;
		for(int j=i+1;j<=n;j++)
			if(fabs(a[j][i])>fabs(a[r][i])) r=j;
		if(r!=i) for(int j=1;j<=n+1;j++) swap(a[r][j],a[i][j]);
		
		for(int k=i+1;k<=n;k++)
		{
			double f=a[k][i]/a[i][i];
			for(int j=i;j<=n+1;j++)
				a[k][j]-=f*a[i][j];
		}
	}
	for(int i=n;i>=1;i--)
	{
		for(int j=i+1;j<=n;j++)
		  a[i][n+1]-=a[j][n+1]*a[i][j];
		a[i][n+1]/=a[i][i];
	}
}
int main()
{
	n=read();
	memset(a,0,sizeof(a));
	for1(i,n) for1(j,n+1) scanf("%lf",&a[i][j]);
	Gauss_elimination();
	for1(i,n)
	{
		printf("%d",(int)(a[i][n+1]+eps));
		if(i!=n) printf(" ");
	}
	return 0;
}

ATTENTION!——

double型数组、f

增广矩阵 注意+1

Category: 数学相关 | Tags: | Read Count: 292

登录 *


loading captcha image...
(输入验证码)
or Ctrl+Enter

Host by is-Programmer.com | Power by Chito 1.3.3 beta | Theme: Aeros 2.0 by TheBuckmaker.com