隐藏页面特效

LA3485 Bridge[(辛普森自适应)微积分]

做此题完全是为了练积分。

【普通求导版】

 Select Code

#include<cstdio> #include<cmath> using namespace std; typedef double real; const real eps=1e-5; int cas,Cas,D,H,B,L; real F(real a,real x){ real a2=a*a; real x2=x*x; return (x*sqrt(a2+x2)+a2*log(fabs(x+sqrt(a2+x2))))/2.0; } real calc(real w,real h){ real a=4.0*h/(w*w); real A=1.0/(2.0*a); return (F(A,w/2.0)-F(A,0))*4.0*a; } int main(){ for(scanf("%d",&cas);cas--;){ scanf("%d%d%d%d",&D,&H,&B,&L); int n=(B+D-1)/D; real D1=(real)B/n; real L1=(real)L/n; real l=0.0,r=H,mid; while(l+eps<r){ mid=(l+r)/2.0; if(calc(D1,mid)<L1) l=mid;else r=mid; } printf("Case %d:\n%.2lf\n",++Cas,H-l); if(cas) putchar('\n'); } return 0; }

【辛普森自适应版】

 Select Code

#include<cstdio> #include<cmath> using namespace std; typedef double real; const real eps=1e-5; int cas,Cas,D,H,B,L; real a; real F(real x){ real a2=a*a; real x2=x*x; return sqrt(1.0+4.0*a2*x2); } real simpson(real a,real b){ real c=(a+b)/2.0; return (F(a)+4.0*F(c)+F(b))*(b-a)/6.0; } real asr(real a,real b,real eps,real A){ real c=(a+b)/2.0; real L=simpson(a,c),R=simpson(c,b); if(fabs(L+R-A)<=15*eps) return L+R+(L+R-A)/15; return asr(a,c,eps/2,L)+asr(c,b,eps/2,R); } real asr(real a,real b,real eps){ return asr(a,b,eps,simpson(a,b)); } real calc(real w,real h){ a=4.0*h/(w*w); return asr(0,w/2.0,eps)*2.0; } int main(){ for(scanf("%d",&cas);cas--;){ scanf("%d%d%d%d",&D,&H,&B,&L); int n=(B+D-1)/D; real D1=(real)B/n; real L1=(real)L/n; real l=0.0,r=H,mid; while(l+eps<r){ mid=(l+r)/2.0; if(calc(D1,mid)<L1) l=mid;else r=mid; } printf("Case %d:\n%.2lf\n",++Cas,H-l); if(cas) putchar('\n'); } return 0; }

__EOF__

本文作者shenben
本文链接https://www.cnblogs.com/shenben/p/6653567.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   神犇(shenben)  阅读(274)  评论(0编辑  收藏  举报
编辑推荐:
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
阅读排行:
· 10年+ .NET Coder 心语 ── 封装的思维:从隐藏、稳定开始理解其本质意义
· 地球OL攻略 —— 某应届生求职总结
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 提示词工程——AI应用必不可少的技术
点击右上角即可分享
微信分享提示