BZOJ2178 圆的面积并 计算几何 辛普森积分
原文链接https://www.cnblogs.com/zhouzhendong/p/BZOJ2178.html
题目传送门 - BZOJ2178
题意
给出 $n(n\leq 1000)$ 个圆,求面积并。
所有圆的圆心坐标和半径都是绝对值不大于 1000 的整数。
题解
自适应辛普森积分模板题。注意先删掉被其他圆包含的圆。
但是 bzoj 大概是加过数据了,网上大部分直接自适应辛普森的代码都 TLE 了。
有一种卡常方法效果很好:
把 x 坐标按照整点划分成 $O(1000)$ 个区间,对于每一个区间,删除不涉及这个区间的所有圆。然后就过了。
代码
#include <bits/stdc++.h>
#define y1 __zzd12323
using namespace std;
const int N=1005;
const double Eps=1e-9,INF=1e9;
double sqr(double x){
return x*x;
}
int n;
struct Circle{
double x,y,r;
}c[N],_c[N];
bool cmpr(Circle a,Circle b){
return a.r>b.r;
}
struct Seg{
double L,R;
}s[N];
bool cmp(Seg a,Seg b){
return a.L<b.L;
}
double F(double x){
int m=0;
for (int i=1;i<=n;i++){
if (fabs(x-c[i].x)>c[i].r-Eps)
continue;
double M=c[i].y,a=sqrt(sqr(c[i].r)-sqr(x-c[i].x));
s[++m].L=M-a,s[m].R=M+a;
}
sort(s+1,s+m+1,cmp);
double res=0.0,mx=-INF;
for (int i=1;i<=m;i++)
if (s[i].L>mx)
res+=s[i].R-s[i].L,mx=s[i].R;
else if (s[i].R>mx)
res+=s[i].R-mx,mx=s[i].R;
return res;
}
double Simpson(double a,double b,double Fa,double Fb,double Fab){
return (b-a)/6*(Fa+Fab*4+Fb);
}
double Asr(double a,double b,double v,double Fa,double Fb,double Fab){
double ab=(a+b)/2;
double LFab=F((a+ab)/2),L=Simpson(a,ab,Fa,Fab,LFab);
double RFab=F((ab+b)/2),R=Simpson(ab,b,Fab,Fb,RFab);
if (fabs(v-L-R)<Eps)
return v;
return Asr(a,ab,L,Fa,Fab,LFab)
+Asr(ab,b,R,Fab,Fb,RFab);
}
int killed[N];
double Dis(double x1,double y1,double x2,double y2){
return sqrt(sqr(x1-x2)+sqr(y1-y2));
}
int main(){
scanf("%d",&n);
double a=INF,b=-INF;
for (int i=1;i<=n;i++){
scanf("%lf%lf%lf",&c[i].x,&c[i].y,&c[i].r);
a=min(a,c[i].x-c[i].r);
b=max(b,c[i].x+c[i].r);
}
sort(c+1,c+n+1,cmpr);
for (int i=1;i<=n;i++)
for (int j=i+1;j<=n;j++)
if (Dis(c[i].x,c[i].y,c[j].x,c[j].y)+c[j].r<c[i].r+Eps)
killed[j]=1;
int _n=0;
for (int i=1;i<=n;i++)
if (!killed[i])
c[++_n]=c[i];
n=_n;
double ans=0;
for (double x=a;x<b;x++){
for (int i=1;i<=n;i++)
_c[i]=c[i];
int _n=n;
n=0;
for (int i=1;i<=_n;i++)
if (c[i].x-c[i].r<i+1&&c[i].x+c[i].r>i)
c[++n]=_c[i];
n=_n;
for (int i=1;i<=n;i++)
c[i]=_c[i];
double Fa=F(x),Fb=F(x+1),Fab=F((x+x+1)/2);
ans+=Asr(x,x+1,Simpson(x,x+1,Fa,Fb,Fab),Fa,Fb,Fab);
}
printf("%.3lf",ans);
return 0;
}
BZOJ2178 圆的面积并 计算几何 辛普森积分的更多相关文章
- 【BZOJ2178】圆的面积并(辛普森积分)
[BZOJ2178]圆的面积并(辛普森积分) 题面 BZOJ 权限题 题解 把\(f(x)\)设为\(x\)和所有圆交的线段的并的和. 然后直接上自适应辛普森积分. 我精度死活一个点过不去,不要在意我 ...
- BZOJ2178 圆的面积并(simpson积分)
板子题.可以转一下坐标防止被卡.精度和常数实在难以平衡. #include<iostream> #include<cstdio> #include<cmath> # ...
- [SPOJ-CIRU]The area of the union of circles/[BZOJ2178]圆的面积并
[SPOJ-CIRU]The area of the union of circles/[BZOJ2178]圆的面积并 题目大意: 求\(n(n\le1000)\)个圆的面积并. 思路: 对于一个\( ...
- BZOJ2178: 圆的面积并(格林公式)
题面 传送门 题解 好神仙-- 先给几个定义 平面单连通区域:设\(D\)是平面内一区域,若属于\(D\)内任一简单闭曲线的内部都属于\(D\),则称\(D\)为单连通区域.通俗地说,单连通区域是没有 ...
- bzoj2178: 圆的面积并
Description 给出N个圆,求其面积并 Input 先给一个数字N ,N< = 1000 接下来是N行是圆的圆心,半径,其绝对值均为小于1000的整数 Output 面积并,保留三位小数 ...
- [BZOJ 2178] 圆的面积并 【Simpson积分】
题目链接:BZOJ - 2178 题目分析 用Simpson积分,将圆按照 x 坐标分成连续的一些段,分别用 Simpson 求. 注意:1)Eps要设成 1e-13 2)要去掉被其他圆包含的圆. ...
- bzoj 2178 圆的面积并【simpson积分】
直接套simpson,f可以直接把圆排序后扫一遍所有圆,这样维护一个区间就可以避免空段. 然而一定要去掉被其他圆完全覆盖的圆,否则会TLE #include<iostream> #incl ...
- BZOJ 2178: 圆的面积并 [辛普森积分 区间并]
2178: 圆的面积并 Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 1740 Solved: 450[Submit][Status][Discus ...
- bzoj 2178 圆的面积并 —— 辛普森积分
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2178 先看到这篇博客:https://www.cnblogs.com/heisenberg- ...
随机推荐
- Find Lines
(Uva 6955可以直接随机,湖大OJ 13348 要优化) 题意:给出 n个点的坐标, 一个 百分数p, 求是否有一条直线上有 n * p /100个点… 随机化算法,但也要优化下……(TLE, ...
- hadoop生态系统主要架构图汇总
1 hadoop1.0时期架构 2 hadoop2.0时期架构 3 hdfs架构 [Active Namenode]:主 Master(只有一个),管理 HDFS 的名称空间,管理数据块映射信息:配置 ...
- Vue -cli 入门 --项目搭建(一)
一. 安装node.js环境. 在node.js官网下载稳定版本(https://nodejs.org/en/) 下载完成后点击安装,安装过程很简单,一直next即可,安装完成会自动添加node及np ...
- 借助dubbo-admin来管理你的服务
1. Github上下载最新的dubbo源码包并解压 2. 修改配置信息(打开 dubbo-admin/src/main/webapp/WEB-INF下的dubbo.properties,修 ...
- 【翻译】关于vertical-align所有你需要知道的
本文是翻译过来的,如果有不对的地方还请指教~,原文链接:Vertical-Align: All You Need To Know 前面一些说明,可以略过不看吧 我经常需要对元素进行垂直方向上的布局. ...
- Java 开源博客 Solo 2.5.0 发布
Java 开源博客 Solo 2.5.0 发布 Solo 是一款一个命令就能搭建好的 Java 开源博客系统,如果你想开个独立博客,请一定不要错过! 2.5.0 版本主要支持了 Markdown/JS ...
- AD9361寄存器配置顺序,循环模式,自收自发
:] cmd_data; :] index; begin case(index) 'h000,8'h00};//set spi -- 'h3df,8'h01};//set init -- 'h037, ...
- uva11426 欧拉函数应用,kuangbin的筛法模板
/* 给定n,对于所有的对(i,j),i<j,求出sum{gcd(i,j)} 有递推式sum[n]=sum[n-1]+f[n] 其中f[n]=gcd(1,n)+gcd(2,n)+gcd(3,n) ...
- cf862d 交互式二分
/* 二分搜索出一个01段或10即可 先用n个0确定1的个数num 然后测试区间[l,mid]是否全是0或全是1 如果是,则l=mid,否则r=mid,直到l+1==r 然后再测试l是1还是r是1 如 ...
- jvm(转)
原:https://blog.csdn.net/luomingkui1109/article/details/72820232 1.JVM简析: 作为一名Java使用者,掌握JVM的体系结构 ...