HDOJ1045解题报告【二分图最大匹配】

题目地址：

　　http://acm.hdu.edu.cn/showproblem.php?pid=1045

题目概述：

　　在图上有一些障碍物，要求在图上放上尽可能多的点使得每一行每一列没有两个点直接相连。

大致思路：

　　类似八皇后问题，而且数据范围很小不清楚搜索加剪枝能不过，有兴趣可以试一下。

　　我的做法是用二分图，建图就在每个点的x坐标和y坐标之间连一条边，障碍物怎么处理呢？事实上在一行里只需要将一个障碍物的两边设成不同的行就好了，列也是同样的道理。

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <ctime>
#include <map>
#include <stack>
#include <queue>
#include <cstring>
#include <algorithm>
using namespace std;

#define sacnf scanf
#define scnaf scanf
#define maxn  10
#define maxm 26
#define inf 1061109567
#define Eps 0.00001
const double PI=acos(-1.0);
#define mod 7
#define MAXNUM 10000
void Swap(int &a,int &b) {int t=a;a=b;b=t;}
int Abs(int x) {return (x<0)?-x:x;}
typedef long long ll;
typedef unsigned int uint;

char G[maxn][maxn];
int rG[maxn][maxn][2];
vector<int> tG[10*maxn];
int l[10*maxn],vis[10*maxn];

bool dfs(int u)
{
    int len=tG[u].size();
    for(int i=0;i<len;i++)
    {
        int v=tG[u][i];
        if(!vis[v])
        {
            vis[v]=1;
            if(l[v]==-1||dfs(l[v]))
            {
                l[v]=u;return true;
            }
        }
    }
    return false;
}

int hungary(int n)
{
    int ans=0;
    memset(l,-1,sizeof(l));
    for(int i=1;i<=n;i++)
    {
        memset(vis,0,sizeof(vis));
        if(dfs(i)) ans++;
    }
    return ans;
}

int main()
{
    //freopen("data.in","r",stdin);
    //freopen("data.out","w",stdout);
    //clock_t st=clock();
    int n;
    while(~scanf("%d",&n))
    {
        if(n==0) break;
        for(int i=1;i<=n;i++)
        {
            getchar();
            for(int j=1;j<=n;j++)
            {
                scanf("%c",&G[i][j]);
            }
        }
        int x=1;bool flag=false,flag2=false;
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                if(G[i][j]=='X')
                {
                    if(flag) x++;
                    flag=false;continue;
                }
                else {flag=true;flag2=true;}
                rG[i][j][0]=x;
            }
            if(flag2) x++;flag=flag2=false;
        }
        int y=1;flag=false;flag2=false;
        for(int j=1;j<=n;j++)
        {
            for(int i=1;i<=n;i++)
            {
                if(G[i][j]=='X')
                {
                    if(flag) y++;
                    flag=false;continue;
                }
                else {flag=true;flag2=true;}
                rG[i][j][1]=y;
            }
            if(flag2) y++;flag=flag2=false;
        }
        int t=x+y;
        for(int i=1;i<=t;i++) tG[i].clear();
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                if(G[i][j]!='X')
                {
                    int tx=rG[i][j][0],ty=rG[i][j][1];
                    tG[tx].push_back(x+ty);
                    tG[x+ty].push_back(tx);
                }
            }
        }
        printf("%d
",hungary(x));
    }
    //clock_t ed=clock();
    //printf("

Time Used : %.5lf Ms.
",(double)(ed-st)/CLOCKS_PER_SEC);
    return 0;
}