Problem 3. area
Input file: area.in
Output file: area.out
Time limit: 1 second
Mr.H最近有些无聊,就在纸上画了n 个开口向上的抛物线,并且这些抛物线与x 轴至多有一个交点。
通过这些函数,我们可以构造一个新的函数:
g(x) = min fi(x)(i=1..n)
现在,Mr.Hu 想问你在x 属于 [L,R] 这个范围内,g(x) 与x 轴围成的面积是多少。
Input
第1 行,2个整数n q,表示抛物线条数和询问次数。
接下来n 行,每行3 个数:a b c,表示fi(x) = ax^2 + bx + c。
接下来q 行,每行两个数:L R,表示一次询问。
Output
输出一行,包含一个数,表示面积,保留三位小数。
Note
. 对于30% 的数据,n = 1。
. 对于50% 的数据,n = 2。
. 对于100% 的数据,1 n 50,L R,其他所有数的绝对值不超过50。
思路:
simpson裸题
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#define inf 1e100
using namespace std;
const int N = 55;
const double eps = 1e-8;
int n, q;
double a[N], b[N], c[N];
double F( double x ) {//高度
double y = inf;
double xx = x * x;
for(int i=0; i<n; i++)
y = min(y, a[i] * xx + b[i] * x + c[i]);
return y;
}
inline double simpson( double a, double b, double fa, double fb, double fc ) {
return (fa + 4*fc + fb) * (b - a) / 6;
}
double adapt( double a, double b, double c, double fa, double fb, double fc ) {
double d = (a+c)/2, e = (c+b)/2;
double fd = F(d), fe = F(e);
double sa = simpson(a,c,fa,fc,fd);//估算左半部分
double sb = simpson(c,b,fc,fb,fe);//估算右半部分
double ss = simpson(a,b,fa,fb,fc);//估算整体
if( fabs(sa+sb-ss)<eps ) return sa+sb;//误差可接受
return adapt(a,c,d,fa,fc,fd) + adapt(c,b,e,fc,fb,fe);//误差太大就继续分
}
int main() {
freopen("area.in", "r", stdin);
freopen("area.out", "w", stdout);
scanf("%d%d", &n, &q);
for(int i = 0; i < n; i++)
scanf("%lf%lf%lf", a + i, b + i, c + i);
for(int i = 0; i < q; i++) {
double L, R;
scanf("%lf%lf", &L, &R);
printf("%.3f\n", adapt(L,R,(L+R)/2.0,F(L),F(R),F((L+R)/2.0)));
}
}