BZOJ_1013_[JSOI2008]_球形空间产生器_(高斯消元)
描述
http://www.lydsy.com/JudgeOnline/problem.php?id=1013
n维空间,给出球上n+1个点的n维坐标,求球心坐标.
提示:给出两个定义:1、 球心:到球面上任意一点距离都相等的点。2、 距离:设两个n为空间上的点A, B
的坐标为(a1, a2, …, an), (b1, b2, …, bn),则AB的距离定义为:dist = sqrt( (a1-b1)^2 + (a2-b2)^2 +
… + (an-bn)^2 )
分析
对于前n个点,第i个点到球心的距离与第i+1个相同,可以列出形式如下的n个方程(x表示球心):
$$\sum_{j=1}^n(a[i][j]^2-a[x][j]^2)=\sum_{j=1}^n(a[i+1][j]^2-a[x][j]^2)$$
我们把它展开,化简,得到:
$$\sum_{j=1}^n2a[x][j](a[i+1][j]-a[i][j])=\sum_{j=1}^n(a[i+1][j]^2-a[i][j]^2)$$
如上的方程中\(a[i][j]与a[i+1][j]\)都是已知的,所以是一个n元一次方程.我们可以列出来n个,问题就转化为了求一个线性方程组的解.
对于这样的问题,我们高斯消元法.其实和我们平常解方程组消元差不多,我大概说一下思路.
我们把方程组列出来,从1~n标号.
先让2~n的方程都加上一定倍数的1号方程,把它们的1号元消去.
再让3~n的方程都加上一定倍数的2号方程,把他们的2号元消去.
...
最后发现方程组变成了一个三角形的样子.n号方程只有一个未知数即n号元.接下来从下到上代入求解.
解出n号元代入n-1号方程,求出n-1号元.
解出n-1号元代入n-2号方程,求出n-2号元.
...
这样就OK了.
但是在实现的时候我们通常用什么什么主元消元法...
意思就是在消去i号元的时候,由于i~n号方程是等价的,所以我们不一定要选i号方程,可以选其他的.
那我们选哪一个呢?我们选i号元的系数绝对值最大的那一个.
这样的话其他方程都会加上一个较小倍数的选中的方程,这样可以减少乘法带来的误差...
差不多就是这样,这东西不复杂,自己感受一下~
#include <bits/stdc++.h>
using namespace std; const int maxn=+;
int n;
double a[maxn][maxn],A[maxn][maxn];
void gause(){
for(int i=;i<=n;i++){
int t=i;
for(int j=i+;j<=n;j++)if(fabs(A[j][i])>fabs(A[t][i])) t=j;
if(t!=i) for(int j=;j<=n+;j++) swap(A[i][j],A[t][j]);
for(int j=i+;j<=n;j++){
double x=A[j][i]/A[i][i];
for(int k=i;k<=n+;k++) A[j][k]-=x*A[i][k];
}
}
for(int i=n;i>=;i--){
for(int j=i+;j<=n;j++) A[i][n+]-=A[j][n+]*A[i][j];
A[i][n+]/=A[i][i];
}
}
int main(){
scanf("%d",&n);
for(int i=;i<=n+;i++)for(int j=;j<=n;j++) scanf("%lf",&a[i][j]);
for(int i=;i<=n;i++){
for(int j=;j<=n;j++) A[i][j]=*(a[i+][j]-a[i][j]);
for(int j=;j<=n;j++) A[i][n+]+=a[i+][j]*a[i+][j]-a[i][j]*a[i][j];
}
gause();
printf("%.3lf",A[][n+]);
for(int i=;i<=n;i++) printf(" %.3lf",A[i][n+]);
return ;
}
1013: [JSOI2008]球形空间产生器sphere
Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 4166 Solved: 2191
[Submit][Status][Discuss]
Description
有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在,你被困在了这个n维球体中,你只知道球
面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧毁这个球形空间产生器。
Input
第一行是一个整数n(1<=N=10)。接下来的n+1行,每行有n个实数,表示球面上一点的n维坐标。每一个实数精确到小数点
后6位,且其绝对值都不超过20000。
Output
有且只有一行,依次给出球心的n维坐标(n个实数),两个实数之间用一个空格隔开。每个实数精确到小数点
后3位。数据保证有解。你的答案必须和标准输出一模一样才能够得分。
Sample Input
0.0 0.0
-1.0 1.0
1.0 0.0
Sample Output
HINT
提示:给出两个定义:1、 球心:到球面上任意一点距离都相等的点。2、 距离:设两个n为空间上的点A, B
的坐标为(a1, a2, …, an), (b1, b2, …, bn),则AB的距离定义为:dist = sqrt( (a1-b1)^2 + (a2-b2)^2 +
… + (an-bn)^2 )
Source
BZOJ_1013_[JSOI2008]_球形空间产生器_(高斯消元)的更多相关文章
- 【BZOJ 1013】【JSOI2008】球形空间产生器sphere 高斯消元基础题
最基础的高斯消元了,然而我把j打成i连WA连跪,考场上再犯这种错误就真的得滚粗了. #include<cmath> #include<cstdio> #include<c ...
- 【BZOJ1013】球形空间产生器(高斯消元)
[BZOJ1013]球形空间产生器(高斯消元) 题面 Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球 面上n+1个点的坐标, ...
- BZOJ 1013: [JSOI2008]球形空间产生器sphere 高斯消元
1013: [JSOI2008]球形空间产生器sphere Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/Judg ...
- [JSOI2008]球形空间产生器 (高斯消元)
[JSOI2008]球形空间产生器 \(solution:\) 非常明显的一道高斯消元.给了你n+1个球上的位置,我们知道球上任何一点到球心的距离是相等,所以我们 可以利用这一个性质.我们用n+1个球 ...
- lydsy1013: [JSOI2008]球形空间产生器sphere 高斯消元
题链:http://www.lydsy.com/JudgeOnline/problem.php?id=1013 1013: [JSOI2008]球形空间产生器sphere 时间限制: 1 Sec 内 ...
- [bzoj1013][JSOI2008][球形空间产生器sphere] (高斯消元)
Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球 面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧 ...
- 洛谷P4035 [JSOI2008]球形空间产生器(高斯消元)
洛谷题目传送门 球啊球 @xzz_233 qaq 高斯消元模板题,关键在于将已知条件转化为方程组. 可以发现题目要求的未知量有\(n\)个,题目却给了我们\(n+1\)个点的坐标,这其中必有玄机. 由 ...
- 【BZOJ1013】[JSOI2008] 球形空间产生器(高斯消元)
点此看题面 大致题意: 给定一个\(n\)维球体上的\(n+1\)个点,请你求出这个球体的圆心的位置. 列出方程 这一看就是一道解方程题. 我们可以设这个球体的圆心的位置为\((x_1,x_2,..x ...
- [luoguP4035] [JSOI2008]球形空间产生器(高斯消元)
传送门 设球心的坐标为未知量 用最后一个点来表示球面到球心的距离,那么它和前n个式子相等 移项乱搞 最后高斯消元 #include <cmath> #include <cstdio& ...
- BZOJ-1013 球形空间产生器sphere 高斯消元+数论推公式
1013: [JSOI2008]球形空间产生器sphere Time Limit: 1 Sec Memory Limit: 162 MB Submit: 3662 Solved: 1910 [Subm ...
随机推荐
- 02_线程的创建和启动_继承Thread方式
[简述] java使用Thread类代表线程,所有的线程都必须是Thread或者其子类的实例. 每个线程的任务就是完成一定的任务,实际上就是执行一段程序流. [创建并启动多线程的步骤(集成Thread ...
- springmvc(五)----异常处理
总结
- SharePoint工作流(workflow)不能自动启动
在定制工作流时,设置了当Item创建或更改时,触发工作流.用系统帐户登录时一直不会触发.这是因为这是SharePoint的安全机制,阻止了在系统帐户登陆时自动启动工作流. 解决方法:使用不是系统账户的 ...
- 安装Debian 7.8 过程,以及该系统的配置过程
最近,我给自己的电脑安装了Debian操作系统,版本7.8 stable.Debian我相信大家都能了解,它是Ubuntu的基础,Ubuntu是基于Debian的unstable开发的.我安装这个系统 ...
- Linux作业控制
在Linux中,利用Shell的作业控制是比较常用的操作,在这一节中我们将探究作业控制相关的操作.为了方便我们查看区分不同的进行,我们编写如下程序,其功能是每间隔2秒输出一次自己的编号. /* ** ...
- 前端资源多个产品整站一键打包&包版本管理(四)—— js&css文件文件打包并生成哈希后缀,自动写入路径、解决资源缓存问题。
问题: 当我们版本更新的时候,我们都要清理缓存的js跟css,如何使得在网页中不需要手动清理呢? 答案: 生成带有哈希后缀的js跟css文件 1.文件路径 路径中的conf.js 是用于放置全局打包的 ...
- Linux下GPIO驱动(五) ----misc_register();
// struct miscdevice { int minor; const char *name; const struct file_operations *fops; struct list_ ...
- 虚拟机添加磁盘LVM分区
参考博客:http://kimjinlsgd.blog.51cto.com/1918030/932210 一.查看磁盘情况 新添加一块磁盘. [root@VMhost /]# fdisk -l Dis ...
- Python 中模块间全局变量的使用上的注意
最近用Python写代码,需要用到模块间的全局变量. 网上四处搜索,发现普遍做法是把全局变量放到一个独立的模块中,使用时,导入此全局变量模块即可. 但是在实际使用过程中发现了些小问题:在使用如下代码导 ...
- OFBiz进阶之HelloWorld(一)创建热部署模块
创建热部署模块 参考文档 https://cwiki.apache.org/confluence/display/OFBIZ/OFBiz+Tutorial+-+A+Beginners+Developm ...