POJ 2194/POJ 2850

题意:

一堆圆，R=1，已知底层各个的圆心坐标，求最上面一个圆心的坐标。

题解：

模拟，尽量少使用反三角函数，不得不用时，尽量用atan2即可。

View Code

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cstring>
 4 #include <cstdlib>
 5 #include <algorithm>
 6 #include <cmath>
 7 
 8 #define N 222
 9 
10 using namespace std;
11 
12 struct PO
13 {
14     double x,y;
15 }p[N][N];
16 
17 int n;
18 
19 inline bool cmp(const PO &a,const PO &b)
20 {
21     return a.x<b.x;
22 }
23 
24 inline void read()
25 {
26     for(int i=1;i<=n;i++)
27     {
28         scanf("%lf",&p[1][i].x);
29         p[1][i].y=1.0;
30     }
31     sort(p[1]+1,p[1]+1+n,cmp);
32 }
33 
34 inline PO rotate(PO a,double hd)
35 {
36     PO c;
37     c.x=a.x*cos(hd)-a.y*sin(hd);
38     c.y=a.x*sin(hd)+a.y*cos(hd);
39     return c;
40 }
41 
42 inline double getlen(PO &a)//向量的模 
43 {
44     return sqrt(a.x*a.x+a.y*a.y);
45 }
46 
47 inline double getdis2(PO &a,PO &b)
48 {
49     return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
50 }
51 
52 inline void go()
53 {
54     int c=1;
55     while(n>1)
56     {
57         c++;
58         for(int i=1;i<n;i++)
59         {
60             double af=atan2(sqrt(4.0-getdis2(p[c-1][i],p[c-1][i+1])*0.25),sqrt(getdis2(p[c-1][i],p[c-1][i+1]))*0.5);
61             //atan2精度较高，其他的反三角函数精度都比较差。。。 
62             PO s;
63             s.y=p[c-1][i+1].y-p[c-1][i].y;
64             s.x=p[c-1][i+1].x-p[c-1][i].x;
65             double k=2.0/getlen(s);
66             s.y*=k; s.x*=k;
67             
68             s=rotate(s,af);
69             p[c][i].y=p[c-1][i].y+s.y;
70             p[c][i].x=p[c-1][i].x+s.x;
71         }
72         n--;
73     }
74     printf("%.4lf %.4lf\n",p[c][1].x,p[c][1].y);
75 }
76 
77 int main()
78 {
79     while(scanf("%d",&n),n) read(),go();
80     return 0;
81 }

依旧在写水题。。。