题目链接

  水比题,把圆方程展开减一下把平方都减掉半径的平方也减掉,高斯消元即可。

  然后我只输出两位小数,爆了两次零。我好菜啊。

  

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cctype>
#include<cstring>
#include<cmath>
#define maxn 20
#define eps 1e-9
using namespace std; double s[maxn*maxn][maxn];
double q[maxn][maxn]; double ans[maxn]; int main(){
int n;
scanf("%d",&n);
for(int i=;i<=n+;++i)
for(int j=;j<=n;++j) scanf("%lf",&q[i][j]);
int m=;
for(int i=;i<=n;++i)
for(int j=i+;j<=n+;++j){
double *c=s[++m];
for(int k=;k<=n;++k){
c[k]=*(q[i][k]-q[j][k]);
c[n+]+=q[i][k]*q[i][k]-q[j][k]*q[j][k];
}
}
//for(int i=1;i<=m;++i,printf("\n"))
// for(int j=1;j<=n+1;++j) printf("%.2lf ",s[i][j]);
for(int i=;i<=n;++i){
int now=i;
for(int j=i+;j<=m;++j)
if(fabs(s[j][i])-fabs(s[now][i])>eps) now=j;
if(now^i) swap(s[i],s[now]);
double ret=s[i][i];
for(int j=i;j<=n+;++j) s[i][j]=s[i][j]/ret;
for(int j=i+;j<=m;++j){
ret=s[j][i];
for(int k=;k<=n+;++k) s[j][k]=s[j][k]-s[i][k]*ret;
}
}
ans[n]=s[n][n+];
//for(int i=1;i<=m;++i,printf("\n"))
// for(int j=1;j<=n+1;++j) printf("%.2lf ",s[i][j]);
for(int i=n-;i;--i){
ans[i]=s[i][n+];
for(int j=i+;j<=n;++j) ans[i]=ans[i]-ans[j]*s[i][j];
}
for(int i=;i<=n;++i) printf("%.3lf ",ans[i]);
return ;
}

最新文章

  1. svn 常用命令总结
  2. Postgresql 数据库维护
  3. 12.我们不是在真空里谈软件工程, 软件要运行在硬件芯片上面, 下面看看一个计算机芯片的发展历史: http://perspectives.mvdirona.com/2014/09/august-21-2014-computer-history-museum-presentation/ http://mvdirona.com/jrh/TalksAndPapers/DileepBhandar
  4. 裴波那契数列 JavaScript 尾递归实现
  5. HDU 5044 (树链剖分+树状数组+点/边改查)
  6. C插入排序
  7. Java网络编程(UDP协议-聊天程序)
  8. 通过判断浏览器的userAgent,用正则来判断是否是ios和Android客户端
  9. 第13章 模版方法模式(Template Method)
  10. 第二篇、vlc-android之源码介绍
  11. 20165225《Java程序设计》第七周学习总结
  12. 给1~3年iOS开发 经验朋友们的一些建议(附BAT面试题)
  13. xcode 5.1打包iOS 7.1应用问题笔记
  14. 【thinkphp5】使用tp5开发api接口 定义全局异常处理
  15. 关于IOS给我的启发
  16. linux可视化桌面安装
  17. website for personal research
  18. yii2美化url
  19. isX字符串方法
  20. C# SuperSocket 消息推送

热门文章

  1. Java 原型模式(克隆模式)
  2. idea spring boot启动项目上面有红色叉
  3. 在ubuntu中用命令安装iperf
  4. Linux学习日记:第一天
  5. linux更新git
  6. 如何挂载一个镜像文件(how to mount an image file)
  7. python生成四位随机数
  8. jmeter操作mysql数据库
  9. hive sql 学习笔记
  10. Python中类的声明,使用,属性,实例属性,计算属性及继承,重写