标签:LG4035 ch int sum JSOI2008 fh getchar 高斯消 register
问题描述
题解
设答案为\((p_1,p_2,p_3,...,p_n)\)
因为是一个球体,令其半径为\(r\),则有
\[\sum_{i=1}^{n}{(a_i-p_i)}^2={\rm dis}^2\]
拆式子可得
\[\sum_{i=1}^{n}a_i^2-2\times\sum_{i=1}^{n}{a_ip_i}=\sum_{i=1}^{n}p_i^2-{\rm dis}^2\]
于是可以构造出新的方程矩阵:
\[f_{i,j}=2 \times (a_{i+1,j}-a_{i,j})\]
\[f_{i,n+1}=\sum_{j=1}^n a_{i+1,j}^2-a_{i,j}^2\]
\(\mathrm{Code}\)
#include<bits/stdc++.h>
using namespace std;
void read(int &x){
x=0;char ch=1;int fh;
while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
if(ch=='-') fh=-1,ch=getchar();
else fh=1;
while(ch>='0'&&ch<='9'){
x=(x<<1)+(x<<3)+ch-'0';
ch=getchar();
}
x*=fh;
}
#define maxn 107
int n;
double a[maxn][maxn],bf[maxn][maxn];
int pla;
int main(){
ios::sync_with_stdio(0);
cin>>n;
for(register int i=1;i<=n+1;i++){
for(register int j=1;j<=n;j++) cin>>bf[i][j];
}
for(register int i=1;i<=n+1;i++){
for(register int j=1;j<=n;j++){
a[i][j]=2*(bf[i+1][j]-bf[i][j]);
a[i][n+1]+=bf[i+1][j]*bf[i+1][j]-bf[i][j]*bf[i][j];
}
}
for(register int i=1;i<=n;i++){
pla=i;
while(pla<=n&&a[pla][i]==0) pla++;
if(pla==n+1){
puts("No Solution");return 0;
}
for(register int j=1;j<=n+1;j++) swap(a[i][j],a[pla][j]);
double tmp=a[i][i];
for(register int j=1;j<=n+1;j++) a[i][j]=a[i][j]/tmp;
for(register int j=1;j<=n;j++){
if(i==j) continue;
double rp=a[j][i];
for(register int k=1;k<=n+1;k++){
a[j][k]=a[j][k]-rp*a[i][k];
}
}
}
for(register int i=1;i<=n;i++){
cout<<fixed<<setprecision(3)<<a[i][n+1]<<" ";
}
return 0;
}
标签:LG4035,ch,int,sum,JSOI2008,fh,getchar,高斯消,register 来源: https://www.cnblogs.com/liubainian/p/11484275.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。