题意:给你一个椭圆,问你一个区间它的积分。
思路:直接套辛普森公式就好。。。
做水题的感觉就是爽!!
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <cctype>
using namespace std;
const double eps = 1e-5;
double D, H, B, L;
double a, b;
double F(double x){
return b * sqrt(a * a - x * x) / a;
}
double simpson(double a, double b){
double c = a + (b - a) / 2;
return (F(a) + 4 * F(c) + F(b)) * (b - a) / 6;
}
double asr(double a, double b, double eps, double A){
double c = a + (b - a) / 2;
double L = simpson(a, c), R = simpson(c, b);
if(fabs(L + R - A) <= eps) return L + R + (L + R - A) / 15;
return asr(a, c, eps / 2, L) + asr(c, b, eps / 2, R);
}
double asr(double a, double b, double eps){
return asr(a, b, eps, simpson(a, b));
}
void solve(){
double l, r;
scanf("%lf%lf%lf%lf", &a, &b, &l, &r);
printf("%.3lf\n", asr(l, r, 1e-5) * 2);
}
int main()
{
int T;
scanf("%d", &T);
while(T--) solve();
return 0;
}