[杂题 贪心] BZOJ 2138 stone

#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;typedef long long ll;inline char nc(){  static char buf[100000],*p1=buf,*p2=buf;  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}inline void read(int &x){  char c=nc(),b=1;  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;}const int N=40005;struct SegMax{  ll T[N<<2],F[N<<2];  void Build(int x,int l,int r,ll *a){    if (l==r) return void(T[x]=a[l]);    int mid=(l+r)>>1;    Build(x<<1,l,mid,a); Build(x<<1|1,mid+1,r,a);    T[x]=max(T[x<<1],T[x<<1|1]);  }  void push(int x){    if (F[x]) T[x<<1]+=F[x],T[x<<1|1]+=F[x],F[x<<1]+=F[x],F[x<<1|1]+=F[x],F[x]=0;  }  void Add(int x,int l,int r,int ql,int qr,ll a){    if (ql<=l && r<=qr){      T[x]+=a; F[x]+=a; return;    }    push(x);    int mid=(l+r)>>1;    if (ql<=mid) Add(x<<1,l,mid,ql,qr,a);    if (qr>mid) Add(x<<1|1,mid+1,r,ql,qr,a);    T[x]=max(T[x<<1],T[x<<1|1]);  }  ll Query(int x,int l,int r,int ql,int qr){    if (ql<=l && r<=qr) return T[x];    push(x); int mid=(l+r)>>1; ll ret=-1LL<<30;    if (ql<=mid) ret=max(ret,Query(x<<1,l,mid,ql,qr));    if (qr>mid) ret=max(ret,Query(x<<1|1,mid+1,r,ql,qr));    return ret;  }}G;struct SegMin{  ll T[N<<2],F[N<<2];  void Build(int x,int l,int r,ll *a){    if (l==r) return void(T[x]=a[l]);    int mid=(l+r)>>1;    Build(x<<1,l,mid,a); Build(x<<1|1,mid+1,r,a);    T[x]=min(T[x<<1],T[x<<1|1]);  }  void push(int x){    if (F[x]) T[x<<1]+=F[x],T[x<<1|1]+=F[x],F[x<<1]+=F[x],F[x<<1|1]+=F[x],F[x]=0;  }  void Add(int x,int l,int r,int ql,int qr,ll a){    if (ql<=l && r<=qr){      T[x]+=a; F[x]+=a; return;    }    push(x);    int mid=(l+r)>>1;    if (ql<=mid) Add(x<<1,l,mid,ql,qr,a);    if (qr>mid) Add(x<<1|1,mid+1,r,ql,qr,a);    T[x]=min(T[x<<1],T[x<<1|1]);  }  ll Query(int x,int l,int r,int ql,int qr){    if (ql<=l && r<=qr) return T[x];    push(x); int mid=(l+r)>>1; ll ret=1LL<<60;    if (ql<=mid) ret=min(ret,Query(x<<1,l,mid,ql,qr));    if (qr>mid) ret=min(ret,Query(x<<1|1,mid+1,r,ql,qr));    return ret;  }}F;int n,m,a[N];ll s[N];int K[N];int main(){  int x,y,z,p,l,r;  freopen("forget.in","r",stdin);  freopen("forget.out","w",stdout);  read(n); read(x); read(y); read(z); read(p);  for (int i=1;i<=n;i++)    a[i]=((ll)(i-x)*(i-x)+(ll)(i-y)*(i-y)+(ll)(i-z)*(i-z))%p,s[i]=s[i-1]+a[i];  F.Build(1,0,n,s); G.Build(1,0,n,s);  read(m); read(K[1]); read(K[2]); read(x); read(y); read(z); read(p);  for (int i=3;i<=m;i++)    K[i]=((ll)x*K[i-1]+(ll)y*K[i-2]+z)%p;  for (int i=1;i<=m;i++){    read(l); read(r);    ll fmin=F.Query(1,0,n,r,n),gmax=G.Query(1,0,n,0,l-1);    ll k=min((ll)K[i],fmin-gmax);    printf("%lld\n",k);    F.Add(1,0,n,r,n,-k);    G.Add(1,0,n,l,n,-k);  }  return 0;}