首页 > 代码库 > 高斯消元-bzoj1013-球形空间产生器
高斯消元-bzoj1013-球形空间产生器
This article is made by Jason-Cow.
Welcome to reprint.
But please post the writer‘s address.
http://www.cnblogs.com/JasonCow/
推方程设n位坐标视为方程组的n个未知数
然后化为一般式。
#include <cmath> #include <cstdio> #include <algorithm> using namespace std; typedef double db; db a[15][15],b[15]; bool Guass(int n); bool work(); int main(){return work();} bool work(){ int n;scanf("%d",&n); for(int i=1;i<=n;i++)scanf("%lf",&b[i]); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++){ db x;scanf("%lf",&x); a[i][j]+=2*(b[j]-x); a[i][n+1]+=(b[j]*b[j])-(x*x); } return Guass(n); } bool Guass(int n){ for(int i=1;i<=n;i++){ db Max=-1e9;int I=i; for(int j=i;j<=n;j++)if(fabs(a[i][j])>Max)Max=fabs(a[j][i]),I=j; if(I!=i)for(int j=1;j<=n+1;j++)swap(a[i][j],a[I][j]); db t=a[i][i]; for(int j=i;j<=n+1;j++)a[i][j]/=t; for(int j=1;j<=n;j++) if(i!=j){ db t=a[j][i]; for(int k=1;k<=n+1;k++)a[j][k]-=t*a[i][k]; } } for(int i=1;i<n;i++)printf("%.3lf ",a[i][n+1]);printf("%.3lf\n",a[n][n+1]); return 0; }
高斯消元-bzoj1013-球形空间产生器
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。