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
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 ",ans1-ans2); 50 } 51 return 0; 52 }