tz大佬的函數&&二分

題目背景

tz大佬從遙遠的喜馬拉雅山挖出了一個N階的函數，爲了研究這個神奇的函數，tzdalao把函數拉到了編程俱樂部。

題目描述

這棵函數是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³-3x² - 3x + 1

思路

首先根據輸入確定函數，然後再對確定的函數求導，二分找出導數爲零的點。

代碼

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define MINX 1e-5
int N;
double l,r,mmp=0,mid,midl=0,midr=0,arr[1005];
int main()
{
    scanf("%d%lf%lf",&N,&l,&r);      //注意：輸入實數而不是整數。
    for(int i=0; i<N+1; i++)
        scanf("%lf",&arr[i]);
        while(r-l>1e-7)      //二分答案。
        {
            mid=(l+r)/2.0;
            mmp=0;
            for(int j=1; j<=N; j++)
            {
                mmp+=(j*arr[N-j]*pow(mid,j-1));     //用求導的方法找零點。
            }
            if(mmp>0)l=mid;
            else if(mmp<0) r=mid;
            if(mmp==0)break;
        }
        printf("%.5lf\n",mid);

return 0;
}