求助 关于高斯消元
  • 板块学术版
  • 楼主ccccc_jy
  • 当前回复3
  • 已保存回复3
  • 发布时间2021/12/23 20:55
  • 上次更新2023/10/28 13:49:11
查看原帖
求助 关于高斯消元
491375
ccccc_jy楼主2021/12/23 20:55

这是代码:

#include<bits/stdc++.h>
using namespace std;
double a[105][105];
int n;
int main(){
  cin>>n;
  for(int i=1;i<=n;i++){
     for(int j=1;j<=n+1;j++){
       scanf("%lf",&a[i][j]);
     }
  }
  for(int i=1;i<=n;i++){
   int flag=1;
   for(int j=i;j<=n;j++){
     if(a[j][i]!=0){
     flag=0;
     for(int k=1;k<=n+1;k++){
       swap(a[i][k],a[j][k]);
     }
    break;
    }
  }
    if(flag==1){
      printf("No Solution\n");
      return 0;
    }
    for(int j=1;j<=n;j++){ //行 
      if(j==i)continue;
      double times=a[i][i]/a[j][i];
      for(int k=i;k<=n+1;k++){
        a[j][k]=a[j][k]*times-a[i][k];
      }
    }
  }
  for(int i=1;i<=n;i++){
  printf("%.2f\n",a[i][n+1]/a[i][i]);
  }
  return 0;
}

第31行开始,我将

double times=a[j][i]/a[i][i];
	for(int k=i;k<=n+1;k++){
	 a[j][k]=a[j][k]-a[i][k]*times;
}//这是100分的

改成了第三十一行的代码代码,惨得十分,望知原因

2021/12/23 20:55
加载中...