HDOJ--2899Strange fuction!!!！二分法（三分）

最新推荐文章于 2021-06-08 20:53:46 发布

y_Cnew

最新推荐文章于 2021-06-08 20:53:46 发布

阅读量377

点赞数

CC 4.0 BY-SA版权

分类专栏：算法+ACM 文章标签： acm 算法二分法三分法

本文链接：https://blog.youkuaiyun.com/y_Cnew/article/details/78608423

算法+ACM 专栏收录该内容

23 篇文章

订阅专栏

该博客介绍了如何找到函数F(x) = 6 * x^7 + 8 * x^6 + 7 * x^3 + 5 * x^2 - y * x在0到100之间的最小值。通过分析函数性质，当y*x对导数产生影响时，可能需要借助二分法或三分法。示例中提供了二种方法的思路，并给出了针对给定y值的最小值计算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

原文链接

Problem Description

Now, here is afuction:
F(x) = 6 * x^7+8*x^6+7*x^3+5*x^2-y*x (0 <= x <=100)
Can you find the minimum value when x is between 0 and 100.

Input

The first line ofthe input contains an integer T(1<=T<=100) which means the number of testcases. Then T lines follow, each line has only one real numbers Y.(0 < Y<1e10)

Output

Just the minimumvalue (accurate up to 4 decimal places),when x is between 0 and 100.

Sample Input

100

200

Sample Output

-74.4291

-178.8534

题目大意：给你一个y，输出这个函数在0-100之间最小值

思路：二分or三分

二分：函数 F(x) = 6 * x^7+8*x^6+7*x^3+5*x^2-y*x(0 <= x <=100)没有y*x的时候是单调递增的

y*x可以影响他的导函数的正负，如果求导<0，100的时候就是最小，否则，找求导为0的点，最小

代码：

#include<stdio.h>
#include<math.h>
double y;
double f(double x)
{
	return 6*pow(x,7.0)+8*pow(x,6.0)+7*pow(x,3.0)+5*pow(x,2.0)-y*x;
}
double daoshu(double x)
{
	return 42*pow(x,6.0)+48*pow(x,5.0)+21*pow(x,2.0)+10*x-y;
}
int main()
{
	int T;
	double l,r;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%lf",&y);
		if(daoshu(100)<=y)
			printf("%.4lf\n",f(100));
		else
		{
			l=0;r=100;
			while(r-l>1e-10)
			{
				double mid=(l+r)/2;
				if(daoshu(mid)<0)
					l=mid;
				else r=mid;
			}
			printf("%.4lf\n",f(r));
		}
	}
	
	return 0;
}

三分法（前提是凸性）第一次使用，半懵半解

代码：

#include<stdio.h>
#include<math.h>
double y;
double f(double x)
{
	return 6*pow(x,7.0)+8*pow(x,6.0)+7*pow(x,3.0)+5*pow(x,2.0)-y*x;
}
int main()
{
	int T;
	double l,r,m,x1,x2;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%lf",&y);
		x1=0;x2=100;
		while(x2-x1>1e-10)
		{
			l=(2*x1+x2)/3;
			r=(x1+2*x2)/3;
			if(f(l)>f(r))	x1=l;
			else			x2=r;
		}
		printf("%.4lf\n",f(r));
	}
	
	return 0;
}