首页 > 其他 > 详细

gym100923C. Por Costel and Bujor (高斯消元)

时间:2019-07-24 09:35:04      阅读:94      评论:0      收藏:0      [点我收藏+]

题意:简化一下 就是解N个 系数矩阵一样 等式右边列矩阵不一样的方程组

题解:系数矩阵一样 为什么我却毫无办法????

   其实只要把等式右边的矩阵都排在后面就好了啊

   就变成解一个N x 2N的方程组了 ...

 

技术分享图片
#include <bits/stdc++.h>
using namespace std;
const double eps = 1e-9;

int n;
double a[205][405];

void gauss()
{
    int now = 1, to;
    for(int i = 1; i <= n; i++)
    {
        to = now;
        for(int j = now; j <= n; j++) {
            if(fabs(a[j][i]) > fabs(a[to][i])) to = j;
        }
        //if(to > n) continue;
        if(fabs(a[to][i]) < eps) continue;
        if(to != now)
            for(int j = 1; j <= n + n; j++) swap(a[to][j], a[now][j]);

        double tmp = a[now][i];
        //for(int j = 1; j <= n + n; j++) a[now][j] /= tmp;
        for(int j = 1; j <= n + n; j++)
            if(j != now)
            {
                tmp = a[j][i];
                for(int k = 1; k <= n + n; k++) a[j][k] -= tmp * a[now][k];
            }
        now++;
    }
}
int main() {
    //freopen("bujor.in","r",stdin);
    //freopen("bujor.out","w",stdout);

    int T;
    scanf("%d", &T);
    while(T--) {
        scanf("%d", &n);
        for(int i = 1; i <= n; i++)
        for(int j = 1; j <= n; j++) scanf("%lf", &a[i][j]);

        for(int j = n + 1; j <= n * 2; j++)
        for(int i = 1; i <= n; i++) {
            if(i + n == j) a[i][j] = 1;
            else a[i][j] = 0;
        }
        gauss();

        for(int i = 1; i <= n; i++) {
            for(int j = n + 1; j <= 2 * n; j++) {
                double tmp = a[i][j] / a[i][i];
                if(fabs(tmp) < eps) tmp = 0;
                if(j != 2 * n) printf("%.9lf ", tmp);
                else printf("%.9lf\n", tmp);
            }
        }
    }
    return 0;
}
View Code

 

gym100923C. Por Costel and Bujor (高斯消元)

原文:https://www.cnblogs.com/lwqq3/p/11235330.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!