題目大意
給你一個的方格,中間有障礙物,給你個機器人以及個出口,機器人開始都在的位置,出口在,開始機器人的前進方向爲豎直向下。
要求你通過在方格中放置若干“NE”,“NW”,“SW”,“SE”四種裝置,要求所有的機器人都從出口出來,判斷是否可行。
其中四種裝置的具體操作如下:
"NW’’ 裝置會把從格子上方走來的機器人的運動方向變成向左,以及把從格子左方走來的機器人的運動方向變成向上,不允許機器人從格子的右方及下方進入。
""NE’’ 裝置會把從格子上方走來的機器人的運動方向變成向右,以及把從格子右方走來的機器人的運動方向變成向上,不允許機器人從格子的左方及下方進入。
""SW’’ 裝置會把從格子下方走來的機器人的運動方向變成向左,以及把從格子左方走來的機器人的運動方向變成向下,不允許機器人從格子的右方及上方進入。
""SE’’ 裝置會把從格子下方走來的機器人的運動方向變成向右,以及把從格子右方走來的機器人的運動方向變成向下,不允許機器人從格子的左方及上方進入。
思路和分析
看到數據範圍其實就能夠想到是最大流問題,因爲數據範圍明顯不是狀壓,而且以前也做過類似的將方格中的每個點拆分建圖的問題,當時卡死在A題上面,沒有早點看這個題,不然應該能做出來。
首先來分析一下問題,觀察一下轉彎裝置的性質,由於只能夠從一邊進一邊出完成一個變向的過程,再根據機器人初始方向一致的條件,不難分析出不存在兩個機器人使用同一個裝置的可能性。所以我們把每個點拆成左右走的和上下走的點就可以很好的解決問題,最後跑一邊最大流,判一下就可。
還有一些顯然的性質,爲了儘量不產生交錯,所以4個方向中可以去掉兩個,維護單調性。
#include <bits/stdc++.h>
#define rep( i , l , r ) for( int i = (l) ; i <= (r) ; ++i )
#define per( i , r , l ) for( int i = (r) ; i >= (l) ; --i )
#define erep( i , u ) for( int i = head[(u)] ; ~i ; i = e[i].nxt )
using namespace std;
const int INF = 0x3f3f3f3f;
const int maxn = 40010;
const int maxm = 1000010;
int dep[maxn] , head[maxn] , _t;
struct edge{
int u , v , nxt;
int c;
} e[maxm];
inline void addedge( int u , int v , int c , int op ){
e[_t].u = u , e[_t].v = v , e[_t].nxt = head[u] , e[_t].c = c , head[u] = _t++;
e[_t].u = v , e[_t].v = u , e[_t].nxt = head[v] , e[_t].c = op ? 0 : c , head[v] = _t++;
}
int N , M , S , T;
int cur[maxn];
bool bfs(){ // 給圖分層
memset( dep , 0x3f , sizeof dep );
memcpy( cur , head , sizeof head );
dep[S] = 0;
queue<int> q;
q.push( S );
while( !q.empty() ){
int u = q.front(); q.pop();
erep( i , u ){
int v = e[i].v;
if( dep[v] > dep[u] + 1 && e[i].c > 0 ){
dep[v] = dep[u] + 1;
q.push( v );
if( v == T ) return 1;
}
}
}
return 0;
}
inline int _min( int a , int b ){
return a > b ? b : a;
}
// recursive version
//bool vis = 0;
//int max_flow = 0;
//int dfs( int u , int flow ){
// int cur_low = 0;
// if( u == T ){
// vis = 1;
// max_flow += flow; // 直接累加
// return flow;
// }
// int used = 0; // 當前已經使用的流量
// for( int i = cur[u] ; ~i ; cur[u] = i = e[i].nxt ){
// int v = e[i].v;
// if( e[i].c > 0 && dep[v] == dep[u] + 1 ){
// if( cur_low = dfs( v , _min( flow - used , e[i].c ) ) ){
// used += cur_low;
// e[i].c -= cur_low;
// e[i^1].c += cur_low;
// if( used == flow ) break; // 流量使用完說明這個點無法增廣了
// }
// }
// }
// return used; // 返回已經使用的流量
//}
//int Dinic(){
// while( bfs() ){
// for( vis = 1 ; vis ; ){
// vis = 0;
// dfs( S , INF );
// }
// }
// return max_flow;
//}
// non-recursive version
int Dinic(){
int res = 0;
int top = 0;
int stack[maxn];
while( bfs() ){
int u = S;
top = 0;
for(;;){
if( u == T ){
int mi = INF , loc;
for( int i = 0 ; i < top ; ++i )
if( mi > e[stack[i]].c ){
loc = i;
mi = e[stack[i]].c;
}
for( int i = 0 ; i < top ; ++i ){
e[stack[i]].c -= mi;
e[stack[i]^1].c += mi;
}
res += mi;
top = loc;
u = e[stack[top]].u;
}
for( int i = cur[u] ; ~i ; cur[u] = i = e[i].nxt )
if( e[i].c > 0 && dep[e[i].v] == dep[u] + 1 ){
break;
}
if( cur[u] != -1 ){
stack[top++] = cur[u];
u = e[cur[u]].v;
}else{
if( top == 0 ) break;
dep[u] = INF;
u = e[stack[--top]].u;
}
}
}
return res;
}
const int MAXN = 111;
int g[MAXN][MAXN];
int m;
inline int id( int x , int y ){
return x * m + y;
}
char s[MAXN];
int main(){
int cas = 0;
scanf("%d" , &cas);
while( cas-- ){
memset( head , 0xff , sizeof head );
_t = 0;
int n , a , b;
scanf("%d%d%d%d" , &n , &m , &a , &b);
char ch;
rep( i , 1 , n ){
scanf("%s" , s);
rep( j , 1 , m ){
g[i][j] = s[j - 1] - '0';
}
}
int x , base = (n + 2) * m;
S = base * 2 + 5 ;
T = S + 1;
rep( i , 1 , a ){
scanf("%d" , &x);
if( g[1][x] != 1 )
addedge( S , id(1,x) , 1 , 1 );
}
rep( i , 1 , b ){
scanf("%d" , &x);
if( g[n][x] != 1 )
addedge( id(n,x) , T , 1 , 1 );
}
rep( i , 1 , n )
rep( j , 1 , m ){
if( g[i][j] == 1 ) continue;
if( i > 1 && g[i - 1][j] == 0 ) addedge( id(i,j) , id(i-1,j) , 1 , 0 );
//if( i < n && g[i + 1][j] == 0 ) addedge( id(i,j) , id(i+1,j) , 1 );
if( j > 1 && g[i][j - 1] == 0 ) addedge( id(i,j) + base , id(i,j-1) + base , 1 , 0 );
//if( j < m && g[i][j + 1] == 0 ) addedge( id(i,j) + base , id(i,j+1) + base , 1 );
addedge( id(i,j) , id(i,j) + base , 1 , 0 );
addedge( id(i,j) + base , id(i,j) , 1 , 0 );
}
puts( Dinic() == a ? "Yes" : "No" );
}
return 0;
}