网络流二十四题之二十二 —— 火星探险问题

火星探险问题

---

Description

火星探险队的登陆舱将在火星表面着陆，登陆舱内有多部障碍物探测车。
登陆舱着陆后，探测车将离开登陆舱向先期到达的传送器方向移动。
探测车在移动中还必须采集岩石标本。
每一块岩石标本由最先遇到它的探测车完成采集。
每块岩石标本只能被采集一次。
岩石标本被采集后，其他探测车可以从原来岩石标本所在处通过。
探测车不能通过有障碍的地面。
本题限定探测车只能从登陆处沿着向南或向东的方向朝传送器移动，而且多个探测车可以在同一时间占据同一位置。
如果某个探测车在到达传送器以前不能继续前进，则该车所采集的岩石标本将全部损失。

用一个PXQ 网格表示登陆舱与传送器之间的位置。登陆舱的位置在($X_1,Y_1$)处，传送器的位置在($X_P ,Y_Q$)处。

给定每个位置的状态，计算探测车的最优移动方案，使到达传送器的探测车的数量最多，而且探测车采集到的岩石标本的数量最多。

---

Input

第 $1$ 行为探测车数，第 $2$ 行为 $P$ 的值，第 $3$ 行为 $Q$ 的值。
接下来的 $Q$ 行是表示登陆舱与传送器之间的位置状态的 PXQ 网格。
用 $3$ 个数字表示火星表面位置的状态：$0$ 表示平坦无障碍，$1$ 表示障碍，$2$ 表示石块。

---

Output

程序运行结束时，将能探索到的岩石标本数目输出。

---

Sample Input

2
10
8
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 1 1 0 0 0
0 0 0 1 0 2 0 0 0 0
1 1 0 1 2 0 0 0 0 1
0 1 0 0 2 0 1 1 0 0
0 1 0 1 0 0 1 1 0 0
0 1 2 0 0 0 0 1 0 0
0 0 0 0 0 0 0 0 0 0

---

Sample Output

2

---

Solution

我们把每个格子与他能够到达的格子连一条边。
如果此格子有岩石标本，那么就再加一个点，再连一条容量为 1 <script type="math/tex" id="MathJax-Element-13">1</script> 的边到这个点，然后从这个点连一条边回到这个格子。
跑一遍最大费用流就行了。

---

Code

[cpp]

1. #include <iostream>  
2. #include <cstdio>  
3. #include <cstring>  
4. #include <queue>  
5.   
6. #define ss 0  
7. #define tt 100000  
8. #define INF 0x3f3f3f3f  
9. #define Min(x,y) ((x)<(y)?(x):(y))  
10. #define PLA(x,y) ((x)*qs+(y))  
11.   
12. #define ANOTHER 1000  
13.   
14. using namespace std;  
15.   
16. int n,cnt,ans,ps,qs;  
17. int nxt[100010],head[100010],data[100010];  
18. int flow[100010],wei[100010];  
19. int dis[100010];  
20. int pre[100010];  
21. int maps[100][100];  
22. bool vis[100][100];  
23. queue<int>q;  
24. bool in_stack[100010];  
25.   
26. void add(int x,int y,int a,int b){    
27.     nxt[cnt]=head[x];data[cnt]=y;wei[cnt]=b;flow[cnt]=a;head[x]=cnt++;        
28.     nxt[cnt]=head[y];data[cnt]=x;wei[cnt]=-b;flow[cnt]=0;head[y]=cnt++;     
29. }    
30.     
31. bool BFS(){    
32.     memset(dis,0x3f,sizeof dis);dis[ss]=0;in_stack[ss]=true;q.push(ss);pre[ss]=pre[tt]=-1;    
33.     while(!q.empty()){    
34.         int now=q.front();    
35.         q.pop();    
36.         in_stack[now]=false;    
37.         for(int i=head[now];i!=-1;i=nxt[i]){    
38.             if(flow[i]!=0&&dis[data[i]]>dis[now]+wei[i]){    
39.                 dis[data[i]]=dis[now]+wei[i];    
40.                 pre[data[i]]=i^1;    
41.                 if(!in_stack[data[i]]){    
42.                     in_stack[data[i]]=true;    
43.                     q.push(data[i]);    
44.                 }    
45.             }    
46.         }    
47.     }    
48.     return pre[tt]!=-1;    
49. }    
50.     
51. void dfs(){    
52.     int Low=INF;    
53.     for(int i=pre[tt];i!=-1;i=pre[data[i]])Low=Min(Low,flow[i^1]);    
54.     for(int i=pre[tt];i!=-1;i=pre[data[i]])flow[i^1]-=Low,flow[i]+=Low;     
55.     ans+=Low*dis[tt];    
56. }  
57.   
58. void Dfs(int p1,int p2){  
59.     vis[p1][p2]=true;  
60.     if(p1+1<=ps&&maps[p1+1][p2]!=1){  
61.         if(!maps[p1+1][p2])add(PLA(p1,p2),PLA(p1+1,p2),n,0);  
62.         if(maps[p1+1][p2]==2){  
63.             add(PLA(p1,p2),PLA(p1+1,p2),n,0);  
64.             add(PLA(p1,p2),PLA(p1+1,p2)+ANOTHER,1,0);  
65.         }  
66.         else add(PLA(p1,p2),PLA(p1+1,p2),n,0);  
67.         if(!vis[p1+1][p2])Dfs(p1+1,p2);  
68.     }  
69.     if(p2+1<=qs&&maps[p1][p2+1]!=1){  
70.         if(!maps[p1][p2+1])add(PLA(p1,p2),PLA(p1,p2+1),n,0);  
71.         if(maps[p1][p2+1]==2){  
72.             add(PLA(p1,p2),PLA(p1,p2+1),n,0);  
73.             add(PLA(p1,p2),PLA(p1,p2+1)+ANOTHER,1,0);  
74.         }  
75.         else add(PLA(p1,p2),PLA(p1,p2+1),n,0);  
76.         if(!vis[p1][p2+1])Dfs(p1,p2+1);  
77.     }  
78. }  
79.   
80. int main(){  
81.   
82.     memset(head,-1,sizeof head);  
83.   
84.     scanf(”%d%d%d”,&n,&qs,&ps);  
85.       
86.     for(int i=1;i<=ps;i++)  
87.         for(int j=1;j<=qs;j++){  
88.             scanf(”%d”,&maps[i][j]);  
89.             add(PLA(i,j)+ANOTHER,PLA(i,j),1,-1);  
90.         }  
91.     Dfs(1,1);  
92.       
93.     add(ss,PLA(1,1),n,0);  
94.     add(PLA(ps,qs),tt,n,0);  
95.   
96.     while(BFS())dfs();  
97.       
98.     if(maps[ps][qs]==2)ans++;  
99.     if(maps[1][1]==2)ans++;  
100.     printf(”%d\n”,-ans);  
101.   
102.     return 0;  
103. }