二分法求函数极大值

题目描述

这棵函数是N次函数，经过观察发现，在闭区间[l,r]内函数存在一个极大值点，请求出这个极大值点x。

输入格式：

第一行，一个正整数N和两个实数l、r，表示闭区间范围。
第二行，N+1个实数，从左到右依次表示函数的系数。

输出格式：

输出x的值，四舍五入保留5位小数。

输入样例#1：

3 -0.9981 0.5
1 -3 -3 1

输出样例#1：

-0.41421

说明

上述样例是n = 3，x^3 -3x^2 - 3x + 1

思路：高中数学知识，题中要求求函数的极大值，所以一定存在极大值，要先对其函数求导，在给定的区间内求出使导数等于0的
近似点即可。我们用一个数组来储存求导之后的系数，再用二分法来使导数的值不断逼近0.

#include<stdio.h>
double arr[100],n;  //数组arr用来储存函数求导之后的系数是多少
double ff(double x) //ff(x)函数表示带入一个函数值x之后导数的值是多少
{
    double t=0;
    double o=1;
    for(int i=n-1;i>=0;i--){
        t+=arr[i]*o;      //用t来求导数的值，导数的系数arr[i]乘以对应的x的次方
        o*=x;
    }
    return t;
}
int main()
{
    double l,r,a,v,t=0,mid;
    scanf("%lf%lf%lf",&n,&l,&r);
    v=n;
    for(int i=0;i<=n;i++){
        scanf("%lf",&a);
        arr[i]=a*v;    //将输入的系数转变为求导之后的系数
        v-=1;
    }
    while(r-l>0.0000001){  //二分法来做
        mid=(r+l)*1.0/2;
        if(ff(mid)>0)l=mid;
        else r=mid;
    }
    printf("%.5lf",mid); //%.5lf有四舍五入的功能，所以这里不需要+0.00005
    return 0;
}

不断的努力，不断的充实自己，加油！！！