王の問題

都市およびnあり、mはチャネルを有しています。道路は市から来た場合はどこへ行く方法はありません、その後、この都市の最初のすべての方法の国に、各都市は唯一の国に属し、当然のリング上の都市は、同じ国に属することができます。Q.は、これらの都市の多くは、少なくとも二つの国に分かれていますどのように。

図面環の存在、および二部グラフには、非循環有向グラフであるのでときに正しく、強いデフレとの最初の点でも非環式有向グラフ再構成できます。そして、二部マッチング。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
#include<vector>
#define INF 0x3f3f3f3f
using namespace std;
const int N=10005;
vector<int> g[N];
stack<int> q;
int n,m,dfs_clock,scc_cnt;
int dfn[N],low[N],sccno[N],link[N];
bool vis[N],map[N][N];
void Init()         
{
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(sccno,0,sizeof(sccno));
    memset(map,false,sizeof(map)); 
    dfs_clock=scc_cnt=0;
    for(int i=0;i<=n;i++)
      g[i].clear();
    while(!q.empty())
      q.pop();
}
void tarjan(int u)       
{
    dfn[u]=low[u]=++dfs_clock;
    q.push(u);
    for(int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(!sccno[v])
          low[u]=min(low[u],dfn[v]);
    } 
    if(low[u]==dfn[u])
    {
        scc_cnt++;
        while(!q.empty())
        {
            int v=q.top();
            q.pop();
            sccno[v]=scc_cnt;
            if(v==u)
              break;
        } 
    }
} 
bool dfs(int u)     
{
    for(int i=1;i<=n;i++)
    {
        if(!vis[i]&&map[u][i])
        {
            vis[i]=true;
            if(!link[i]||dfs(link[i]))
            {
                link[i]=u;
                return true;
            }
        }
    }
    return false;
}
int hungary()         
{
    int ans=0;
    memset(link,0,sizeof(link));
    for(int i=1;i<=n;i++)
    {
        memset(vis,false,sizeof(vis));
        if(dfs(i))
          ans++;
    }
    return scc_cnt-ans;
}
int main()
{
    int t,u,v;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        Init();
        while(m--)
        {
            scanf("%d%d",&u,&v);
            g[u].push_back(v);
        }
        for(int i=1;i<=n;i++)
          if(!dfn[i])
            tarjan(i);
        for(int u=1;u<=n;u++)          
          for(int i=0;i<g[u].size();i++)
          {
              int v=g[u][i];
              if(sccno[u]!=sccno[v])
                map[sccno[u]][sccno[v]]=1;
          } 
        printf("%d\n",hungary());
    }
    return 0;
}
公開された44元の記事 ウォン称賛35 ビュー765

おすすめ

転載: blog.csdn.net/huangziguang/article/details/104482839