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以内的植物提供水。(如图)

bnu 4067 美丽的花环-LMLPHP

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

bnu 4067 美丽的花环-LMLPHP

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 ;

 }