首页 > 代码库 > poj Muddy Fields

poj Muddy Fields

                            Muddy Fields

 

原题去我创的专题里找,在目录首页。

题目:

   给出N*M矩阵,其中*表示泥土,.表示小草。要你用最少的木板把泥土覆盖。木板长度不限,但是只能水平和竖直。

   行列式二分匹配配。是poj3041建图的加强版。

算法:

  由poj3041的建图可以得到灵感,即把木板当作顶点,X轴,Y轴当作边。求解最小顶点覆盖。

3041的激光束不存在障碍物的问题,而这题却有不能压坏小草的限制。

   我们可以换个角度思考,既然不能压坏小草,那就相当于该激光束到达这里必须停止了。所以,我们就以小草为界分块的进行建图分析。每当遇到了小草就到做是行或者是列到头了。这样一来就是普通的二分图了。

 

T_T这一题的数据是25*25是上限 。因为,最坏情况下是草和泥土相间出现。

 

还有就是以后做题的时候,特别是矩阵输入一定要注意输入行列问题!!!!!!!!!

说多了都是泪啊!!!!

 

 

/*
    把对行列覆盖的木板当作点,则每个格子变成了边,
    则问题转换成了最小顶点覆盖问题

*/

#include <iostream>
#include <algorithm>
#include <vector>
#include <cstdio>
#include <cstring>
using namespace std;

const int MAXN = 50 + 10;
const int MAXV = 1500;
vector<int> G[MAXV];

char mp[MAXN][MAXN];
int numx[MAXN][MAXN],numy[MAXN][MAXN];
int match[MAXV];
bool used[MAXV];
int N,M,V;

void init(){
    for(int i = 0;i <= N*M;++i)
        G[i].clear();
}

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

void solve(){
    int res = 0;
    memset(match,-1,sizeof(match));
    for(int i = 1;i <= V;++i){
        memset(used,0,sizeof(used));
        if(dfs(i)) res++;
    }

    printf("%d\n",res);
}

int main()
{
   // freopen("Input.txt","r",stdin);

    while(~scanf("%d%d",&N,&M)){
         init();

         for(int i = 0;i < N;++i){
            scanf("%s",mp[i]);
         }

         int cnt = 0;
         memset(numx,0,sizeof(numx));
         memset(numy,0,sizeof(numy));

         for(int i = 0;i < N;++i)
            for(int j = 0;j < M;++j)
                 if(mp[i][j] == '*'){
                     ++cnt;
                     while(j < M&&mp[i][j] == '*'){
                        numx[i][j] = cnt;
                        ++j;
                     }
                 }

         V = cnt;
         for(int j = 0;j < M;++j)
            for(int i = 0;i < N;++i)
                if(mp[i][j] == '*'){
                    ++cnt;
                    while(i < N&&mp[i][j] == '*'){
                        numy[i][j] = cnt;
                        ++i;
                    }
                }

         for(int i = 0;i < N;++i){
            for(int j = 0;j < M;++j){
                if(numx[i][j] && numy[i][j]){
                    G[numx[i][j]].push_back(numy[i][j]);
                }
            }
         }

         solve();

    }
    return 0;
}


 

poj Muddy Fields