simpson
使用二次函数拟合复杂的连续函数求积分
对于(l,r)拟合的积分为(r-l)*(f(l)+4*f((l+r)/2)+f(r))/6
___________________________
BZOJ2178
- #include <cstdio>
- #include <cmath>
- #include <algorithm>
- #define LDB double
- using namespace std;
- const LDB eps=1e-;
- LDB ans=;
- struct data{
- LDB x,y,r,va;
- }a[],tmp[];
- int st,en,n,del[];
- int mycomp (const data&a,const data&b){
- return(a.va<b.va);
- }
- LDB dis(data a,data b){
- return(sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)));
- }
- LDB getf(LDB xx){
- int cnt=;
- for (int i=st;i<=en;i++)
- if (fabs(a[i].x-xx)<a[i].r-eps){
- tmp[++cnt]=a[i];tmp[cnt].va=a[i].y-sqrt(a[i].r*a[i].r-(a[i].x-xx)*(a[i].x-xx));
- }
- sort(tmp+,tmp+cnt+,mycomp);
- LDB ret=;int po;
- for (int i=;i<=cnt;i=po+){
- po=i;LDB last=*tmp[i].y-tmp[i].va;
- while (po<cnt&&tmp[po+].va<last){
- po++;
- last=max(last,*tmp[po].y-tmp[po].va);
- }
- ret+=(last-tmp[i].va);
- }
- return(ret);
- }
- LDB cal(LDB l,LDB r,LDB fl,LDB fmid,LDB fr){
- return((r-l)*(fl+*fmid+fr)/);
- }
- void simpson(LDB l,LDB r,LDB fl,LDB fmid,LDB fr){
- LDB mid=(l+r)/;
- LDB mid1=(l+mid)/,mid2=(mid+r)/;
- LDB f1=getf(mid1),f2=getf(mid2);
- LDB num1=cal(l,r,fl,fmid,fr);
- LDB num2=cal(l,mid,fl,f1,fmid)+cal(mid,r,fmid,f2,fr);
- if (fabs(num1-num2)<eps) {ans+=num1;return;}
- simpson(l,mid,fl,f1,fmid);simpson(mid,r,fmid,f2,fr);
- }
- void work(){
- int po;LDB last;
- for (int i=;i<=n;i=po+){
- po=i,last=a[i].x+a[i].r;
- while (po<n&&a[po+].x-a[po+].r<last+eps){
- po++;
- last=max(last,a[po].x+a[po].r);
- }
- st=i;en=po;
- simpson(a[i].x-a[i].r,last,getf(a[i].x-a[i].r),getf((a[i].x-a[i].r+last)/),getf(last));
- }
- }
- int main(){
- scanf("%d",&n);
- for (int i=;i<=n;i++) scanf("%lf%lf%lf",&tmp[i].x,&tmp[i].y,&tmp[i].r);
- for (int i=;i<=n;i++)
- for (int j=;j<=n;j++)
- if (i!=j)
- if (dis(tmp[i],tmp[j])<tmp[j].r-tmp[i].r) {del[i]=;break;}
- int cnt=;
- for (int i=;i<=n;i++)
- if (!del[i]) a[++cnt]=tmp[i],a[cnt].va=a[cnt].x-a[cnt].r;
- n=cnt;
- sort(a+,a+n+,mycomp);
- work();
- printf("%.3lf\n",ans);
- }
simpson的更多相关文章
- 【BZOJ-1502】月下柠檬树 计算几何 + 自适应Simpson积分
1502: [NOI2005]月下柠檬树 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1017 Solved: 562[Submit][Status] ...
- 农场阳光 (simpson)
计算若干个圆与一个矩形的面积并 simpson公式 ans = ( f[l] + f[r] + 4 * f[mid] ) * (r-l) / 6 uses math; type arr=record ...
- BZOJ 1502 月下柠檬树(simpson积分)
题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=1502 题意:给出如下一棵分层的树,给出每层的高度和每个面的半径.光线是平行的,与地面夹角 ...
- UVa 10465 Homer Simpson (枚举)
10465 - Homer Simpson Time limit: 3.000 seconds http://uva.onlinejudge.org/index.php?option=com_onli ...
- 10465 - Homer Simpson 贪心
Homer Simpson Time Limit: 3 seconds Memory Limit: 32 MB Home ...
- [BZOJ 2178] 圆的面积并 【Simpson积分】
题目链接:BZOJ - 2178 题目分析 用Simpson积分,将圆按照 x 坐标分成连续的一些段,分别用 Simpson 求. 注意:1)Eps要设成 1e-13 2)要去掉被其他圆包含的圆. ...
- [BZOJ 1502] [NOI2005] 月下柠檬树 【Simpson积分】
题目链接: BZOJ - 1502 题目分析 这是我做的第一道 Simpson 积分的题目.Simpson 积分是一种用 (fl + 4*fmid + fr) / 6 * (r - l) 来拟合 fl ...
- UVA 10465 Homer Simpson(dp + 完全背包)
Problem C: Homer Simpson Time Limit: 3 seconds Memory Limit: 32 MB Homer Simpson, a very smart guy, ...
- hdu 1724 Ellipse simpson积分
/* hdu 1724 Ellipse simpson积分 求椭圆的部分面积 simpson积分法 http://zh.wikipedia.org/zh-tw/%E8%BE%9B%E6%99%AE%E ...
随机推荐
- iOS: 在iPhone和Apple Watch之间共享数据: App Groups
我们可以在iPhone和Apple Watch间通过app groups来共享数据.方法如下: 首先要在dev center添加一个新的 app group: 接下来创建一个新的single view ...
- 利用rebase来压缩多次提交
我们可以用Git merge –squash来将分支中多次提交合并到master后,只保留一次提交历史.但是有些提交到github远程仓库中的commit信息如何合并呢? 历史记录 首先我们查看一下m ...
- Navisworks Api Quantification
Quantification 国外有的叫定量 我们国内一些施工方叫工程量. 通过TakeOff API的开发者有机会获得更多的数据和数据可通过图形用户界面. 1 添加Navisworks的Api ...
- 设计模式(十四)模板方法模式(Template Pattern)
一.引言 提到模板,大家肯定不免想到生活中的“简历模板”.“论文模板”.“Word中模版文件”等,在现实生活中,模板的概念就是——有一个规定的格式,然后每个人都可以根据自己的需求或情况去更新它,例如简 ...
- 【夯实PHP基础】php开发时遇到白页的调试方法
本文地址 分享提纲: 1. 设置报错报错级别,显示报错 2. 白页的可能原因 1.[设置报错报错级别,显示报错] php开发时,访问地址也对,但就是不出来页面,显示的是 白的页面,所以就可 ...
- ASP.NET MVC搭建项目后台UI框架—1、后台主框架
目录 ASP.NET MVC搭建项目后台UI框架—1.后台主框架 ASP.NET MVC搭建项目后台UI框架—2.菜单特效 ASP.NET MVC搭建项目后台UI框架—3.面板折叠和展开 ASP.NE ...
- dubbo连接zookeeper注册中心因为断网导致线程无限等待问题【转】
最近维护的系统切换了网络环境,由联通换成了电信网络,因为某些过滤规则导致系统连不上zookeeper服务器(应用系统机器在深圳,网络为电信线路,zookeeper服务器在北京,网络为联通线路),因为我 ...
- java静态修饰符static的使用
class Person { private String name; private int age; /* * 假设每个Person对象的国籍都一样, * 那么每次调用都要赋值就会不合理. * 使 ...
- iOS 对模型对象进行归档
归档是指一种形式的序列化,专门编写用于保存数据的任何对象都应该支持归档.使用对模型对象进行归档的技术可以轻松将复杂的对象写入文件,然后再从中读取它们. 只要在类中实现的每个属性都是标量或者都是遵循NS ...
- ArcGIS Engine开发之地图导出
关于地图导出的方法有很多,但是核心技术就那么一点.下面是从项目实战中总结的一部分地图导出的方法:(以全域导出和区域导出为例) 1.由于地图导出用到的函数和方法容易重复,增加了工作量故首先将其进行封装成 ...