编程珠玑8.4节讲扫描算法,我看了半天都没看明白,最后自己写了一遍,终于搞懂了,把它记下来,以免今后忘了。
首先,书上的算法是这样写的:
1
2
3
4
5
|
maxsofar
=
0
maxendinghere
=
0
for
i
=
[
0
n)
maxendinghere
=
max
(maxendinghere
+
x[i],
0
)
maxsofar
=
max
(maxsofar, maxendinghere)
|
这个max so far和max ending here究竟表示什么呢?
先来看我写的程序的输出:
输入:{31, -41, 59, 26, -53, 58, 97, -93, -23, 84}
输出:
其中Accu就是maxsofar,而Max就是maxending。
原理是这样的:
我们把成为最大字串和的字串SubSeq看成两部分:
前一部分为整个字串和贡献的值要大于0。
后一部分不能减少前一部分的值。
比如:{59, 26, -53}是前一部分,而{58, 97}是后一半。
好,现在再来看算法:
每次将x[i]加到Accu上,看Accu是不是减小了;
如果是,就先暂时不要x[i],直到Accu比它达到过的最大值还要大,这时就可以把暂时不要的的所有x[..]加到字串中,因为它们对字串的贡献>0。
而Max记录的就是Accu达到过的最大值。
废话不多说了,直接上代码:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
#include <stdio.h>
#include <limits.h>
#include <math.h>
#include <string.h>
int
seq_start = 0;
int
seq_end = 0;
int
max_sum(
int
x[],
int
n)
{
int
max = INT_MIN, temp_sum = 0;
bool
hasNewStart =
true
;
printf
(
"No Accu Max SPos EPos\tSubSeq\n"
);
for
(
int
i = 0; i < n; i++)
{
temp_sum += x[i];
if
(temp_sum > max)
{
max = temp_sum;
seq_end = i;
if
(hasNewStart)
{
seq_start = i;
hasNewStart =
false
;
}
}
else
if
(temp_sum < 0)
{
temp_sum = 0;
hasNewStart =
true
;
}
char
* tmp_str = arr2str((
int
*) &(x[seq_start]), seq_end - seq_start + 1);
printf
(
"%2d %4d %4d %3d %3d\t%s \n"
, i, temp_sum, max, seq_start, seq_end, tmp_str);
delete
[] tmp_str;
}
return
max;
}
|