HDU 4973 A simple simulation problem.

最新推荐文章于 2017-08-16 10:04:46 发布

原创最新推荐文章于 2017-08-16 10:04:46 发布 · 544 阅读

0 ·

CC 4.0 BY-SA版权

线段树专栏收录该内容

3 篇文章

订阅专栏

本文探讨了线段树在解决区间更新和查询问题时的优化策略，通过具体实例展示了如何在多校赛中有效利用线段树进行复杂操作，并在时间限制下成功解决问题。

一道线段树多校的时候有点想歪了最后时间又不够结果没写出~~。。。后来改了20分钟就过了

#include <stdio.h>
#include <math.h>
#include <iostream>
#include <string.h>
#include <algorithm>
#include <vector>
#include <map>
using namespace std;
int n, m;
struct SEG
{
    int l, r;
    __int64 va, maxv, col;
} seg[50005<<2];
void build(int l,int r,int rt)
{
    seg[rt].l=l, seg[rt].r=r;
    seg[rt].col=1;
    seg[rt].maxv=1;
    if(l==r){
        seg[rt].va=1;
        return;
    }
    int mid=(l+r)>>1;
    build(l,mid,rt<<1);
    build(mid+1,r,rt<<1|1);
    seg[rt].va=seg[rt<<1].va+seg[rt<<1|1].va;
}
void push_down(int rt)
{
    if(seg[rt].col!=1){
        __int64 col=seg[rt].col;
        int ls, rs;
        ls=rt<<1, rs=rt<<1|1;
        seg[ls].va*=col;
        seg[rs].va*=col;
        seg[ls].maxv*=col;
        seg[rs].maxv*=col;
        seg[ls].col*=col, seg[rs].col*=col;
        seg[rt].col=1;
    }
}
void push_up(int rt)
{
    seg[rt].va=seg[rt<<1].va+seg[rt<<1|1].va;
    seg[rt].maxv=max(seg[rt<<1].maxv,seg[rt<<1|1].maxv);
}
void update(__int64 L,__int64 R,int rt)
{
    if(seg[rt].l==seg[rt].r)
    {
        seg[rt].va+=R-L+1;
        seg[rt].maxv+=R-L+1;
        return;
    }
    if(L==1&&seg[rt].va==R){
        seg[rt].col=seg[rt].col*2;
        seg[rt].maxv=seg[rt].maxv*2;
        seg[rt].va=seg[rt].va*2;
        return;
    }
    push_down(rt);
    int ls=rt<<1, rs=rt<<1|1;
    __int64 tmp=seg[ls].va;
    if(R<=tmp) update(L,R,ls);
    else if(L>tmp) update(L-tmp,R-tmp,rs);
    else{
        update(L,tmp,ls);
        update(1,R-tmp,rs);
    }
    push_up(rt);
}

__int64 qmax(__int64 L,__int64 R,int rt)
{
    if(seg[rt].l==seg[rt].r)
    {
        return R-L+1;
    }
    if(L==1&&seg[rt].va==R){
        return seg[rt].maxv;
    }
    push_down(rt);
    __int64 ret=1;
    int ls=rt<<1, rs=rt<<1|1;
    if(R<=seg[ls].va) ret=qmax(L,R,ls);
    else if(L>seg[ls].va) ret=qmax(L-seg[ls].va,R-seg[ls].va,rs);
    else{
        ret=qmax(L,seg[ls].va,ls);
        ret=max(ret,qmax(1,R-seg[ls].va,rs));
    }
    return ret;
}
int main()
{
    int T, i, j, l, r, ri=0;
    __int64 L, R;
    __int64 maxv, lsum, rsum;
    char op[11];
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        build(1,n,1);
        printf("Case #%d:\n",++ri);
        for(i=1;i<=m;i++){
            scanf("%s%I64d%I64d",op,&L,&R);
            if(op[0]=='D'){
                update(L,R,1);
            }else{
                printf("%I64d\n",qmax(L,R,1));
            }
        }
    }
    return 0;
}