poj2226Muddy Fields（二分图的最小点覆盖）

http://poj.org/problem?id=2226

模板抄错，结果错了一下午。将连续横行上的*和连续纵桁上点分别看做两个点集，连接两集合中点的边就是图中的点。即求最小顶点覆盖最多的边。最小点覆盖=最大匹配。

有篇图片的讲解，挺好，很形象。http://ip96cns.blog.163.com/blog/static/170095192201117465473/

#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int map[1000][1000];
int r,c,link[1000],vis[1000];
int find(int x,int n)//二分匹配模板
{
    int i,j;
    for(i = 1; i <= n ; i++)
    {
        if(map[x][i]&&!vis[i])
        {
            vis[i] = 1;
            if(link[i]==0||find(link[i],n))
            {
                link[i] = x;
                return 1;
            }
        }
    }
    return 0;
}
int main()
{
    int i,j,k,a[100][100],b[100][100];
    char s[100][100];
    while(cin>>r>>c)
    {
        memset(link,0,sizeof(link));

        memset(map,0,sizeof(map));
        k = 0;
        int sum = 0;
        int ii=1,jj=1;
        for(i = 1 ; i <= r ; i++)
        {
            getchar();
            for(j = 1; j <= c ; j++)
            {
                scanf("%c", &s[i][j]);
            }
        }
        for(i = 1; i <= r ; i++)
        for(j = 1; j <= c ; j++)//建立X集合中的点
        {
            if(j==1&&i!=1&&s[i-1][c]=='*')
            ii++;
            if(s[i][j]=='*')
            {
                a[i][j] = ii;
            }
            if(j!=c&&(s[i][j]=='*'&&s[i][j+1]!='*'))
                ii++;
        }
        for(i = 1; i <= c ; i++)//建立Y集合中的点
        for(j = 1; j <= r ; j++)
        {
            if((j==1&&i!=1&&s[r][i-1]=='*'))
            jj++;
            if(s[j][i]=='*')
            b[j][i] = jj;
            if(j!=r&&(s[j][i]=='*'&&s[j+1][i]!='*'))
            {
                jj++;
            }
        }
        for(i = 1; i <= r ; i++)
        for(j = 1; j <= c ; j++)//建立连接两集合点的边，即图中的"*"
        {
            if(s[i][j]=='*')
            {
                map[a[i][j]][b[i][j]] = 1;
            }
        }
        for(i = 1; i <= ii ; i++)
        {
            memset(vis,0,sizeof(vis));
            if(find(i,jj))//寻找最大匹配
            {
                sum++;
            }
        }
        cout<<sum<<endl;
    }
    return 0;
}