zoukankan      html  css  js  c++  java
  • 高斯消元

    #include <bits/stdc++.h>
    using namespace std;
    int n,pl;
    double a[1001][1001];
    int main() {
        cin>>n;
        for(int i=1;i<=n;i++) for(int j=1;j<=n+1;j++) cin>>a[i][j];
        for(int i=1;i<=n;i++) {
            pl=i;
            while(a[pl][i]==0&&pl<=n) pl++;                                    
            if(pl==n+1) { cout<<"No Solution"; return 0; }    
            for(int j=1;j<=n+1;j++) swap(a[i][j],a[pl][j]);
            double k=a[i][i];
            for(int j=1;j<=n+1;j++) a[i][j]=a[i][j]/k;
            for(int j=1;j<=n;j++) {
                if(i!=j){
                    double ki=a[j][i];
                    for(int m=1;m<=n+1;m++) a[j][m]=a[j][m]-ki*a[i][m];
                }
            }
        }
        for(int i=1;i<=n;i++) printf("%.2lf
    ",a[i][n+1]);///x1 x2 x3...
        return 0;
    }
  • 相关阅读:
    夜生活
    LOST
    TMD
    玩AJAX要先知道这些
    水漫金陵
    家装
    JAVA里面的小名堂
    朋友值多少钱
    猪八戒
    冬枣
  • 原文地址:https://www.cnblogs.com/iwomeng/p/11535158.html
Copyright © 2011-2022 走看看