bzoj 1043 [HAOI2008]下落的圆盘——圆的周长
2024-08-26 04:11:42
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1043
算每个圆被它后面的圆盖住了多少圆弧即可。注意判断这个圆完全被后面盖住的情况。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define db double
using namespace std;
const int N=; const db pi=acos(-),pi2=pi*;
int n;db x[N],y[N],r[N],rf[N],ans;
struct Node{
db l,r;
Node(db a=,db b=):l(a),r(b) {}
bool operator< (const Node &b)const
{return l<b.l;}
}a[N<<];
db Sqr(db x){return x*x;}
db dMx(db a,db b){return a>b?a:b;}
db Fix(db a){while(a>pi2)a-=pi2;while(a<)a+=pi2;return a;}
int main()
{
scanf("%d",&n);
for(int i=;i<=n;i++)
scanf("%lf%lf%lf",&r[i],&x[i],&y[i]),rf[i]=r[i]*r[i];
for(int i=;i<=n;i++)
{
int cnt=;bool flag=;
for(int j=i+;j<=n;j++)
{
db d=Sqr(x[i]-x[j])+Sqr(y[i]-y[j]),sd=sqrt(d);
if(sd>=r[i]+r[j]||r[j]+sd<=r[i])continue;
if(r[i]+sd<=r[j]){flag=;break;}
db fx=atan2(y[j]-y[i],x[j]-x[i]);
db k=acos((rf[i]+d-rf[j])/(*r[i]*sd));
db l=Fix(fx-k),r=Fix(fx+k);//Fix
if(l<r)a[++cnt]=Node(l,r);
else a[++cnt]=Node(,r),a[++cnt]=Node(l,pi2);
}
if(flag)continue;///
sort(a+,a+cnt+);
db ret=pi2;
for(int x=,y;x<=cnt;x=y)
{
db R=a[x].r;
for(y=x+;y<=cnt&&a[y].l<=R;y++)R=dMx(R,a[y].r);
ret-=R-a[x].l;
}
ans+=r[i]*ret;
}
printf("%.3f\n",ans);
return ;
}
最新文章
- Linux Distribution / ROM
- 习题: codevs 2492 上帝造题的七分钟2 解题报告
- sql server 创建只读帐号
- select 通过jq赋值
- the philosophy behind of the design of the STL
- [stm32] USART USART1收发功能工程
- Django项目深入理解
- 锋利的jQuery-2--一个显示和隐藏的例子,主要看写法
- NeHe OpenGL教程 第四课:旋转
- java中的生产者和消费者的问题
- 获取新浪天气api显示天气情况(转)
- javascript中的字典
- Orcale 之基本术语二
- 【python】ftp连接,主被动,调试等级
- 目标检测(一)RCNN--Rich feature hierarchies for accurate object detection and semantic segmentation(v5)
- Matlab:五点差分方法求解椭圆方程非导数边值问题
- 2016年总结 - Java程序员
- opencv学习之路(5)、鼠标和滑动条操作
- Twitter开发2
- VSC KeyNote