Sicily.1012. Stacking Cylinders

最新推荐文章于 2017-06-30 18:47:00 发布

豪运SYSU

最新推荐文章于 2017-06-30 18:47:00 发布

阅读量542

点赞数

CC 4.0 BY-SA版权

分类专栏： Sicily

本文链接：https://blog.csdn.net/xiehaoyun2012/article/details/8438285

Sicily 专栏收录该内容

99 篇文章

订阅专栏

/*1012. Stacking Cylinders
思路:求相切球的球心，依次求得每层的球心。

*/

#include <stdlib.h>
#include <iostream>
#include <stdio.h>
#include <math.h>
#include <algorithm>
using namespace std;

struct Center{
   double x;
   double y;       
};

bool cmp(Center a, Center b){
   return a.x < b.x;     
}
//求C与A B两个球相切的第三个球的x，y坐标。 通过勾股定理得到 
void Cal(double x1,double y1,double x2,double y2,double& x3,double& y3){
    //AB两球心的距离 
    double dis=sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
    //三个球心构成的等腰三角形的中线长 
    double len=sqrt((2*2) - (dis/2)*(dis/2));
    double s,c;
    //中线与垂直线的夹角的sin值 
    s=(y2-y1)/dis;
    //中线与垂直线的夹角的cos值 
    c=(x2-x1)/dis;
    
    //球心C的x坐标 = AB球心连线的中点x坐标 - 多出的长度 
    x3=(x1+x2)/2-len*s;
    //球心C的y坐标 = AB球心连线的中点y坐标 + 少的长度 
    y3=(y1+y2)/2+len*c;
}

int main()
{
    int n;
    Center center[10];
    n = -1;
    while(n!=0){
        scanf("%d", &n);
        //cin >> n;
        if(n == 0)
          break;
        for(int i=0; i<n; i++){
            //cin >> center[i].x;
            scanf("%lf", ¢er[i].x);
            center[i].y = 1.0;  
        }        
         sort(center, center+n, cmp); 
         while(n>1){
          int k = 0;
          for( ; k<n-1; k++){
             Cal(center[k].x, center[k].y, center[k+1].x, center[k+1].y, center[k].x, center

[k].y);
          } 
          n--;                               
        }
        printf("%.4lf %.4lf\n", center[0].x, center[0].y);
        //cout << center[0].x << " " << center[0].y << endl;  
    }
    
    system("pause");
    return 0;
}