高斯消元模板

最新推荐文章于 2022-02-16 09:37:32 发布

扩展的灰

最新推荐文章于 2022-02-16 09:37:32 发布

阅读量273

点赞数

分类专栏：数学 ----线性代数 OI 文章标签：高斯消元扩展的灰

本文链接：https://blog.csdn.net/JacaJava/article/details/81873550

版权

OI 同时被 3 个专栏收录

321 篇文章 1 订阅

订阅专栏

数学

105 篇文章 0 订阅

订阅专栏

----线性代数

13 篇文章 0 订阅

订阅专栏

本文介绍了通过高斯消元方法解决模板题Bzoj1013，探讨在dp题目中应用的不便，并展示了如何建立球心坐标与距离方程之间的线性方程组，利用向量简化计算。

摘要由CSDN通过智能技术生成

高斯消元还是很好理解的，于是就找了个模板题做了下:Bzoj1013
不过好像用在那些dp题上面还是很不方便，感觉还是自己太菜了
这个题就很简单了，设出球心坐标 $O(O_1,O_2...O_n)$
写出距离方程式
$\sum_{i=1}^n(X_{0,i}-O_i)^2=\sum_{i=1}^n(X_{1,i}-O_i)^2=...=\sum_{i=1}^n(X_{n,i}-O_i)^2$
两边展开就得到了n个线性方程组，消元即可
偷懒用了向量，其实应该学一学stl的valarray的

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define db double
using namespace std;
int n; db x[20][20],k[20],A[20];
struct vec{
    db v[20];
    inline db& operator[] (int j){ return v[j]; }
    inline vec operator- (vec b){
        vec c;
        for(int i=0;i<20;++i) c[i]=v[i]-b[i];
        return c;
    }
    inline vec operator* (db j){
        vec c;
        for(int i=0;i<20;++i) c[i]=v[i]*j;
        return c;
    }
} s[20];
int main(){
    scanf("%d",&n);
    for(int i=0;i<=n;++i)
        for(int j=1;j<=n;++j) scanf("%lf",x[i]+j);
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n;++j){
            s[i][j]=2*(x[i][j]-x[0][j]);
            k[i]+=-x[0][j]*x[0][j]+x[i][j]*x[i][j];
        }
    }
    for(int i=1;i<=n;++i){
        for(int j=i;j<=n;++j) if(s[j][i]){ swap(s[j],s[i]); break; }
        for(int j=i+1;j<=n;++j){
            k[j]-=k[i]*(s[j][i]/s[i][i]); 
            s[j]=s[j]-s[i]*(s[j][i]/s[i][i]);
        }
    }
    for(int i=n;i;--i){
        for(int j=n;j>i;--j) k[i]-=A[j]*s[i][j];
        A[i]=k[i]/s[i][i];
    }
    for(int i=1;i<=n;++i) printf("%.3lf ",A[i]);
}