BZOJ5326 : [Jsoi2017]博弈

将所有物品按照$b$的选择顺序排序,则先手在任意前$i$个物品中最多只能拿走$lceilfrac{i}{2} ceil$个物品。

将每个物品的价值设为$a+b$,那么答案为先手拿走的价值和减去所有物品的$b$之和,目标是最大化先手拿走的价值和。

如果不考虑修改,则满足拟阵,可以贪心选取,修改时最优解最多发生一处变动。

每个修改可以看作是先删除一个物品,再插入一个物品。

设$f[i]$表示$[1,i]$里选择的物品数减去$lceilfrac{i}{2} ceil$,则方案合法的条件为$max(f[1,n])leq 0$。

对于插入物品$x$:

  • 如果$x$可以直接加入最优解,那么直接加入。
  • 否则加入$x$会导致$f$中出现$1$,那么需要找到一个最优解中的价值最小的物品$y$,满足$y$在$f$从左往右第一个$1$之前,要么不选$x$,要么将$y$替换成$x$。

对于删除物品$x$:

  • 如果$x$不在最优解中,那么直接删除。
  • 否则删除$x$后可以继续多加入一个物品,那么需要找到一个不在最优解中的价值最大的物品$y$,满足$y$到$n$的所有$f$都是负数。

以上所有操作都可以用线段树维护,时间复杂度$O((n+m)log n)$。

#include<cstdio>
#include<algorithm>
using namespace std;
typedef pair<int,int>P;
const int N=100010,M=262150,inf=~0U>>1,BUF=10000000;
char Buf[BUF],*buf=Buf;long long ans;
int n,m,i,x,y,pos[N];bool used[N];
P ma[M],mi[M];int f[M],tag[M];
struct E{int a,b,p;}e[N];
inline bool cmp(const E&a,const E&b){
  if(a.b!=b.b)return a.b>b.b;
  return a.p<b.p;
}
inline void read(int&a){for(a=0;*buf<48;buf++);while(*buf>47)a=a*10+*buf++-48;}
inline void up(int x){
  ma[x]=max(ma[x<<1],ma[x<<1|1]);
  mi[x]=min(mi[x<<1],mi[x<<1|1]);
}
void build(int x,int a,int b){
  if(a==b){
    ma[x]=P(-inf,a);
    mi[x]=P(inf,a);
    f[x]=-(a+1)/2;
    return;
  }
  int mid=(a+b)>>1;
  build(x<<1,a,mid),build(x<<1|1,mid+1,b),up(x);
  f[x]=max(f[x<<1],f[x<<1|1]);
}
inline void tag1(int x,int p){f[x]+=p;tag[x]+=p;}
inline void pb(int x){if(tag[x])tag1(x<<1,tag[x]),tag1(x<<1|1,tag[x]),tag[x]=0;}
void maketag(int x,int a,int b,int c,int d,int p){
  if(c<=a&&b<=d){tag1(x,p);return;}
  pb(x);
  int mid=(a+b)>>1;
  if(c<=mid)maketag(x<<1,a,mid,c,d,p);
  if(d>mid)maketag(x<<1|1,mid+1,b,c,d,p);
  f[x]=max(f[x<<1],f[x<<1|1]);
}
void setused(int x,int a,int b,int c,int p){
  if(a==b){
    used[a]=1;
    ma[x]=P(-inf,a);
    mi[x]=P(p,a);
    return;
  }
  int mid=(a+b)>>1;
  if(c<=mid)setused(x<<1,a,mid,c,p);else setused(x<<1|1,mid+1,b,c,p);
  up(x);
}
void setunused(int x,int a,int b,int c,int p){
  if(a==b){
    used[a]=0;
    ma[x]=P(p,a);
    mi[x]=P(inf,a);
    return;
  }
  int mid=(a+b)>>1;
  if(c<=mid)setunused(x<<1,a,mid,c,p);else setunused(x<<1|1,mid+1,b,c,p);
  up(x);
}
int askf(int x,int a,int b,int c,int d){
  if(c<=a&&b<=d)return f[x];
  pb(x);
  int mid=(a+b)>>1,t=-inf;
  if(c<=mid)t=askf(x<<1,a,mid,c,d);
  if(d>mid)t=max(t,askf(x<<1|1,mid+1,b,c,d));
  return t;
}
P askmi(int x,int a,int b,int c,int d){
  if(c<=a&&b<=d)return mi[x];
  int mid=(a+b)>>1;P t(inf,0);
  if(c<=mid)t=askmi(x<<1,a,mid,c,d);
  if(d>mid)t=min(t,askmi(x<<1|1,mid+1,b,c,d));
  return t;
}
P askma(int x,int a,int b,int c,int d){
  if(c<=a&&b<=d)return ma[x];
  int mid=(a+b)>>1;P t(-inf,0);
  if(c<=mid)t=askma(x<<1,a,mid,c,d);
  if(d>mid)t=max(t,askma(x<<1|1,mid+1,b,c,d));
  return t;
}
inline int findfirstpos(){
  int x=1,a=1,b=n,mid;
  while(a<b){
    pb(x);
    mid=(a+b)>>1;
    if(f[x<<1]>0)x=x<<1,b=mid;else x=x<<1|1,a=mid+1;
  }
  return a;
}
inline int findlastneg(){
  int x=1,a=1,b=n,mid,ret=n+1;
  while(a<b){
    pb(x);
    mid=(a+b)>>1;
    if(f[x<<1|1]<0){
      x=x<<1;
      b=mid;
      ret=mid+1;
    }else{
      x=x<<1|1;
      a=mid+1;
    }
  }
  if(f[x]<0)ret=a;
  return ret;
}
inline void additem(int x){
  int val=e[x].a+e[x].b,tmp=askf(1,1,n,x,n);
  maketag(1,1,n,x,n,1);
  if(tmp<0){
    setused(1,1,n,x,val);
    ans+=val;
    return;
  }
  int o=findfirstpos();
  maketag(1,1,n,x,n,-1);
  P y=askmi(1,1,n,1,o);
  if(y.first>=val){
    setunused(1,1,n,x,val);
    return;
  }
  ans+=val-y.first;
  setunused(1,1,n,y.second,y.first);
  maketag(1,1,n,y.second,n,-1);
  setused(1,1,n,x,val);
  maketag(1,1,n,x,n,1);
}
inline void delitem(int x){
  if(!used[x])return;
  ans-=e[x].a+e[x].b;
  setunused(1,1,n,x,-inf);
  maketag(1,1,n,x,n,-1);
  int o=findlastneg();
  if(o>n)return;
  P y=askma(1,1,n,o,n);
  if(y.first<=0)return;
  ans+=y.first;
  setused(1,1,n,y.second,y.first);
  maketag(1,1,n,y.second,n,1);
}
int main(){
  fread(Buf,1,BUF,stdin);read(n);
  for(i=1;i<=n;i++)read(e[i].a);
  for(i=1;i<=n;i++)read(e[i].b),e[i].p=i,ans-=e[i].b;
  sort(e+1,e+n+1,cmp);
  for(i=1;i<=n;i++)pos[e[i].p]=i;
  build(1,1,n);
  for(i=1;i<=n;i++)additem(i);
  printf("%lld
",ans);
  read(m);
  while(m--){
    read(x),read(y);
    x=pos[x];
    delitem(x);
    e[x].a=y;
    additem(x);
    printf("%lld
",ans);
  }
  return 0;
}

  

原文地址:https://www.cnblogs.com/clrs97/p/10360779.html