bnu 4067 美丽的花环

http://www.bnuoj.com/bnuoj/problem_show.php?pid=4067

美丽的花环

1000ms
1000ms
65536KB
Special Judge
 
64-bit integer IO format: %lld      Java class name: Main
Submit Status
Font Size:  

 学校的草坪上最近种植了一些漂亮的花卉，所有的花围成了一个环形(内径为r,外径为R,0 <r < R) 。原来这片地上有一个用于喷灌的喷头。这个喷头可以为半径K以内的植物提供水。(如图) 

现在，HK请你帮忙计算一下，花构成的环形当中有多大面积的可以由喷头提供灌溉。 

 

Input

 输入数据只有两行。 
格式为(中间使用空格隔开) 

x1 y1 R r 
x2 y2 K 

花环中心在(x1,y1) ，外径为R,内径为r。 
喷头的位置为(x2,y2),覆盖范围为K。 
其中，x1,y1,R,r,x2,y2,K均为实数。

 

Output

 输出题目描述中所求面积(保留2位小数)。

 

Sample Input

5 5 20 10
2 2 20
 

Sample Output

773.09
 

Hint

π的值请取：3.1415926或acos(-1.0) 

 

 

 1 /*
 2   两个圆的相交面积
 3 */
 4 
 5 
 6 #include<iostream>
 7 #include<cstdio>
 8 #include<cstring>
 9 #include<cmath>
10 using namespace std;
11 const double PI=3.141592653;
12 struct point
13 {
14     double x,y;
15     double r;
16 };
17 double dis(point a,point b)
18 {
19     return sqrt( (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y) );
20 }
21 double fun(point a,point b)  //求两圆的相交面积函数。
22 {
23     double d=dis(a,b);
24     if(d>=a.r+b.r) return 0;
25     double r=(a.r>b.r?b.r:a.r);
26     if( d<=fabs(a.r-b.r) ) return PI*r*r;
27     double A1=acos( (a.r*a.r+d*d-b.r*b.r)/2/a.r/d );
28     double A2=acos( (b.r*b.r+d*d-a.r*a.r)/2/b.r/d );
29     double res=A1*a.r*a.r + A2*b.r*b.r;
30     res-=sin(A1)*a.r*d;
31     return res;
32 }
33 int main()
34 {
35     point a,b;
36     double x1,y1,RR,rr;
37     double x2,y2,K;
38     while(scanf("%lf%lf%lf%lf",&x1,&y1,&RR,&rr)>0)
39     {
40         scanf("%lf%lf%lf",&x2,&y2,&K);
41         a.x=x1;a.y=y1;a.r=RR;
42         b.x=x2;b.y=y2;b.r=K;
43         double ans1=fun(a,b);
44 
45         a.x=x1;a.y=y1;a.r=rr;
46         b.x=x2;b.y=y2;b.r=K;
47         double ans2=fun(a,b);
48 
49         printf("%.2lf\n",ans1-ans2);
50     }
51     return 0;
52 }

 

posted @ 2013-10-01 18:39  芷水  阅读(147)  评论(0)  编辑  收藏  举报