BZOJ 2178 Simpson积分
思路:
我发现能用Simpson积分水的题 好像都是裸题诶233333
//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
#define pr pair<double,double>
const int N=;
int n,ban[N];double inf=1e100,l=inf,r=-inf;pr p[N];
struct Point{int x,y;};
double dis(Point a,Point b){
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
struct Circle{
Point p;int r;
pr f(double x){
if(r<=abs(p.x-x))return pr(,);
double t=sqrt(r*r-(p.x-x)*(p.x-x));
return pr(1.0*p.y-t,1.0*p.y+t);
}
}c[N];
double Cut(double x){
double ret=,last=-inf;int cnt=;
for(int i=;i<=n;i++){
p[++cnt]=c[i].f(x);
if(p[cnt]==pr(,))cnt--;
}
sort(p+,p++cnt);
for(int i=;i<=cnt;i++){
if(p[i].first>last)ret+=p[i].second-p[i].first,last=p[i].second;
else if(p[i].second>last)ret+=p[i].second-last,last=p[i].second;
}return ret;
}
double Simpson(double l,double r,double mid,double Cl,double Cr,double Cm){
double tCl=Cut((l+mid)/),tCr=Cut((mid+r)/);
double ans=(r-l)*(Cl+Cr+*Cm)/,lans=(mid-l)*(Cl+Cm+*tCl)/,rans=(r-mid)*(Cr+Cm+*tCr)/;
if(abs(lans+rans-ans)<1e-)return ans;
else return Simpson(l,mid,(l+mid)/,Cl,Cm,tCl)+Simpson(mid,r,(mid+r)/,Cm,Cr,tCr); }
int main(){
scanf("%d",&n);
for(int i=;i<=n;i++){
scanf("%d%d%d",&c[i].p.x,&c[i].p.y,&c[i].r);
l=min(l,(double)c[i].p.x-c[i].r);
r=max(r,(double)c[i].p.x+c[i].r);
}
for(int i=;i<=n;i++){
if(ban[i])continue;
for(int j=i+;j<=n;j++){
if(ban[j])continue;
if(dis(c[i].p,c[j].p)<=c[i].r-c[j].r)ban[j]=;
}
}
for(int i=;i<=n;i++)if(ban[i])swap(ban[i],ban[n]),swap(c[i--],c[n--]);
printf("%.3lf\n",Simpson(l,r,(l+r)/,,,Cut((l+r)/)));
}
BZOJ 2178 Simpson积分的更多相关文章
- BZOJ 1845 Simpson积分
思路: Simpson积分直接上 限制一下递归深度+精度就好了 (难以理解为什么这么多人写扫描线) //By SiriusRen #include <bits/stdc++.h> usi ...
- [BZOJ 2178] 圆的面积并 【Simpson积分】
题目链接:BZOJ - 2178 题目分析 用Simpson积分,将圆按照 x 坐标分成连续的一些段,分别用 Simpson 求. 注意:1)Eps要设成 1e-13 2)要去掉被其他圆包含的圆. ...
- BZOJ 2178 圆的面积并 ——Simpson积分
[题目分析] 史上最良心样例,史上最难调样例. Simpson积分硬上. 听说用long double 精度1e-10才能过. 但是double+1e-6居然过了. [代码] #include < ...
- [BZOJ 1502] [NOI2005] 月下柠檬树 【Simpson积分】
题目链接: BZOJ - 1502 题目分析 这是我做的第一道 Simpson 积分的题目.Simpson 积分是一种用 (fl + 4*fmid + fr) / 6 * (r - l) 来拟合 fl ...
- 自适应Simpson积分
自适应Simpson积分 作用 如标题所示,这玩意就是当你不会微积分的时候来求积分的. 总所周知,积分的定义就是函数的某一段与坐标轴之间的面积. 那么,自适应Simpson积分就是一种可以再某些精度下 ...
- hdu 1724 Ellipse simpson积分
/* hdu 1724 Ellipse simpson积分 求椭圆的部分面积 simpson积分法 http://zh.wikipedia.org/zh-tw/%E8%BE%9B%E6%99%AE%E ...
- 【BZOJ1502】[NOI2005]月下柠檬树 Simpson积分
[BZOJ1502][NOI2005]月下柠檬树 Description 李哲非常非常喜欢柠檬树,特别是在静静的夜晚,当天空中有一弯明月温柔地照亮地面上的景物时,他必会悠闲地坐在他亲手植下的那棵柠檬树 ...
- 【bzoj1502】[NOI2005]月下柠檬树 自适应Simpson积分
题目描述 李哲非常非常喜欢柠檬树,特别是在静静的夜晚,当天空中有一弯明月温柔地照亮地面上的景物时,他必会悠闲地坐在他亲手植下的那棵柠檬树旁,独自思索着人生的哲理.李哲是一个喜爱思考的孩子,当他看到在月 ...
- SPOJ CIRU The area of the union of circles ——Simpson积分
[题目分析] 圆的面积并. 直接Simpson积分,(但是有计算几何的解法,留着flag). simpson积分,如果圆出现了不连续的情况,是很容易出事情的.(脑补一下) 但是没有什么办法,本来就是一 ...
随机推荐
- 洛谷 3833 SHOI 2012 魔法树
[题解] 树链剖分模板题.. #include<cstdio> #include<algorithm> #include<queue> #define N 5000 ...
- saltstack(三) grains、pillar的使用
一,grains grains: 这个跟puppet的facter功能一样.主要负责采集客户端一些基本信息, 这个也完全可以自定义,可以在客户端自定义,然后自动汇报上来:也可以从服务器端定义然后推下去 ...
- enum 的使用
public enum Color { RED(), GREEN(), BLANK(), YELLO(); // 成员变量 private String name; private int index ...
- 如何通过js在子页面调用父页面元素的click事件
//获取父页面的某个元素var node = window.parent.document.getElementById("btnReturn");//调用该元素的Click事件 ...
- Ural 1114 Boxes
Boxes Time Limit: 600ms Memory Limit: 16384KB This problem will be judged on Ural. Original ID: 1114 ...
- Oracle 关于几个随机函数sys_guid、dbms_random.random、dbms_random.value(取随机的结果集)
sys_guid():SYS_GUID (),是Oracle 8i 后提供的函数.SYS_GUID产生并返回一个全球唯一的标识符(原始值)由16个字节组成.更适合多个数据库数据集成时使用(--源自百度 ...
- rm -rf & node
rm -rf & node rm -rf $ rm -rf mydir https://www.computerhope.com/issues/ch000798.htm https://sta ...
- jquery ajax获取json并解析,获取的json是object对象格式
首先我们使用的是ajax方式,推荐一个学习网址: http://blog.csdn.net/shiyaru1314/article/details/51065410 这个博主写的特别好.现在来看我们的 ...
- [bzoj3131]淘金[sdoi2013][数位DP]
求出每个数i可以被转移到的数目$f[i]$,则点$(i,j)$中的金子数目为$f[i]*f[j]$,我们就可以用优先队列求解前$k$大. 首先所有的积数目在$10^4$左右,可以先Dfs搜索出所有的数 ...
- SQLalchemy 查询总结
#简单查询 print(session.query(User).all()) print(session.query(User.name, User.fullname).all()) print(se ...