Codeforces Round #673 (Div. 1) D.Graph and Queries

D.Graph and Queries
题意:
一开始n个点m条边的无向图,这个无向图不一定是联通的,每个点都有一个点权然后有q个操作,操作1给出一个x,然后询问x的所在的连通图中权值最大的点的权值,然后将这个点的点权改为0,操作2给出一个x,将第x条边断开。
思路:
可以离线,所以可以倒着合并每条断开的边,每次合并断开的边就有可能合并两个集合,因此可以用链表代表两个集合,合并时直接用splic函数o(1)合并,所以总结就是倒着来两遍,第一遍用链表合并,使得每次断开一条边得到的两个集合是相邻的,这样就可以进行线段树区间操作了,第二遍就可以维护当前的询问该询问哪个区间,然后正着来一遍询问和修改即可,反正询问的区间是固定不变的了。
ps:可能还有一个克鲁斯重构树+dfs序的解法,待补。
重构树思路(无代码,纯口嗨):克鲁斯卡尔重构树裸题,没有删除的边的边权都为0,然后倒着来,删除的边的边权依次递增赋值,然后维护一个dfs序,这样正着跑一遍询问即可,操作1就是单纯的区间修改和查询了。

#include<iostream>
#include<cstdio>
#include<list>
#include<algorithm>
using namespace std;
const int MAX_N=501000;
struct skt{
    int op,x,l,r;
}qy[MAX_N];
struct edges{
    int x,y;
}edge[MAX_N];
struct kt{
    int l,r;
}d[MAX_N];//记录区间范围 
int b[MAX_N],vb[MAX_N];
int p[MAX_N];
bool vis[MAX_N],viss[MAX_N];
int pre[MAX_N];
int find(int x){
    if(x!=pre[x])
    pre[x]=find(pre[x]);
    return pre[x];
}
list<int>li[201000];
struct node{
     int l,r,max;
}a[201000*4];
int pos[MAX_N];
void update(int k){
     a[k].max=max(a[k<<1].max,a[k<<1|1].max);
}
void build(int k,int l,int r){
     a[k].l=l;a[k].r=r;
     if(l==r){
          a[k].max=p[b[l]];
          return;
     }
     int mid=(l+r)>>1;
     build(k<<1,l,mid);
     build(k<<1|1,mid+1,r);
     update(k);
}
void change(int k,int x,int y){
     if(a[k].l==a[k].r){
          a[k].max=y;
          return;
     }
     int mid=(a[k].l+a[k].r)>>1;
     if(x<=mid)
     change(k<<1,x,y);
     else
     change(k<<1|1,x,y);
     update(k);
}
int query(int k,int l,int r){
    if(a[k].l>=l&&a[k].r<=r)
     return a[k].max;
     int mid=(a[k].l+a[k].r)>>1;
     int x=0;
     if(r>mid)
     x=max(x,query(k<<1|1,l,r));
     if(l<=mid)
     x=max(x,query(k<<1,l,r));
     return x;
}
int main(void){
    int n,m,q,i;
    scanf("%d%d%d",&n,&m,&q);
    for(i=1;i<=n;i++){
        scanf("%d",&p[i]);
        pos[p[i]]=i;
    }
    for(i=1;i<=m;i++)
    scanf("%d%d",&edge[i].x,&edge[i].y);
    for(i=1;i<=q;i++){
        scanf("%d%d",&qy[i].op,&qy[i].x);
        if(qy[i].op==2)
        vis[qy[i].x]=true;
    }
    for(i=1;i<=n;i++)
    li[i].insert(li[i].end(),i),pre[i]=i;
    for(i=1;i<=m;i++){
        if(vis[i])
        continue;
        int fx=find(edge[i].x);
        int fy=find(edge[i].y);
        if(fx!=fy){
            pre[fx]=fy;
            li[fy].splice(li[fy].end(),li[fx]);
            //li[fy].insert(li[fy].end(),li[fx].begin(),li[fx].end());
            //li[fx].clear();
        }
    }
    for(i=q;i>=1;i--){
        if(qy[i].op==2){
            int fx=find(edge[qy[i].x].x);
            int fy=find(edge[qy[i].x].y);
            if(fx!=fy){
                pre[fx]=fy;
                li[fy].splice(li[fy].end(),li[fx]);
                //li[fy].insert(li[fy].end(),li[fx].begin(),li[fx].end());
                //li[fx].clear();
            }
        }
    }
    int cnt=0;
    for(i=1;i<=n;i++){
        int fx=find(i);
        if(viss[fx])
        continue;
        viss[fx]=true;
        list<int>::iterator it=li[fx].begin();
        while(it!=li[fx].end()){
            b[++cnt]=*it;
            vb[b[cnt]]=cnt;
            it++;
        }
    }
    build(1,1,n);
    for(i=1;i<=n;i++)
    pre[i]=i,d[i].l=d[i].r=vb[i];
    for(i=1;i<=m;i++){
        if(vis[i])
        continue;
        int fx=find(edge[i].x);
        int fy=find(edge[i].y);
        if(fx!=fy){
            pre[fx]=fy;
            d[fy].l=min(d[fx].l,d[fy].l);
            d[fy].r=max(d[fx].r,d[fy].r);
        }
    }
    for(i=q;i>=1;i--){
        if(qy[i].op==2){
            int fx=find(edge[qy[i].x].x);
            int fy=find(edge[qy[i].x].y);
            if(fx!=fy){
                pre[fx]=fy;
                d[fy].l=min(d[fx].l,d[fy].l);
                d[fy].r=max(d[fx].r,d[fy].r);
            }
        }
        else if(qy[i].op==1){
            int fx=find(qy[i].x);
            qy[i].l=d[fx].l;
            qy[i].r=d[fx].r;
            //cout<<i<<" "<<fx<<" "<<d[fx].l<<" "<<d[fx].r<<" !!!!\n";
        }
    }
    for(i=1;i<=q;i++){
        if(qy[i].op==1){
            int ans=query(1,qy[i].l,qy[i].r);
            printf("%d\n",ans);
            if(ans){
                int x=pos[ans];
                change(1,vb[x],0);
            }
        }
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值