[COGS2175][SDOI2012]象棋-二分图最大权匹配

象棋

【题目描述】

小云和小南两姐妹从小喜欢下象棋，现在作为象棋高手的她们，已经不屑于玩平常的象棋了，于是她们便开始用棋盘和棋子玩各种各样的新游戏。
今天天气晴朗，阳光明媚，她们将在n×m的棋盘上进行游戏。
棋盘上有k颗棋子和若干有障碍格子，令棋盘左上角格子坐标为(1,1)，右下角格子坐标为(n,m)，参数a、b规定了所有棋子的走法：在(x,y)的棋子下一步能走到
(x+a,y+b),(x+a,y−b),(x–a,y+b),(x–a,y–b),(x+b,y+a),(x+b,y−a),(x–b,y+a),(x–b,y–a)这八个格子中的一个，棋子任何时候不能跃出棋盘或走到有障碍的格子上。
这k颗棋子是相同的，小云和小南的目标是用最少步数把所有棋子移动到特定格子，要求移动过程中不能出现多颗棋子同时在某一格的情况。
她们已经想出步数较少方案，但无法确定这是否为最少步数，所以向作为程序员的你求助。

【输入格式】

第一行五个空格隔开的整数n、m、k、a以及b；
接下来n行，每行为长度m的字符串，描述棋盘，‘.’表示没有障碍的格子，‘*’表示有障碍的格子；
接下来k行，每行两个整数x和y，分别表示k颗棋子的初始位置；
接下来k行，每行两个整数x和y，分别表示k颗棋子的目标位置。

【输出格式】

一个整数，为把所有棋子移动到’t’位置的最少步数，数据保证有解。

【样例输入】

1 8 2 2 0
…….*
1 1
1 3
1 5
1 7

【样例输出】

4

【提示】

样例说明
一可行方案如下：第二颗棋子向右跳两步，随后第一颗棋子向右跳两步，共4步。值得注意的是，第一颗棋子向右跳三步，随后第二颗棋子向右跳一步的方案尽管能把棋子都移动到目标位置，但途中两颗棋子曾经同时在(1,3)，违反了规则，所以不能选用此方案。
数据范围
其中20%的数据，n×m≤20；
另外10%的数据，n=1；
对于100%的数据，n、m≤100，k≤500。

---

学习了一发KM算法……
意外的比费用流好写？

---

思路:
可以发现给出的条件中的“不能出现多颗棋子同时在某一格的情况”毫无意义。
因为出现此情况的方案都可以做到不修改路径而避免相交发生。

那么剩下的就简单了~
首先，用BFS预处理出每个起节点到每个终点所需的步数。
然后，可以发现这就是个二分图最大权匹配~
裸费用流可得到90pts，采用神奇的KM算法可得到满分。

至于KM算法，就是在匈牙利的基础上加了一个可行顶标，来不断调整进行匈牙利算法的子图的大小，从而一步步得到最优解。
具体自行百度吧~

这里给出$O(n^3)$的KM算法版本:

#include<bits/stdc++.h>
using namespace std;

inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0' || '9'<ch)ch=getchar();
    while('0'<=ch && ch<='9')x=x*10+(ch^48),ch=getchar();
    return x;
}

inline bool chkmin(int &a,int b){if(a>b){a=b;return 1;}return 0;}
inline bool chkmax(int &a,int b){if(a<b){a=b;return 1;}return 0;}

typedef pair<int,int> pr;
const int N=109;
const int K=509;
const int M=K*K;

int n,m,k,a,b,Inf;
int sx[K],sy[K],tx[K],ty[K],dis[N][N];
int dx[]={1,1,-1,-1,1,1,-1,-1};
int dy[]={1,-1,1,-1,1,-1,1,-1};
int f[K][K],lx[K],ly[K],rec[K],match[K];
bool visx[K],visy[K];
char g[N][N];
pr q[N*N];

inline bool in(int x,int y)
{
    return 1<=x && x<=n && 1<=y && y<=m;
}

inline void bfs(int stx,int sty,int id)
{
    memset(dis,127,sizeof(dis));
    Inf=dis[stx][sty];
    dis[stx][sty]=0;
    q[1]=pr(stx,sty);

    for(int l=1,r=1;l<=r;l++)
    {
        int cx=q[l].first,cy=q[l].second;
        for(int j=0;j<=7;j++)
        {
            int tx=cx+dx[j],ty=cy+dy[j];
            if(in(tx,ty) && g[tx][ty]=='.' && dis[tx][ty]==Inf)
            {
                dis[tx][ty]=dis[cx][cy]+1;
                q[++r]=pr(tx,ty);
            }
        }
    }

    lx[id]=-1e9;
    for(int i=1;i<=k;i++)
    {
        f[id][i]=-dis[tx[i]][ty[i]];
        chkmax(lx[id],f[id][i]);
    }
}

inline bool find(int x)
{
    if(visx[x])return 0;
    visx[x]=1;
    for(int i=1;i<=k;i++)
        if(!visy[i] && f[x][i]==lx[x]+ly[i])
        {
            visy[i]=1;
            if(!match[i] || find(match[i]))
            {
                match[i]=x;
                return 1;
            }
        }
        else if(!visy[i])
            chkmin(rec[i],lx[x]+ly[i]-f[x][i]);
    return 0;
}

int main()
{
    freopen("chessc.in","r",stdin);
    freopen("chessc.out","w",stdout);

    n=read();m=read();
    k=read();a=read();b=read();
    for(int i=0;i<=3;i++)
        dx[i]*=a,dy[i]*=b;
    for(int i=4;i<=7;i++)
        dx[i]*=b,dy[i]*=a;
    for(int i=1;i<=n;i++)
        scanf("%s",g[i]+1);
    for(int i=1;i<=k;i++)
        sx[i]=read(),sy[i]=read();
    for(int i=1;i<=k;i++)
        tx[i]=read(),ty[i]=read();
    for(int i=1;i<=k;i++)
        bfs(sx[i],sy[i],i);

    for(int v=1;v<=k;v++)
    {
        memset(visx,0,sizeof(bool)*(k+2));
        memset(visy,0,sizeof(bool)*(k+2));
        memset(rec,127,sizeof(int)*(k+2));
        while(!find(v))
        {
            int d=1e9;
            for(int i=1;i<=k;i++)
                chkmin(d,rec[i]);
            for(int i=1;i<=k;i++)
            {
                if(visx[i])visx[i]=0,lx[i]-=d;
                if(visy[i])visy[i]=0,ly[i]+=d;
            }
        }
    }

    int ans=0;
    for(int i=1;i<=k;i++)
        ans-=f[match[i]][i];
    printf("%d\n",ans);
    return 0;
}