洛谷 P4525 & P4526 [模板] 自适应辛普森积分

题目：https://www.luogu.org/problemnew/show/P4525

https://www.luogu.org/problemnew/show/P4526

学习辛普森积分：https://blog.csdn.net/VictoryCzt/article/details/80660113

代码如下：

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<cmath>

#define eps 1e-8

using namespace std;

typedef double db;

db a,b,c,d,L,R;

db f(db x){return (c*x+d)/(a*x+b);}

db simp(db a,db b){return (b-a)/*(f(a)+*f((a+b)/)+f(b));}

db asr(db l,db r,db lst)

{

  db mid=(l+r)/;

  db ls=simp(l,mid),rs=simp(mid,r);

  if(fabs(ls+rs-lst)<eps)return ls+rs;

  else return asr(l,mid,ls)+asr(mid,r,rs);

}

db asme(db l,db r){return asr(l,r,simp(l,r));}

int main()

{

  scanf("%lf%lf%lf%lf%lf%lf",&a,&b,&c,&d,&L,&R);

  printf("%.6f\n",asme(L,R));

  return ;

}

P4525

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<cmath>

#define eps 1e-8

using namespace std;

typedef double db;

db a;

db f(db x){return pow(x,a/x-x);}

db simp(db l,db r){return (r-l)/*(f(l)+*f((l+r)/)+f(r));}

db asr(db l,db r,db lst)

{

  db mid=(l+r)/;

  db ls=simp(l,mid),rs=simp(mid,r);

  if(fabs(ls+rs-lst)<eps)return ls+rs;

  return asr(l,mid,ls)+asr(mid,r,rs);

}

db asme(db l,db r){return asr(l,r,simp(l,r));}

int main()

{

  scanf("%lf",&a);

  if(a<){puts("orz"); return ;}

  printf("%.5f\n",asme(eps,));

  return ;

}

P4526