HDU 1281 棋盘游戏

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1281

对于这道题，很容易使人想起1045那道题，但区别在于，1045那道题数据量比较小，可以暴力枚举所有方案，复杂度O(2^(n*m))。而这里的n*m高达10000，显然不能暴力出所有方案，那么我们就要寻求更高效的算法。

对于这道题，首先，我们要算出最多能够放置多少个“车”。

假设车放在(x,y)点，很显然，第x行的其它列，以及第y列的其它行都不能放置车。

先考虑行不动，比如放在x0行，设x0能放置小车的列有y0,y1...yn，比如放在y0列，那么可以说x0匹配到了y0。放在y1列，那么可以说x0匹配到了y1。

X集合代表行的集合，Y集合代表列的集合。每个可以放小车的点(x,y)就表示x和y联接的一条可行边。

很显然，最多放置多少个小车这个问题，就是选尽量多的可放点，也就是找尽量多的x和y的匹配边，就是求X和Y集合的最大匹配。

比如对于第一个样例

3 3 4
1 2
1 3
2 1
2 2

那么它建完图对应的二分图如下图所示

可求出其中的一个最大匹配(最大匹配方案可能不止一种)，例如这里选的是(1,2)和(2,1)点放车。

求出最多放几个小车后，那来算重要点就比较简单了。

枚举删除每个可以放小车的点，每次再计算出最大匹配，如果此时的最大匹配数不等于先前那个，那说明这个点就是重要点，因为没有它，最大匹配数会发生变化。

二分图算法为O(n^2)再加上枚举的复杂度，所以最后的复杂度为O(k*n^2)，显然可以过掉题目要求的复杂度。

至于建图，其实只要用邻接矩阵就行了。

#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;
const int maxn=105;
vector<int>G[maxn*maxn];
struct node{int x,y;}p[maxn*maxn];
int girl[maxn*maxn],vis[maxn*maxn],Map[maxn][maxn],link[maxn][maxn];
int n,m,k,cnt1,cnt2,row[maxn][maxn],col[maxn][maxn];
bool Find(int x)
{
    for(int i=1;i<=m;i++)
    {
        if(Map[x][i]&&!vis[i])
        {
            vis[i]=1;
            if(girl[i]==0||Find(girl[i]))
            {
                girl[i]=x;
                return true;
            }
        }
    }
    return false;
}
int match()//匈牙利算法求二分图最大匹配
{
    memset(girl,0,sizeof(girl));
    int ans=0;
    for(int i=1;i<=n;i++)
    {
        memset(vis,0,sizeof(vis));
        if(Find(i)) ans++;
    }
    return ans;
}
int main()
{
    int Case=1;
    while(~scanf("%d%d%d",&n,&m,&k))
    {
        memset(Map,0,sizeof(Map));//初始化邻接矩阵
        for(int i=0;i<k;i++)
        {
            scanf("%d%d",&p[i].x,&p[i].y);
            Map[p[i].x][p[i].y]=1;//可放点置为1
        }
        int Max=match(),ans=0;//先求出最大匹配
        for(int i=0;i<k;i++)//枚举删点
        {
            Map[p[i].x][p[i].y]=0;//删点
            int t=match();//求删点后的最大匹配
            if(t<Max) ans++;
            Map[p[i].x][p[i].y]=1;//恢复
        }
        printf("Board %d have %d important blanks for %d chessmen.\n",Case++,ans,Max);
    }
    return 0;
}