社区讨论

Wa 7 9 11 12 13 80pts

P2455[SDOI2006] 线性方程组参与者 1已保存回复 0

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
0 条
当前快照
1 份
快照标识符
@mhjgwbhr
此快照首次捕获于
2025/11/04 02:22
4 个月前
此快照最后确认于
2025/11/04 02:22
4 个月前
查看原帖
CPP
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n;
double a[1010][1010];
void swaper(int i,int j){
	for(int k=1;k<=n+1;k++){
		swap(a[i][k],a[j][k]);
	}
}
void simpler(int i,double d){
	for(int k=1;k<=n+1;k++){
		a[i][k]=a[i][k]*1.0/d;
	}
}
void minuser(int i,double d,int j){
	for(int k=1;k<=n+1;k++){
		a[j][k]-=a[i][k]*1.0*d;
	}
}
int gause(){
	for(int i=1;i<=n;i++){
		//cout<<i;
		int mx=i;
		for(int j=1;j<i;j++){
			if(fabs(a[j][j])<1e-10&&fabs(a[j][i]>=1e-10)){
				if(fabs(a[j][i])>fabs(a[mx][i])){
					mx=j;
				}
			}
		}
		for(int j=i+1;j<=n;j++){
			if(fabs(a[j][i])>fabs(a[mx][i])){
				mx=j;
			}
		}
		//cout<<mx;
		swaper(i,mx);
		//cout<<'y';
		if(fabs(a[i][i])<1e-10){
			return 0;
		}
		simpler(i,a[i][i]);
		for(int j=1;j<=n;j++){
			if(j==i){
				continue;
			}
			minuser(i,a[j][i]/a[i][i],j);
		}
	}
	return 1;
}
int main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n+1;j++){
			cin>>a[i][j];
		}
	}
	if(!gause()){
		int ans=1;
		for(int i=n;i>=1;i--){
			if(fabs(a[i][i])<1e-10){
				int ok=0;
				for(int j=1;j<=n;j++){
					if(fabs(a[i][j])>=1e-10){
						ok=1;
						break;
					}
				}
				if(!ok&&fabs(a[i][n+1])>=1e-10){
					cout<<-1;
					return 0;
				}
			}
			
		}
		for(int i=n;i>=1;i--){
			if(fabs(a[i][i])<1e-10){
				int ok=0;
				for(int j=1;j<=n;j++){
					if(fabs(a[i][j])>=1e-10){
						ok=1;
						break;
					}
				}
				if(fabs(a[i][n+1])<1e-10){
					cout<<0;
					return 0;
				}
			}
		}
	}
	else{
		for(int i=1;i<=n;i++){
			printf("x%d=%.2lf\n",i,a[i][n+1]);
		}
	}
	return 0;
}

回复

0 条回复,欢迎继续交流。

正在加载回复...