BZOJ3747 POI2015Kinoman(线段树)

时间:2023-03-09 19:58:59
BZOJ3747 POI2015Kinoman(线段树)

  考虑固定左端点,求出该情况下能获得的最大值。于是每次可以在某数第一次出现的位置加上其价值,第二次出现的位置减掉其价值,查询前缀最大值就可以了。每次移动左端点在线段树上更新即可。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int read()
{
int x=,f=;char c=getchar();
while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
return x*f;
}
#define N 1000010
#define ll long long
int n,m,a[N],b[N],p[N],nxt[N];
int L[N<<],R[N<<];
ll mx[N<<],sum[N<<],ans;
void build(int k,int l,int r)
{
L[k]=l,R[k]=r;
if (l==r) return;
int mid=l+r>>;
build(k<<,l,mid);
build(k<<|,mid+,r);
}
void up(int k)
{
sum[k]=sum[k<<]+sum[k<<|];
mx[k]=max(mx[k<<],sum[k<<]+mx[k<<|]);
}
void modify(int k,int p,int x)
{
if (L[k]==R[k]) {mx[k]=sum[k]=x;return;}
int mid=L[k]+R[k]>>;
if (p<=mid) modify(k<<,p,x);
else modify(k<<|,p,x);
up(k);
}
ll query(int k,int l,int r)
{
if (L[k]==l&&R[k]==r) return mx[k];
int mid=L[k]+R[k]>>;
if (r<=mid) return query(k<<,l,r);
else if (l>mid) return query(k<<|,l,r);
else return max(query(k<<,l,mid),sum[k<<]+query(k<<|,mid+,r));
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("bzoj3747.in","r",stdin);
freopen("bzoj3747.out","w",stdout);
const char LL[]="I64d\n";
#else
const char LL[]="%lld\n";
#endif
n=read(),m=read();
for (int i=;i<=n;i++)
{
a[i]=read();
nxt[p[a[i]]]=i,p[a[i]]=i;
}
for (int i=;i<=n;i++) if (!nxt[i]) nxt[i]=n+;nxt[n+]=n+;
for (int i=;i<=m;i++) b[i]=read();
build(,,n+);
memset(p,,sizeof(p));
for (int i=;i<=n;i++)
if (!p[a[i]]) p[a[i]]=,modify(,i,b[a[i]]),modify(,nxt[i],-b[a[i]]);
for (int i=;i<=n;i++)
{
ans=max(ans,query(,i,n));
modify(,i,),modify(,nxt[i],b[a[i]]),modify(,nxt[nxt[i]],-b[a[i]]);
}
cout<<ans;
return ;
}