三分法(洛谷3382 【模板】三分法)
阿新 • • 發佈:2017-06-03
printf log 含義 三分 tps ans 區間 bits int
如題,給出一個N次函數,保證在範圍[l,r]內存在一點x,使得[l,x]上單調增,[x,r]上單調減。試求出x的值。
輸入格式:
第一行一次包含一個正整數N和兩個實數l、r,含義如題目描述所示。
第二行包含N+1個實數,從高到低依次表示該N次函數各項的系數。
輸出格式:
輸出為一行,包含一個實數,即為x的值。四舍五入保留5位小數。
輸入輸出樣例
輸入樣例#1:
3 -0.9981 0.5 1 -3 -3 1輸出樣例#1:
-0.41421
說明
時空限制:50ms,128M
數據規模:
對於100%的數據:7<=N<=13
樣例說明:
如圖所示,紅色段即為該函數f(x)=x^3-3x^2-3x+1在區間[-0.9981,0.5]上的圖像。
當x=-0.41421時圖像位於最高點,故此時函數在[l,x]上單調增,[x,r]上單調減,故x=-0.41421,輸出-0.41421。
就是三分啦,思路看代碼~
#include<bits/stdc++.h> using namespace std; int n; double l,r; double s[20000]; double solve(double x) { double ans=0; for(int i=1;i<=n;i++) ans=(ans+s[i])*x; return ans+s[n+1]; } int main() { cin>>n; cin>>l>>r; for(int i=1;i<=n+1;i++) cin>>s[i]; while(abs(l-r)>=0.0000001) { double m1=l+(r-l)/3,m2=l+(r-l)/3*2; if(solve(m1)<solve(m2))l=m1; else r=m2; } printf("%.5f",l); return 0; }
三分法(洛谷3382 【模板】三分法)