【算法和数据结构】三分法

1. 三分解决的问题

开门见山的说，三分法就是解决最值问题的方法。所以它能解决以下的问题
（两种说法等价）：

在定义域上，函数的导函数是单调的。
函数是单峰的。

原理很简单——拿出所求区间 [ l , r ] 中的两个值进行比较，一般选用中间值（mid = ( l + r ) / 2 ）与中间值和右侧（r）的中间值（midmid = ( mid + r ) / 2 ）进行比较：

若mid更大（小）则抛弃（midmid，r ] 这个区间
若midmid更大（小）则抛弃 [ l，mid）这个区间。

然后在剩下的区间中查找即可。
如下图：

道理很好理解，对于上图中的情况直接选择[ L，midmid ]区间继续计算即可。
下面给出一道例题（HDU 2899）

Problem Description
Now, here is a fuction:
$$F(x)=6\ast x^7+8\ast x^6+7\ast x^3+5\ast x^2-y\ast x(0<=x<=100)$$
Can you find the minimum value when x is between 0 and 100.

Input
The first line of the input contains an integer T(1<=T<=100) which means the number of test cases. Then T lines follow, each line has only one real numbers Y.(0 < Y <1e10)

Output
Just the minimum value (accurate up to 4 decimal places),when x is between 0 and 100.

#include<cstdio>
#include<cmath>
#include<iostream>
using namespace std;
double y;
double fx(double x)
{
	return 6 * pow(x, 7.0) + 8 * pow(x, 6.0) + 7 * pow(x, 3.0) + 5 * pow(x, 2.0) - y * x;
}
double dfx(double x)
{
	return 42 * pow(x, 6) + 48 * pow(x, 5) + 21 * pow(x, 2) + 10 * x - y;
}
double sanfen(double l, double r)
{
	double mid = (l + r) / 2;
	double midmid = (mid + r) / 2;
	if (r - l < 1e-6) return l;
	if (fx(mid) < fx(midmid))
	{
		return sanfen(l, midmid);
	}
	else
	{
		return sanfen(mid, r);
	}
}
int main()
{
	int T;
	scanf("%d", &T);
	while (T--)
	{
		cin >> y;
		if (dfx(0) > 0)//当时写的多余操作（恒增的函数直接返回f(0)）
		{
			printf("%.4f\n", fx(0));
		}
		else if (dfx(100) < 0)//（恒减的函数直接返回f(0)）
		{
			printf("%.4f\n", fx(100));
		}
		else
		{
			printf("%.4f\n", fx(sanfen(0,100)));
		}
	}
}