题意:http://www.lydsy.com/JudgeOnline/problem.php?id=2178

sol  :是谁.......是谁往题里下毒......

   辛普森积分,每次判断左边+右边是否-lastans是否<eps,递归处理

   公式:f'(l~r)=(r-l)*(f(l)+f(r)+4*f(mid))/6

   f(i)表示x=i时与圆相交部分的和,然后直接上辛普森积分就可以水过了

   P.S.此题有坑QAQ

     eps设成1e-13才能过,1e-12会WA

     需要剪枝:将包含在其他圆内的圆删去

   P.S.此题有毒.....不知道为什么排序时写cmp就TLE了,然而重载运算符就AC了也是很迷啊QAQ

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#define eps 1e-13
using namespace std;
const int Mx=;
int n,tot,st,to;
bool vis[Mx];
double L[Mx],R[Mx];
struct Node { double x,y,r; } str[Mx];
struct Line { double l,r; } line[Mx];
double dis(Node a,Node b) { return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); }
bool cmp1(Node a,Node b) { return a.r<b.r; }
bool cmp2(Node a,Node b) { return a.x-a.r<b.x-b.r; }
bool operator < (Line a,Line b) { return a.l<b.l; }
double cal(double fl,double fmid,double fr,double r_l) { return r_l*(fl+fr+*fmid)/; }
double get_f(double x)
{
double ans=; tot=;
for(int i=;i<=n;i++)
if(x>=L[i]&&x<=R[i])
{
double len=sqrt(str[i].r*str[i].r-(str[i].x-x)*(str[i].x-x));
line[++tot].l=str[i].y-len; line[tot].r=str[i].y+len;
}
sort(line+,line+tot+);
for(int i=,j;i<=tot;i++)
{
double minl=line[i].l,minr=line[i].r;
for(j=i+;j<=tot;j++)
{
if(line[j].l>minr) break;
if(line[j].r>minr) minr=line[j].r;
}
ans+=minr-minl; i=j-;
}
return ans;
}
double simpson(double l,double r,double mid,double fl,double fr,double fmid,double lastans)
{
double mid1=(l+mid)/,mid2=(mid+r)/,fmid1=get_f(mid1),fmid2=get_f(mid2);
double ans1=cal(fl,fmid1,fmid,mid-l),ans2=cal(fmid,fmid2,fr,r-mid);
if(fabs(ans1+ans2-lastans)<eps) return ans1+ans2;
return simpson(l,mid,mid1,fl,fmid,fmid1,ans1)+simpson(mid,r,mid2,fmid,fr,fmid2,ans2);
}
void solve()
{
double ans=0.00;
for(int i=,j;i<=n;i++)
{
st=i; double minl=L[i],minr=R[i],minmid;
for(j=i+;j<=n;j++)
{
if(L[j]>minr) break;
if(R[j]>minr) minr=R[j];
}
to=j-,i=to,minmid=(minl+minr)/;
double fl=get_f(minl),fr=get_f(minr),fmid=get_f(minmid),lstans=cal(fl,fmid,fr,minr-minl);
ans+=simpson(minl,minr,minmid,fl,fr,fmid,lstans);
}
printf("%.3lf",ans);
}
int main()
{
scanf("%d",&n);
for(int i=;i<=n;i++) scanf("%lf%lf%lf",&str[i].x,&str[i].y,&str[i].r);
sort(str+,str++n,cmp1);
//剪枝
for(int i=;i<=n;i++)
for(int j=i+;j<=n;j++)
if(fabs(str[j].r-str[i].r)>=dis(str[j],str[i]))
{
vis[i]=; break;
}
for(int i=,tmp=,mx=n;i<=mx;i++) if(!vis[i]) str[++tmp]=str[i],n=tmp;
//
sort(str+,str++n,cmp2);
for(int i=;i<=n;i++) L[i]=str[i].x-str[i].r,R[i]=str[i].x+str[i].r;
solve();
return ;
}

bzoj2178:圆的面积并的更多相关文章

  1. [SPOJ-CIRU]The area of the union of circles/[BZOJ2178]圆的面积并

    [SPOJ-CIRU]The area of the union of circles/[BZOJ2178]圆的面积并 题目大意: 求\(n(n\le1000)\)个圆的面积并. 思路: 对于一个\( ...

  2. BZOJ2178: 圆的面积并(格林公式)

    题面 传送门 题解 好神仙-- 先给几个定义 平面单连通区域:设\(D\)是平面内一区域,若属于\(D\)内任一简单闭曲线的内部都属于\(D\),则称\(D\)为单连通区域.通俗地说,单连通区域是没有 ...

  3. BZOJ2178 圆的面积并 计算几何 辛普森积分

    原文链接https://www.cnblogs.com/zhouzhendong/p/BZOJ2178.html 题目传送门 - BZOJ2178 题意 给出 $n(n\leq 1000)$ 个圆,求 ...

  4. bzoj2178: 圆的面积并

    Description 给出N个圆,求其面积并 Input 先给一个数字N ,N< = 1000 接下来是N行是圆的圆心,半径,其绝对值均为小于1000的整数 Output 面积并,保留三位小数 ...

  5. BZOJ2178 圆的面积并(simpson积分)

    板子题.可以转一下坐标防止被卡.精度和常数实在难以平衡. #include<iostream> #include<cstdio> #include<cmath> # ...

  6. 【BZOJ2178】圆的面积并(辛普森积分)

    [BZOJ2178]圆的面积并(辛普森积分) 题面 BZOJ 权限题 题解 把\(f(x)\)设为\(x\)和所有圆交的线段的并的和. 然后直接上自适应辛普森积分. 我精度死活一个点过不去,不要在意我 ...

  7. 【题解】CIRU - The area of the union of circles [SP8073] \ 圆的面积并 [Bzoj2178]

    [题解]CIRU - The area of the union of circles [SP8073] \ 圆的面积并 [Bzoj2178] 传送门: \(\text{CIRU - The area ...

  8. JAVA求圆的面积

    import java.text.DecimalFormat;import java.util.Scanner; public class TheAreaOfCircle { public stati ...

  9. c语言求平面上2个坐标点的直线距离、求俩坐标直线距离作为半径的圆的面积、递归、菲波那次数列、explode

    #include <stdio.h> #include <math.h> #include <string.h> char explode( char * str ...

  10. 【BZOJ】2178: 圆的面积并

    http://www.lydsy.com/JudgeOnline/problem.php?id=2178 题意:给出n<=1000个圆,求这些圆的面积并 #include <cstdio& ...

随机推荐

  1. 返回固定数据的web服务器

    import socket def handle_client(socket_con): """ 接收来自客户端的请求,并接收请求报文,解析,返回 "" ...

  2. FastDFS文件管理系统

    一.FastDFS介绍 FastDFS 是一个开源的高性能分布式文件系统(DFS). 它的主要功能包括:文件存储,文件同步和文件访问,以及高容量和负载平衡.主要解决了海量数据存储问题,特别适合以中小文 ...

  3. 前端vue项目部署到tomcat,一刷新报错404解决方法

    公司前端写的后台部署到tomcat webapps目录下后,无法进行刷新,一刷新就会报错404,自动跳的404页面.在网上查了下,官方说是HTML5 History 模式引发的问题,但是解决方案中,并 ...

  4. 虚拟机VMware安装Kali Linux

    本文讲解如何在虚拟机上安装Kali Linux,希望对大家有所帮助. 准备:一台电脑,VMware(VMware安装教程) 一.下载系统镜像文件 1.首先下载系统镜像,进入kali官网,在Downlo ...

  5. PHP.19-验证码生成

    生成验证码 思路:先定义验证码函数getCode() //绘制验证码 $num = 4; //字符长度 getCode($num, 2); 1.创建画布,分配颜色 imagecreatetruecol ...

  6. 14 Django的用户认证组件

    用户认证 auth模块 from django.contrib import auth django.contrib.auth中提供了许多方法,这里主要介绍其中的三个: 1.1 .authentica ...

  7. springmvc上传图片并显示--支持多图片上传

    实现上传图片功能在Springmvc中很好实现.现在我将会展现完整例子. 开始需要在pom.xml加入几个jar,分别是: <dependency> <groupId>comm ...

  8. 《数据结构与算法分析:C语言描述》复习——第九章“图论”——拓扑排序

    2014.07.04 17:23 简介: 我们考虑一种特殊的图: 1. 有向图 2. 只有一个连通分量 3. 不存在环 那么这样的图里,必然可以找到一种排序方式,来确定谁在谁的“前面”. 简单的来说可 ...

  9. Java线程的两种实现形式

    一.创建线程的第一种方式:继承Thread类 class Demo extends Thread{ @Override public void run() { super.run(); for(int ...

  10. 油田(DFS)

    //DFS:油田问题 #include <iostream> using namespace std; ][]; int n,m; //一个网格的8个方向 ][] = {{-,-},{-, ...