题目描述 Description
给定平面上两个圆,第i个圆的圆心为(xi,yi),半径为ri。求两个圆公共部分的面积
输入描述 Input Description
给出6个实数 x1,y1,r1,x2,y2,r2 其中 r1,r2>0
输出描述 Output Description
输出公共部分面积,结果保留3位小数。
样例输入 Sample Input
20.0 30.0 15.0 40.0 30.0 30.0
样例输出 Sample Output
608.366
题解
很简单,当做一个模板吧!
#include<cstdio>
#include<iostream>
#include<math.h>
#include<cstring>
#include<cstdlib>
using namespace std;
const double pi=3.14159265;
double r1,r2,d/*圆心距*/,alf/*极角*/,bt1,bt2/*两圆心角的一半*/;
struct dian
{double x,y;} p[10];/*我是保存点的,不过一般是保存向量的*/
double lenth(double a,double b,double c,double e)/*求线段长度,或者叫向量的模*/
{return sqrt((a-c)*(a-c)+(b-e)*(b-e));}
double pointl(double jj,double zj)/*向量旋转*/
{
double m=jj+zj,n=jj-zj;
double p1x=p[1].x+cos(m)*r1, p1y=p[1].y+sin(m)*r1,
p2x=p[1].x+cos(n)*r1, p2y=p[1].y+sin(n)*r1;
//printf("%.3lf %.3lf %.3lf %.3lf\n",p1x,p1y,p2x,p2y);
return lenth(p1x,p1y,p2x,p2y);
}
void chuli()
{
if(p[2].x==p[1].x) alf=90.0;
else alf=atan((p[2].y-p[1].y)/(p[2].x-p[1].x));//算圆心连线向量AB的极角
bt1=acos((r1*r1+d*d-r2*r2)/(2*r1*d));
bt2=acos((r2*r2+d*d-r1*r1)/(2*r2*d));
//printf("%.3lf %.3lf %.3lf\n",alf,bt1,bt2);
double l=pointl(alf,bt1),jd1=bt1*2,jd2=bt2*2;//求交点之间长度。
//if(jd1<0) jd1=0-jd1; if(jd2<0) jd2=0-jd2;
printf("%.3lf",r1*r1*jd1/2+r2*r2*jd2/2-(l*d)/2);//求面积
}
int main()
{
scanf("%lf%lf%lf%lf%lf%lf",&p[1].x,&p[1].y,&r1,&p[2].x,&p[2].y,&r2);
d=lenth(p[1].x,p[1].y,p[2].x,p[2].y);
//printf("%.3lf\n",d);
if(d>=r1+r2) printf("0.000\n");
else if(d<=max(r1,r2)-min(r1,r2))
{if(r1>r2) printf("%.3lf",pi*r2*r2);
else printf("%.3lf",pi*r1*r1);
}
/*以上是特判*/
else chuli();
return 0;
}