[BZOJ1502]月下柠檬树(自适应辛普森积分)

1502: [NOI2005]月下柠檬树

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 1387  Solved: 739
[Submit][Status][Discuss]

Description

李哲非常非常喜欢柠檬树,特别是在静静的夜晚,当天空中有一弯明月温柔地照亮地面上的景物时,他必会悠闲地
坐在他亲手植下的那棵柠檬树旁,独自思索着人生的哲理。李哲是一个喜爱思考的孩子,当他看到在月光的照射下
柠檬树投在地面上的影子是如此的清晰,马上想到了一个问题:树影的面积是多大呢?李哲知道,直接测量面积是
很难的,他想用几何的方法算,因为他对这棵柠檬树的形状了解得非常清楚,而且想好了简化的方法。李哲将整棵
柠檬树分成了n 层,由下向上依次将层编号为1,2,…,n。从第1到n-1 层,每层都是一个圆台型,第n 层(最上面一
层)是圆锥型。对于圆台型,其上下底面都是水平的圆。对于相邻的两个圆台,上层的下底面和下层的上底面重合
。第n 层(最上面一层)圆锥的底面就是第n-1 层圆台的上底面。所有的底面的圆心(包括树顶)处在同一条与地面垂
直的直线上。李哲知道每一层的高度为h1,h2,…,hn,第1 层圆台的下底面距地面的高度为h0,以及每层的下底面
的圆的半径r1,r2,…,rn。李哲用熟知的方法测出了月亮的光线与地面的夹角为alpha。
为了便于计算,假设月亮的光线是平行光,且地面是水平的,在计算时忽略树干所产生的影子。
李哲当然会算了,但是他希望你也来练练手

Input

第1行包含一个整数n和一个实数alpha,表示柠檬树的层数和月亮的光线与地面夹角(单位为弧度)。
第2行包含n+1个实数h0,h1,h2,…,hn,表示树离地的高度和每层的高度。
第3行包含n个实数r1,r2,…,rn,表示柠檬树每层下底面的圆的半径。
上述输入文件中的数据,同一行相邻的两个数之间用一个空格分隔。
输入的所有实数的小数点后可能包含1至10位有效数字。
1≤n≤500,0.3<alpha<π/2,0<hi≤100,0<ri≤100

Output

输出1个实数,表示树影的面积。四舍五入保留两位小数。

Sample Input

2 0.7853981633
10.0 10.00 10.00
4.00 5.00

Sample Output

171.97

HINT

Source

[Submit][Status][Discuss]

这个题的正解是各种分类讨论求面积,但是可以用辛普森积分骗分。
关于辛普森积分的推导可以看 https://www.zhihu.com/question/47728235
其实这个东西很好卡掉,但是对于很多数据还是可以较为精确地出解的。

首先发现要求的东西是轴对称的,这意味着我们可以只求x坐标上面的函数区域没。这个图像的解析式还是很难求,但是对于每个横坐标上对应的函数值可以比较轻松地求出,这个时候就可以用Simpson积分了。
普通辛普森积分根本不可能出解,因为仅仅一个二次函数不可能拟合这么复杂的图像。但是要分成很多份的话时间复杂度又过高,这个时候就需要引入自适应:如果要对[a,b]积分,可以先对[a,b]拟合,再对[a,(a+b)/2]和[(a+b)/2,b]拟合,如果两个相差很小则直接退出,否则继续递归。
个人认为这个方法在这里可以成功很大程度是因为图像中有很多直线围城的区域,而自适应Simpson积分在直线积分方面的表现应该是很不错的。

具体做法看 https://www.cnblogs.com/DaD3zZ-Beyonder/p/5676841.html
上面讲的很清楚了,直接作为模板即可。

 1 #include<cmath>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #define rep(i,l,r) for (int i=l; i<=r; i++)
 5 typedef double db;
 6 using namespace std;
 7 
 8 const int N=1010;
 9 double eps=1e-5,inf=1e12;
10 int n,num;
11 db alpha;
12 struct P{ db x,y; P(db X=0,db Y=0){ x=X; y=Y; } };
13 struct Ci{ db r; P c; Ci (P C=(P(0,0)),db R=0):c(C),r(R){}; }C[N];
14 
15 struct Li{
16     P s,t; db k,b;
17     Li (P S=P(0,0),P T=P(0,0)){
18         s=S,t=T; if (s.x>t.x) swap(s,t);
19         k=(s.y-t.y)/(s.x-t.x); b=s.y-k*s.x;
20     }
21     db f(db x){ return k*x+b; }
22 }l[N];
23 
24 int dcmp(db x){ if (fabs(x)<eps) return 0; return (x<0) ? -1 : 1; } 
25 db F(db x){
26     db res=0;
27     rep(i,1,n){
28         db d=fabs(x-C[i].c.x);
29         if (dcmp(d-C[i].r)>0) continue;
30         db len=2*sqrt(C[i].r*C[i].r-d*d);
31         res=max(res,len);
32     }
33     rep(i,1,num) if (x>l[i].s.x && x<=l[i].t.x) res=max(res,2*l[i].f(x));
34     return res;
35 }
36 
37 db calc(db l,db r){ db mid=(l+r)/2; return (F(l)+F(r)+F(mid)*4)*(r-l)/6; }
38 db Simpson(db l,db r,db now){
39     db mid=(l+r)/2,x=calc(l,mid),y=calc(mid,r);
40     if (!dcmp(now-x-y)) return now; else return Simpson(l,mid,x)+Simpson(mid,r,y);
41 }
42 
43 void solve(){
44     db L=inf,R=-inf;
45     rep(i,1,n+1) L=min(L,C[i].c.x-C[i].r),R=max(R,C[i].c.x+C[i].r);
46     rep(i,1,n){
47         db d=C[i+1].c.x-C[i].c.x;
48         if (dcmp(d-fabs(C[i].r-C[i+1].r))<0) continue;
49         db sina=(C[i].r-C[i+1].r)/d,cosa=sqrt(1-sina*sina);
50         l[++num]=Li(P(C[i].c.x+C[i].r*sina,C[i].r*cosa),P(C[i+1].c.x+C[i+1].r*sina,C[i+1].r*cosa));
51     }
52     printf("%.2lf\n",Simpson(L,R,calc(L,R)));
53 }
54 
55 int main(){
56     scanf("%d%lf",&n,&alpha); db h,r;
57     rep(i,1,n+1)
58         scanf("%lf",&h),C[i]=Ci((P((h/tan(alpha))+C[i-1].c.x,0)),0);
59     rep(i,1,n) scanf("%lf",&r),C[i].r=r;
60     solve();
61     return 0;
62 }

 

转载于:https://www.cnblogs.com/HocRiser/p/8884527.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/398210.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

擎天出口退税软件_新版申报系统退税申报注意事项

前言&#xff1a;随着税务信息化建设“后金三”时代的不断深入,2019年年底国家税务总局启动出口退税管理系统整合(进入金税三期)项目&#xff0c;经过广东、大连两个地区试点&#xff0c;即将于今年年底前在全国完成金税三期审核系统的上线工作&#xff0c;至此&#xff0c;出口…

词性标注

4.10 词性标注 词性用来描写叙述一个词在上下文中的作用。比如描写叙述一个概念的词叫做名词&#xff0c;在下文引用这个名词的词叫做代词。有的词性常常会出现一些新的词&#xff0c;比如名词&#xff0c;这种词性叫做开放式词性。另外一些词性中的词比較固定&#xff0c;比如…

MVCC浅析(转)

在并发读写数据库时&#xff0c;读操作可能会不一致的数据&#xff08;脏读&#xff09;。为了避免这种情况&#xff0c;需要实现数据库的并发访问控制&#xff0c;最简单的方式就是加锁访问。由于&#xff0c;加锁会将读写操作串行化&#xff0c;所以不会出现不一致的状态。但…

关于单链表,二叉树,图,查找和排序的软件编程

课程名称&#xff1a;计算机软件 使用软件&#xff1a;devcpp 注意&#xff1a;这里列出了关于单链表&#xff0c;二叉树&#xff0c;图&#xff0c;查找和排序的编程&#xff0c;全部程序由博主一人编写&#xff0c;会有瑕疵&#xff0c;谨慎使用。 1.单链表 要求&#xff1a;…

pcie 的function_PCIe扫盲——BDF与配置空间

前面的文章中介绍过&#xff0c;每一个PCIe设备可以只有一个功能(Function)&#xff0c;即Fun0。也可以拥有最多8个功能&#xff0c;即多功能设备(Multi-Fun)。不管这个PCIe设备拥有多少个功能&#xff0c;其每一个功能都有一个唯一独立的配置空间(Configuration Space)与之对应…

微信二维码扫描下载APK

前几天给客户制作的app需要上线&#xff0c;生成二维码扫描进行下载&#xff0c;把生成好的apk挂在服务器端&#xff0c;将地址复制下来&#xff0c;通过草料二维码&#xff08;http://cli.im/&#xff09;生成一个二维码它需要一个应用宝ID。上传到应用宝&#xff0c;光审批就…

mysql 字段 as_mysql 字段as详解及实例代码

mysql 字段使用as在mysql中&#xff0c;select查询可以使用AS关键字为查询的字段起一个别名&#xff0c;该别名用作表达式的列名&#xff0c;并且别名可以在GROUP BY&#xff0c;ORDER BY或HAVING等语句中使用。例如&#xff1a;SELECT CONCAT(last_name,, ,first_name) AS ful…

dubbo接口快速测试技巧

在分布式系统的开发中&#xff0c;用到了dubbozookeeper技术&#xff0c;最近遇到一个问题&#xff0c;产品上线后&#xff0c;我负责的模块出了问题&#xff0c;某个bean中某个字段的值一直为null&#xff0c;而这个bean是我调用注册在zookeeper上的一个服务查询到的&#xff…

跳出多重循环 JS

关于如何跳出多重循环loop: //循环标记for(var key in jsonObj){for(var i0;i<jsonObj[key].length;i){if(jsonObj[key][i].password123456){break loop;//跳出双重循环 }}}转载于:https://blog.51cto.com/9381188/1790457

webservice引用spring的bean

1 <jaxws:endpoint address"/test/webservice" implementor"#testBean" /> 这行代码里面的#号表示webservice自动装配到spring的bean。转载于:https://www.cnblogs.com/xiluhua/p/4472544.html

mysql 5.0 数据库_mysql5.0常用命令

MySQL常用操作基本操作&#xff0c;以下都是MySQL5.0下测试通过首先说明下&#xff0c;记住在每个命令结束时加上&#xff1b;(分号)1.导出整个数据库mysqldump -u 用户名 -p --default-character-setlatin1 数据库名 > 导出的文件名(数据库默认编码是latin1)mysqldump -u w…

IISASP.NET 站点IP跳转到域名

前言&#xff1a;先到微软的 https://www.iis.net/downloads/microsoft/url-rewrite 下载URL Rewrite 目标&#xff1a;输入ip跳转到域名所在的网站 比如58的115.159.231.173 跳转到https://passport.58.com/login 先看下58的例子 我们在地址栏输入ip之后 箭头指向的地方是跳转…

java基础五

继承 1. 继承含义 在面向对象编程中&#xff0c;可以通过扩展一个已有的类&#xff0c;并继承该类的属性和行为&#xff0c;来创建一个新的类&#xff0c;这种方式称为继承&#xff08;inheritance&#xff09;. 2. 继承的优点 A&#xff0e;代码的可重用性 B&#xff0e;子类可…

Mac OS X 下Java开发环境配置

首先自己又装的jdk1.8版本目录&#xff1a;/Library/Java/JavaVirtualMachines/jdk1.8.0_91.jdk/Contents/Home待添加。。。转载于:https://blog.51cto.com/jiyanle/1790537

mysql emoji表情_让MySQL支持Emoji表情 mysql 5.6

最近在做微信相关的项目&#xff0c;其中MySQL 要存储emoji表情&#xff0c;因此发现我们常用的utf8 字符集根本无法存储表情。网上有不少替代方案。本人还是采用了修改MySQL字符集的方案简单快捷。首先将我们数据库默认字符集由utf8 更改为utf8mb4&#xff0c;对应的表默认字符…

高效开发iOS系列 -- 那些不为人知的KVC

我的简书地址&#xff1a;http://www.jianshu.com/p/a6a0abac1c4a valueForKeyPath 本篇来解说一下那些不为人知&#xff0c;也常常被忽略掉&#xff0c;而且非常有用的KVC干货小技巧 获取数组里的,最大、最小、平均、求和 NSArray *array ["1","3",2,9.5…

::selection

改变浏览器文字选中背景 ::selection CSS的伪类选择器 只有一小部分CSS属性可以用于::selection选择器&#xff1a; color, background-color, cursor, outline, text-decoration, text-emphasis-color和text-shadow。 要特别注意的是&#xff0c;background-image会如同其他属…

B00013 字符串哈希函数

哈希算法将任意长度的二进制值映射为较短的固定长度的二进制值&#xff0c;这个小的二进制值称为哈希值。 哈希函数用途广泛&#xff0c;这个程序给出了绝大多数常用的哈希函数。源程序来自&#xff1a;哈希算法_百度百科。 程序员可以根据自己的需要取用这些代码。 需要注意的…

*hdu5632Rikka with Array

$n \leq 10^300$&#xff0c;问所有$i<j$且$f_i>f_j$的$(i,j),1 \leq i \leq n,1 \leq j \leq n$数量。对某个数取模。 $f(a,b,0/1/2/3)$--前$a$位&#xff0c;两个数相差$b$个1&#xff0c;$b$可以为负数&#xff0c;然后四种情况&#xff1a;$i<j<n$&#xff0c;…

java正则表达式提取字符串中的数字

http://stackoverflow.com/questions/2367381/extract-numbers-from-a-string-java Pattern p Pattern.compile(\\d); 使用这个会有空字符串出现 Pattern p Pattern.compile("-?\\d"); 这是OK,转载于:https://www.cnblogs.com/pengxinglove/p/4476634.html