[NOIP2010]引水入城

 1 #include <cstdio>
 2 #include <iostream>
 3 #include <algorithm>
 4 #include <cmath>
 5 #include <cctype>
 6 #include <queue>
 7 using namespace std;
 8 #if defined DEBUG
 9 FILE *in = fopen("test","r");
10 #define out stdout
11 #else
12 FILE *in = fopen("flow.in","r");
13 FILE *out = fopen("flow.out","w");
14 #endif
15 inline void getint(int &x){
16     int c = fgetc(in);
17     while(!isdigit(c))c = fgetc(in);
18     x = c - ‘0‘;
19     while(isdigit(c = fgetc(in)))x = x * 10 - ‘0‘ + c;
20 }
21 /*==================================================*/
22 const int maxn = 500 + 2;
23 typedef pair<int, int> RNG;//Range
24 inline bool cmp(const RNG &a, const RNG &b){
25     if(!(a.first | a.second))return 0;
26     if(!(b.first | b.second))return 1;
27     if(a.first == b.first)return a.second > b.second;
28     return a.first < b.first;
29 }
30 inline RNG operator + (const RNG &a, const RNG &b){//Union
31     RNG x, y;
32     int l, r;
33     if(cmp(b, a))x = b, y = a;
34     else x = a, y = b;
35     if(!(y.first | y.second))return x;
36     l = x.first, r = max(x.second, y.second);
37     return RNG(l, r);
38 }
39 RNG rng[maxn][maxn]; //‘Wet‘ Ranges (caused by water in G[i][j])
40 int N, M, G[maxn][maxn];
41 bool known[maxn][maxn] = {0}, wet[maxn] = {0};
42 inline void dfs(int x, int y){ //Dynamic Programming
43     RNG ans = RNG(0, 0);
44     if(x == N-1)wet[y] = 1, ans = RNG(y, y+1);
45     if(x && G[x-1][y] < G[x][y]){
46         if(!known[x-1][y])dfs(x-1, y);
47         ans = ans + rng[x-1][y];
48     }
49     if(y && G[x][y-1] < G[x][y]){
50         if(!known[x][y-1])dfs(x, y-1);
51         ans = ans + rng[x][y-1];
52     }
53     if(x < N-1 && G[x+1][y] < G[x][y]){
54         if(!known[x+1][y])dfs(x+1, y);
55         ans = ans + rng[x+1][y];
56     }
57     if(y < M-1 && G[x][y+1] < G[x][y]){
58         if(!known[x][y+1])dfs(x, y+1);
59         ans = ans + rng[x][y+1];
60     }
61     known[x][y] = 1;
62     rng[x][y] = ans;
63 }
64 inline int cover(){ //Greedy Algorithm - Line Coverage
65     sort(rng[0], rng[0] + M, cmp);
66     int ans = 1, curr, maxr, i;
67     curr = rng[0][0].second;
68     for(i = 1;i < M && (rng[0][i].first | rng[0][i].second);){
69         maxr = curr;
70         while((rng[0][i].first|rng[0][i].second)&&rng[0][i].first <= curr){
71             if(rng[0][i].second > maxr)maxr = rng[0][i].second;
72             ++i;
73         }
74         if(maxr > curr)++ans, curr = maxr;
75     }
76     return ans;
77 }
78 inline void work(){
79     getint(N), getint(M);
80     int i, j, dry = 0;
81     for(i = 0;i < N;++i)for(j = 0;j < M;++j)
82         getint(G[i][j]);
83     for(j = 0;j < M;++j)
84         dfs(0, j);
85     for(j = 0;j < M;++j)
86         if(!wet[j])++dry;
87     if(dry) fprintf(out, "0\n%d\n", dry);
88     else fprintf(out, "1\n%d\n", cover());
89 }
90 int main(){
91     work();
92     return 0;
93 }