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;
}