这道题是一个集合题,很简单的,只需要把题理解清楚,然后利用数学知识求解;
首先题上,需要注意,每个切割部分不会相交,并且切掉之后就被丢掉了,意思就是大圆的外部所对应的周长会减去,而加上小圆的内部弧长;
数学知识:
正弦(余弦)定理,弧长公式,周长公式,两点间距离公式,圆与圆的位置关系;
复习一下圆的位置关系吧(d为两个圆心的距离,计算时避免误差就可以采用平方比较大小);
几何分析:
然后这道题就可以先把有效的点给存下来,无效的点直接抛掉;合法的条件是(|R-r|<=d<|R+r|)
所以这道几何数学题就出来了,我还是第二次遇见利用数学知识求解,我之前写过一道斜率题,结果发现数学是真的有用!!
也只有比赛才能用的到实践上吧!
AC代码:
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long LL;
#define PI acos(-1)
struct Point{
LL x,y,r;
}p[200];
LL Dis(LL x,LL y){
return x*x+y*y;
}
double Angle(LL x,LL y,LL R,LL r){
LL dis=Dis(x,y);
double angle=acos((dis+r*r-R*R)/(2*sqrt(dis)*r));
return angle;
}
bool D(LL R,LL r,LL x,LL y){
LL t=Dis(x,y);
return abs(R-r)*abs(R-r)<=t&&t<(r+R)*(r+R);//说明相交
}
int main(){
LL T,x,y,m,r,R,k;
cin>>T;
while(T--){
cin>>m>>R;
k=0;
for(LL i=0;i<m;i++){
cin>>x>>y>>r;
if(D(R,r,x,y)){//有效点
p[k].x=x;p[k].y=y;p[k++].r=r;
}
}
double sum=2*PI*R;
double Angle1,Angle2,BP,IL;
for(LL i=0;i<k;i++){
Angle1=2*Angle(p[i].x,p[i].y,R,p[i].r);
Angle2=2*Angle(p[i].x,p[i].y,p[i].r,R);
BP=R*Angle2;
IL=p[i].r*Angle1;
sum+=IL;
sum-=BP;
}
printf("%.20lf\n",sum);
}
return 0;
}