poj2031(prime()double printf 不能用%lf）

来源：互联网发布：千牛软件何用编辑：程序博客网时间：2024/05/22 19:26

这里写图片描述

大致题意：

就是给出三维坐标系上的一些球的球心坐标和其半径，搭建通路，使得他们能够相互连通。如果两个球有重叠的部分则算为已连通，无需再搭桥。求搭建通路的最小费用（费用就是边权，就是两个球面之间的距离）。

解题思路：

不要被三维吓到了，其实就是图论的最小生成树问题

球心坐标和半径是用来求两点之间的边权的，求出边权后，把球看做点，用邻接矩阵存储这个无向图，再求最小生成树，非常简单的水题。

把球A和球B看做无向图图的两个结点，那么

边权 = AB球面距离 = A球心到B球心的距离 – A球半径 – B球半径。
心塞….最后被一个输出卡了一下午….double 输出用printf输出时 %lf不标准…可能不会四舍五入…..wrong了一下午…要么printf(“%.3f”,…)
或是 cout .

#include<iostream>#include<iomanip>#include<cstdio>#include<cstring>#include<cmath>#define eps 1e-8using namespace std;double M[220][220],dist[220];int book[220];struct Node{    double x,y,z,r;}node[220];int n;int EPS(double k){  if(fabs(k)<eps)    return 0;   return k>0?1:-1;}double dis(Node a,Node b){    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)+(a.z-b.z)*(a.z-b.z))-a.r-b.r;}double prime(int cur){    memset(book,0,sizeof(book));    book[cur]=1;    for(int i=1;i<=n;i++)        dist[i]=M[cur][i];    int index;double sum=0.0;    for(int i=2;i<=n;i++)    {        double Min=999999.00000;        for(int j=1;j<=n;j++)        {            if(!book[j]&&dist[j]<Min)            {                Min=dist[j];                index=j;            }        }        book[index]=1;        sum+=Min;        for(int j=1;j<=n;j++)        {            if(!book[j]&&M[index][j]<dist[j])                dist[j]=M[index][j];        }    }    return sum;}int main(){    while(~scanf("%d",&n)&&n)    {        for(int i=1;i<=n;i++)            scanf("%lf %lf %lf %lf",&node[i].x,&node[i].y,&node[i].z,&node[i].r);        memset(M,0,sizeof(M));        for(int i=1;i<=n;i++)        {            for(int j=1;j<=n;j++)            {                double d=dis(node[i],node[j]);                if(EPS(d)<=0)                    M[i][j]=M[j][i]=0;                else                    M[i][j]=M[j][i]=d;            }        }       //cout<<fixed<<setprecision(3)<<prime(1)<<endl;       printf("%.3f",prime());    }}

0 0