Description
初始时有\(n(n\leq10^5)\)个孤立的点,依次向图中加入\(m(m\leq3\times10^5)\)条带权无向边。使得图中每个点的度数均为奇数的边集是合法的,其权值定义为集合中的最大边权。每次加入边后,询问权值最小的合法边集的权值,不存在合法边集时输出\(-1\)。
Solution
存在合法边集 \(\Leftrightarrow\) 每个连通块的大小均为偶数。如果某连通块大小为奇数,那么该块的总度数是奇数,但一条无向边会提供两个度数,所以不存在合法边集。如果大小为偶数,那么可以按照如下的方式构造合法边集:
求出该连通块的一棵生成树。由深到浅观察每个点\(u\),如果\(u\)与\(\{ch_u\}\)的边中有偶数个在边集中,那么将边\((u,fa[u])\)加入边集。
做到根的时候,因为除了根以外有奇数个度数为奇数的点,但总度数为偶数,所以根的度数必然是奇数。
该方法将每个偶数连通块划分为了若干个偶数连通块。对于划分后的每个连通块,其最小生成树就是其权值最小的合法边集。
那么现在我们可以进行操作了。首先开一个set<pair<int,int>>
记录lct中的边,first
记录边权,second
记录边的编号。初始用\(cnt\)表示剩余奇数连通块的个数,在其减至\(0\)前一直输出\(-1\)。
当存在合法边集后,加入边\((u,v)\)时,若\(u,v\)已经连通则维护最小生成树,否则在lct中连接\((u,v)\),同时维护set
。
如果该边被加入lct,那么尝试获得更优答案。从set
中最大的边开始,检查能否移除这条边;若能移除则移除并检查次大,不能则break
。因为如果当前的最大边不能移除,答案就不会变的更优了。最后输出set
中最大的边权。
时间复杂度\(O(mlog(n+m))\)。
Code
//Pastoral Oddities
#include <cstdio>
#include <set>
#include <algorithm>
using namespace std;
typedef pair<int,int> Ipair;
inline char gc()
{
static char now[1<<16],*S,*T;
if(S==T) {T=(S=now)+fread(now,1,1<<16,stdin); if(S==T) return EOF;}
return *S++;
}
inline int read()
{
int x=0; char ch=gc();
while(ch<‘0‘||‘9‘<ch) ch=gc();
while(‘0‘<=ch&&ch<=‘9‘) x=x*10+ch-‘0‘,ch=gc();
return x;
}
int const N=4e5+10;
int const INF=0x3F3F3F3F;
int n,m;
struct edge{int u,v;} ed[N];
set<Ipair> st;
int fa[N],ch[N][2],siz[N],isiz[N];
int val[N],maxL[N]; bool rev[N];
int wh(int p) {return p==ch[fa[p]][1];}
int notRt(int p) {return p==ch[fa[p]][wh(p)];}
void update(int p)
{
int ch0=ch[p][0],ch1=ch[p][1];
siz[p]=siz[ch0]+1+siz[ch1]+isiz[p];
int maxP=max(val[p],max(val[maxL[ch0]],val[maxL[ch1]]));
if(val[p]==maxP) maxL[p]=p;
else if(val[maxL[ch0]]==maxP) maxL[p]=maxL[ch0];
else if(val[maxL[ch1]]==maxP) maxL[p]=maxL[ch1];
}
void pushdw(int p) {if(rev[p]) swap(ch[p][0],ch[p][1]),rev[ch[p][0]]^=1,rev[ch[p][1]]^=1,rev[p]=false;}
void rotate(int p)
{
int q=fa[p],r=fa[q],w=wh(p);
fa[p]=r; if(notRt(q)) ch[r][wh(q)]=p;
fa[ch[q][w]=ch[p][w^1]]=q;
fa[ch[p][w^1]=q]=p;
update(q),update(p);
}
void pushRt(int p) {if(notRt(p)) pushRt(fa[p]); pushdw(p);}
void splay(int p)
{
pushRt(p);
for(int q=fa[p];notRt(p);rotate(p),q=fa[p]) if(notRt(q)) rotate(wh(p)^wh(q)?p:q);
}
void access(int p) {for(int q=0;p;q=p,p=fa[p]) splay(p),isiz[p]+=siz[ch[p][1]]-siz[q],ch[p][1]=q,update(p);}
void makeRt(int p) {access(p),splay(p),rev[p]^=1;}
void path(int p,int q) {makeRt(p),access(q),splay(q);}
void link(int i)
{
int p=ed[i].u,q=ed[i].v; makeRt(p),makeRt(q);
fa[p]=n+i,isiz[n+i]+=siz[p]; update(n+i);
fa[n+i]=q,isiz[q]+=siz[n+i]; update(q);
}
void clear(int p) {fa[p]=ch[p][0]=ch[p][1]=0,isiz[p]=0; update(p);}
void cut(int i)
{
int p=ed[i].u,q=ed[i].v; path(p,q);
fa[p]=ch[p][1]=0; ch[q][0]=0;
clear(n+i); update(p),update(q);
}
int find(int p) {access(p),splay(p); while(ch[p][0]) p=ch[p][0]; return p;}
int odd(int p) {return siz[p]+1>>1&1;}
int main()
{
n=read(),m=read();
for(int i=1;i<=n;i++) val[i]=-INF,update(i);
for(int i=1,cnt=n;i<=m;i++)
{
int u=read(),v=read(),w=read();
ed[i].u=u,ed[i].v=v; val[n+i]=w,maxL[n+i]=n+i;
if(find(u)==find(v))
{
path(u,v); int id=maxL[v]-n,w0=val[id+n];
if(w<w0) cut(id),link(i),st.erase(make_pair(w0,id));
else {printf("%d\n",cnt?-1:st.rbegin()->first); continue;}
}
else
{
makeRt(u),cnt-=odd(u); makeRt(v),cnt-=odd(v);
link(i),cnt+=odd(v);
}
st.insert(make_pair(val[i+n],i));
if(cnt) {puts("-1"); continue;}
while(true)
{
int id=st.rbegin()->second; u=ed[id].u,v=ed[id].v;
path(u,v); if(odd(u)) break;
cut(id); st.erase(*st.rbegin());
}
printf("%d\n",st.rbegin()->first);
}
return 0;
}
P.S.
第一次用set
,不大懂...%%%zhx