wikioi 两圆的交

给定平面上两个圆,第i个圆的圆心为(xi,yi),半径为ri。求两个圆公共部分的面积

给出6个实数 x1,y1,r1,x2,y2,r2     其中 r1,r2>0

输出公共部分面积,结果保留3位小数。

20.0 30.0 15.0 40.0 30.0 30.0

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;
}




  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值