【博弈+二分图匹配】[NOI2011]兔兔与蛋蛋游戏

题目描述

Description

这里写图片描述

Input

输入的第一行包含两个正整数 n、m。 接下来 n行描述初始棋盘。其中第i 行包含 m个字符,每个字符都是大写英文字母"X"、大写英文字母"O"或点号"."之一,分别表示对应的棋盘格中有黑色棋子、有白色棋子和没有棋子。其中点号"."恰好出现一次。 接下来一行包含一个整数 k(1≤k≤1000) ,表示兔兔和蛋蛋各进行了k次操作。 接下来 2k行描述一局游戏的过程。其中第 2i – 1行是兔兔的第 i 次操作(编号为i的操作) ,第2i行是蛋蛋的第i次操作。每个操作使用两个整数x,y来描述,表示将第x行第y列中的棋子移进空格中。 输入保证整个棋盘中只有一个格子没有棋子, 游戏过程中兔兔和蛋蛋的每个操作都是合法的,且最后蛋蛋获胜。

Output

输出文件的第一行包含一个整数r,表示兔兔犯错误的总次数。 接下来r 行按递增的顺序给出兔兔“犯错误”的操作编号。其中第 i 行包含一个整数ai表示兔兔第i 个犯错误的操作是他在游戏中的第 ai次操作。 1 ≤n≤ 40, 1 ≤m≤ 40

Sample Input

样例一:
1 6
XO.OXO
1
1 2
1 1
样例二:
3 3
XOX
O.O
XOX
4
2 3
1 3
1 2
1 1
2 1
3 1
3 2
3 3
样例三:
4 4
OOXX
OXXO
OO.O
XXXO
2
3 2
2 2
1 2
1 3

Sample Output


样例一:
1
1
样例二:
0
样例三:
2
1
2

样例1对应图一中的游戏过程
样例2对应图三中的游戏过程

HINT



这里写图片描述

Source


分析

75分做法

我们把操作看做是空格在移动。
我们做这道题最先想到的就是,操作是否可能成环?
答案是否定的。我们假设环长为 n ,显然n为偶数,因为环中的操作总是成对出现的,即有向上就要有向下,有向左就有向右。假设操作为 A1,A2,A3An ,显然操作 A1 An 移入空格的是同一个棋子,但是操作的人却不一样,显然是不可能的,所以操作不可能成环。
既然不能成环,而且数据还这么小,似乎搜索很好写。敲一发,75分到手。

满分做法

我们发现,移动路径上,相邻两个点的颜色总是不一样的,这然我们想到了二分图。
我们在相邻而且颜色不相同的两个点之间连边,由于空格能够走到白色格子,我们不妨把空格看做黑色。
先手在当前点能够获胜的条件是从当前点出发,能够找到一条路径长度是奇数而且先手一定能够使路径长度为奇数。
由于是二分图,那我们求个最大匹配试试。
我们发现,匈牙利算法增广时的交错轨不就是一条长度为奇数的路径吗。
那么我们来思考一下两种情况。

  1. 如果当前点不一定最大匹配上,那么它的邻接点一定在最大匹配上。 我们假设它的邻接点不一定在最大匹配上:

    1. 存在邻接点不在时,当前点也不在。那么久可以增广了,和假设是最大匹配不符。
    2. 如果邻接点和当前点一定有一个在,那么邻接点和当前点一定还通过一个有偶数条边的交错轨相连,加上这条边就是一个奇环了,不是二分图。
  2. 如果当前点一定在最大匹配上,则一定获胜
    首先,当前点一定在一条增广路径上,即在一条交错轨上。这条交错轨一定是奇数的,而且被这个点分成了一边是奇数,一边是偶数。但是奇数那一边也有可能有一条偶数的交错轨,即可能有岔路。
    我们来分析一下,岔路可能有两种情况。
    这里写图片描述
    显然,这种情况你自己不去走那边就好了。

    1. 出现在别人的节点
      这里写图片描述
      如果朝奇数条交错轨的方向走,自己的节点到别人的节点这条边应该在最大匹配中(图中蓝色的边),但是这条边显然可以被橙色的边代替,和当前点一定在最大匹配上的假设不符。
      所以,最终一定走的是奇数条边,使后手无法行动。 也可以这么看,先手一定可以有边可走(沿着匹配边走),但是后手不一定有边可走。

由此我们得出结论,如果先手所在的点在一定在最大匹配上,则先手有必胜策略,否则,后手有必胜策略。
每次操作时,我们就检查当前空格所在位置是不是一定在最大匹配上,然后将这个位置删掉,将空格的位置设置为这次操作的位置。
那么我们怎么看当前点是不是一定在最大匹配上呢?

  • 如果当前点本来就没有匹配的点,显然不在。
  • 如果当前点所匹配的点在删去当前点后,能够继续增广,就说明当前点不一定在最大匹配上,否则一定在。

代码

#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
#define MAXN 40
#define MAXM 40
#define MAXK 1000
queue<int>q;
int n,m,cx[MAXN*MAXM+10],dir[4][2]={{1,0},{-1,0},{0,1},{0,-1}},bx,by,cnt,ans[MAXK+10],k;
char s[MAXN+10][MAXN+10];
bool vis[MAXN*MAXM+10],ban[MAXN*MAXM+10];
struct node{
    int v;
    node *next;
}*adj[MAXN*MAXN+10],edge[MAXN*MAXN*4+10],*ecnt=edge;
inline void addedge(int u,int v){
    node *p=++ecnt;
    p->v=v;
    p->next=adj[u];
    adj[u]=p;
}
void Read(int &x){
    static char c;
    while(c=getchar(),c!=EOF)
        if(c>='0'&&c<='9'){
            x=c-'0';
            while(c=getchar(),c>='0'&&c<='9')
                x=x*10+c-'0';
            ungetc(c,stdin);
            return;
        }
}
void read(){
    Read(n),Read(m);
    int i,j;
    for(i=1;i<=n;i++)
        scanf("%s",s[i]+1);
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
            if(s[i][j]=='O')
                s[i][j]=0;
            else if(s[i][j]=='X')
                s[i][j]=1;
            else
                s[i][j]=1,bx=i,by=j;
}
inline int Get_id(int i,int j){
    return (i-1)*m+j;
}
void bfs(){
    q.push(Get_id(bx,by));
    int u,x,y,tx,ty,v,d;
    while(!q.empty()){
        u=q.front();
        q.pop();
        x=(u+m-1)/m;
        y=u-(x-1)*m;
        for(d=0;d<4;d++){
            tx=x+dir[d][0];
            ty=y+dir[d][1];
            if(tx&&ty&&tx<=n&&ty<=m&&s[x][y]^s[tx][ty]){
                v=Get_id(tx,ty);
                addedge(u,v);
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
}
void print(){
    printf("%d\n",cnt);
    int i;
    for(i=1;i<=cnt;i++)
        printf("%d\n",ans[i]);
}
bool dfs(int u){
    for(node *p=adj[u];p;p=p->next){
        if(!vis[p->v]&&!ban[p->v]){
            vis[p->v]=1;
            if(!cx[p->v]||dfs(cx[p->v])){
                cx[p->v]=u;
                cx[u]=p->v;
                return 1;
            }
        }
    }
    return 0;
}
void solve(){
    int i,u,mt;
    for(i=n*m;i;i--){
        if(!cx[i]){
            memset(vis,0,sizeof vis);
            dfs(i);
        }
    }
    bool r1,r2;
    Read(k);
    for(i=1;i<=k;i++){
        u=Get_id(bx,by);
        ban[u]=1;
        if(cx[u]){
            mt=cx[u];
            cx[u]=cx[mt]=0;
            memset(vis,0,sizeof vis);
            r1=!dfs(mt);
        }
        else
            r1=0;
        Read(bx),Read(by);
        u=Get_id(bx,by);
        ban[u]=1;
        if(cx[u]){
            mt=cx[u];
            cx[u]=cx[mt]=0;
            memset(vis,0,sizeof vis);
            r2=!dfs(mt);
        }
        else
            r2=0;
        if(r1&&r2)
            ans[++cnt]=i;
        Read(bx),Read(by);
    }
}
int main()
{
    read();
    bfs();
    solve();
    print();
}
  • 0
    点赞
  • 0
    收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
©️2022 CSDN 皮肤主题:编程工作室 设计师:CSDN官方博客 返回首页
评论

打赏作者

outer_form

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值