南昌网络赛-I. Max answer （RMQ+单调栈）

题目链接：https://nanti.jisuanke.com/t/38228
题目大意：
给一个n个数的数组，找到其中连续几个数，使得这几个数的和乘上这几个数中的最小值的结果最大。
和poj2796很像，不过poj2796的数列只有正数，这道题有正有负，所以需要分类讨论。
poj2796：http://poj.org/problem?id=2796

本题需要的基础知识：
RMQ：https://blog.csdn.net/lr7682/article/details/89715286
单调栈：https://www.cnblogs.com/tham/p/8038828.html

解题思路：
先将题目简化一下，假设数列中只有正数时，
我们可以先用一个sum数组记录一下前缀和。
然后用单调栈找到以第i个数为最小值的一个最大的区间（l,r），复杂度为O(n)。
然后遍历数组，找到最大的一个$a_i$乘上区间和的结果就行了。
比如数组4,1,3,5,2，以第i个数为最小值的左端点分别为0,0,2,3,1，右端点分别为0,4,3,3,4。
所以以第i个数字为最小值的区间就分别是（0,0）（0,4）（2,3）（3,3）（1,4）。
计算max{$a_i$*（sum[r]-sum[l-1]）}，得到最大值为36。

如果数组中有负数，就要分类讨论了：
如果$a_i$>0，那么在区间(l,r)中应该找到一个最大的sum[a]-sum[b-1] (l≤a≤i≤b≤r)
使得$a_i$(sum[a]-sum[b-1])最大。
如果 $a_i$<0，那么在区间(l,r)中应该找到一个最小的sum[a]-sum[b-1] (l≤a≤i≤b≤r)
使得$a_i$(sum[a]-sum[b-1])最大。（我们期望sum[a]-sum[b-1]可以是一个绝对值很大的负数）
因为要在一个已知的区间里找到sum[a]-sum[b-1]的最大值/最小值，而前缀和又不会被修改，因此想到RMQ的方法计算。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <stack>
#define N 500005
#define logn 25
#define INF 1e17
#define LL long long
#define CL(a,b) memset(a,b,sizeof(a))

using namespace std;

int n;
LL logg[N],a[N],l[N],r[N];
LL sum[N],stmax[N][logn],stmin[N][logn];
stack<int>indexst;

void init()
{
    logg[0]=-1;
    for(int i=1;i<=n;i++)
        logg[i]=logg[i>>1]+1;
    for(int i=1;i<=n;i++)
        sum[i]=sum[i-1]+(LL)(a[i]);
    for(int i=1;i<=n;i++)
        stmax[i][0]=stmin[i][0]=sum[i];
    for(int j=1;j<=logn;j++)
        for(int i=0;i+(1<<j)-1<=n;i++)
        {
            stmax[i][j]=max(stmax[i][j-1],stmax[i+(1<<(j-1))][j-1]);
            stmin[i][j]=min(stmin[i][j-1],stmin[i+(1<<(j-1))][j-1]);
        }
    for(int i=1;i<=n;i++)
    {
        while(!indexst.empty()&&a[i]<=a[indexst.top()])
            indexst.pop();
        if(indexst.empty())
            l[i]=1;
        else
            l[i]=indexst.top()+1;
        indexst.push(i);
    }
    indexst=stack<int>();
    for(int i=n;i>=1;i--)
    {
        while(!indexst.empty()&&a[i]<=a[indexst.top()])
            indexst.pop();
        if(indexst.empty())
            r[i]=n;
        else
            r[i]=indexst.top()-1;
        indexst.push(i);
    }
}

LL querymax(int l,int r)
{
    int k=logg[r-l+1];
    return max(stmax[l][k],stmax[r-(1<<k)+1][k]);
}

LL querymin(int l,int r)
{
    int k=logg[r-l+1];
    return min(stmin[l][k],stmin[r-(1<<k)+1][k]);
}

int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%lld",&a[i]);
    init();
    LL ans=-INF;
    for(int i=1;i<=n;i++)
    {
        if(a[i]<0)
        {
            LL rsum=querymin(i,r[i]);
            LL lsum=querymax(l[i]-1,i-1);
            ans=max(ans,(rsum-lsum)*(LL)a[i]);
        }
        else if(a[i]>0)
        {
            LL rsum=querymax(i,r[i]);
            LL lsum=querymin(l[i]-1,i-1);
            ans=max(ans,(rsum-lsum)*(LL)a[i]);
        }
        else
            ans=ans>0?ans:0;
    }
    printf("%lld\n",ans);
    return 0;
}