【算法】高斯消元

【题解】

建矩阵;

for i

  找到同列绝对值最大数字;

  交换;

  for k(行) j(列)(倒序)   除法;

for i(倒序)

  for j 减去已知元素

  除到右边;

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
const double eps=1e-6;
const int maxn=15;
double a[maxn][maxn],f[maxn];
int n;
void gauss()
{int r;for(int i=1;i<=n;i++){r=i;for(int j=i+1;j<=n;j++)if(fabs(a[j][i])>fabs(a[r][i]))r=j;if(r!=i)for(int j=1;j<=n+1;j++)swap(a[r][j],a[i][j]);for(int k=i+1;k<=n;k++)for(int j=n+1;j>=i;j--)a[k][j]-=a[k][i]/a[i][i]*a[i][j];      }for(int i=n;i>=1;i--){for(int j=i+1;j<=n;j++)a[i][n+1]-=a[j][n+1]*a[i][j];a[i][n+1]/=a[i][i];}}
int main()
{scanf("%d",&n);for(int i=1;i<=n;i++)scanf("%lf",&f[i]);double t;for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){scanf("%lf",&t);a[i][j]=2*(t-f[j]);a[i][n+1]+=t*t-f[j]*f[j];}}gauss();for(int i=1;i<n;i++)printf("%.3lf ",a[i][n+1]);printf("%.3lf",a[n][n+1]);return 0;
}

View Code

转载于:https://www.cnblogs.com/onioncyc/p/6611073.html

【BZOJ】1013 [JSOI2008]球形空间产生器sphere相关推荐

  1. BZOJ 1013: [JSOI2008]球形空间产生器sphere 高斯消元

    1013: [JSOI2008]球形空间产生器sphere Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/Judg ...

  2. BZOJ 1013: [JSOI2008]球形空间产生器sphere( 高斯消元 )

    可以得到N条N元一次方程, 高斯消元就OK了..时间复杂度O(N3) ----------------------------------------------------------------- ...

  3. bzoj 1013: [JSOI2008]球形空间产生器sphere(高斯消元)

    1013: [JSOI2008]球形空间产生器sphere Time Limit: 1 Sec  Memory Limit: 162 MB Submit: 5598  Solved: 2943 [Su ...

  4. BZOJ 1013 JSOI2008 球形空间产生器sphere 高斯消元

    题目大意:给定n维空间下的n+1个点,求这n个点所在的球面的球心 以前尝试了非常久的模拟退火0.0 至今仍未AC 0.0 今天挖粪涂墙怒学了高斯消元-- 我们设球心为X(x1,x2,...,xn) 如 ...

  5. 【BZOJ】1013: [JSOI2008]球形空间产生器sphere(高斯消元)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1013 只要列出方程组就能套高斯来解了. 显然距离相等,所以开不开平方都无所谓. b表示圆心,可列 s ...

  6. BZOJ.1013.[JSOI2008]球形空间产生器(高斯消元)

    题目链接 HDU3571 //824kb 40ms //HDU3571弱化版 跟那个一比这个太水了,练模板吧. //列出$n+1$个二次方程后两两相减,就都是一次方程了. #include <c ...

  7. 【高斯消元】[JSOI2008]球形空间产生器sphere

    [JSOI2008]球形空间产生器sphere 省选难度的板子 + 普及思维 = 紫题 (没毛病 日推天天推紫题再见了同学们我自请退群 为什么要我一个三维生物想象n维空间呢 以二维举例(糊一波我的手稿 ...

  8. bzoj1013 [JSOI2008]球形空间产生器sphere

    bzoj1013 [JSOI2008]球形空间产生器sphere 原题链接 题解 对\(k\in[2,n]\)列方程 \[\sum_{i=1}^{n}x_{i,0}(2x_{i,k}-2x_{i,k- ...

  9. [2018.12.18]BZOJ1013 [JSOI2008]球形空间产生器sphere

    设中心点坐标\(B=(b_1,b_2,b_3,...,b_n)\),球面上的点的坐标A=\((a_1,a_2,a_3,...,a_n)\) 则 \(dist_{A,B}=\sqrt{(a_1-b_1) ...

最新文章

  1. Python打包程序
  2. ftp服务器账户密码修改,ftp服务器如何设置账户密码
  3. Windows XP中防火墙后如何实现共享(图)
  4. 使用go来实现类似erlang otp里面的gen_server功能
  5. 制作五彩纸屑转场动效_何时以及如何将五彩纸屑添加到产品UI
  6. C#LeetCode刷题之#169-求众数(Majority Element)
  7. 使用 ale.js 制作一个小而美的表格编辑器(1)
  8. Spring Boot基础学习笔记08:Spring Boot整合Redis
  9. [独家放送]Unity2020规划预览,可视化编程又双叒叕来了!
  10. Git 历史记录内容对比
  11. The proxy server is refusing connections 怎么解决呢?
  12. Pandas高级教程之:时间处理
  13. 吃鸡 反作弊服务器未运行,绝地求生反作弊服务(BE)修复工具
  14. 从选课系统看软件开发周期
  15. 当年轻人在说“夸克真香”,他们在说什么?
  16. matlab波形转换,matlab波形图怎么转换为矢量图并导出?
  17. Oracle不走索引hint,oracle不走hint原因1:依据hint会出现错误结果
  18. php的设计模式和mvc,mvc设计模式是什么
  19. 计算机主板diy,DIY电脑配件立即选——主板篇
  20. MSSQL 服务器主体 无法在当前安全上下文下访问数据库,错误:916

热门文章

  1. 中医点滴 1 --- 冬季外感喉咙痛 寒性荨麻疹
  2. 【Kaidi安装问题】install_mkl.sh报错,没有数字签名
  3. 超级计算机 任务提交,vasp在超算中心的任务提交 - 第一原理 - 小木虫 - 学术 科研 互动社区...
  4. 昌邑机器人_上下料机器人昌邑机器人生产工厂
  5. 查看回调几个选项含义_C 盘总是莫名变大?更改这个文件位置至少腾出几个 G !...
  6. 如何在简历中添加自己的CSDN博客链接
  7. Linux 帧缓存 数据,linux--帧缓存技术
  8. linux 更改wp版本号,代码实现移除 WordPress 版本号
  9. 【RocketMQ工作原理】offset管理
  10. 【带你重拾Redis】Redis事务