弱菜的各种模板

Posted cx233666

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了弱菜的各种模板相关的知识,希望对你有一定的参考价值。

高斯-约当消元(这里以JSOI2008球形空间产生器)

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define gt getchar()
#define ll long long
#define db double
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
inline int in()
{
    int k=0;char ch=gt;
    while(ch<'-')ch=gt;
    while(ch>'-')k=k*10+ch-'0',ch=gt;
    return k;
}
db a[15][15];
int p[15],o[15];
const db eps=1e-3;
int main()
{
    int n=in()+1;
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<n;++j)
            scanf("%lf",a[i]+j),a[i][n+1]+=a[i][j]*a[i][j],a[i][j]*=2;
        a[i][n]=1;
    }   
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<=n;++j)if(!o[j]&&fabs(a[j][i])>eps){p[i]=j;break;}
        double t=a[p[i]][i];o[p[i]]=1;
        for(int j=1;j<=n+1;++j)a[p[i]][j]/=t;
        for(int k=1;k<=n;++k)
        {
            if(p[i]==k)continue;
            t=a[k][i];
            for(int j=1;j<=n+1;++j)a[k][j],a[k][j]-=a[p[i]][j]*t;
        }
    }
    for(int i=1;i<n;++i)printf("%.3lf ",a[p[i]][n+1]);puts("");
    return 0;
}

以上是关于弱菜的各种模板的主要内容,如果未能解决你的问题,请参考以下文章

2016-6-6

VSCode自定义代码片段——.vue文件的模板

VSCode自定义代码片段1——vue主模板

VSCode自定义代码片段2——.vue文件的模板

VSCode自定义代码片段(vue主模板)

Eclipse 中的通用代码片段或模板