zoj 2403 || poj 2194 || poj 2850 Stacking Cylinders
买一赠二啊。。。。T T 。。。poj这两题居然完全一样，就输入方式不一样，情何以堪！什么质量啊。。
给你最后一排球的坐标，求最顶层的球的坐标。
我直接一排一排模拟的，无奈用的三角函数，精度死活达不到，WA死了。。。搜题解，学到一个向量旋转的方法，http://dumbear.com/blog/?p=143 可以看下这个。总之呢，能不用精度的就不要用了。
还有一种方法就是，因为题目已经告诉你了，最顶的球的x坐标一定是最左和最右球坐标的中点，http://fghtech.blogbus.com/logs/56860901.html，这个做法，好聪明，直接推算出来y坐标。
这道题就学到了，精度是个很纠结的东西。。。以后能不用三角函数就不用。
#include <queue> #include <stack> #include <math.h> #include <stdio.h> #include <stdlib.h> #include <iostream> #include <limits.h> #include <string.h> #include <algorithm> using namespace std; const int MAX = 15; const double eps = 1e-6; bool xy(double x,double y) { return x < y - eps;} // x < y struct point{ double x,y;}; point p[MAX]; point t[MAX]; int n; double disp2p(point a,point b) { return sqrt((a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y)); } point intersection(point u1,point u2,point v1,point v2) { point ans = u1; double t = ((u1.x - v1.x)*(v1.y - v2.y) - (u1.y - v1.y)*(v1.x - v2.x))/ ((u1.x - u2.x)*(v1.y - v2.y) - (u1.y - u2.y)*(v1.x - v2.x)); ans.x += (u2.x - u1.x)*t; ans.y += (u2.y - u1.y)*t; return ans; } point Whirl(double cosl, double sinl, point a, point b) { b.x -= a.x; b.y -= a.y; point c; c.x = b.x * cosl - b.y * sinl + a.x; c.y = b.x * sinl + b.y * cosl + a.y; return c; } point find(point a,point b) { double d = disp2p(a,b); double cosl = d / 4.0; double sinl = sqrt(1.0 - cosl*cosl); point a1, b1; b1 = Whirl(cosl,sinl,a,b); //底边线段ab 绕a旋转b->b1 a1 = Whirl(cosl,-sinl,b,a); //底边线段ab绕b选择a->a1 return intersection(a,b1,a1,b); } bool cmp(point a,point b) { return xy(a.x,b.x); } int main() { while( ~scanf("%d",&n) && n ) { for(int i=0; i<n; i++) { scanf("%lf",&p[i].x); p[i].y = 1.0; } sort(p,p+n,cmp); memcpy(t,p,sizeof(p)); for(int k=n-1; k>=1; k--) { for(int i=1; i<=k; i++) t[i-1] = find(p[i-1],p[i]); for(int i=0; i<k; i++) p[i] = t[i]; } printf("%.4lf %.4lf/n",t[0].x,t[0].y); } return 0; }