POJ 2194 2850 Stacking Cylinders

很简单，直接用高中的解析几何知识，将两个与两个圆都想切的圆的圆心求出来就可以了，在加个两重循环找到最上面的点。

#include <stdio.h>
#include <math.h>
#include<algorithm>
using namespace std;
double num[20][20], numy[20][20];

int n;

void getpoint( double x1, double y1, double x2, double y2 , double & x, double & y )
{
	double a, b, c, d, e, f;
	a= x2-x1;
	b= y2-y1;
	c= (x1+x2)/2.0;
	d= (y1+y2)/2.0;
	e= (y1-y2)/2.0;
	f= (x1-x2)/2.0;
	y=sqrt( (4-e*e-f*f)*a*a/(a*a+b*b) ) + d;
	x= ( b*( d-y )/a )+c;
}

int main()
{
	int t;
	while(scanf("%d", &t) && t)
	{
		int i, j;
		for ( i = 0; i < t ; i++ )
		{
			scanf("%lf", num[t-1]+i);
			numy[t-1][i]=1.0;
		}
		sort(num[t-1],num[t-1]+t);
		for ( j = t-1; j > 0 ; j-- )
		{
			for ( i = 0; i < t; i++ )
			{
				getpoint(num[j][i], numy[j][i], num[j][i+1], numy[j][i+1], num[j-1][i], numy[j-1][i]);
			}
		}
		printf("%.4f %.4f\n",num[0][0], numy[0][0]);
	}
	return 0;
}