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) 

 /*
   两个圆的相交面积
 */

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cmath>
 using namespace std;
 const double PI=3.141592653;
 struct point
 {
     double x,y;
     double r;
 };
 double dis(point a,point b)
 {
     return sqrt( (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y) );
 }
 double fun(point a,point b)  //求两圆的相交面积函数。
 {
     double d=dis(a,b);
     if(d>=a.r+b.r) return ;
     double r=(a.r>b.r?b.r:a.r);
     if( d<=fabs(a.r-b.r) ) return PI*r*r;
     double A1=acos( (a.r*a.r+d*d-b.r*b.r)//a.r/d );
     double A2=acos( (b.r*b.r+d*d-a.r*a.r)//b.r/d );
     double res=A1*a.r*a.r + A2*b.r*b.r;
     res-=sin(A1)*a.r*d;
     return res;
 }
 int main()
 {
     point a,b;
     double x1,y1,RR,rr;
     double x2,y2,K;
     while(scanf("%lf%lf%lf%lf",&x1,&y1,&RR,&rr)>)
     {
         scanf("%lf%lf%lf",&x2,&y2,&K);
         a.x=x1;a.y=y1;a.r=RR;
         b.x=x2;b.y=y2;b.r=K;
         double ans1=fun(a,b);

         a.x=x1;a.y=y1;a.r=rr;
         b.x=x2;b.y=y2;b.r=K;
         double ans2=fun(a,b);

         printf("%.2lf\n",ans1-ans2);
     }
     return ;
 }