[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算法…… 
意外的比费用流好写?


思路: 

  本题的难点是“移动过程中不能出现多颗棋子同时在某一格的情况”。

  事实上,可以忽略此条件,因为棋子是相同的,我们可以用合法的等效方案替代一棋子越过另一棋子的情况:A、B、C三格,A能在一步走到B,B也能在一步走到C。

在A的棋子需要走到存在棋子的B,接着走到C。此情形我们可以看成在B的棋子先走到C,接着在A的棋子走到B。

  BFS预处理出每个初始位置走到每个终止位置的最少步数。

  把初始位置抽象成二部图的左部,终止位置抽象成二部图的右部,左右之间边权为最少步数。

  那么次二部图的完备匹配对应着一种方案,匹配的边权和对应最少总步数。

  可用最佳匹配解决。

#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
",ans);
    return 0;
}

(转载自:https://blog.csdn.net/zlttttt/article/details/79394475

原文地址:https://www.cnblogs.com/zbsy-wwx/p/11680650.html