模拟赛 怨灵退治 题解（Hall定理+线段树）

题意：
有 n 群怨灵排成一排，燐每秒钟会选择一段区间，消灭至多 k 只怨灵。
如果怨灵数量不足 k，则会消灭尽量多的怨灵。
燐作为一只有特点的猫，它选择的区间是不会相互包含的。它想要知道它每秒最多能消灭多少怨灵。
要求：在之前每次都消灭尽量多的怨灵的情况下，求第 i 秒最多能消灭的怨灵的数量。

首先，这题可以用网络流做部分分。
考虑如何判断是否可行：
有一种显然的二分图匹配：把每个询问放在X部，怨灵放在Y部。
然后，把询问，怨灵分别拆点，进行区间连边，做匹配，如果有完美匹配，则可行。

但是，如果直接而匹配，会超时。我们考虑使用Hall定理进行判断。
Hall定理，就是说，在X部上取任意子集，与其相连的Y部结点都不小于集合大小。
如果我们要取子集，会变得更慢。

首先，把区间排序。
通过分析，我们知道：只要判断连续区间即可。

在X不连续时：
1、X部中一个点拆出来的点选子集，不如都选，因为都选不会使Y部大小增加。、
2、类似的，若X部中选的点在Y部中连续，那么X一定要连续，因为中间增加不会使Y部大小增加（区间是不会相互包含的）。
3、若X部中选的点在Y部中不连续，那可以拆成多个，拆的前后是等价的。

所以，我们只要快速判断连续区间就行。
记怨灵数量的前缀和为 sum[i]，区间消灭的怨灵数量的前缀和为 sumval[i]，区间左右端点为 l[i]，r[i]。
那么我们要求对于任意 \(1\leq L\leq R\leq m\),满足：

\(sumval[R]-sumval[L-1]\leq sum[r[R]]-sum[l[L]-1]\)
即\(sumval[R]-sum[r[R]]\leq sumval[L-1]-sum[l[L]-1]\)。

由于sum不变，只有sumval会边，所以当区间 i 被加入时：
若L-1>=i，R>=i，那么大小关系不变。
若L-1<i，R<i，那么大小关系不变。
只有对L-1<i，R>=i，即L<=i,R>=i 的情况产生影响。
所以区间 i 被加入时设结果为val[i]，那么对于任意 \(1\leq L\leq i\leq R\leq m\)，要求 val[i]满足：

\(sumval[R]+val[i]-sum[r[R]]\leq sumval[L-1]-sum[l[L]-1]\)
即\(val[i]\leq sumval[L-1]-sum[l[L]-1]]-(sumval[R]-sum[r[R]])\)。

那么，找到最小的 \(sumval[L-1]-sum[l[L]-1]\)，最大的 \(sumval[R]-sum[r[R]]\)，相减就是 val[i]的值。(当然,如果 k 小于这个值的话,val[i]就是 k)
使用线段树，支持区间加，区间最值，即可实现上述操作。
时间复杂度：\(O(mlog_2 m)\)。
感觉这个思路非常巧妙。

代码：

#include <stdio.h>
#include <stdlib.h>
#define inf 2000000000
int sz[100010],px[100010],wz[100010];
int xl[100010],xr[100010],xk[100010];
int cmp(const void*a,const void*b)
{
    return xl[*(int*)a]-xl[*(int*)b];
}
int z1[400010],z2[400010],su[100010];
int l1[400010],l2[400010];
void pushup(int i)
{
    z1[i]=z1[(i<<1)|1];
    z2[i]=z2[(i<<1)|1];
    if(z1[i<<1]<z1[i])
        z1[i]=z1[i<<1];
    if(z2[i<<1]>z2[i])
        z2[i]=z2[i<<1];
}
void pushdown(int i)
{
    l1[i<<1]+=l1[i];
    l1[(i<<1)|1]+=l1[i];
    z1[i<<1]+=l1[i];
    z1[(i<<1)|1]+=l1[i];
    l1[i]=0;
    l2[i<<1]+=l2[i];
    l2[(i<<1)|1]+=l2[i];
    z2[i<<1]+=l2[i];
    z2[(i<<1)|1]+=l2[i];
    l2[i]=0;
}
void jianshu(int i,int l,int r)
{
    if(l+1==r)
    {
        z1[i]=-su[xl[px[l]]-1];
        z2[i]=-su[xr[px[l]]];
        return;
    }
    int m=(l+r)>>1;
    jianshu(i<<1,l,m);
    jianshu((i<<1)|1,m,r);
    pushup(i);
}
void xiugai1(int i,int l,int r,int L,int R,int x)
{
    if(R<=l||r<=L)
        return;
    if(L<=l&&r<=R)
    {
        z1[i]+=x;
        l1[i]+=x;
        return;
    }
    int m=(l+r)>>1;
    pushdown(i);
    xiugai1(i<<1,l,m,L,R,x);
    xiugai1((i<<1)|1,m,r,L,R,x);
    pushup(i);
}
void xiugai2(int i,int l,int r,int L,int R,int x)
{
    if(R<=l||r<=L)
        return;
    if(L<=l&&r<=R)
    {
        z2[i]+=x;
        l2[i]+=x;
        return;
    }
    int m=(l+r)>>1;
    pushdown(i);
    xiugai2(i<<1,l,m,L,R,x);
    xiugai2((i<<1)|1,m,r,L,R,x);
    pushup(i);
}
void xiugai(int n,int j,int x)
{
    xiugai1(1,1,n+1,j+1,n+1,x);
    xiugai2(1,1,n+1,j,n+1,x);
}
int getz1(int i,int l,int r,int L,int R)
{
    if(R<=l||r<=L)
        return inf;
    if(L<=l&&r<=R)
        return z1[i];
    int m=(l+r)>>1;
    pushdown(i);
    int rt1=getz1(i<<1,l,m,L,R),rt2=getz1((i<<1)|1,m,r,L,R);
    if(rt1<rt2)return rt1;
    else return rt2;
}
int getz2(int i,int l,int r,int L,int R)
{
    if(R<=l||r<=L)
        return -inf;
    if(L<=l&&r<=R)
        return z2[i];
    int m=(l+r)>>1;
    pushdown(i);
    int rt1=getz2(i<<1,l,m,L,R),rt2=getz2((i<<1)|1,m,r,L,R);
    if(rt1>rt2)return rt1;
    else return rt2;
}
int getans(int n,int x)
{
    int x1=getz1(1,1,n+1,1,x+1);
    int x2=getz2(1,1,n+1,x,n+1);
    return x1-x2;
}
int main()
{
    freopen("rin.in","r",stdin);
    freopen("rin.out","w",stdout);
    int n,m;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&sz[i]);
        su[i]=su[i-1]+sz[i];
    }
    scanf("%d",&m);
    if(m==0)return 0;
    for(int i=1;i<=m;i++)
        scanf("%d%d%d",&xl[i],&xr[i],&xk[i]);
    for(int i=1;i<=m;i++)
        px[i]=i;
    qsort(px+1,m,sizeof(int),cmp);
    for(int i=1;i<=m;i++)
        wz[px[i]]=i;
    jianshu(1,1,m+1);
    for(int i=1;i<=m;i++)
    {
        int jg=getans(m,wz[i]);
        if(xk[i]<jg)
            jg=xk[i];
        xiugai(m,wz[i],jg);
        printf("%d\n",jg);
    }
    return 0;
}