题面

有一个

n

n

n 维球,给定

n

+

1

n+1

n+1 个在球面上的点,求球心坐标。

n

10

n\leq 10

n≤10 。

题解

好久以前的题了,昨天首 A 。


n

n

n 太小了!明明可以开 100!

看到题解里有两种做法:

  • 模拟退火
  • 距离半径法高斯消元

后面那种,以前一直搞不懂,不知道怎么解二次的。于是很久都没做出来。

现在,我有一种更好的方法:向量高斯消元法。

我们知道,对于球上两个点

A

,

B

A,B

A,B 和球心

P

P

P ,满足:

P

P

P 到

A

B

AB

AB 中点的连线与

A

B

AB

AB 垂直。

于是,我们就可以列出这样的等式:

A

B

(

A

P

+

B

P

)

=

0

\stackrel{\longrightarrow}{AB}\cdot(\stackrel{\longrightarrow}{AP}+\stackrel{\longrightarrow}{BP})=0\\

AB⟶⋅(AP⟶+BP⟶)=0

(

O

B

O

A

)

(

2

O

P

(

O

B

+

O

A

)

)

=

0

(

O

B

O

A

)

2

O

P

=

(

O

B

O

A

)

(

O

B

+

O

A

)

=

O

B

2

O

A

2

(\stackrel{\longrightarrow}{OB}-\stackrel{\longrightarrow}{OA})\cdot(2\stackrel{\longrightarrow}{OP}-(\stackrel{\longrightarrow}{OB}+\stackrel{\longrightarrow}{OA}))=0\\ (\stackrel{\longrightarrow}{OB}-\stackrel{\longrightarrow}{OA})\cdot2\stackrel{\longrightarrow}{OP}=(\stackrel{\longrightarrow}{OB}-\stackrel{\longrightarrow}{OA})\cdot(\stackrel{\longrightarrow}{OB}+\stackrel{\longrightarrow}{OA})=\stackrel{\longrightarrow}{OB}^2-\stackrel{\longrightarrow}{OA}^2

(OB⟶−OA⟶)⋅(2OP⟶−(OB⟶+OA⟶))=0(OB⟶−OA⟶)⋅2OP⟶=(OB⟶−OA⟶)⋅(OB⟶+OA⟶)=OB⟶2−OA⟶2

从坐标来看,就是

i

=

1

n

2

(

X

B

[

i

]

X

A

[

i

]

)

X

P

[

i

]

=

O

B

2

O

A

2

\sum_{i=1}^n2(X_B[i]-X_A[i])\cdot \underline{~X_P[i]~}=|OB|^2-|OA|^2

i=1∑n​2(XB​[i]−XA​[i])⋅ XP​[i] ​=∣OB∣2−∣OA∣2

其中画横线的就是方程未知量,右边是两个距离的平方相减。

我们把

n

+

1

n+1

n+1 个点中每两个相邻的点列一个这样的方程。

下一步拉板子。

CODE

一定有解,所以不考虑这么多了。

#include<set>
#include<map>
#include<queue>
#include<stack>
#include<vector>
#include<bitset>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 205
#define LL long long
#define DB double
#define lowbit(x) (-(x) & (x))
#define ENDL putchar('\n')
#define FI first
#define SE second
#define SQ 450
#define eps 1e-9
LL read() {
LL f=1,x=0;int s = getchar(); if(s<0) return -1;
while(s < '0' || s > '9') {if(s=='-')f = -f;s = getchar();}
while(s >= '0' && s <= '9') {x=x*10+(s^48);s = getchar();}
return f*x;
}
void putpos(LL x) {if(!x)return ;putpos(x/10);putchar('0'+(x%10));}
void putnum(LL x) {
if(!x) {putchar('0');return ;}
if(x<0) {putchar('-');x = -x;}
return putpos(x);
}
void AIput(LL x,int c) {putnum(x);putchar(c);} int MOD = 1;
int n,m,s,o,k;
DB pt[105][105];
DB a[105][105];
DB Abs(DB x) {return x<0 ? -x:x;} void Gauss(int n) {
for(int i = 1;i <= n;i ++) {
for(int j = i+1;j <= n;j ++) {
if(Abs(a[j][i]) > Abs(a[i][i])) {
swap(a[j],a[i]); break;
}
}
for(int j = i+1;j <= n;j ++) {
DB nm = a[j][i] / a[i][i];
for(int k = n+1;k >= i;k --) {
a[j][k] -= a[i][k] * nm;
}
}
}
for(int i = n;i > 0;i --) {
for(int j = n;j > i;j --) {
a[i][n+1] -= a[i][j] * a[j][n+1];
}
a[i][n+1] /= a[i][i];
}return ;
}
int main() {
n = read();
for(int i = 1;i <= n+1;i ++) {
for(int j = 1;j <= n;j ++) {
scanf("%lf",&pt[i][j]);
}
}
for(int i = 1;i <= n;i ++) {
for(int j = 1;j <= n;j ++) {
a[i][j] = 2.0*(pt[i+1][j]-pt[i][j]);
a[i][n+1] += pt[i+1][j]*pt[i+1][j] - pt[i][j]*pt[i][j];
}
}
Gauss(n);
for(int i = 1;i <= n;i ++) {
printf("%.3f",a[i][n+1]);
putchar(i==n ? '\n':' ');
}
return 0;
}

P4035 [JSOI2008]球形空间产生器 (向量,高斯消元)的更多相关文章

  1. 洛谷P4035 [JSOI2008]球形空间产生器(高斯消元)

    洛谷题目传送门 球啊球 @xzz_233 qaq 高斯消元模板题,关键在于将已知条件转化为方程组. 可以发现题目要求的未知量有\(n\)个,题目却给了我们\(n+1\)个点的坐标,这其中必有玄机. 由 ...

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

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

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

    [JSOI2008]球形空间产生器 \(solution:\) 非常明显的一道高斯消元.给了你n+1个球上的位置,我们知道球上任何一点到球心的距离是相等,所以我们 可以利用这一个性质.我们用n+1个球 ...

  4. lydsy1013: [JSOI2008]球形空间产生器sphere 高斯消元

    题链:http://www.lydsy.com/JudgeOnline/problem.php?id=1013 1013: [JSOI2008]球形空间产生器sphere 时间限制: 1 Sec  内 ...

  5. [bzoj1013][JSOI2008][球形空间产生器sphere] (高斯消元)

    Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球 面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧 ...

  6. 【BZOJ1013】[JSOI2008] 球形空间产生器(高斯消元)

    点此看题面 大致题意: 给定一个\(n\)维球体上的\(n+1\)个点,请你求出这个球体的圆心的位置. 列出方程 这一看就是一道解方程题. 我们可以设这个球体的圆心的位置为\((x_1,x_2,..x ...

  7. [luoguP4035] [JSOI2008]球形空间产生器(高斯消元)

    传送门 设球心的坐标为未知量 用最后一个点来表示球面到球心的距离,那么它和前n个式子相等 移项乱搞 最后高斯消元 #include <cmath> #include <cstdio& ...

  8. 【BZOJ1013】球形空间产生器(高斯消元)

    [BZOJ1013]球形空间产生器(高斯消元) 题面 Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球 面上n+1个点的坐标, ...

  9. BZOJ_1013_[JSOI2008]_球形空间产生器_(高斯消元)

    描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1013 n维空间,给出球上n+1个点的n维坐标,求球心坐标. 提示:给出两个定义:1. 球心:到 ...

  10. 【BZOJ 1013】【JSOI2008】球形空间产生器sphere 高斯消元基础题

    最基础的高斯消元了,然而我把j打成i连WA连跪,考场上再犯这种错误就真的得滚粗了. #include<cmath> #include<cstdio> #include<c ...

随机推荐

  1. 【转载】解决k8s中的长连接负载均衡问题

    原文链接:一流铲屎官二流程序员[解决k8s中的长连接负载均衡问题] 长连接与短连接: 简介 长连接是指在一个TCP连接上可以连续发送多个数据包,在TCP连接保持期间,如果没有数据包发送,需要双方发检测 ...

  2. 想学嵌入式?要不一起玩 Arduino 吧

    作者:HelloGitHub-Anthony 这里是 HelloGitHub 推出的<讲解开源项目>系列,本期介绍的是如何用开源硬件开发平台 Arduino,自己动手做一个温湿度显示器. ...

  3. # 【由浅入深_打牢基础】WEB缓存投毒(上)

    image-20220611092344882 [由浅入深_打牢基础]WEB缓存投毒(上) 1. 什么是WEB缓存投毒 简单的来说,就是利用缓存将有害的HTTP响应提供给用户 什么是缓存,这里借用Bu ...

  4. 针对elementUI 中InfiniteScroll按需引入的一点注意事项

    大家为了节省空间,常常进行按需引入来节省空间,这里我给大家来介绍一下element中按需引入无限滚动指令注意的事项. 针对前面element 按需引入的一些配置这里就不再详细介绍了. 那么这里讲的是在 ...

  5. Python|range函数用法完全解读

    写在前面的一些过场话: 迭代器是 23 种设计模式中最常用的一种(之一),在 Python 中随处可见它的身影,我们经常用到它,但是却不一定意识到它的存在.在关于迭代器的系列文章中(链接见文末),我至 ...

  6. 零基础学Java(1)初识Java程序

    前言 就国内来说,Java毫无疑问是后端语言中的No.1没有之一,所以今天我们也来0基础学习Java!!! Java的好处(针对测试工程师) 面试加分->涨薪 大多数公司服务端用的都是Java, ...

  7. C语言求100以内的和的4种方式

    C语言的一个很经典的例子,帮助熟练运行几个循环的写法 * 方法一(do---while语句) #include main () { int i,sum=0; do { sum=sum+i; i++; ...

  8. 10道不得不会的JavaEE面试题

    10道不得不会的 JavaEE 面试题 我是 JavaPub,专注于面试.副业,技术人的成长记录. 以下是 JavaEE 面试题,相信大家都会有种及眼熟又陌生的感觉.看过可能在短暂的面试后又马上忘记了 ...

  9. Django WEB开发 - Django 3.0 Tutorial入门

    1. Django 官方网站 https://docs.djangoproject.com/zh-hans/3.0/intro/tutorial01/ 2. Model View Controller ...

  10. 【cartographer_ros】四: 发布和订阅里程计odom信息

    上一节介绍了激光雷达Scan传感数据的订阅和发布. 本节会介绍里程计Odom数据的发布和订阅.里程计在cartographer中主要用于前端位置预估和后端优化. 官方文档: http://wiki.r ...