bnu 4067 美丽的花环

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

美丽的花环

Time Limit: 1000ms

Case Time Limit: 1000ms

Memory Limit: 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 }