zoukankan      html  css  js  c++  java
  • P4169 [Violet]天使玩偶/SJY摆棋子 cdq分治

      

     按照时间分治  只要对四个方向cdq一遍即可

    // luogu-judger-enable-o2
    // luogu-judger-enable-o2
    #include<bits/stdc++.h>
    using namespace std;
    #define rep(i,a,b) for(int i=(a);i<=(b);i++)
    #define repp(i,a,b) for(int i=(a);i>=(b);--i)
    #define ll long long
    #define see(x) (cerr<<(#x)<<'='<<(x)<<endl)
    #define inf 0x3f3f3f3f
    #define CLR(A,v)  memset(A,v,sizeof A)
    //////////////////////////////////
    const int N=1e7+10;
    int n,m,S,x,y,x2,y2,z,t[N],op,tim,cnt,cntq,ans[N],lmax,rmax,tot,maxx,maxy;
    void add(int x,int v){for(;x<=N;x+=x&-x)t[x]=max(t[x],v);}
    int qmax(int x){int ans=0;for(;x;x-=x&-x)ans=max(ans,t[x]);return ans;}
    void clean(int x){for(;x<=N;x+=x&-x)t[x]=0;}
    struct node
    {
        int x,y,t,op,id;
    }s[N],temp[N],a[N];
    bool cmp(node a,node b){return a.x<b.x||a.x==b.x&&a.y<b.y;}
    void cdq(int l,int r)
    {
        if(l==r)return;int mid=(l+r)>>1;
        cdq(l,mid);cdq(mid+1,r);
        sort(s+l,s+mid+1,cmp);sort(s+mid+1,s+r+1,cmp);
        int j=l;
        rep(i,mid+1,r)if(s[i].op==2)
        {
            while(j<=mid&&s[j].x<=s[i].x)
            {
                if(s[j].op==1)add(s[j].y,s[j].x+s[j].y);j++;
            }
            int tmp=qmax(s[i].y);
            if(tmp)
            ans[s[i].id]=min(ans[s[i].id],s[i].x+s[i].y-tmp);
        }
        while(--j>=l)
        {
            if(s[j].op==1)clean(s[j].y);
        }
    }
    bool cmpt(node a,node b){return a.t<b.t;}
    void init()
    {
        lmax=rmax=tot=0;
        rep(i,1,cnt)if(s[i].op==2)lmax=max(lmax,s[i].x),rmax=max(rmax,s[i].y);
        rep(i,1,cnt)
        {
            if(s[i].x<=lmax&&s[i].y<=rmax)s[++tot]=s[i];
        }
    }
    
    int main()
    {
        scanf("%d%d",&n,&m);
        rep(i,1,n)
        {
            scanf("%d%d",&x,&y);x++;y++;
            s[++cnt]=(node){x,y,cnt,1,0};
        }
        while(m--)
        {
            scanf("%d%d%d",&op,&x,&y);x++;y++;
            if(op==1)
            {   
                s[++cnt]=(node){x,y,cnt,1,0};
            }
            else 
            {
                s[++cnt]=(node){x,y,cnt,2,++cntq};
            }
        }
        maxx=maxy=1e6+2;
        rep(i,1,cntq)ans[i]=1e9;
    
        rep(i,1,cnt)a[i]=s[i];
        init();
        cdq(1,tot);
    
        rep(i,1,cnt)s[i]=a[i],s[i].x=maxx-s[i].x;
        init();
        cdq(1,tot);
    
        rep(i,1,cnt)s[i]=a[i],s[i].y=maxy-s[i].y;
        init();
        cdq(1,tot);
    
        rep(i,1,cnt)s[i]=a[i],s[i].x=maxx-s[i].x,s[i].y=maxy-s[i].y;
        init();
        cdq(1,tot);
    
        rep(i,1,cntq)
        printf("%d
    ",ans[i]);
        return 0;
    }
    View Code
  • 相关阅读:
    网络安全等级保护系统定级流程与示例
    分析设计之类图
    分析设计中用例图、类图与时序图关系
    工业网络安全产品应用场景
    网络安全拟态防御技术
    浅说:网络空间拟态防御是个什么鬼?
    邬江兴院士:鲁棒控制与内生安全
    邬江兴院士:工业互联网安全&拟态防御
    邬江兴院士:多模态智慧网络与内生安全
    Samba CVE-2017-7494验证实验
  • 原文地址:https://www.cnblogs.com/bxd123/p/11454784.html
Copyright © 2011-2022 走看看