[HDU3861]The King's Problem（Tarjan缩点+二分图匹配）

题目描述

传送门
题意：将一个图划分成若干个部分，要求：1、可以互达的点在一个部分中 2、在一个部分中的点至少单向到达。

题解

题意比较坑。
我们发现互达的点一定在一个部分中，很容易想到Tarjan缩点。但是不互达的点也有可能在一个部分中，所以缩点之后就变成了一个有向图的最小路径覆盖问题，用二分图匹配就可以解决。
刚开始写了dinic发现TLE，不情愿地改成了二分图匹配。

代码

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;

const int max_n=5e3+5;
const int max_m=1e5+5;
const int max_e=max_m*2;
const int INF=1e9;

int T,n,m,x,y,cnt,N,temp,ans;
int tot,point[max_n],nxt[max_e],v[max_e];
int Dfn[max_n],Low[max_n],belong[max_n],strack[max_n];
bool vis[max_n];
int girl[max_n];
struct hp{int x,y;}e[max_m];

inline void clear(){
    cnt=N=temp=ans=0;
    for (int i=1;i<=n;++i) Dfn[i]=Low[i]=belong[i]=strack[i]=vis[i]=girl[i]=0;
}
inline void add(int x,int y){++tot;nxt[tot]=point[x];point[x]=tot;v[tot]=y;}
inline void tarjan(int now){
    Dfn[now]=Low[now]=++N; vis[now]=true; strack[++temp]=now;
    for (int i=point[now];i;i=nxt[i])
      if (!Dfn[v[i]]){
        tarjan(v[i]);
        Low[now]=min(Low[now],Low[v[i]]);
      }
      else if (vis[v[i]]) Low[now]=min(Low[now],Dfn[v[i]]);
    if (Dfn[now]==Low[now]){
        ++cnt;
        while (strack[temp]!=now) belong[strack[temp]]=cnt,vis[strack[temp]]=false,--temp;
        belong[strack[temp]]=cnt,vis[strack[temp]]=false,--temp;
    }
}

inline bool find(int x){
    for (int i=point[x];i;i=nxt[i])
      if (!vis[v[i]]){
        vis[v[i]]=true;
        if (girl[v[i]]==0||find(girl[v[i]])){
            girl[v[i]]=x;
            return true;
        }
      }
    return false;
}

int main(){
    scanf("%d",&T);
    while (T--){
        tot=0; memset(point,0,sizeof(point)); memset(nxt,0,sizeof(nxt)); memset(v,0,sizeof(v));
        scanf("%d%d",&n,&m);
        clear();
        for (int i=1;i<=m;++i){
            scanf("%d%d",&x,&y),add(x,y);
            e[i].x=x; e[i].y=y;
        }
        for (int i=1;i<=n;++i) if (!Dfn[i]) tarjan(i);
        tot=0; memset(point,0,sizeof(point)); memset(nxt,0,sizeof(nxt)); memset(v,0,sizeof(v));
        for (int i=1;i<=m;++i)
          if (belong[e[i].x]!=belong[e[i].y])
            add(belong[e[i].x],belong[e[i].y]);
        N=cnt*2;
        for (int i=1;i<=cnt;++i){
            memset(vis,0,sizeof(vis));
            if (find(i)) ans++;
        }
        printf("%d\n",cnt-ans);
    }
}

总结

情况允许的时候也要考虑写二分图，虽然网络流写的熟练，但是匈牙利省空间常数小。