| ||||||||||
| Online Judge | Problem Set | Authors | Online Contests | User | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Web Board Home Page F.A.Qs Statistical Charts | Current Contest Past Contests Scheduled Contests Award Contest | |||||||||
火候欠佳啊。wa了好久。。。#include <iostream>
#include <math.h>
using namespace std;
#define eps 1e-7
struct node
{
double x,y,z;
};
node mm[50];
node ppx[20];
double dist(node aa,node bb)
{
return sqrt((aa.x-bb.x)*(aa.x-bb.x)+(aa.y-bb.y)*(aa.y-bb.y)+(aa.z-bb.z)*(aa.z-bb.z));
}
double cal(node p,int k)
{
double dis;
dis=0;
for(int i=0;i<k;i++)
{
dis=max(dis,dist(p,mm[i]));
}
return dis;
}
void work(int p,int k,double now)
{
for(int i=0;i<p;i++)
{
ppx[i].x=(rand()%10000+1)/10000.0*100;
ppx[i].y=(rand()%10000+1)/10000.0*100;
ppx[i].z=(rand()%10000+1)/10000.0*100;
}
node d;
double dis,len;
int t;
while(now>eps)
{
for(int i=0;i<p;i++)
{
dis=0;
for(int j=0;j<k;j++)
{
len=dist(ppx[i],mm[j]);
if(dis<len)
{
dis=len;
t=j;
}
}
ppx[i].x+=(mm[t].x-ppx[i].x)/dis*now;
ppx[i].y+=(mm[t].y-ppx[i].y)/dis*now;
ppx[i].z+=(mm[t].z-ppx[i].z)/dis*now;
}
now*=0.99;
}
}
int main ()
{
// freopen("in.txt","r",stdin);
int n;
double dis,now;
node p;
while(scanf("%d",&n)&&n)
{
for(int i=0;i<n;i++)
{
scanf("%lf%lf%lf",&mm[i].x,&mm[i].y,&mm[i].z);
}
work(10,n,150);
dis=1e50;
for(int i=0;i<10;i++)
{
dis=min(dis,cal(ppx[i],n));
}
printf("%.5lf\n",dis);
}
}
Followed by: Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator