【模板】三分法

题目描述

如题，给出一个N次函数，保证在范围[l,r]内存在一点x，使得[l,x]上单调增，[x,r]上单调减。试求出x的值。

输入输出格式

输入格式：

 

第一行一次包含一个正整数N和两个实数l、r，含义如题目描述所示。

第二行包含N+1个实数，从高到低依次表示该N次函数各项的系数。

 

输出格式：

 

输出为一行，包含一个实数，即为x的值。四舍五入保留5位小数。

说明

时空限制：50ms,128M

数据规模：

对于100%的数据：7<=N<=13

思路：

这是一道三分的模板题，在这里我讲一下三分

三分是什么？？：

三分是一种优化的暴力，从他的名字就可以知道，这是一种很类似于二分的方法，通过一系列操作时复杂度由O（n）降到O（logn）

三分能干什么？

求单峰多次函数最值位置（近似值）

如何实现三分？（这个是重点）

我们先看一幅图（picture by luogu）

 

 现在我们的已知条件有峰所在的区间和函数表达式，这就很好办

首先我们取区间中点

然后取离中点很近的两个点（距离必须小于等于要求精度），并求出它们的函数值

这时候通过判断函数值的大小就可以知道此处是递增还是递减

如果递增，那点就在左区间，峰在其右，区间左端点就变成区间中点

反之亦然

一直缩小区间大小

到满足精度时，答案就呼之欲出

一点小优化：

高中的同学应该都知道有个东西叫秦九韶算法

它可以将a0+a1*x^1+a2*x^2……an^xn这个高复杂度的东西

优化为an（x+an-1（x……a1（x+a0）））这个线性的式子

复杂度下降了一个log

很好使

见代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#define rii register int i
using namespace std;
int n;
double xs[15],l,r,final,mid;
double eps=0.000001;
double dxs(double z)
{
	double ans=0;
	for(rii=n;i>=0;i--)
	{
		ans=ans*z+xs[i];
	}
	return ans; 
}
int main()
{
	cin>>n>>l>>r;
	for(rii=n;i>=0;i--)
	{
		cin>>xs[i];
	}
	eps=0.000001;
	while(1)
	{
		if(fabs(r-l)<eps)
		{
			final=r;
			break;
		}
		mid=(l+r)/2;
		if(dxs(mid+eps)>dxs(mid-eps))
		{
			l=mid;
		}
		else
		{
			r=mid;
		}
	}
	printf("%.5lf",final);
}