POJ 2226 Muddy Fields 最小点覆盖+加建图（好题）

题目链接

---

题目一看就是最小点覆盖，这道题与POJ - 3041 算是一类题，但是3041算是一道十分裸的，因为删除的是整行或者整列，所以图其实是现成的，但是本题的难点就在如何建图。

思路：首先还是尽量去建立行与列之间的关系去构造出一条边，因为我们不能将 ※ 以外的点给连接起来，所以二分图的左右的含义不再是简单的行和列，我们需要将行分块以及列分块，就比如离散的思想一样，将每一行分成几块，列也是相同。

建图的代码：

memset(x,0,sizeof(x));///第i行的x[i]个点
memset(y,0,sizeof(y));///第j列的y[i]个点
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
            {
                if(mp[i][j]=='*')
                {
                    if(j==0||mp[i][j-1]=='.') x[i]++;///第i行的x[i]个点
                    if(i==0||mp[i-1][j]=='.') y[j]++;///地j列的y[i]个点
                    add(i*m+x[i],j*n+y[j]);///在i行与第j列建一条边,加到邻接表里
                }
            }

或者我们可以这样就建图，看起来还是比较直观的，横着遍历一遍用r[i][j]的行的号，如果一个点是※，将他右边能延伸的点都与赋成相同的值rx，如果碰到换行或者点的时候rx++；同样对列这样来一遍。

int rx=0,cx=0,c[105][105],r[105][105];//rx行的离散号，cx列的离散号，c,r数组来装
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
            {
                if(mp[i][j]=='*')
                {
                   rx++;
                   while(j<m&&mp[i][j]=='*')
                       r[i][j++]=rx;
                   j--;
                }
            }
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
            {
                if(mp[i][j]=='*')
                {
                    cx++;
                    while(j<m&&mp[i][j]=='*')
                        c[i++][j]=cx;
                    i--;
                }
            }
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
                if(mp[i][j]=='*')
                     if(mp[i][j]=='*') add(r[i][j],c[i][j]);

这样的话，加进去的边我们确定的是如果有两个※号相连，这两条边一定有一个共同的覆盖点，如果是行相邻，就是rx相等，如果列相邻就是cx相等，反正我们就给转化为正常的二分图了。

现在 最小点覆盖=最大匹配。

匈牙利算法就求的最小数量了。

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<stack>
#include<vector>
#include<map>
#include<queue>
#define myself i,l,r
#define lson i<<1
#define rson i<<1|1
#define Lson i<<1,l,mid
#define Rson i<<1|1,mid+1,r
#define half (l+r)/2
#define inff 0x3f3f3f3f
#define lowbit(x) x&(-x)
#define PI 3.14159265358979323846
#define min4(a,b,c,d) min(min(a,b),min(c,d))
#define min3(x,y,z) min(min(x,y),min(y,z))
#define pii make_pair
#define pr pair<int,int>
const int dir[4][2]= {0,-1,-1,0,0,1,1,0};
typedef long long ll;
const ll inFF=9223372036854775807;
typedef unsigned long long ull;
using namespace std;
const int maxn=1e5+5;
char a[105][105];
int used[maxn],girl[maxn];
int n,m;
int x[maxn],y[maxn];
int head[maxn],sign;
struct node
{
    int to,p;
}edge[maxn<<1];
void add(int u,int v)
{
    edge[sign]=node{v,head[u]};
    head[u]=sign++;
}
void init()
{
    memset(head,-1,sizeof(head));
    sign=0;
}
bool find(int u)
{
    for(int i=head[u];~i;i=edge[i].p)
    {
        int v=edge[i].to;
        if(!used[v])
        {
            used[v]=1;
            if(girl[v]==-1||find(girl[v]))
            {
                girl[v]=u;
                return true;
            }
        }
    }
    return false;
}
int hungry()
{
    memset(girl,-1,sizeof(girl));
    int ans=0;
    for(int i=0;i<=n*m;i++)
    {
        memset(used,0,sizeof(used));
        if(find(i)) ans++;
    }
    return ans;
}
int main()
{
    while(cin>>n>>m)
    {
        init();
        for(int i=0;i<n;i++)
            scanf("%s",a[i]);
        memset(x,0,sizeof(x));
        memset(y,0,sizeof(y));
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
            {
                if(a[i][j]=='*')
                {
                    if(j==0||a[i][j-1]=='.') x[i]++;///第i行的x[i]个点
                    if(i==0||a[i-1][j]=='.') y[j]++;///地j列的y[i]个点
                    add(i*m+x[i],j*n+y[j]);///在i行与第j列建一条边
                }
            }
        printf("%d\n",hungry());
    }
}