bzoj2150: 部落战争(匈牙利)

2150: 部落战争

题目：传送门 

---

题解：

　　 辣鸡数据..毁我AC率

　　 先说做法，很容易就可以看出是二分图匹配的最小路径覆盖（可能是之前不久刚做过类似的题）

　　 一开始还傻逼逼的去直接连边然后准备跑floyd...肯定是做祭祀做傻了

　　 二分图嘛，将每个点拆成两个集合再连啊...

　　 然后最小路径覆盖=总点数-最大匹配数

　　 强烈吐槽！个人习惯用数组存点编号...50*50的数据范围我开55*55...结果WA n次

　　 然后...改成56*56...AC。。。ORZ

 　　 

---

代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
struct node
{
    int x,y,next;
}a[25000];int len,last[25000];
void ins(int x,int y)
{
    len++;a[len].x=x;a[len].y=y;
    a[len].next=last[x];last[x]=len;
}
char s[55][55];
int n,m,R,C,t;
int match[25000],chw[25000];
bool findmuniu(int x)
{
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(chw[y]!=t)
        {
            chw[y]=t;
            if(match[y]==0 || findmuniu(match[y]))
            {
                match[y]=x;
                return true;
            }
        }
    }
    return false;
}
int d[56][56];//辣鸡数据
int main()
{
    scanf("%d%d%d%d",&n,&m,&R,&C);
    for(int i=1;i<=n;i++)scanf("%s",s[i]+1);
    int ss=0;memset(d,-1,sizeof(d));
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)d[i][j]=++ss;
    const int dx[5]={0,R,R,C,C};
    const int dy[5]={0,-C,C,-R,R};
    int sum=0;len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(s[i][j]=='.')
            {
                sum++;
                for(int k=1;k<=4;k++)
                {
                    int x=i+dx[k],y=j+dy[k];
                    if(d[x][y]!=-1 && s[x][y]!='x')
                        ins(d[i][j],d[x][y]+n*m);
                }
            }
    memset(match,0,sizeof(match));
    memset(chw,0,sizeof(chw));
    int ans=0;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(s[i][j]=='.')
            {
                t++;
                if(findmuniu(d[i][j]))ans++;
            }
    printf("%d
",sum-ans);
    return 0;
}