本解法启发自http://blog.csdn.net/pickless/article/details/12034365, 以下文字部分亦转自该文
O(n^2)的算法很容易想到:
找寻一个点j,将原来的price[0..n-1]分割为price[0..j]和price[j..n-1],分别求两段的最大profit。
进行优化:
对于点j+1,求price[0..j+1]的最大profit时,很多工作是重复的,在求price[0..j]的最大profit中已经做过了。
类似于Best Time to Buy and Sell Stock,可以在O(1)的时间从price[0..j]推出price[0..j+1]的最大profit。
但是如何从price[j..n-1]推出price[j+1..n-1]?反过来思考,我们可以用O(1)的时间由price[j+1..n-1]推出price[j..n-1]。
最终算法:
数组l[i]记录了price[0..i]的最大profit,
数组r[i]记录了price[i..n]的最大profit。
已知l[i],求l[i+1]是简单的,同样已知r[i],求r[i-1]也很容易。
最后,我们再用O(n)的时间找出最大的l[i]+r[i],即为题目所求。
#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
class Solution
{
public:
//注释部分为 O(n^2)的解法
/*
int subProfit(vector<int> &p, int lhs, int rhs)
{
int len = rhs - lhs + 1;
int minv = p[lhs], profit = 0;
for(int i = lhs + 1; i <= rhs; ++i)
{
if(profit < p[i] - minv)
profit = p[i] - minv;
if(p[i] < minv)
minv = p[i];
}
return profit;
}
int maxProfit(vector<int> &prices)
{
int len = prices.size();
if(len < 2)
return 0;
int maxv = 0;
for(int i = 1; i < len - 1; ++i)
{
int tmp = subProfit(prices, 0, i) + subProfit(prices, i+1, len-1);
if(maxv < tmp)
maxv = tmp;
}
return maxv;
}*/
int maxProfit(vector<int> &prices)
{
int len = prices.size();
if(len < 2)
return 0;
int *l = new int[len]; //用以记录 prices[0 ~ i] 最大的收益
int *r = new int[len]; //用以记录 prices[i ~ len -1] 最大的收益
memset(l, 0, sizeof(int) * len);
memset(r, 0, sizeof(int) * len);
int minv = prices[0];
for(int i = 1; i < len; ++i)
{
l[i] = max(prices[i] - minv, l[i-1]);
minv = min(prices[i], minv);
}
int maxv = prices[len - 1];
for(int j = len - 2; j >= 0; --j)
{
r[j] = max(maxv - prices[j], r[j + 1]);
maxv = max(prices[j], maxv);
}
int re = 0;
for(int i = 0; i < len; ++i)
{
re = re > l[i] + r[i] ? re : l[i] + r[i];
}
delete[] l;
delete[] r;
return re;
}
};
int main(void)
{
/*
int a[] = {
3, 7, 8, 2, 9, 6, 1, 10, 5
};*/
//int a[] = {
// 4, 2, 1
// };
//
int a[] = {
1, 2, 4
};
int len = sizeof(a) / sizeof(int);
vector<int> v(a, a + len);
Solution s;
cout << s.maxProfit(v) << endl;
return 0;
}