洛谷 P3382 【模板】三分法

题目描述

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

输入格式

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

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

输出格式

输出为一行，包含一个实数，即为 xx 的值。若你的答案与标准答案的相对或绝对误差不超过 10^{-5} 则算正确。

输入输出样例

输入 #1

3 -0.9981 0.5
1 -3 -3 1

输出 #1

-0.41421

说明/提示

对于 100% 的数据，6≤N≤13，函数系数均在 [−100,100] 内且至多 15 位小数，∣l∣,∣r∣≤10 且至多 15 位小数。l≤r。

【样例解释】

如图所示，红色段即为该函数 f(x) = x^3 - 3 x^2 - 3x + 1 在区间 [-0.9981, 0.5]上的图像。

当 x = -0.41421时图像位于最高点，故此时函数在 [l, x] 上单调增，[x, r]上单调减，故 x = -0.41421，输出 -0.41421。

---

对于一道给定了函数和L,R两点时第一反应是二分答案，然后代入函数判断答案是否正确，然而此题没有“答案”，只能判断代入产生的值是否是极值点，因此用该点对左右两点精度要求内的最小范围进行判断，是否是极值即可。

然而朴素算法计算一元n次多项式需要经过(n+1)*n/2次乘法和n次加法，不仅代码长而且时间复杂度高，因此引用秦九韶算法。具体证明见百度百科 秦九韶算法_百度百科

double Fqjs(double x){//x即f(x)中的代数x
    double sum=0;
    for(int i=0;i<=n;i++){
        sum=sum*x+a[i];
    }
    return sum;
}

这里需要注意的是，秦九韶算法是自内而外展开的，从最高项的系数开始调用，然后依次递减。

故方法中的for循环体i的方向由输入决定。

---

const double precision = 1e-6;
while(fabs(l-r)>precision){
    mid=(l+r)/2;
    if(Fqjs(mid-precision)> Fqjs(mid+precision))
        r=mid;
    else
        l=mid;
}
cout<<mid;

接下来就是二分，在while(1){2}，1来判断精度是否达标，2来逼近答案

虽然看起来是二分实际上是三分，因为判断时用到的是

Fqjs(mid-precision)> Fqjs(mid+precision)

故mid1=mid-precision mid2=mid+precision

即三分。

---

完整代码

#include <bits/stdc++.h>
using namespace std;

int n;
double l,r,a[15],mid;
const double precision=1e-6;

double Fqjs(double x){
    double sum=0;
    for(int i=0;i<=n;i++){
        sum=sum*x+a[i];
    }
    return sum;
}

int main(){
    cin>>n>>l>>r;
    for(int i=0;i<=n;i++)
        cin>>a[i];

    while(fabs(l-r)>precision){
        mid=(l+r)/2;
        if(Fqjs(mid-precision)> Fqjs(mid+precision))r=mid;           
        else l=mid;
    }

    cout<<mid;
    return 0;
}

这道题似乎还有其他解法（？）

禁止在麻瓜面前使用魔法%%%%%%%%%%%orz orz orz