BZOJ 1502 月下柠檬树(simpson积分)
题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=1502
题意:给出如下一棵分层的树,给出每层的高度和每个面的半径。光线是平行的,与地面夹角alpha。求树在地面上投影的面积。
思路:样例的投影如下,所有面积就是这些的面积交,包括一些圆、等腰梯形。首先要计算出等腰梯形,也就是两个圆公切线与圆的交点。之后就是simpson公式:
simpson(L,R)表示图形在x=L,x=R与x轴围成图形的面积。若simpson(L,mid)+simpon(mid,R)=simpson(L,R),则直接返回simpson(L,R),否则递归计算rsimpson(L,R)。计算时我们调用的是rsimpson(L,R)。
struct point
{
double x,y;
void print()
{
printf("%lf %lf\n",x,y);
}
};
int sgn(double x)
{
if(x>EPS) return 1;
if(x<-EPS) return -1;
return 0;
}
point a[N],s[N],e[N];
int n;
double ang;
void cal(point &s,point &e,point a,point b)
{
if(sgn(a.y-b.y)==0)
{
s=a; e=b;
return;
}
double x0=a.x-a.y*(b.x-a.x)/(b.y-a.y);
double SinAng=a.y/(a.x-x0);
s.x=a.x-a.y*SinAng;
s.y=sqrt(sqr(a.y)-sqr(s.x-a.x));
e.x=b.x-b.y*SinAng;
e.y=sqrt(sqr(b.y)-sqr(e.x-b.x));
}
double f(double x)
{
double y=0;
int i;
for(i=1;i<=n+1;i++) if(fabs(x-a[i].x)<=a[i].y) y=max(y,sqrt(sqr(a[i].y)-sqr(x-a[i].x)));
for(i=1;i<=n;i++) if(a[i+1].x-a[i].x-fabs(a[i+1].y-a[i].y)>EPS&&s[i].x<=x&&x<=e[i].x)
{
y=max(y,s[i].y+(e[i].y-s[i].y)*(x-s[i].x)/(e[i].x-s[i].x));
}
return y;
}
double DFS1(double L,double R)
{
return (f(L)+4*f((L+R)/2.0)+f(R))*(R-L)/6.0;
}
double DFS(double L,double R)
{
double mid=(L+R)/2.0;
double x=DFS1(L,R),y=DFS1(L,mid),z=DFS1(mid,R);
if(fabs(x-y-z)<EPS) return y+z;
return DFS(L,mid)+DFS(mid,R);
}
int main()
{
RD(n); RD(ang); ang=1.0/tan(ang);
int i;
double h=0;
FOR1(i,n+1) RD(a[i].x),h+=a[i].x,a[i].x=h*ang;
FOR1(i,n) RD(a[i].y); a[i+1].y=0;
double L=a[1].x,R=a[n+1].x;
FOR1(i,n)
{
L=min(L,a[i].x-a[i].y);
R=max(R,a[i].x+a[i].y);
if(a[i+1].x-a[i].x-fabs(a[i+1].y-a[i].y)>EPS)
{
cal(s[i],e[i],a[i],a[i+1]);
}
}
PR(2*DFS(L,R));
}
BZOJ 1502 月下柠檬树(simpson积分)的更多相关文章
- BZOJ 1502 月下柠檬树(simpson积分)
题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=1502 题意:给出如下一棵分层的树,给出每层的高度和每个面的半径.光线是平行的,与地面夹角 ...
- 【BZOJ1502】[NOI2005]月下柠檬树 Simpson积分
[BZOJ1502][NOI2005]月下柠檬树 Description 李哲非常非常喜欢柠檬树,特别是在静静的夜晚,当天空中有一弯明月温柔地照亮地面上的景物时,他必会悠闲地坐在他亲手植下的那棵柠檬树 ...
- bzoj 1502 月下柠檬树【Simpson积分】
投影到地面之后,会发现圆形在平行光下面积和形状是不会变的,也就是所要求的图形是若干个圆和把相邻两个圆连起来的公切线所组成的. 公切线和圆间距瞎求一下就行,注意要去掉被完全覆盖的圆 然后simpson即 ...
- 【BZOJ1502】【NOI2005】月下柠檬树 simpson 积分
特别提醒:eps至少要5e-6 首先我们来研究下平行光对投影的影响. 一个二维的图形,若它与光屏平行,那么不论平行光与光屏的夹角为多少,所得图形与原图形全等的(只是位置会有影响) 通过这么一分析,我们 ...
- BZOJ 1502: [NOI2005]月下柠檬树 [辛普森积分 解析几何 圆]
1502: [NOI2005]月下柠檬树 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1070 Solved: 596[Submit][Status] ...
- [NOI2005]月下柠檬树(计算几何+积分)
题目描述 李哲非常非常喜欢柠檬树,特别是在静静的夜晚,当天空中有一弯明月温柔 地照亮地面上的景物时,他必会悠闲地坐在他亲手植下的那棵柠檬树旁,独自思 索着人生的哲理. 李哲是一个喜爱思考的孩子,当他看 ...
- 【BZOJ-1502】月下柠檬树 计算几何 + 自适应Simpson积分
1502: [NOI2005]月下柠檬树 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1017 Solved: 562[Submit][Status] ...
- 【bzoj 1502】月下柠檬树
月下柠檬树 题意 求n个圆与他们的公切线的定积分. 解法 求出圆的公切线就可以了. 特别坑的一点 : 最两端的圆,有可能会被其他的圆所包含,所以要重新求一下最左端与最右端. 比较坑的一点 : 精度要设 ...
- [BZOJ1502]月下柠檬树(自适应辛普森积分)
1502: [NOI2005]月下柠檬树 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1387 Solved: 739[Submit][Status] ...
随机推荐
- zoj 2760 How Many Shortest Path 最大流
题目链接:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1760 Given a weighted directed graph ...
- NYOJ-86 找球号(一)AC 分类: NYOJ 2014-02-02 10:45 160人阅读 评论(0) 收藏
NO.1 单纯的傻傻的代码: #include<stdio.h> long long num[100000005]={0}; int main(){ int n, m, k; scanf( ...
- Swift-4-数组和字典
// Playground - noun: a place where people can play import UIKit // 数组 字典 // 集合的可变性 赋值给var的集合是可变的mut ...
- On Explainability of Deep Neural Networks
On Explainability of Deep Neural Networks « Learning F# Functional Data Structures and Algorithms is ...
- WSDL
WSDL 彻底理解webservice SOAP WSDL WSDL 详解 http://www.cnblogs.com/hujian/p/3494064.html http://www.cnblog ...
- 解决iptables和vsftpd设置的问题
解决iptables和vsftpd设置的问题 博客分类: linux/centos/ubuntu 防火墙J#工作 解决iptables和vsftpd设置的问题 修改 vi /etc/sysconfig ...
- Ruby探针的基本实现原理
李哲 - MAY 13, 2015 语言本身 Ruby语言支持语法级别的系统,框架,甚至语言本身的方法复写,一般叫做元编程(meta programming), 此基础之上还有一些术语为mixin,方 ...
- javascript 注意事项
1.submit input 标签的 type里的submit是提交表单的按钮 提交时应注意 为 标签加上 “name=""”字段和form标签,才能达到想要的效果 2.正则表达 ...
- 【面试题041】和为s的两个数字VS和为s的连续正数序列
[面试题041]和为s的两个数字VS和为s的连续正数序列 题目一: 输入一个递增排序的数组和一个数字s,在数组中查找两个数,使得它们的和正好是s.如果有多对数字的和等于s,输出任意一对即可. ...
- iOS后台运行
http://www.cocoachina.com/bbs/read.php?tid=149564 文一 我从苹果文档中得知,一般的应用在进入后台的时候可以获取一定时间来运行相关任务,也就是说可以在后 ...