从零开始刷HDOJ(3)【HDOJ2899 - Strange fuction 】

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

Edward-Tsui

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

阅读量327

点赞数

CC 4.0 BY-SA版权

分类专栏： C++ 算法 HDOJ 从零开始刷HDOJ 水题数学文章标签： C++ ACM HDOJ 算法

本文链接：https://blog.youkuaiyun.com/edward00324258/article/details/72614431

C++ 同时被 3 个专栏收录

21 篇文章

订阅专栏

算法

18 篇文章

订阅专栏

水题

12 篇文章

订阅专栏

本文详细解析了HDOJ2899 Strangefunction问题，通过数学方法证明了该函数在[0,100]区间内为单峰函数，并给出了使用二分法寻找最小值的具体实现过程。

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

从零开始刷HDOJ(3)【HDOJ2899 - Strange function 】

题面

Strange fuction

Time limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 7321 Accepted Submission(s): 5057

Problem Description

Now, here is a function:

F (x) = 6 x 7 + 8 x 6 + 7 x 3 + 5 x 2 - y x (0 \leq x \leq 100)

$F(x) = 6x^7+8x^6+7x^3+5x^2-yx (0 ≤ x ≤100)$

Can you find the minimum value when x is between 0 and 100.

Input

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

Output

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

Sample Input

2
100
200

Sample Output

-74.4291
-178.8534

提交传送门：Submit

翻译

就是把求一下当 $y$ 等于你输入的值时 $F(x)$ 在[0, 100]上的零点。

思路

其实这个函数在[0, 100]上是单峰的。

Proof:

$\because (kx^\mu)' = k\mu (x)^{(\mu - 1)}$

$\therefore F'(x) = 42x^6 + 48x^5 + 21x^2 + 10x – y$

$\because y=ax^6在[0, 100]上单调增，y=ax^5在[0, 100]上单调增，y=ax^2在[0, 100]上单调增，y=ax在[0, 100]上单调增$

$\therefore F'(x)在[0, 100]上单调增$

既然我们可以证明 $F(x)$ 的导函数 $F'(x)$ 在定义域单调增，那么我们就可以在定义域中二分零点（具体方法参考高中数学必修一），找到的这个点就是我们要求的取值点了。最后输出 $F(ans)$ 就好了。

代码

#include <cmath>
#include <iostream>

// F(x) = 6x7 + 8x6 + 7x3 + 5x2 – xy 

double y;

inline double f(const double x)
{
    return 6 * std::pow(x, 7) + 8 * std::pow(x, 6) + 7 * std::pow(x, 3) + 5 * std::pow(x, 2) - y * x;
}

inline double F(const double x)
{
    return 42 * std::pow(x, 6) + 48 * std::pow(x, 5) + 21 * std::pow(x, 2) + 10 * x;
}

int main(int argc, char ** argv)
{
    std::ios_base::sync_with_stdio(false);
    std::cout.setf(std::ios_base::fixed);
    std::cout.precision(4);
    double ans;
    double l = 0.0, r = 100.0;
    int times = 1000;
    int T;
    std::cin >> T;
    while (T--)
    {
        std::cin >> y;
        l = 0.0;
        r = 100.0;
        times = 10000;
        double eps = 1e-10;
        while (times--)
        {
            double mid = (l + r) / 2.0;
            if (F(mid)>y)
                r = mid;
            else
                l = mid;
        }
        std::cout << f(l) << std::endl;
    }
    std::cin.get();
    std::cin.get();
    return 0;
}