2018-2019 ACM-ICPC, Asia Nanjing Regional Contest D Country Meow(최소 오버레이)
제목:
3차원 공간은 n개의 점을 주고 공간에서 하나의 점에서 이 점까지의 최장 거리, 출력 거리를 구한다.
사고방식: 판자 문제
코드:
#include
using namespace std;
typedef long long ll;
const double eps=1e-10;
struct Tpoint{
double x,y,z;
};
Tpoint operator +(Tpoint a,Tpoint b){
Tpoint c;
c.x=a.x+b.x;
c.y=a.y+b.y;
c.z=a.z+b.z;
return c;
}
Tpoint operator -(Tpoint a,Tpoint b){
Tpoint c;
c.x=a.x-b.x;
c.y=a.y-b.y;
c.z=a.z-b.z;
return c;
}
Tpoint operator /(Tpoint a,double k){
Tpoint c;
c.x=a.x/k;
c.y=a.y/k;
c.z=a.z/k;
return c;
}
Tpoint operator *(Tpoint a,double k){
Tpoint c;
c.x=a.x*k;
c.y=a.y*k;
c.z=a.z*k;
return c;
}
int npoint,nouter;
Tpoint pt[200000],outer[4],res;
double radius,tmp;
double dist2(Tpoint a,Tpoint b){
return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)+(a.z-b.z)*(a.z-b.z);
}
double dot(Tpoint a,Tpoint b){
return a.x*b.x+a.y*b.y+a.z*b.z;
}
void ball(){
Tpoint q[3];
double m[3][3],sol[3],l[3],det;
int i,j;
res.x=res.y=res.z=radius=0;
switch(nouter){
case 1:
res=outer[0];break;
case 2:
res=(outer[0]+outer[1])/2;
radius=dist2(res,outer[0]);
break;
case 3:
for(i=0;i<2;i++) q[i]=outer[i+1]-outer[0];
for(i=0;i<2;i++)
for(j=0;j<2;j++)
m[i][j]=dot(q[i],q[j])*2;
for(i=0;i<2;i++) sol[i]=dot(q[i],q[i]);
if(fabs(det=m[0][0]*m[1][1]-m[0][1]*m[1][0])eps){
outer[nouter++]=pt[i];
minball(i);
--nouter;
if(i>0){
Tpoint Tt=pt[i];
memmove(&pt[1],&pt[0],sizeof(Tpoint)*i);
pt[0]=Tt;
}
}
}
int main()
{
cin>>npoint;
for(int i=0;i>pt[i].x>>pt[i].y>>pt[i].z;
random_shuffle(pt,pt+npoint);
radius=-1;
for(int i=0;ieps)
nouter=1,outer[0]=pt[i],minball(i);
printf("%.5lf
",sqrt(radius));
return 0;
}
이 내용에 흥미가 있습니까?
현재 기사가 여러분의 문제를 해결하지 못하는 경우 AI 엔진은 머신러닝 분석(스마트 모델이 방금 만들어져 부정확한 경우가 있을 수 있음)을 통해 가장 유사한 기사를 추천합니다:
다양한 언어의 JSONJSON은 Javascript 표기법을 사용하여 데이터 구조를 레이아웃하는 데이터 형식입니다. 그러나 Javascript가 코드에서 이러한 구조를 나타낼 수 있는 유일한 언어는 아닙니다. 저는 일반적으로 '객체'{}...
텍스트를 자유롭게 공유하거나 복사할 수 있습니다.하지만 이 문서의 URL은 참조 URL로 남겨 두십시오.
CC BY-SA 2.5, CC BY-SA 3.0 및 CC BY-SA 4.0에 따라 라이센스가 부여됩니다.