hdu 1054 Strategic Game 匈牙利算法 || Hopcroft-Carp算法

原题链接：点击打开链接

题意：在节点上放置哨兵，若某点上已经放置了 则与它相连的点不用再放置，给出各点连接情况 求最少放置哨兵的数目---->最少覆盖点

思路：将各点拆分为两个点 分别位于两个点集（存图 注意），得到对称的二分图，最大匹配数增大1倍，最小覆盖数同样 所以结果要除以2；

匈牙利算法：

code：

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
#include<vector>
#include<stack>
using namespace std;

const int MAXN=1505;
int linker[MAXN];
int sum,n;
int vis[1505];
vector<int>map[MAXN];

bool dfs(int s)
{
    for(int u=0;u<map[s].size();u++)
    {
        int t=map[s][u];
        if(!vis[t])
        {
            vis[t]=1;
            if(linker[t]==-1 || dfs(linker[t]))
            {
                linker[t]=s;
                return true;
            }
        }
    }
    return false;
}
int hungary()
{
    int ans=0;
    memset(linker,-1,sizeof(linker));
    for(int i=0;i<n;i++)
    {
        memset(vis,0,sizeof(vis));
        if(dfs(i))
            ans++;

    }
   // cout<<"ans= "<<ans<<endl;
    return ans;
}
int main()
{//
   // freopen("i.txt","r",stdin);
    int u,k,v;
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<MAXN;i++)
           map[i].clear();

        for(int i=0;i<n;i++)
        {
            scanf("%d:(%d)",&u,&k);
            while(k--)
            {
                scanf("%d",&v);
                map[u].push_back(v);
                map[v].push_back(u);
            }
        }
        //sum=nn=0;
        printf("%d\n",hungary()/2);
    }
    return 0;
}

Hopcroft-Carp的算法，我不会：

code：

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<iostream>
#include<queue>
using namespace std;

const int MAXN=1505;
const int INF=1<<28;
int g[MAXN][MAXN],Mx[MAXN],My[MAXN],Nx,Ny;
int dx[MAXN],dy[MAXN],dis;
bool vst[MAXN];
bool searchP()
{
    queue<int>Q;
    dis=INF;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=0;i<Nx;i++)
        if(Mx[i]==-1)
        {
            Q.push(i);
            dx[i]=0;
        }
    while(!Q.empty())
    {
        int u=Q.front();
        Q.pop();
        if(dx[u]>dis)  break;
        for(int v=0;v<Ny;v++)
            if(g[u][v]&&dy[v]==-1)
            {
                dy[v]=dx[u]+1;
                if(My[v]==-1)  dis=dy[v];
                else
                {
                    dx[My[v]]=dy[v]+1;
                    Q.push(My[v]);
                }
            }
    }
    return dis!=INF;
}
bool DFS(int u)
{
    for(int v=0;v<Ny;v++)
       if(!vst[v]&&g[u][v]&&dy[v]==dx[u]+1)
       {
           vst[v]=1;
           if(My[v]!=-1&&dy[v]==dis) continue;
           if(My[v]==-1||DFS(My[v]))
           {
               My[v]=u;
               Mx[u]=v;
               return 1;
           }
       }
    return 0;
}
int MaxMatch()
{
    int res=0;
    memset(Mx,-1,sizeof(Mx));
    memset(My,-1,sizeof(My));
    while(searchP())
    {
        memset(vst,0,sizeof(vst));
        for(int i=0;i<Nx;i++)
          if(Mx[i]==-1&&DFS(i))  res++;
    }
    return res;
}

int main()
{
    int n;
    int u,k,v;
    while(scanf("%d",&n)!=EOF)
    {
        memset(g,0,sizeof(g));
        for(int i=0;i<n;i++)
        {
            scanf("%d:(%d)",&u,&k);
            while(k--)
            {
                scanf("%d",&v);
                g[u][v]=1;
                g[v][u]=1;
            }
        }
        Nx=Ny=n;
        int ans=MaxMatch();
        printf("%d\n",ans/2);
    }
    return 0;
}