编程珠玑 求解最大字串和

 

编程珠玑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达到过的最大值。

更多信息请查看 java进阶网 http://www.javady.com

废话不多说了，直接上代码：

?

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; }