本篇章大部分出自于: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]
例题:
解析:
求6x7+8x6+7x3+5x2-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));
}