[BZOJ3205][APIO2013]Robot(斯坦纳树)

3205: [Apio2013]机器人

Time Limit: 15 Sec  Memory Limit: 128 MB
Submit: 1007  Solved: 240
[Submit][Status][Discuss]

Description

VRI（Voltron 机器人学会）的工程师建造了 n个机器人。任意两个兼容的机器人站在同一个格子时可以合并为一个复合机器人。我们把机器人用 1至 n编号（n ≤ 9）。如果两个机器人的编号是连续的，那么它们是兼容的，可以合并成一个复合机器人。最初这   n   个机器人各自都只有唯一的编号。而一个由两个或以上的机器人合并构成的复合机器人拥有两个编号，分别是构成它的所有机器人中最小和最大的编号。例如， 2号机器人只可以与 1号或 3号机器人合并。若 2号机器人与 3号机器人合并，可构成编号为 2-3的复合机器人。如果编号为 2-3的复合机器人与编号为 4-6的复合机器人合并，可构成编号为 2-6的复合机器人。当所有机器人合并以后则构成 1-n复合机器人。工程师把这 n个机器人放在了一个封闭的房间中，房间四周均是墙。该房间被划分成 w     h    个方格。有些方格有障碍物，机器人不可经过或停留；其余方格允许多个机器人停留，同时允许机器人经过。任何时候一个机器人只占用一个方格。初始时刻，所有 机器人均在不同的方格中。这些原始的机器人不会自发地移动。它们只有被工程师沿   x轴或 y轴推动后，才会沿推动的方向不断向前直线移动，直至碰到障碍物或墙停止移动。停止移动后，它会扫描当前的格子是否存在可以与它合并的机器人，如果有，则 合并并继续检查，直至不能再合并为止。工程师只能沿水平向左、水平向右、竖直向上、竖直向下四个方向推动机器人，并且，在机器人尚未停止移动时，不允许推 动其它机器人，因此任何时刻，房间中都只能有一个机器人移动,为了帮助机器人转向，工程师在一些格子中放置了转向器。具体地说，转向器分为顺时针转向器 （右转器）和逆时针转向器（左转器），顺时针转向器可以使到达该格子的机器人沿顺时针方向转向   90_；逆时针转向器可以使到达该格子的机器人沿逆时针方向转向 90_。现在，我们将告诉你初始时刻房间内的信息。请你计算工程师最少共计需要推动机器人多少次，才能把所有的 n个机器人全部合并（如果可能的话）。

Input

你的程序必须从标准输入读入。输入的第 1行包含 3个整数 n、w和 h，用空格隔开。输入文件中接下来的 h行描述初始时刻房间内的信息，每行包含w个字符。这w* h 字符中每一个表示房间中的一个格子，意义如下：
 
‘ 1’至‘9’：表示该方格中有一个机器人，编号为这个数字；
‘ x’：表示该方格有障碍物；
 
‘ A’：表示该方格中有一个逆时针转向器；
 
‘ C’：表示该方格中有一个顺时针转向器；
‘ .’：表示该方格为空地。

Output

你的程序必须输出到标准输出。输出仅一个整数，表示最少需要推动的次数。
若不能使所有机器人全部合并，输出-1。

Sample Input

4 10 5
1.........
AA...x4...
..A..x....
2....x....
..C.3.A...

Sample Output

5

HINT

第一步：向右推动 3 号机器人，当它碰到转向器后会向上继续移动，直至碰到墙壁停止移动。第二步：向上推动 4 号机器人，当它碰到墙壁后停止移动，与3 号机器人合并，构成  3-4 号机器人 第三步：向上推动 2 号机器人，当它碰到转向器后会向左移动，由于左侧为墙壁，故停留在原地。第四步：向右推动  2 号机器人，由于它在一个转向器上，故它会向上移动，直至碰到墙壁停止移动，与  1 号机器人合并，构成 1-2 号机器人。第五步：向左推动  3-4 号机器人，当它碰到墙壁后停止移动，与 1-2 号机器人合并，构成  1-4 号机器人。

n ≤ 9，w ≤ 500且   h ≤ 500

Source

[Submit][Status][Discuss]

都说是比较明显的斯坦纳树，但我并没有看出来，只是状态转移比较像而已，和连通性并没有什么关系（或许是合并的过程比较像最小生成树）。

http://wyfcyx.is-programmer.com/posts/75949.html

不好写，卡空间又卡时间。注意能用栈的不要用队列，栈的常数小不少。

不知道为什么按照lych的代码写的，改得几乎一样了，常数还要大一半。

还有那些用std::queue过的真心服气。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
typedef long long ll;
using namespace std;

const int N=505,M=N*N+5,inf=0x3f3f3f3f;
const int dx[4]={-1,0,1,0},dy[4]={0,1,0,-1};
char ch;
struct P{ int x,y; }go[N][N][4],q1[M],q2[M],pos[N];
int n,m,s,ans,tim,a[N][N],f[10][10][N][N],inq[N][N],vis[N][N][4],c[M],val[M];

int add(int x){ x++; return (x>=M) ? x-M : x; }
void inc(int &x){ x=add(x); }
bool jud(int x,int y){ return x>=1 && x<=n && y>=1 && y<=m && a[x][y]!=12; }

P dfs(int x,int y,int d){
    if (vis[x][y][d]==tim) return (P){-1,-1};
    if (go[x][y][d].x || go[x][y][d].y) return go[x][y][d];
    vis[x][y][d]=tim; int t=d;
    if (a[x][y]==10) t=(t+3)%4;
    if (a[x][y]==11) t=(t+1)%4;
    int xx=x+dx[t],yy=y+dy[t];
    if (!jud(xx,yy)) return go[x][y][d]=(P){x,y};
    return go[x][y][d]=dfs(xx,yy,t);
}

void spfa(int l,int r){
    memset(inq,0,sizeof(inq)); P ss;
    int top=0,st=0,ed=0,mn=inf,mx=-inf;
    rep(x,1,n) rep(y,1,m) if (f[l][r][x][y]!=inf)
        q1[++top]=(P){x,y},val[top]=f[l][r][x][y],mn=min(mn,val[top]),mx=max(mx,val[top]);
    
    rep(i,0,mx) c[i]=0;
    rep(i,1,top) c[val[i]]++;
    rep(i,mn+1,mx) c[i]+=c[i-1];
    for (int i=top; i; i--) q2[c[val[i]]--]=q1[i];
    rep(i,1,top) q1[i]=q2[top-i+1];
    
    while (top || st!=ed){
        int now=add(st);
        if (st==ed || (top && f[l][r][q1[top].x][q1[top].y]<f[l][r][q2[now].x][q2[now].y])) ss=q1[top--]; else ss=q2[st=now];
        int x=ss.x,y=ss.y; inq[x][y]=0;
        rep(d,0,3){
            P pp=go[x][y][d];
            if (pp.x==x && pp.y==y) continue;
            if (~pp.x && f[l][r][pp.x][pp.y]>f[l][r][x][y]+1){
                f[l][r][pp.x][pp.y]=f[l][r][x][y]+1;
                if (!inq[pp.x][pp.y]) inq[pp.x][pp.y]=1,inc(ed),q2[ed]=pp;
            }
        }
    }
}

int main(){
    freopen("bzoj3205.in","r",stdin);
    freopen("bzoj3205.out","w",stdout);
    scanf("%d%d%d",&s,&m,&n); ans=inf;
    memset(f,0x3f,sizeof(f));
    rep(i,1,n) rep(j,1,m){
        scanf(" %c",&ch);
        if (ch>='1' && ch<='9'){
            int t=ch-'0';
            a[i][j]=t; pos[ch-'0']=(P){i,j}; f[t][t][i][j]=0;
        }
        if (ch=='A') a[i][j]=10;
        if (ch=='C') a[i][j]=11;
        if (ch=='x') a[i][j]=12;
    }
    rep(i,1,n) rep(j,1,m) if (a[i][j]!=12) rep(k,0,3) if (a[i][j]!=12) tim++,go[i][j][k]=dfs(i,j,k);
    rep(i,1,s) f[i][i][pos[i].x][pos[i].y]=0;
    rep(l,1,s) rep(i,1,s-l+1){
        int j=i+l-1;
        rep(k,i,j-1) rep(x,1,n) rep(y,1,m)
            f[i][j][x][y]=min(f[i][j][x][y],f[i][k][x][y]+f[k+1][j][x][y]);
        spfa(i,j);
    }
    rep(x,1,n) rep(y,1,m) ans=min(ans,f[1][s][x][y]);
    printf("%d
",(ans<inf)?ans:-1);
    return 0;
}