暑训记录7.16

早上

早上先补下昨晚晚训的c，然后去写网络流吧

Skyscraper Gym - 102220H

题目大意

给出一个数组，每个位置有值ai，我们每次操作可以将任意l到r区间每个值减一。
有两个事件：
1.将区间l到r的所有数加1
2.输入将l到r区间每个值减为0的最小操作数

题目思路

我们考虑第二个事件：
对于l到r区间，第一个位置需要ai次操作，清除第一个位置后，第二个位置需要的操作数就是max(0,a[1]-a[2])……第i个位置就是max(0,a[i-1]-a[i])。那么每个位置的值就可以表示为max(0,a[i-1]-a[i]),而a[i-1]-a[i]可以看作a[i]的差分数组b[i]的第i项。那么我们对于这个区间的答案就是l到r的差分数组正值的和。
我们可以用线段树记录a[i]的差分数组b[i]，和答案需要的数组值c[i]（差分数组正值不变，负值为0），这样我们要输出l到r的最小操作数时，只需要求得a[l]的值（即差分数组的b[1]到b[l]的和）加上l+1到r的差分数组正值总和就好了（即c[l+1]到c[r]的和），这样我们就完成了事件2。
对于事件1，我们可以利用差分更新区间加的性质，只改变b[l]和b[r+1]就好了

ac代码

const int maxn = 2e6+10;
const int inf = 1e9;
const ll llinf =1e18+10;
const ll mod = 1e9+7;

int n,m;
ll a[maxn];
ll b[maxn];
ll t[maxn][10];
ll lz[maxn][10];

void build(int rt,int l,int r)
{
    t[rt][2]=0;
    t[rt][1]=0;
    if(l==r)
    {
        t[rt][2]=b[l];
        t[rt][1]=max((ll)0,b[l]);
        return;
    }
    int mid=(l+r)>>1;
    build(lson,l,mid);
    build(rson,mid+1,r);
    t[rt][1]=t[lson][1]+t[rson][1];
    t[rt][2]=t[lson][2]+t[rson][2];
}



void upd(int rt,int l,int r,int p,ll v)
{

    if(l==r)
    {
        t[rt][2]+=v;
        t[rt][1]=max(t[rt][2],(ll)0);
        return;
    }
    int mid=(l+r)>>1;
    if(mid>=p)
        upd(lson,l,mid,p,v);
    else
        upd(rson,mid+1,r,p,v);
    t[rt][1]=t[lson][1]+t[rson][1];
    t[rt][2]=t[lson][2]+t[rson][2];
}

ll query(int rt,int l,int r,int L,int R)
{
    if(L<=l&&r<=R)
    {
        return t[rt][1];
    }
    int mid=(l+r)>>1;

    ll ans=0;
    if(L<=mid)
        ans+=query(lson,l,mid,L,R);
    if(R>mid)
        ans+=query(rson,mid+1,r,L,R);
    return ans;
}

ll query2(int rt,int l,int r,int L,int R)
{
    if(L<=l&&r<=R)
    {
        return t[rt][2];
    }
    int mid=(l+r)>>1;

    ll ans=0;
    if(L<=mid)
        ans+=query2(lson,l,mid,L,R);
    if(R>mid)
        ans+=query2(rson,mid+1,r,L,R);
    return ans;
}



int main()
{
    int _,cas=0;
    scanf("%d",&_);
    while(_--)
    {

        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
        }
        b[1]=a[1];
        for(int i=2;i<=n;i++)b[i]=a[i]-a[i-1];
        build(1,1,n);
        while(m--)
        {
            int t,l,r,x;
            scanf("%d%d%d",&t,&l,&r);
            if(t==1)
            {
                scanf("%d",&x);
                upd(1,1,n,l,x);
                if(r!=n)
                {
                    upd(1,1,n,r+1,-1*x);
                }

            }else
            {
                ll ans=query2(1,1,n,1,l);
                //printf("%lld\n",ans);
                if(l!=r)
                    ans+=query(1,1,n,l+1,r);
                printf("%lld\n",ans);
            }
        }
    }
}

写代码40分钟，调bug半个小时，t2，再改bug半小时，我们线段树真是太好玩辣！

网络流24题目前施工进度：3/24

下午

网络流24题目前施工进度：5/24

晚上

晚训b题数论果断跳过，c题简单结论题但是边的长度我以为都是1，然后卡了一个多小时。。。
网络流24题目前施工进度：6/24