【MT原创题 丘比特的烦恼】Hopcroft-Karp算法、tarjan

题目H8OJ上有，这里就不罗嗦了

题目不难，简单讲一下解法

先求个最大匹配，再把无向边定向，匹配边由x集->y集，非匹配边由y集->x集，这样定向了之后，任意一条从x集出发y集结束的路径就对应了一条增广路

求出这个图的强连通分量，如果一条边的两个端点处于同一个分量内，那这条边其实没什么影响力，删了这条边原图仍然存在完备匹配

而如果这条边的两个端点不在一个分量内，那这条边如果是匹配边那就一定会被选，如果不是匹配边那就一定不会被选

可是，由于这个题的数据规模比较重口味，N<=200000; M<=600000，传统的匈牙利算法已经无法胜任了，我们需要更NB的求最大匹配的算法，这个算法就是HK算法

这个算法名字看起来很唬人，实际上没什么

它的思想跟dinic有点像，大概就是每次先广搜出一个距离标号，然后匈牙利的增广过程就只找标号等于自己加一的点，一次标号多次增广

可以证明增广次数是O(sqrt(n))的，每次增广的复杂度是O(m)，总的复杂度就只有O(sqrt(n)*m)

一个细节：

由于要保证广搜出来的路也是交错轨，所以广搜时x集的点只扩展y集的已盖点，而y集的点只扩展它对应的那个匹配点，但距离标号还是要更新的

哎呀，感觉这里有点说不清，参考代码吧

program syj; var n,m,t,i,j,k,e,st,ed:longint; ok:boolean; a,h,d,q,low,f:array[0..400005]of longint; next,point:array[0..1200005]of longint; aa,bb:array[0..600005]of longint; v:array[0..400005]of boolean; u:array[0..1200005]of boolean; function find(i:longint):boolean; var j,k:longint; begin j:=h[i]; while j<>0 do begin k:=point[j]; if (d[k]=d[i]+1)and not v[k] then begin v[k]:=true; if (a[k]=0)or find(a[k]) then begin a[k]:=i;a[i]:=k; exit(true); end; end; j:=next[j]; end; find:=false; end; procedure dfs(i:longint); var j,k:longint; begin inc(t); d[i]:=t;low[i]:=t; inc(ed); q[ed]:=i; v[i]:=true; j:=h[i]; while j<>0 do begin if u[j] then begin j:=next[j];continue; end; u[j]:=true; k:=point[j]; if d[k]=0 then begin dfs(k); if low[k]<low[i] then low[i]:=low[k]; end else if v[k]and(d[k]<low[i]) then low[i]:=d[k]; // if low[k]>d[i] then g[j>>1]:=1+ord(i>k); j:=next[j]; end; if d[i]=low[i] then while q[ed+1]<>i do begin v[q[ed]]:=false;f[q[ed]]:=i; dec(ed); end; end; begin assign(input,'cupid.in');reset(input); assign(output,'cupid.out');rewrite(output); readln(n,m); e:=1; for m:=1 to m do begin readln(j,k); aa[m]:=j;bb[m]:=k; inc(e); next[e]:=h[j];point[e]:=k+n;h[j]:=e; inc(e); next[e]:=h[k+n];point[e]:=j;h[k+n]:=e; end; ok:=true; while ok do begin st:=0;ed:=0; fillchar(d,sizeof(d),$FF); for i:=1 to n do if a[i]=0 then begin inc(ed);q[ed]:=i;d[i]:=0; end; while st<ed do begin inc(st);i:=q[st]; if i<=n then begin j:=h[i]; while j<>0 do begin k:=point[j]; if d[k]=-1 then begin d[k]:=d[i]+1; if a[k]>0 then begin inc(ed);q[ed]:=k; end; end; j:=next[j]; end; end else begin inc(ed);q[ed]:=a[i];d[a[i]]:=d[i]+1; end; end; ok:=false; fillchar(v,sizeof(v),0); for i:=1 to n do begin if a[i]=0 then ok:=find(i)or ok; end; end; // for i:=1 to n+n do if a[i]=0 then writeln(i); for i:=1 to n+n do begin j:=h[i]; while j<>0 do begin if (point[j]<>a[i])xor(i>n) then u[j]:=true; j:=next[j]; end; end; ed:=0; t:=0; fillchar(v,sizeof(v),0); fillchar(d,sizeof(d),0); for i:=1 to n do if d[i]=0 then dfs(i); for i:=1 to m do if f[aa[i]]=f[bb[i]+n] then write(0) else if a[aa[i]]=bb[i]+n then write(1) else write(2); writeln; close(input);close(output); end.