zoj 1654

很好的一道题目。
题意：给定一个 M * N 大小的方格，地图有3中方格，墙、草地、空地。他的老板希望Robert能在地图中放置尽可能多的机器人。每个机器人可以向四个方向开火，激光可以穿透草地，但不能穿透墙壁。（http://www.cnblogs.com/g0feng/archive/2012/11/10/2763771.html）
求出最大可以放置的机器人数。
思路：
建图方式：矩阵横竖块建二分图，将问题转化为最大匹配问题。
这题是《图论算法理论、实现及应用》书中的例题，有详细的分析P361

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
// Accepted  1654 C++ 50ms 6840K 
#define MAXN 51
int xx, yy, x[MAXN*MAXN], y[MAXN*MAXN], xs[MAXN][MAXN], ys[MAXN][MAXN];
char Graph[MAXN][MAXN];
int vis[MAXN*MAXN];
bool map[MAXN*MAXN][MAXN*MAXN];
int row, col;
bool dfs( int u )
{
    int v;
    for( v = 1; v <= yy; v++ )
    {
        if( map[u][v] && !vis[v] )
        {
            vis[v] = 1;
            if( !y[v] || dfs( y[v] ) )
            {
                x[u] = v;
                y[v] = u;
                return true;
            }
        }
    }
    return false;
}
void MaxMatch()
{
    int u, ans = 0;
    memset( x, 0, sizeof(x) );
    memset( y, 0, sizeof(y) );
    for( u = 1; u <= xx; u++ )
    {
        if( !x[u] )
        {
            memset( vis, 0, sizeof(vis) );
            if( dfs( u ) )ans++;
        }
    }
    printf( "%d\n", ans );
}
int main()
{
    int cas, i, j, k, tot, flag;
    scanf( "%d", &cas );
    for( k = 1; k <= cas; k++ )
    {
        memset( xs, 0, sizeof(xs) );
        memset( ys, 0, sizeof(ys) );
        scanf( "%d%d", &row, &col );
        for( i = 0; i < row; i++ )
            scanf( "%s", Graph[i] );
        tot = 0;
        for( i = 0; i < row; i++ )
        {
            flag = 0;
            for( j = 0; j < col; j++ )
            {
                if( Graph[i][j] == 'o' )
                {
                    if( !flag )tot++;
                    xs[i][j] = tot;    flag = 1;
                }
                else if( Graph[i][j] == '#' )
                    flag = 0;
            }
        }
        xx = tot;
        tot = 0;
        for( j = 0; j < col; j++ )
        {
            flag = 0;
            for( i = 0; i < row; i++ )
            {
                if( Graph[i][j] == 'o' )
                {
                    if( !flag )tot++;
                    ys[i][j] = tot;    flag = 1;
                }
                else if( Graph[i][j] == '#' )
                    flag = 0;
            }
        }
        yy = tot;
        memset( map, 0, sizeof(map) );
        for( i = 0; i < row; i++ )
            for( j = 0; j < col; j++ )
                if( xs[i][j] )map[xs[i][j]][ys[i][j]] = 1;
        printf( "Case :%d\n", k );
        MaxMatch();
    }
}