pku3020 Antenna Placement

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

二分图匹配，匈牙利算法(dfs+邻接矩阵)

建图：相邻点双向连接，建无向图

结果（最少需要的“圈”数） ==（总点数 - 最大匹配数）/*没有匹配成功的点需要自己单独一个“圈”*/  +  (最大匹配数 / 2) /*匹配成功的点，2个点共用一个"圈"*/

化简为：结果 == 总点数 - （最大匹配数 / 2）

#include <stdio.h>
#include <string.h>
#define N 432

int n, m;
int map[N][N], flag[N], girl[N];

int find(int x)
{
    int i;
    for(i=1; i<=n*m; i++)
    {
        if(!flag[i] && map[x][i])
        {
            flag[i] = 1;
            if(girl[i]==-1 || find(girl[i]))
            {
                girl[i] = x;
                return 1;
            }
        }
    }
    return 0;
}

int main()
{
    int t, i, j, k, p, q, temp, sum, count;
    int dir[4][2] = {1,0,0,1,-1,0,0,-1},map1[43][12];
    char c;
    scanf("%d", &t);
    while(t-- && scanf("%d%d%*c", &n, &m))
    {
        count = 0;
        for(i=1; i<=n; i++)
        {
            for(j=1; j<=m; j++)
            {
                scanf("%c", &c);
                map1[i][j] = c-'o'? (++count,1): 0;
            }
            getchar();
        }
        temp = n*m;
        for(i=1; i<=temp; i++)
        {
            girl[i] = -1;
            for(j=1; j<=temp; j++)
            {
                map[i][j] = 0;
            }
        }
        for(i=1; i<=n; i++)
        {
            for(j=1; j<=m; j++)
            {
                if(map1[i][j])
                {
                    temp = m*(i-1) + j;
                    //printf("%d: ", temp);
                    for(k=0; k<4; k++)
                    {
                        p = i + dir[k][0];
                        q = j + dir[k][1];
                        if(p<1 || p>n || q<1 || q>m)
                        {
                            continue;
                        }
                        if(map1[p][q])
                        {
                            //printf("(%d, %d) ", p, q);
                            map[temp][m*(p-1)+q] = 1;
                        }
                    }
                    //printf("\n");
                }
            }
        }
        temp = n*m;
        sum = 0;
        for(i=1; i<=temp; i++)
        {
            for(j=1; j<=temp; j++)
            {
                //printf("%d ", map[i][j]);
                flag[j] = 0;
            }
            //printf("\n");
            sum += find(i);
        }
        printf("%d\n", count-sum/2);
    }
    return 0;
}