每日一题 2019/4/11

啊，明天就要去武汉了，这个怕是要咕几天了，回来再继续吧。

今天学一下三分。

好像也没啥好学的，现在学的东西都是很好理解但是不容易想到或者推出来公式，就把板子整理一下吧，以前板子的常数不怎么优秀。

//求凸函数最大值，答案取l
while(r - l > eps){
	double mid1 = (2 * l + r) / 3, mid2 = (l + 2 * r) / 3;
	if(f(mid1) < f(mid2)) l = mid1;
	else r = mid2;
}


//求凹函数最小值，答案取r
while(r - l > eps){
	double mid1 = (2 * l + r) / 3, mid2 = (l + 2 * r) / 3;
	if(f(mid1) < f(mid2)) r = mid2;
	else l = mid1;
}

然后一个例题：P3382 【模板】三分法

int n; 
double l, r;
double a[15];

double f(double x){
	double p = 0;
	for(int i = n; i >= 0; i--) p += a[i] * pow(x, i);
	return p;
}

int main()
{
	scanf("%d %lf %lf", &n, &l, &r);
	Rep(i, n, 0) scanf("%lf", a + i);
	while(r - l > eps){
		double mid1 = (2 * l + r) / 3, mid2 = (l + 2 * r) / 3;
		if(f(mid1) < f(mid2)) l = mid1;
		else r = mid2;
	}
	printf("%.5f\n", r);
}