1. 程式人生 > >三分法(洛谷3382 【模板】三分法)

三分法(洛谷3382 【模板】三分法)

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 【模板】三分法)