洛谷P2455 线性方程组做题记录

 P2455 [SDOI2006] 线性方程组 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

这里只谈心得，下帖正解代码

#include <iostream>
#include <cmath>
#include <cstdio>
using namespace std;
double a[100][100];
int n;
bool eq(double x,double y){
	return fabs(x-y)<1e-9;
}
void see(){
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n+1;j++){
			cout<<a[i][j]<<" ";
		}
		cout<<endl;
	}
	cout<<endl;
}
int p=0;
int main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n+1;j++) cin>>a[i][j];
	}
	int nwline=1;
	for(int k=1;k<=n;k++){
		int maxi=nwline;
		for(int i=nwline+1;i<=n;i++){
			if(fabs(a[i][k]>fabs(a[maxi][k]))) maxi=i;
		}
		if(eq(a[maxi][k],0)){
				//see();
//			swap(a[maxi],a[n+1-++p]);
			continue;
		}
		swap(a[nwline],a[maxi]);
		
		for(int i=1;i<=n;i++){
			if(i==nwline) continue;
			if(!eq(a[nwline][k],0)){
    				double div=a[i][k]/a[nwline][k];	
    			for(int j=k;j<=n+1;j++){
    				a[i][j]-=a[nwline][j]*div;
    			}
			}
		}
		nwline++;
	}
	if(nwline<=n){
		while(nwline<=n){
			if(!eq(a[nwline][n+1],0)){
				cout<<-1;
				return 0;
			}
			nwline++;
			}
		cout<<0<<endl;	
	}else{
		for(int i=1;i<=n;i++)
         if(eq(a[i][n+1]/a[i][i],0)) cout<<"x"<<i<<"=0.00"<<endl;
			else printf("x%d=%.4lf\n", i, a[i][n+1]/a[i][i]);
	}
	return 0;
}

在网上找了很多题解发现已经被评论区hack掉了，最终写出这一个对的。这题的难点在于如何判断无穷解和有解。

回顾一下手算线性方程组的过程：把行列式化为行简化阶梯式，然后看是否出现了美丽的阶梯型，如果有一行全是零（出现了线性相关），且其他列正常，就是无穷解。如果有一行前面都是0，最后一个数字不是0，那么出现0=x的情况，也就是无解。

大多数被hack掉的原因是直接在循环中，把循环变量当成了当前处理行，导致当出现第一列都是0这样的数据时，第一行不会被处理，导致没有化为行简化阶梯型。

还有一个易错点就是误差传播：

浮点数小数据除以大数据比大数据除以小数据误差更小：

**大数据除以小数据**

当大数据除以小数据时，结果通常会得到一个非常大的值。这是因为除法运算本质上是将被除数（大数据）分成相等的份数（除数）。当除数很小时，份数就会变得非常多，导致结果非常大。

例如，如果我们用 1000000 除以 0.0001，结果将是 10000000000。

**小数据除以大数据**

当小数据除以大数据时，结果通常会得到一个非常小的值。这是因为除法运算本质上是将被除数（小数据）分成相等的份数（除数）。当除数很大时，份数就会变得非常小，导致结果非常小。

例如，如果我们用 0.0001 除以 1000000，结果将是 0.0000001。

**浮点数计算中的注意事项**

在浮点数计算中，当除以非常小的数字时，可能会出现精度问题。这是因为浮点数只能表示有限精度的数字，当除数非常小时，结果可能会四舍五入为零。

为了避免这种情况，可以使用以下技巧：

* 使用双精度或更高精度的浮点数类型。
* 在除法之前将除数乘以一个常数，使其变大。
* 使用舍入到最近偶数的舍入模式，以减少舍入误差。

**示例代码（C++）：**

```cpp
#include <iostream>

using namespace std;

int main() {
    double big_num = 1000000;
    double small_num = 0.0001;

    // 大数据除以小数据
    double result1 = big_num / small_num;
    cout << "大数据除以小数据：" << result1 << endl;

    // 小数据除以大数据
    double result2 = small_num / big_num;
    cout << "小数据除以大数据：" << result2 << endl;

    return 0;
}
```

**输出：**

```
大数据除以小数据：10000000000
小数据除以大数据：0.0000001
```

最终导致结果差之毫厘，失之千里：

改良版本如代码展示：要把当前列最大数字所在行与当前处理行交换一下，这样永远都是小数字除以大数字。