博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
洛谷P4526 【模板】自适应辛普森法2(Simpson法)
阅读量:6428 次
发布时间:2019-06-23

本文共 903 字,大约阅读时间需要 3 分钟。

题面

题解

据说这函数在\(x>15\)的时候趋近于\(0\)

据说当且仅当\(a<0\)时积分发散

所以直接套自适应\(simpson\)吧……

//minamoto#include
#define R register#define fp(i,a,b) for(R int i=(a),I=(b)+1;i
I;--i)#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)using namespace std;const double eps=1e-7;double a;inline double f(const R double &x){return pow(x,a/x-x);}inline double simpson(const R double &l,const R double &r){ return (r-l)*(f(l)+f(r)+4.0*f((l+r)/2))/6;}double query(double l,double r,double eps,double res){ double mid=(l+r)/2,ql=simpson(l,mid),qr=simpson(mid,r); if(fabs(ql+qr-res)<=15*eps)return ql+qr+(ql+qr-res)/15; return query(l,mid,eps/2,ql)+query(mid,r,eps/2,qr);}int main(){// freopen("testdata.in","r",stdin); scanf("%lf",&a); a<0?puts("orz"):printf("%.5lf\n",query(eps,20,eps,simpson(eps,20))); return 0;}

转载于:https://www.cnblogs.com/bztMinamoto/p/10516296.html

你可能感兴趣的文章
关于ASP.NET MVC4 Web API简单总结
查看>>
BGP最新的AS号:4-byte-as 转换为十进制及AS号兼容性
查看>>
Windows2008server R2 组策略批量更改本地管理员密码
查看>>
ubutnu安装geany
查看>>
webservice 之 Java CXF实战效果 RS WS(一)
查看>>
我的友情链接
查看>>
Repository 与 DAO
查看>>
Zabbix监控Windows主机
查看>>
IBM x3850 RAID5数据恢复方案及过程
查看>>
移动计算领域五大机遇:交通运输优势待挖掘
查看>>
如何把win7 旗舰版升级到sp1最新版本
查看>>
Software Enginering-------using git
查看>>
浅谈IP地址-1
查看>>
我的友情链接
查看>>
C#中的线程池使用(一)
查看>>
利用Windows Server Backup功能备份活动目录
查看>>
RAC维护手记08-ASM磁盘组信息查看常用命令
查看>>
实验08 磁盘和文件系统管理
查看>>
我的友情链接
查看>>
我的友情链接
查看>>