bzoj 2482: [Spoj1557] Can you answer these queries II

Description

给定n个元素的序列。

给出m个询问：求l[i]~r[i]的最大子段和（可选空子段）。

这个最大子段和有点特殊：一个数字在一段中出现了两次只算一次。

比如：1，2，3，2，2，2出现了3次，但只算一次，于是这个序列的和是1+2+3=6。

Input

第一行一个数n。

第二行n个数，为给定的序列，这些数的绝对值小于等于100000。

第三行一个数m。

接下来m行，每行两个数，l[i]，r[i]。

Output

M行，每行一个数，为每个询问的答案。

Sample Input

9

4 -2 -2 3 -1 -4 2 2 -6

3

1 2

1 5

4 9

Sample Output

4

5

3

HINT

【数据说明】

30%：1 <= n, m <= 100

100%：1 <= n, m <= 100000

Source

spoj gss2

枚举右端点。我们用线段树维护以这个点开始到当前点的ans和过去最大ans
然后离线回答询问即可

#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
using namespace std;
struct tree
{
int l,r;
long long m,vm;
long long tag,vtag;
}tr[800001];
struct que
{
int l,r;
int p;
long long ans;
}ask[100001];
inline bool cmp1(que x,que y)
{
if(x.r<y.r)
return true;
return false;
}
inline bool cmp2(que x,que y)
{
if(x.p<y.p)
return true;
return false;
}
inline void up(int p)
{
tr[p].m=max(tr[p*2].m,tr[p*2+1].m);
tr[p].vm=max(tr[p*2].vm,tr[p*2+1].vm);
}
inline void push(int p)
{
tr[p*2].vm=max(tr[p*2].vm,tr[p*2].m+tr[p].vtag);
tr[p*2+1].vm=max(tr[p*2+1].vm,tr[p*2+1].m+tr[p].vtag);
tr[p*2].m+=tr[p].tag;
tr[p*2+1].m+=tr[p].tag;

tr[p*2].vtag=max(tr[p*2].vtag,tr[p*2].tag+tr[p].vtag);
tr[p*2+1].vtag=max(tr[p*2+1].vtag,tr[p*2+1].tag+tr[p].vtag);
tr[p*2].tag+=tr[p].tag;
tr[p*2+1].tag+=tr[p].tag;
tr[p].tag=0;
tr[p].vtag=0;
}
inline void build(int p,int l,int r)
{
tr[p].l=l;
tr[p].r=r;
if(l!=r)
{
int mid=(l+r)/2;
build(p*2,l,mid);
build(p*2+1,mid+1,r);
up(p);
}
}
inline void add(int p,int l,int r,long long x)
{
if(l>r)
return ;
/* if(l==0)
l=1;*/
if(l<=tr[p].l&&tr[p].r<=r)
{
tr[p].m+=x;
tr[p].tag+=x;
tr[p].vm=max(tr[p].vm,tr[p].m);
tr[p].vtag=max(tr[p].vtag,tr[p].tag);
}
else
{
push(p);
int mid=(tr[p].l+tr[p].r)/2;
if(l<=mid)
add(p*2,l,r,x);
if(r>mid)
add(p*2+1,l,r,x);
up(p);
}
}
inline long long askx(int p,int l,int r)
{
if(l<=tr[p].l&&tr[p].r<=r)
return tr[p].vm;
else
{
push(p);
int mid=(tr[p].l+tr[p].r)/2;
long long ans=0;
if(l<=mid)
ans=max(ans,askx(p*2,l,r));
if(r>mid)
ans=max(ans,askx(p*2+1,l,r));
return ans;
}
}
int a[1000001];
int pre[1000001],la[1000001];
int main()
{
//	 freopen("sequence.in","r",stdin);
//	 freopen("sequence.out","w",stdout);
int n;
scanf("%d",&n);
int i;
for(i=1;i<=n;i++)
scanf("%d",&a[i]);
for(i=1;i<=n;i++)
{
pre[i]=la[a[i]];
la[a[i]]=i;
}
build(1,1,n);
int m;
scanf("%d",&m);
for(i=1;i<=m;i++)
{
scanf("%d%d",&ask[i].l,&ask[i].r);
ask[i].p=i;
}
sort(ask+1,ask+1+m,cmp1);
int d=1;
for(i=1;i<=n;i++)
{
add(1,pre[i]+1,i,a[i]);
while(d<=m&&ask[d].r==i)
{
ask[d].ans=askx(1,ask[d].l,ask[d].r);
d++;
}
}
sort(ask+1,ask+1+m,cmp2);
for(i=1;i<=m;i++)
printf("%lld\n",ask[i].ans);
return 0;
}