题目链接: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
a_i
ai乘上区间和的结果就行了。
比如数组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
a_i
ai*(sum[r]-sum[l-1])},得到最大值为36。
如果数组中有负数,就要分类讨论了:
如果
a
i
a_i
ai>0,那么在区间(l,r)中应该找到一个最大的sum[a]-sum[b-1] (l≤a≤i≤b≤r)
使得
a
i
a_i
ai(sum[a]-sum[b-1])最大。
如果
a
i
a_i
ai<0,那么在区间(l,r)中应该找到一个最小的sum[a]-sum[b-1] (l≤a≤i≤b≤r)
使得
a
i
a_i
ai(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;
}