三分法讲解---------------------------模板

本篇章大部分出自于:https://blog.csdn.net/u012469987/article/details/50897291

先减后增有极小值，若lm<rm 且(lm对应的函数值<rm对应的函数值)那么有极小点
在[left,rm] 反之在[lm,right].

lm和rm的取值分别为[left,right] 的1/3 和 2/3；
lm=l+(r-l)/3;
rm=r-(r-l)/3;

先增后减有极大值，若lm<rm 则极大在[lm,right] 反之在[left,rm]

例题:

解析:
求6x⁷+8x⁶+7x³+5x²-y*x
y>0且x∈[0,100]的极小值
该函数的图像

所以直接套模板

#include<iostream>
#include<cstdio>
using namespace std;
int t;
double y;
double val(double x)
{
	return 6*x*x*x*x*x*x*x+8*x*x*x*x*x*x+7*x*x*x+5*x*x-y*x;
}
double  slove(double l,double r)
{
	double eps=1e-7;
	while(l+eps<r)
	{
		double lmid=l+(r-l)/3;
		double rmid=r-(r-l)/3;
		if(val(lmid)<val(rmid)){
			r=rmid;
		}
		else
		{
			l=lmid;
		}
		
	}
	return val(l);
}
int main()
{
	cin>>t;
	while(t--)
	{
		cin>>y;
		printf("%.4lf\n",slove(0,100.0));	
	}
}

题目

解析:
三分求极小值模板

#include<bits/stdc++.h>
using namespace std;
double a,b,c,x,y;
double val(double X)
{
	
    return sqrt((X-x)*(X-x)+(a*X*X+b*X+c-y)*(a*X*X+b*X+c-y));
}
double slove(double l,double r)
{
	double eps=1e-7;
	while(l+eps<r)
	{
		double lmid=l+(r-l)/3;
		double rmid=r-(r-l)/3;
		if(val(lmid)<val(rmid))
		{
			r=rmid;
		}
		else
		{
			l=lmid;
		}
	}
	return val(l);
	
}
int main()
{
	 while(cin>>a>>b>>c>>x>>y){
        printf("%.3f\n", slove(-200.0,200.0));
    }
//	printf("%.3lf\n",slove(-200.0,200.0));
}