4185 Oil Skimming 最大匹配 奇偶建图

题目大意：

　　统计相邻（上下左右）的‘#’的对数。

解法：

　　与题目hdu1507 Uncle Tom's Inherited Land*类似，需要用奇偶建图。就是行+列为奇数的作为X集合，偶尔作为Y集合，都是‘#’就连边。最后求最大匹配。

　　数据有点大，直接建图会出错（我试过）。可以按照‘#’出现的顺序给顶点编号。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
#include<algorithm>
using namespace std;
const int N=36005,INF=0x3f3f3f3f;
int cx[N],cy[N],dx[N],dy[N];
bool bmask[N];
vector<int > bmap[N];
int nx,ny,dis,ans;
bool searchpath()
{
    queue<int> q;
    dis=INF;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=1;i<=nx;i++)
    {
        if(cx[i]==-1){ q.push(i); dx[i]=0; }
        while(!q.empty())
        {
            int u=q.front(); q.pop();
            if(dx[u]>dis) break;
            for(int i=0;i<bmap[u].size();i++)
            {
                int v=bmap[u][i];
                if(dy[v]==-1)
                {
                    dy[v]= dx[u] + 1;
                    if(cy[v]==-1) dis=dy[v];
                    else
                    {
                        dx[cy[v]]= dy[v]+1;
                        q.push(cy[v]);
                    }
                }
            }
        }
    }
    return dis!=INF;
}
int findpath(int u)
{
    for(int i=0;i<bmap[u].size();i++)
    {
        int v=bmap[u][i];
        if(!bmask[v]&&dy[v]==dx[u]+1)
        {
            bmask[v]=1;
            if(cy[v]!=-1&&dy[v]==dis) continue;
            if(cy[v]==-1||findpath(cy[v]))
            {
                cy[v]=u; cx[u]=v;
                return 1;
            }
        }
    }
    return 0;
}
void maxmatch()
{
    ans=0;
    memset(cx,-1,sizeof(cx));
    memset(cy,-1,sizeof(cy));
    while(searchpath())
    {
        memset(bmask,0,sizeof(bmask));
        for(int i=1;i<=nx;i++)
            if(cx[i]==-1) ans+=findpath(i);
    }
}
void init()
{
    for(int i=0;i<=N;i++) bmap[i].clear();
}
char s[605][605];
int Map[605][605];
int main()
{
    //freopen("test.txt","r",stdin);
    int cas,n,i,j,k=1,a,b,t=1;
    char ch;
    scanf("%d",&cas);
    while(cas--)
    {
        scanf("%d",&n);
        t=0;
        for(i=1;i<=n;i++){
            getchar();
            for(j=1;j<=n;j++)
            {
                scanf("%c",&s[i][j]);
                if(s[i][j]=='#') Map[i][j]=++t;
            }
        }
        init();
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=n;j++)
            {
                if(s[i][j]=='#'&&((i+j)%2))
                {
                    a=Map[i][j];
                    if(i>1&&s[i-1][j]=='#'){
                        b=Map[i-1][j];
                        bmap[a].push_back(b);
                    }
                    if(i<n&&s[i+1][j]=='#'){
                        b=Map[i+1][j];
                        bmap[a].push_back(b);
                    }
                    if(j>1&&s[i][j-1]=='#'){
                        b=Map[i][j-1];
                        bmap[a].push_back(b);
                    }
                    if(j<n&&s[i][j+1]=='#'){
                        b=Map[i][j+1];
                        bmap[a].push_back(b);
                    }
                }
            }
        }
        nx=t;
        maxmatch();
        printf("Case %d: %d
",k++,ans);
    }
    return 0;
}