BestCoder18 1002.Math Problem(hdu 5105) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5105

题目意思:给出一个6个实数:a, b, c, d, l, r。通过在[l, r]中取数 x,使得这个函数 f(x)= |a∗x3+b∗x2+c∗x+d| 最大。

我一开始做的时候,很天真的认为数据量这么小,一个一个试,暴力搜就肯定能得到答案啦。但是一个很严重的问题是,x 没有说是实数还是整数,所以枚举根本不可行。

  于是就联想到应该使用高中求一元三次方程的方法来做。我当时是直接从 f(l), f(r),f(x1),f(x2)(x1、x2代表方程两个根)中选择最大的来输出,但是没有考虑是否在[l, r]范围内,竟然过了,那个开心啊。可想而知最后被 hack 了,连终评都到达不了,泪~~~

  后来认认真真复习了具体求解方法,于是很正规地做:对ax^3 + bx^2 + c*x + d求导,得到一元二次方程 3ax^2 + 2bx + c,用delta(那个高中经常用到的三角形Δ)讨论根的情况,画出二次函数图象,根据单调性再画出原函数大致走向,再讨论根是否在区间[l, r]上来选择应该输出哪个f()。前前后后做了31次。那 30 次当中,知道自己很多的不足:

(1)竟然用 abs 来求浮点数绝对值,实际上是用 fabs!

(2)浮点数跟 0 比较是不可行的,要设定一个可接受的误差 eps,例如为1e-9,如果 < eps 就代表是 < 0 的。

(3)在做 delta 运算之前,漏了 a == 0(对应代码中 fabs(a) < eps) 的判断,这关乎到根存不存在的问题: 如果 b < eps(表示 b == 0),就代表无解啦,否则就有一个根 x = (-c) / (2b)

(4)最最致命的一点,就是按传统解题步骤做!今天在课上终于想明白为什么这样中规中矩做是不可行的,两天的努力没有白费了。试想下,这样做的情况非常多!

  举个例子吧,假设我们已经知道 a != 0 啦,很自然地进行delta讨论啦,假设讨论到 delta > 0(其实delta = 0可以归为该类情况) 的情况,也就是有x1, x2两个根啦。然后根据在 x 坐标轴以上表示原函数递增,以下原函数递减,最后对应到原函数(ax^3 + bx^2 + c*x + d)画出大致走向,再根据 x 的位置来判断是否应该计算f(x)(包含f(x1), f(x2)),最后选出f(x), f(l),f(r)最大的那个输出。真的是复杂到不能再复杂啦。

我当时就在想为什么是错的,因为原函数不是单纯的 ax^3 + bx^2 + c*x + d,而是这个东西的绝对值!绝对值意味着这个东西的最小值或者最大值就是答案。而我常规做的,是非绝对值的情形,所以势必会漏了不容易觉察的情况。

那么最简单,最靠谱,最正确的方法就是,除了那种 a == 0 && b == 0 没根(不用算x)的情况,其他都需要算f(x),然后选出最大的,前提是 x 在 [l, r]之间。

  这里补充说明一点,为什么我说的进行到 delta 判断中,可以将 delta < 0 与 delta >= 0 一并讨论,因为始终只会算出 x 在 [l, r]的 f(x),也就是不在或者不存在(对应delta < 0)的话就不计算(函数是返回-1),那么delta < 0 其实就是选择f(l)、 f(r)较大的一个输出。至于delta = 0,虽然只有一个根,即x1, x2相等,也不影响合并。

辛苦读者看我这篇长篇大论,不过通过这题真的学到很多,多想一个为什么一定会有所进步的!^_^ !

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cstdlib>
 4 #include <cstring>
 5 #include <cmath>
 6 using namespace std;
 7
 8 const double eps = 1e-9;
 9 double a, b, c, d, l, r;
10
11 double get_max(double x, double y)
12 {
13     if (x > y)
14         return x;
15     return y;
16 }
17
18 double cal(double x)
19 {
20     if (x >= l && x <= r)               // 前提是在[l, r] 之间
21         return fabs(a*x*x*x + b*x*x + c*x + d);
22     return -1;
23 }
24
25 int main()
26 {
27     while (scanf("%lf%lf%lf%lf%lf%lf", &a, &b, &c, &d, &l, &r) != EOF)
28     {
29         double fl = cal(l);
30         double fr = cal(r);
31         double ans = get_max(fl, fr);
32
33         if (fabs(a) < eps)        // a == 0
34         {
35             if (fabs(b) < eps)    // b == 0,无解
36                 printf("%.2lf\n", ans);
37             else
38             {
39                 double x = -c/(2*b);    // 导数为一次函数,原函数为二次函数
40                 printf("%.2lf\n", get_max(ans, cal(x)));
41             }
42         }
43         else      // delta 环节
44         {
45             double delta = 4*b*b - 12*a*c;
46             delta = sqrt(delta);
47             double x1 = (-2*b - delta) / (6*a);
48             double x2 = (-2*b + delta) / (6*a);
49             double tans = get_max(cal(x1), cal(x2));
50             printf("%.2lf\n", get_max(ans, tans));
51         }
52     }
53     return 0;
54 }
时间: 2024-10-17 20:18:50

BestCoder18 1002.Math Problem(hdu 5105) 解题报告的相关文章

BestCoder12 1002.Help him(hdu 5059) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5059 题目意思:就是输入一行不多于 100 的字符串(除了'\n' 和 '\r' 的任意字符),问是否是合法的整数,如果是的话问是否在[a, b] 范围内,是则输出 YES,否则输出 NO 合法的整数:(1)非负整数:只有数字,没有前导0 (2)负数:负号后面只能跟着数字,负号前没有任何字符 首先这条题感觉不是出得太好,不过都是硬着头皮学人家做啦.反正一些很变态的数据可能还是过不了,但却AC的. 模

BestCoder3 1002 BestCoder Sequence(hdu 4908) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4908 题目意思:给出 一个从1~N 的排列你和指定这个排列中的一个中位数m,从这个排列中找出长度为奇数,中位数是m的子序列有多少个. 我的做法被discuss 中的测试数据一下子就否定了. 这个是别人的做法,暂时留下来,有些地方还没真正弄懂,应该是缺了部分的知识没有学到... 留着先: (1)http://blog.csdn.net/hcbbt/article/details/38377815 (2

hdu 2112 HDU Today 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2112 题目意思:又是求最短路的,不过结合埋字符串来考查. 受之前1004 Let the Balloon Rise 学到用map 的解法做之后,有点蠢蠢欲动,当时见到要用字典树做有点吓坏了(之前看过下,非一般难理解,所以暂时放下了).于是,死就死吧,硬住头皮用map做,反反复复修改终于过了. 首先是memory limit exceeded,因为无读清题目意思,直接开10000 * 10000的数组

BestCoder17 1002.Select(hdu 5101) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5101 题目意思:给出 n 个 classes 和 Dudu 的 IQ(为k),每个classes 都有相应的人,每个人又有相应的IQ.现在要求从这些classes挑选两个人,满足IQ之和 > k,不过要满足这两个人不能来自同一个class的. 解题思路不难想出,直接所有人两两之和 > k - 同一个class 两两之和 > k  就是答案了. 不过很容易超时!!!! 用二分就可以过了.二分有

BestCoder16 1002.Revenge of LIS II(hdu 5087) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5087 题目意思:找出第二个最长递增子序列,输出长度.就是说,假如序列为 1 1 2,第二长递增子序列是1 2(下标为2 3),而第一长递增子序列也是(下标为 1 3). 我一开始天真的以为,还是利用求最长递增子序列的算法啦,第二长不就是对dp 数组sort完后(从小到大)的dp[cnt-1] 啦,接着就呵呵啦----= = 题解说,要加多一个 dp 数组,以便对当前下标值为 i 的数 a[i] 为结

BestCoder22 1002.NPY and arithmetic progression(hdu 5143) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5143 题目意思:给出 1, 2, 3, 4 的数量,分别为a1, a2, a3, a4,问是否在每个数只使用一次的前提下,分成若干部分,每部分数的长度 >= 3且满足是等差数列.可以的话输出 Yes ,否则 No . 比赛的时候过了pretest,那个开心啊--然后跟 XX 兽讨论了之后,才发现自己理解错了题目= = 要用到深搜,问组合嘛--组合就是有可能是(1, 2, 3, 4).(1, 2, 3

BestCoder18 1001.Alexandra and Prime Numbers(hdu 5108) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5108 题目意思:给出一个数正整数 N,N <= 1e9,现在需要找出一个最少的正整数 M,使得 N/M 是素数.如果找不到就输出0. 一开始有想过将所有 <= 1e9 的素数求出来的,不过绝对超时就放弃了:然后就开始从题目中挖掘简便的处理方法.受到求素数的方法启发,枚举的因子 i 如果在 i * i <= N 之内都没有找到符合条件的素数,那么那些 > N 的因子就更不可能了.于是时间

BestCoder15 1002.Instruction(hdu 5083) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5083 题目意思:如果给出 instruction 就需要输出对应的 16-bit binary code,给出16-bit binary code 就需要输出对应的instruction. 由于不会截取的技巧,代码量非常可观 = =,所以说,一直很讨厌做模拟题!!! 留下这代码,纪念一个代码还是不够精简的自己!!!内存和时间还能接受,也比较容易理解,不过好多重复代码= =.以下这个可以代码可以忽略,

BestCoder20 1002.lines (hdu 5124) 解题报告

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5124 题目意思:给出 n 条线段,每条线段用两个整数描述,对于第 i 条线段:xi,yi 表示该条线段的左端点和右端点.设 A 表示最多线段覆盖的点(当然这个 A 可以有多个啦,但这个无关紧要).现在需要求的是 A 被多少条线段覆盖. 我是不会做啦.......一开始还看不懂题目= = 以下是按照题解做的, 题解中的最大区间和,实际上就是求最大连续子序列的和. 1 #include <iostrea