Блог пользователя Binary_ToothLess

Автор Binary_ToothLess, 10 лет назад, По-английски

Don't be afraid this is not very difficult....

For new programmer, this is very essential... We know a little bit for that we don't answer many problem as like "largest sum contiguous sub array" that's really very easy to determine largest sum if we know the following algorithm.....

Initialize: max_so_far = 0

max_ending_here = 0

Loop for each element of the array

(a) max_ending_here = max_ending_here + a[i]

(b) if(max_ending_here < 0) max_ending_here = 0

(c) if(max_so_far < max_ending_here) max_so_far = max_ending_here

return max_so_far

Explanation: Simple idea of the Kadane's algorithm is to look for all positive contiguous segments of the array (max_ending_here is used for this). And keep track of maximum sum contiguous segment among all positive segments (max_so_far is used for this). Each time we get a positive sum compare it with max_so_far and update max_so_far if it is greater than max_so_far

Lets take the example:

{-2, -3, 4, -1, -2, 1, 5, -3}

max_so_far = max_ending_here = 0

for i=0, a[0] = -2 max_ending_here = max_ending_here + (-2) Set max_ending_here = 0 because max_ending_here < 0

for i=1, a[1] = -3 max_ending_here = max_ending_here + (-3) Set max_ending_here = 0 because max_ending_here < 0

for i=2, a[2] = 4 max_ending_here = max_ending_here + (4) max_ending_here = 4 max_so_far is updated to 4 because max_ending_here greater than max_so_far which was 0 till now

for i=3, a[3] = -1 max_ending_here = max_ending_here + (-1) max_ending_here = 3

for i=4, a[4] = -2 max_ending_here = max_ending_here + (-2) max_ending_here = 1

for i=5, a[5] = 1 max_ending_here = max_ending_here + (1) max_ending_here = 2

for i=6, a[6] = 5 max_ending_here = max_ending_here + (5) max_ending_here = 7 max_so_far is updated to 7 because max_ending_here is greater than max_so_far

for i=7, a[7] = -3 max_ending_here = max_ending_here + (-3) max_ending_here = 4

Code :

for(int i = 0; i < N; i++){

max_ending_here += Number[i];

        if(max_ending_here > max_so_far)max_so_far = max_ending_here;

        if(max_ending_here < 0)max_ending_here = 0;

    }

So the complexity for this code is O(n)...

Happy coding...:)

  • Проголосовать: нравится
  • +9
  • Проголосовать: не нравится

»
10 лет назад, # |
  Проголосовать: нравится +1 Проголосовать: не нравится

First time I saw this problem I didn't understand very well the idea. Then I saw it was an optimization in memory use. I would like to share the logic behind this optimization.

dp[i] -> max value from first position to i

For each position:

dp[i] = max(dp[i - 1] + Number[i], 0) // This is max_ending_here

maxsingleelement = max(maxsingleelement, Number[i]) // Check case all negative values.

Then we will get our solution taking the max dp[i]

For each position again:

sol = max(sol, dp[i]) // This is max_so_far

This implementation has O(n) in time complexity and memory usage.

Then we can use only two variables in the solution you gave us. The complexity will be the same but memory usage will be constant .

By the way, what happen when all the elements are negative values? This is a possible case. Based in this case I like to use a variable (maxsingleelement) that stores the max value from each Number[i]. At the end, if this value is lower than 0 then I return maxsingleelement otherwise return sol (max_so_far)

»
10 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

just to fine-tune ur implementation a bit.

if the input array is [ - 1,  - 2,  - 2,  - 1,  - 3] (i.e. all negative elements), then ur code would return 0.
to have it return  - 1 (i.e. the least negative element) instead, u can initialize max_so_far to  - ∞ instead of 0.

»
9 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

test case 1: For input array {-2, -3, 4, -1, -2, 1, 5, -3}, the max sum should be 7, with sub-arrary {4, -1, -2, 1, 5} test case 2: For input array {-2, -3, 4, -1, -2, 1, 5, -3, 4}, the max sum should be 8, with sub-array {4, -1, -2, 1, 5, -3, 4}.

I tried using the algorithm above, the result of test case 2 still 7 instead 8. Any issue on it?

»
9 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

How to implement the same for a 2D array?

  • »
    »
    8 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    Compute the prefix sums for all rows, then iterate over all pair of colums and for each pair l, r run Kadane's algorithm using prefix[i][r] - prefix[i][l - 1] as array elements.

»
8 лет назад, # |
  Проголосовать: нравится +10 Проголосовать: не нравится

I find this version slightly easier to understand. Sum for some contiguous sub array can be obtained using partial sums, by subtracting sum at the beginning from sum at the end. To obtain the largest sum you want to subtract as little as possible. So it is necessary to maintain minimum sum that can be subtracted.

int min_sum = 0, max_sum = 0; // max_sum = Number[0] to not allow empty sub array
int sum = 0;
for (int i=0; i < N; i++)
{
    sum += Number[i];
    if (sum - min_sum > max_sum) max_sum = sum - min_sum;
    if (sum < min_sum) min_sum = sum;
}
»
8 лет назад, # |
Rev. 2   Проголосовать: нравится 0 Проголосовать: не нравится

I think this two pointer approach works with O(n) :)

Code:
  • »
    »
    8 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    you can't apply two pointer because function is not monotonic. There can be negative numbers in the input. just try with input given in the post. Your code outputs 6 whereas ans is 7.

»
8 лет назад, # |
  Проголосовать: нравится +31 Проголосовать: не нравится

Three years later, I still can't believe this algorithm has a name.

  • »
    »
    3 года назад, # ^ |
      Проголосовать: нравится +19 Проголосовать: не нравится

    It's interesting to know the history behind this problem. We should understand that this all happened during early years of the applied computer science... For example, it includes such now-famous scientists as Michael Shamos and Jon Bentley, who spent 1 week (!) on optimizing an O(n^2) (!) solution for this problem down to O(n log n), using a divide-and-conquer algorithm (!). They also believed for some period of time that it's the best possible solution, given that other similar problems, mostly related to sorting, can't be solved faster than O(n log n). Only some time afterwards Jay Kadane, a statistician, accidentally learned about this problem and the colleagues' solutions and suggested the now-well-known linear-time algorithm.

    • »
      »
      »
      3 года назад, # ^ |
      Rev. 2   Проголосовать: нравится 0 Проголосовать: не нравится

      It’s actually pretty cool to know that people would historically come up with D&C algorithms to solve the problem before this. You would think D&C would not have existed at all at the time this algorithm was novel.

»
6 лет назад, # |
Rev. 4   Проголосовать: нравится -16 Проголосовать: не нравится

Nice one!

»
5 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

Any bug in it?

int kadane(int A[], int n){

  int ans = INT_MIN, curr = 0;

  for(int i = 0; i < n; ++i){
    curr = curr + A[i];
    ans = max(ans, curr);
    if(curr < 0)
      curr = 0;
  }
}

»
5 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

This algorithm can also be derived from the brute force algorithm using Bird-Meertens formalism

»
4 года назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

Can we get all the elements of the maximum subarray by kadane's algorithm??

»
4 года назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится