zoj1654

最新推荐文章于 2022-03-31 11:30:30 发布

原创最新推荐文章于 2022-03-31 11:30:30 发布 · 743 阅读

0 ·

CC 4.0 BY-SA版权

每天一题同时被 2 个专栏收录

69 篇文章

订阅专栏

acm

13 篇文章

订阅专栏

本文介绍了一种基于最大二分匹配算法的应用实例——解决皇家卫士问题。通过将地图上的横向线段与纵向线段转化为集合A和B，并利用匈牙利算法找到A集合与B集合间的最大匹配数，从而求解出可以放置的最大数量的皇家卫士。

参考《算法艺术与信息学竞赛》P334 皇家卫士

把每个横向线段当作集合A 每个纵向线段当作集合B

如果是空地则A，B中的点有一条边，最大二分匹配

如o***

就对应了一条横向线段A1，和4条纵向线段B1，B2，B3，B4

其中A1和B1的交点是空地，所以他们之间有一条边，所以每个边都会对应一个点，这个点为一个横向线段和纵向线段的交点

#include<iostream>
#include<stdio.h>
#include<string.h>
#define MAX 100
using namespace std;
int n,p;
int map[MAX*26][MAX*26];
bool visit[MAX*26];
int match[MAX*26];
bool find(int k)
{
    int i,j;
    for(i=1;i<=n;i++)
    {
        if(map[k][i]&&visit[i]==false)
        {
            visit[i]=true;
            if(match[i]==0||find(match[i]))
            {
                match[i]=k;
                return true;
            }
        }
    }
    return false;
}
int hungary() 
{
    memset(match,0,sizeof(match));
    int sum=0;
    for(int i=1;i<=p;i++)
    {
        memset(visit,0,sizeof(visit));
        if(find(i))
          sum++;
    }
    return sum;
}
struct node
{
    int x;
    int y;
}point[MAX][MAX];
int main()
{
    int i,j,k,t;
    int r,c;
    char str[MAX][MAX]; 
    scanf("%d",&t);
    for(k=1;k<=t;k++)
    {
        memset(map,0,sizeof(map));
        n=0;
        p=0;
        scanf("%d%d",&r,&c);
        for(i=0;i<r;i++)
        {
            scanf("%s",&str[i]);
        }
        for(i=0;i<r;i++)
        {
            for(j=0;j<c;j++)
            {
                if(str[i][j]=='#')continue;
                if(j!=0)
                {
                    if(str[i][j-1]=='#')
                      point[i][j].x=++n;
                    else
                      point[i][j].x=point[i][j-1].x;
                }
                else
                {
                    point[i][j].x=++n;              
                }
                if(i!=0)
                {
                    if(str[i-1][j]=='#')
                      point[i][j].y=++p;
                    else
                      point[i][j].y=point[i-1][j].y;
                }
                else
                {
                    point[i][j].y=++p;
                }
                if(str[i][j]=='o')
                {
                    map[point[i][j].y][point[i][j].x]=1;
                }

            }
        }
        printf("Case :%d\n",k);
        printf("%d\n",hungary());
    }
}