poj1796 Feel Good

原题
题目大意：
给出一个长度为n（n<100000）的序列，求出一个子序列，使得这个序列中的最小值乘以这个序列的和的值最大。
算法分析：
维护一个单调递增栈，每次进栈新元素与栈顶元素比较，小于栈顶元素就出栈，维护l,r两个变量，表示以num为最小值的最大区间。

一个长度为n的序列，对于一个区间中的value为这个区间的最小值乘上这个区间的所有数的和。
这就是找每一个点，找以其为最小值的区间。

显然直接暴力法搜索必然超时，可以开一个单调递增栈，对于每一个新元素a[i]，当前值小于等于栈顶元素，表示以栈顶元素为最小值的区间的右端是i，栈顶元素a[j]出栈。

由于是单调递增栈，新的栈顶元素a[k]必然小于a[j]，故需要更新新的栈顶元素的右端值，新元素的左端的值也要更新为a[j]的左端。一直重复操作直到栈顶元素大于新元素，这样这个栈的每一个元素出栈时其左右区间必然是刚好最大的以该元素为最小值的区间。

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<stack>
using namespace std;
const int maxn=1e5+7;
int n,ansl,ansr;
struct node{
	int l,r,num;
}s[maxn];
int a[maxn];
long long sum[maxn];//sum数组做前缀和
long long ans;
void update(int top)//某一元素出栈时比较函数
{
     if((sum[s[top].r] - sum[s[top].l-1])*s[top].num>ans)
     {
        ans=(sum[s[top].r]-sum[s[top].l-1])*s[top].num ;
        ansl=s[top].l;
        ansr=s[top].r;
     }
     if(top>0)//栈不为空
     s[top-1].r=s[top].r;
}
int main()
{
	while(scanf("%d",&n)!=EOF)
	{
		ans=-1000000;
		int top=-1;
	    sum[0]=0;
	    for(int i=1;i<=n;i++)
	    {
		   scanf("%d",&a[i]);
		   sum[i]=sum[i-1]+a[i] ;
	    }
	    for(int i=1;i<=n;i++)
	    {
    		node temp={i,i,a[i]};//形成新元素
    		while(top!=-1&&s[top].num>=a[i])//栈不为空，维护单增站
            {
                 update(top) ;
                 temp.l=s[top].l ;//因为top值大于a[i],所以top的左区间一定是新元素的左区间
                 top--;            }
            s[++top].l=temp.l ;
            s[top].r=temp.r;
            s[top].num=temp.num ;
        }
        while(top!=-1)//如果栈内还有元素
        {
            update(top);
            top--;
        }
        printf("%lld\n",ans) ;
        printf("%d %d\n",ansl,ansr) ;
 	}
    return 0;
}
/*int n,ansl,ansr;
long long num,ans;
struct node{
	int l,r,val;
	long long sum;
}a[maxn];
stack<node> h;
void update(node x)
{
	num=(a[x.r].sum-a[x.l-1].sum)*x.val;
	if(ans<num)
	{
       ans=num;
	   ansl=x.l;
	   ansr=x.r;
	}
	if(!h.empty())
	{
		node y1,y2;
		y1=h.top(),h.pop();
		y2=h.top(),h.pop();
		y2.r=y1.r;
		h.push(y2),h.push(y1);
	}
	return ;
}
void in(node x)
{
    node temp=h.top();
	while(x.val<temp.val&&!h.empty())
	{
		temp=h.top();
		update(temp);
	    x.l=temp.l;
    	h.pop();
	}
	h.push(x);
	return ;
}
int main()
{
	scanf("%d",&n);
	ans=-10000;
	a[0].sum=0;
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&a[i].val);
		a[i].l=i,a[i].r=i;
		a[i].sum=a[i-1].sum+a[i].val;
	}
	for(int i=1;i<=n;i++)
	{
		in(a[i]);	
	}
	printf("%lld\n%d %d",ans,ansl,ansr);
	return 0;
}*/