zoukankan      html  css  js  c++  java
  • 1013: [JSOI2008]球形空间产生器sphere

    #include<iostream>
    #include<cstdio>
    #include<cmath>
    #include<cstring>
    #include<algorithm>
    using namespace std;
    
    const int N=15;
    
    int n;
    double a[N][N];
    
    int main()
    {
        scanf("%d",&n);
        ++n;
        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+1;++j)
    //            printf("%.2lf ",a[i][j]);
    //        puts("");
    //    }
        for(int i=1;i<=n;++i)
        {
            int pos=i;
            for(int j=i;j<=n;++j)
            {
                if(fabs(a[pos][i])<fabs(a[j][i]))
                    pos=j;
            }
            if(pos!=i)
                swap(a[i],a[pos]);
            for(int j=i+1;j<=n;++j)
            {
                double ratio=a[j][i]/a[i][i];
                for(int k=i;k<=n+1;++k)
                {
                    a[j][k]-=a[i][k]*ratio;
                }
            }
        }
        for(int i=n;i;--i)
        {
            for(int j=i+1;j<=n;++j)
            {
                a[i][n+1]-=a[i][j]*a[j][n+1];
            }
            a[i][n+1]/=a[i][i];
        }
        for(int i=1;i<n;++i)
            printf("%.3lf ",a[i][n+1]);
        return 0;
    }
  • 相关阅读:
    反向映射和写时复制
    内存分析的好blog
    minicom使用
    tee的妙用
    网络带宽
    mem analyse
    linux 应用层常用API/命令
    ubuntu 库依赖问题
    Python基础学习笔记(一:hello world)
    第7章 取消与关闭
  • 原文地址:https://www.cnblogs.com/lovewhy/p/9633727.html
Copyright © 2011-2022 走看看