P4035 [JSOI2008]球形空间产生器

最新推荐文章于 2023-03-12 16:03:07 发布

原创最新推荐文章于 2023-03-12 16:03:07 发布 · 197 阅读

0 ·

CC 4.0 BY-SA版权

数学专栏收录该内容

7 篇文章

订阅专栏

本文介绍了一个使用C++实现的高斯消元法求解线性方程组的代码示例。该算法通过矩阵操作逐步消元未知数，最终求得方程组的解。代码中详细展示了如何初始化系数矩阵和常数项，执行前向消元和后向替代的过程。

之前自己写过的PDF版的题解：

AC Code：

#include<cstdio>
#include<iostream>
#define ri register int
using namespace std;

const int MAXN=15;
int n;
double fst[MAXN][MAXN],k[MAXN],a[MAXN][MAXN],b[MAXN],x[MAXN];

void Gauss()
{
    for(ri i=1;i<=n;i++)
        for(ri j=i+1;j<=n;j++)
        {
            double chu=a[j][i]/a[i][i];
            for(ri k=i;k<=n;k++)	a[j][k]-=a[i][k]*chu;
            b[j]-=chu*b[i];
        }
    for(ri i=n;i>=1;i--)
    {
        double ans=0;
        for(ri j=i+1;j<=n;j++)	 ans+=x[j]*a[i][j];
        x[i]=(b[i]-ans)/a[i][i];
    }
}

int main()
{
    scanf("%d",&n);
    for(ri i=1;i<=n+1;i++)
        for(ri j=1;j<=n;j++)	scanf("%lf",&fst[i][j]);
    for(ri i=1;i<=n+1;i++)
        for(ri j=1;j<=n;j++)	k[i]+=fst[i][j]*fst[i][j];
    for(ri i=1;i<=n;i++)
    {
        for(ri j=1;j<=n;j++)	a[i][j]=(fst[i][j]-fst[i+1][j])*2;
        b[i]=k[i]-k[i+1];
    }
    Gauss();
    for(ri i=1;i<=n;i++)	printf("%.3f ",x[i]);
    return 0;
}