jzxx2860引水入城

题目描述

在一个遥远的国度，一侧是风景秀美的湖泊，另一侧则是漫无边际的沙漠。该国的行政区划十分特殊，刚好构成一个 N行 M 列的矩形，如上图所示，其中每个格子都代表一座城市，每座城市都有一个海拔高度。
为了使居民们都尽可能饮用到清澈的湖水，现在要在某些城市建造水利设施。水利设施有两种，分别为蓄水厂和输水站。蓄水厂的功能是利用水泵将湖泊中的水抽取到所在城市的蓄水池中。因此，只有与湖泊毗邻的第 1 行的城市可以建造蓄水厂。而输水站的功能则是通过输水管线利用高度落差，将湖水从高处向低处输送。故一座城市能建造输水站的前提，是存在比它海拔更高且拥有公共边的相邻城市，已经建有水利设施。
由于第 N 行的城市靠近沙漠，是该国的干旱区，所以要求其中的每座城市都建有水利设施。那么，这个要求能否满足呢？如果能，请计算最少建造几个蓄水厂；如果不能，求干旱区中不可能建有水利设施的城市数目。

输入
输入文件的每行中两个数之间用一个空格隔开。
输入的第一行是两个正整数 N和 M，表示矩形的规模。
接下来 N行，每行 M 个正整数，依次代表每座城市的海拔高度。

输出
输出有两行。如果能满足要求，输出的第一行是整数 1，第二行是一个整数，代表最少建造几个蓄水厂；如果不能满足要求，输出的第一行是整数 0，第二行是一个整数，代表有几座干旱区中的城市不可能建有水利设施。

样例
输入1
2 5
9 1 5 4 3 1
8 7 6 1 2
输出1
1
输入2
3 6
8 4 5 6 4 4
7 3 4 3 3 3
3 2 2 1 1 2
输出2
1
3

提示
【样例 1 说明】
只需要在海拔为 9 的那座城市中建造蓄水厂，即可满足要求。

【样例 2 说明】

上图中，在 3 个粗线框出的城市中建造蓄水厂，可以满足要求。以这 3个蓄水厂为源头在干旱区中建造的输水站分别用3种颜色标出。当然，建造方法可能不唯一。

传送门

满分代码：
const h:array[1..4,1..2] of longint=((1,0),(-1,0),(0,-1),(0,1));
var col,a:array[0..501,0..501] of longint;f:array[0..500] of longint;c:array[0..500,1..2] of longint;d:array[0..500000] of record x,y:longint;end;n,m,i,j,color:longint;
function min(a,b:longint):longint;
begin if a<b then exit(a) else exit(b); end;
procedure Judge;
var i,j:longint;procedure bfs(x,y:integer);var i,t,f:longint;begint:=1;f:=1;d[t].x:=x;d[t].y:=y;col[x,y]:=color;repeatfor i:=1 to 4 doif col[d[f].x+h[i,1],d[f].y+h[i,2]]=0 thenif a[d[f].x+h[i,1],d[f].y+h[i,2]]<a[d[f].x,d[f].y] then begininc(t);d[t].x:=d[f].x+h[i,1]; d[t].y:=d[f].y+h[i,2];col[d[t].x,d[t].y]:=color;end;inc(f);until f>t;end;   {end bfs}
begin  {judge}for i:=0 to n+1 do begina[i,0]:=maxlongint;a[i,m+1]:=maxlongint;end;for i:=0 to m+1 do begina[0,i]:=maxlongint;a[n+1,i]:=maxlongint;end;for color:=1 to m do bfs(1,color);j:=0;for i:=1 to m doif col[n,i]=0 then inc(j);if j>0 then beginwriteln(0);writeln(j);halt;end;
end;
procedure floodfill;procedure bfs(x,y:integer);var i,t,f:longint;begint:=1;f:=1;d[t].x:=x;d[t].y:=y;col[x,y]:=color;repeatfor i:=1 to 4 doif col[d[f].x+h[i,1],d[f].y+h[i,2]]=0 thenif a[d[f].x+h[i,1],d[f].y+h[i,2]]>a[d[f].x,d[f].y] then begininc(t);d[t].x:=d[f].x+h[i,1]; d[t].y:=d[f].y+h[i,2];col[d[t].x,d[t].y]:=color;end;inc(f);until f>t;end;
var i,j:longint;
beginfor i:=0 to n+1 do begina[i,0]:=0; a[i,m+1]:=0;end;for i:=0 to m+1 do begina[0,i]:=0; a[n+1,i]:=0;end;fillchar(col,sizeof(col),0);for color:=1 to m doif col[n,color]=0 thenbfs(n,color);for i:=1 to m doc[i,1]:=col[1,i];fillchar(col,sizeof(col),0);for color:=m downto 1 doif col[n,color]=0 then bfs(n,color);for i:=1 to m do c[i,2]:=col[1,i];
end;
procedure DP;
var i,j:longint;
begin
f[0]:=0;
for i:=1 to m do
beginf[i]:=maxint;for j:=1 to m doif (c[j,2]>=i)and(c[j,1]<=i) then
f[i]:=min(f[i],f[c[j,1]-1]+1);
end;
writeln(1);
writeln(f[m]);
end;
BEGIN   {main}read(n,m);for i:=1 to n dofor j:=1 to m doread(a[i,j]);judge;floodfill;DP;
END.