UVa 12265 - Selling Land

链接：

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=3417

题意：

输入一个n*m（1≤n，m≤1000）矩阵，每个格子可能是空地，也可能是沼泽。对于每个空地格子，
求出以它为右下角的空地矩形的最大周长，然后统计每个周长出现了多少次。

分析：

扫描法 + 单调栈。
按照从上到下的顺序处理每一行，在每一行中从左到右处理每个格子（以下称为“当前格”），
找出以该格子为右下角的最大周长矩形（以下简称最优矩形）。只要找到了以每个
格子为右下角的最优矩形，本题就可以得到解决。
假定“当前格”已经固定，则只需要再确定一个左上角，就可以得到一个矩形。
下面用c表示第c列，h表示第c列的空地高度。

首先从上到下枚举“当前行”，然后从左到右枚举“当前列”。
在移动“当前列”的过程中，保存若干个(c,h)，按照c从小到大排列成有序表，
则h也是从小到大排列，并且h－c也是从小到大排列。则可以在O(1)时间
内求出每个当前格对应的最优矩形（因为最后一个矩形就是最优的），然后根据需要从右到
左删除一些矩形（也可能不删除），并且可能会把最右边的矩形变矮。然后，当且仅当新矩
形的h－c比它左边的矩形大时，加到表的最右边。由于添加和删除都在表的最右端，用一个
栈来实现即可。

代码：

 #include <cstdio>
 #include <cstring>

 const int UP =  + ;

 struct NODE {
     int c, h;
 } K[UP]; //栈

 int hei[UP], ans[UP*];
 char s[UP][UP];

 int main(){
     int T, n, m;
     scanf("%d", &T);
     while(T--){
         scanf("%d%d", &n, &m);
         for(int i = ; i < n; i++) scanf("%s", s[i]);

         memset(hei, , sizeof(hei));
         memset(ans, , sizeof(ans));
         for(int r = ; r < n; r++){
             int p = ; //栈顶指针
             for(int c = ; c < m; c++){
                 if(s[r][c] == '#') p = hei[c] = ;
                 else{
                     hei[c]++;
                     NODE nd = (NODE){c, hei[c]};
                     if(!p) K[++p] = nd;
                     else{
                         while(p && nd.h <= K[p].h) nd.c = K[p--].c;
                         if(!p || nd.h - nd.c > K[p].h - K[p].c) K[++p] = nd;
                     }
                     ans[K[p].h + c - K[p].c + ]++;
                 }
             }
         }

         for(int i = ; i <= m + n; i++) if(ans[i]){
             printf("%d x %d\n", ans[i], i * );
         }
     }
     return ;
 }