[NOI2005]月下柠檬树
题意
<body> <center> <div style="width:90%; text-align:left"> <img src="image/logo.png"> </div> <table width="96%"> <tbody><tr align="center" class="hd" valign="top"> <th><a href="faqs.php">F.A.Qs</a></th> <th><a href="./">Home</a></th> <th><a href="./bbs.php">Discuss</a></th> <th><a href="./problemset.php">ProblemSet</a></th> <th><a href="./status.php">Status</a></th> <th><a href="./ranklist.php">Ranklist</a></th> <th><a href="./contest.php">Contest</a></th> <th><a href="http://begin.lydsy.com/JudgeOnline">入门OJ</a></th> <th><a href="./modifypage.php"><b>ModifyUser</b></a> <a href="userinfo.php?user=autoint"> <font color="red">autoint</font></a></th><th><a href="logout.php">Logout</a></th> <th><a href="./donation.php"><font color="red">捐赠本站</font></a></th> </tr> </tbody></table> </center> <title>Problem 1502. -- [NOI2005]月下柠檬树</title><center><h2>1502: [NOI2005]月下柠檬树</h2><span class="green">Time Limit: </span>5 Sec <span class="green">Memory Limit: </span>64 MB<br><span class="green">Submit: </span>1534 <span class="green">Solved: </span>813<br>[<a href="submitpage.php?id=1502">Submit</a>][<a href="problemstatus.php?id=1502">Status</a>][<a href="bbs.php?id=1502">Discuss</a>]</center><h2>Description</h2><div class="content"><div>李哲非常非常喜欢柠檬树，特别是在静静的夜晚，当天空中有一弯明月温柔地照亮地面上的景物时，他必会悠闲地</div> <div>坐在他亲手植下的那棵柠檬树旁，独自思索着人生的哲理。李哲是一个喜爱思考的孩子，当他看到在月光的照射下</div> <div>柠檬树投在地面上的影子是如此的清晰，马上想到了一个问题：树影的面积是多大呢？李哲知道，直接测量面积是</div> <div>很难的，他想用几何的方法算，因为他对这棵柠檬树的形状了解得非常清楚，而且想好了简化的方法。李哲将整棵</div> <div>柠檬树分成了n 层，由下向上依次将层编号为1,2,…,n。从第1到n-1 层，每层都是一个圆台型，第n 层(最上面一</div> <div>层)是圆锥型。对于圆台型，其上下底面都是水平的圆。对于相邻的两个圆台，上层的下底面和下层的上底面重合</div> <div>。第n 层(最上面一层)圆锥的底面就是第n-1 层圆台的上底面。所有的底面的圆心(包括树顶)处在同一条与地面垂</div> <div>直的直线上。李哲知道每一层的高度为h1,h2,…,hn，第1 层圆台的下底面距地面的高度为h0，以及每层的下底面</div> <div>的圆的半径r1,r2,…,rn。李哲用熟知的方法测出了月亮的光线与地面的夹角为alpha。</div> <div><img src="https://www.lydsy.com/JudgeOnline/upload/201802/111.png" width="407" height="217" alt=""></div> <div></div> <div>为了便于计算，假设月亮的光线是平行光，且地面是水平的，在计算时忽略树干所产生的影子。</div> <div>李哲当然会算了，但是他希望你也来练练手</div></div><h2>Input</h2><div class="content"><div> <div>第1行包含一个整数n和一个实数alpha，表示柠檬树的层数和月亮的光线与地面夹角(单位为弧度)。</div> <div>第2行包含n+1个实数h0,h1,h2,…,hn，表示树离地的高度和每层的高度。</div> <div>第3行包含n个实数r1,r2,…,rn，表示柠檬树每层下底面的圆的半径。</div> <div>上述输入文件中的数据，同一行相邻的两个数之间用一个空格分隔。</div> <div>输入的所有实数的小数点后可能包含1至10位有效数字。</div> <div>1≤n≤500，0.3<alpha<π/2，0<hi≤100，0<ri≤100</div> </div> <p></p></div><h2>Output</h2><div class="content"><p>输出1个实数，表示树影的面积。四舍五入保留两位小数。</p></div><h2>Sample Input</h2> <div class="content"><span class="sampledata">2 0.7853981633<br> 10.0 10.00 10.00<br> 4.00 5.00</span></div><h2>Sample Output</h2> <div class="content"><span class="sampledata">171.97</span></div><h2>HINT</h2> <div class="content"><p></p></div><h2>Source</h2> <div class="content"><p><a href="problemset.php?search="></a></p></div><center>[<a href="submitpage.php?id=1502">Submit</a>][<a href="problemstatus.php?id=1502">Status</a>][<a href="bbs.php?id=1502">Discuss</a>]</center><br> <a href="./"><span class="red">HOME</span></a> <a href="javascript:history.go(-1)"><span class="red">Back</span></a>分析
参照GFY的题解。
首先，做这题需要知道一点：一个圆从任意一个角度投影都永远是一个圆。
我们可以画出一个简图如下：
如图，这棵树倒影之后，有图中两个圆心p1,p2,他们的横坐标即为这颗树上他们原先的高度乘以cotΘ，而他们的半价却不会变化，因为月光是平行光，所以在圆面与地面平行时，两点间距离不会变化。
如图，倒影最终是圆和他们之间的公切线构成的图形，最右边的点可以看做是半径为eps的圆。之后，可以利用simpson积分公式计算，simpson(l,r)=(f(l)+f(r)+4f(mid))(r-l)/6,若是精度差距大可以继续递归，注意：本题的eps要1e-6以下才能过。
时间复杂度？不知道
代码
#include<bits/stdc++.h>
#define rg register
#define il inline
#define co const
typedef long long ll;
co double eps=1e-6;
co int N=501;
struct Point{
double x,y;
}s[N],e[N],a[N],b[N];
int n;
il double sqr(double x) {return x*x;}
void cal(Point&s,Point&e,Point a,Point b){
if(fabs(a.y-b.y)<eps) return s=a,e=b,void();
double k=(b.y-a.y)/(b.x-a.x);
s.x=a.x-k*a.y,s.y=sqrt(sqr(a.y)-sqr(a.x-s.x));
e.x=b.x-k*b.y,e.y=sqrt(sqr(b.y)-sqr(b.x-e.x));
// s=(Point){a.x-k*a.y,sqrt(sqr(a.y)-sqr(a.x-s.x))};
// e=(Point){b.x-k*b.y,sqrt(sqr(b.y)-sqr(b.x-e.x))};
}
double f(double x){
double y=0;
for(int i=1;i<=n+1;++i)if(fabs(x-a[i].x)<=a[i].y)
y=std::max(y,sqrt(sqr(a[i].y)-sqr(x-a[i].x)));
for(int i=1;i<=n;++i)if(a[i+1].x-a[i].x-fabs(a[i].y-a[i+1].y)>eps&&s[i].x<=x&&x<=e[i].x)
y=std::max(y,s[i].y+(x-s[i].x)*(e[i].y-s[i].y)/(e[i].x-s[i].x));
return y;
}
double work(double L,double R){
double mid=(L+R)/2;
return (f(L)+f(R)+f(mid)*4)*(R-L)/6;
}
double simpson(double L,double R){
double mid=(L+R)/2;
double x1=work(L,mid),x2=work(mid,R),x3=work(L,R);
if(fabs(x1+x2-x3)<eps) return x1+x2;
else return simpson(L,mid)+simpson(mid,R);
}
int main(){
// freopen(".in","r",stdin),freopen(".out","w",stdout);
double theta;
scanf("%d%lf",&n,&theta);
theta=1/tan(theta);
double h=0;
for(int i=1;i<=n+1;++i){
scanf("%lf",&a[i].x);
h+=a[i].x,a[i].x=h*theta;
}
for(int i=1;i<=n;++i) scanf("%lf",&a[i].y);
a[n+1].y=a[n+2].y=0;
double L=a[1].x,R=a[n+1].x;
for(int i=1;i<=n;++i){
L=std::min(L,a[i].x-a[i].y),R=std::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) // not in
cal(s[i],e[i],a[i],a[i+1]);
}
printf("%.2lf\n",2*simpson(L,R));
return 0;
}
静渊以有谋，疏通而知事。
