HDU1003 Max Sum

题目地址

看算法导论看的有点害怕了，做做题找找成就感。然而。。学了这么长时间算法，这样一道dp竟然还不能直接AC。。。

拿到这个题可以看出是一道dp题目，然而我开始考虑的是备忘录d[i][j]存储从第i个到第j个数的和，d[i][j] = d[i][j-1]+a[j]，然后每次判断是否更新最大值，基本就跟遍历一样了（但是效率要比普通遍历要高一些，毕竟有备忘录），代码如下：

#include<iostream>
#include<vector>
using namespace std;
short int a[100005];
vector<int> d[100005];  //d[i][j]表示第i个数开头的总共j+1个数的和，i>=0
int main()
{
    short T;
    short n;
    cin>>T;
    short t=0;
    while(t!=T)
    {
        for(int i=0; i<100005; i++)
            vector<int> ().swap(d[i]);
        long long s = 0;
        int be=1, en=1;
        cin>>n;
        for(int i=0; i<n; i++)
        {
            cin>>a[i];
        }
        for(int i=0; i<n; i++)
        {
            d[i].push_back(a[i]);
            if(a[i]>s)
            {
                s=a[i];
                be=i+1;
                en=i+1;
            }
        }

        for(int i=2; i<=n; i++)  //个数
        {
            for(int j=0; j<=n-i; j++)  //从第几个开始
           {
                if(d[j][i-2]+a[j+i-1]>=s)
                {
                    if(d[j][i-2]+a[j+i-1]==s)
                    {
                        if(j+1<be)     //题目要求如果存在多个结果，则输出第一个
                        {
                            s=d[j][i-2]+a[j+i-1];
                            be=j+1;
                            en=be+i-1;
                        }
                    }
                    else
                    {
                        s=d[j][i-2]+a[j+i-1];
                        be=j+1;
                        en=be+i-1;
                    }

                }
                d[j].push_back(d[j][i-2]+a[j+i-1]);
           }
        }
        t++;

        cout<<"Case "<<t<<":"<<endl;
        cout<<s<<" "<<be<<" "<<en<<"\n"<<endl;
    }
    return 0;
}

嗯，这样当然过不了。。。。Memory Limit Exceeded

毕竟该问题的最终答案毕竟是以序列中的某个数结尾的，所以，我们将各个数作结尾时序列和求出来，然后找出以哪个数结尾时值最大，即为最终解。

所以，这个问题中，我们可以把以某个数结尾时的最大值看作一个状态。

比如序列a,b,c. f(a)表示以a结尾时序列和的最大值。
那么只需将f(a),f(b),f(c)依次求出，找出最大值即可。

上面的例子中f(b)其实是依赖于f(a)的，如果f(a)>=0,那么f(b)=f(a)+b，否则：f(b)=b。因为f(b)必须是以b结尾并且连续的。这样状态转移方程就出来了：f(i+1)=(f(i)+a[i]>=a[i])?(f(i)+a[i]):a[i]
代码：

//问题的最终解不一定是以哪个数结尾的，也就是说找出以各个数结尾的子序列的最大值，问题也就解决了。（分解为各个子问题）
//从解决一个问题到解决下一个问题时，有两种选择，要么重新组织一个序列，要么将当前数作为序列结尾。
#include<iostream>
using namespace std;
int main()
{
    int T;
    cin>>T;
    int t = 1, TT = T;
    while(T--)
    {
        int n,x,left,right,leftRes,rightRes;
        long long res,maxRes;
        cin>>n;
        cin>>res;  //只有一个数时，以该数结尾最大值就是该数。
        maxRes = res;
        left = right = leftRes= rightRes = 1;
        n--;
        while(n--)   //res表示以上一个数结尾的最大值
        {
            cin>>x;
            if(res+x>=x)
            {
                res+=x;
                right++;
            }
            else     //重新开始新的序列
            {
                res = x;
                right++;
                left=right;
            }
            if(res>maxRes)
            {
                maxRes = res;
                leftRes = left;
                rightRes = right;
            }
        }
        cout<<"Case "<<t<<":"<<endl;
        cout<<maxRes<<" "<<leftRes<<" "<<rightRes<<endl;
        if(t!=TT)   cout<<endl;
        t++;
    }
    return 0;
}