寻找图中的数字8

这是一道来自UASCO的一道有趣的题目（2013 US Open, Gold），给出一个图：

...............
...............
...*******.....
.*....*.......*
.*......*....*.
....*..........
...*...****....
...............
..**.*..*..*...
...*...**.*....
*..*...*.......
...............
.....*..*......
.........*.....
...............

求出图中美观程度最大的一个数字8，要求不能含有*字符：

1  ..88888888888..
2  ..8.........8..
3  ..8*******..8..
4  .*8...*.....8.*
5  .*8.....*...8*.
6  ..8.*.......8..
7  ..8*...****.8..
8  .88888888888888
9  .8**.*..*..*..8
10 .8.*...**.*...8
11 *8.*...*......8
12 .8............8
13 .8...*..*.....8
14 .8.......*....8
15 .88888888888888

所谓的数字8， 就是由上下两个矩形组成，并且下面矩形的顶部要包含（或重叠）上面矩形的底部的图形。美观程度就是两个矩形面积的乘积。如上图，上面的矩形是6 * 9 = 54，下面是12 * 6 = 72，所以美观程度就是65 * 72 = 3888。

首先，枚举上下两个矩形的分界点i，如上图就是第8行。因为上下两个矩形有限制条件，不好计算，所以先确定下面矩形的左右两条边所在的列，如下图，即left和right两条线。然后，上下两个矩形的面积就可以分开算了。当确定了i，left，right三条线后，下面矩形的最大面积就好求了（用一条扫描线即可），而上面的就要动点脑筋了。

设f（i, left, right）为当上面矩形的底部在线i上，左边的边在left的右边，右边的边在right的左边，的最大面积。如图，蓝色矩形是合法的。

那么有f（i, left, right） = max（f（i, left + 1, right）， f（i, left, right - 1）， 上面矩形的底部在i上且左右两边在线left、right的最大面积（算法和下面的矩形一样））。

这时，问题就完美解决了，时间复杂度是三次方的。实现的时候，空间要优化成平方级别的，这就需要以right - left的大小为第一枚举顺序，然后用滚动数组。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int N = 307;

// 滚动数组 
int n, f[N][N], g[N][N];
// 读入数据 
char d[N][N];
// con[i][j]: 点i, j，往右连续size个地方都是字符'.' 
bool con[N][N];

int main() {
    freopen("eight.in", "r", stdin);
    freopen("eight.out", "w", stdout);
    
    scanf("%d\n", &n);
    for (int i = 0; i < n; i ++) {
        for (int j = 0; j < n; j ++) {
            scanf("%c", &d[i][j]);
            if (j) con[i][j - 1] = (d[i][j - 1] == '.' && d[i][j] == '.');
        }
        scanf("\n");
    }
    
    long long res = 0LL;
    memset(f, 0, sizeof(f));
    
    // right - left 的长度 
    for (int s = 2; s < n; s ++) {
        memcpy(g, f, sizeof(f));
        
        for (int j = 0; j + s < n; j ++) {
            int top = n;
            // 更新con 
            for (int i = n - 1; i >= 0; i --) {
                con[i][j] &= (d[i][j + s] == '.');
                if (con[i][j]) top = i;
            }
            
            // line为扫描线
            // top能使上矩形面积最大的顶部所在的行
            // 同理，bottom为能使下矩形面积最大的底部所在的行 
            int line = top, bottom = top;
            for (int i = top; i + 4 < n; i ++) {
                // 枚举时同时更新top 
                if (d[i][j] == '*' || d[i][j + s] == '*') top = i + 1;
                
                // bottom需要line配合更新，bottom只增不减，不会提高复杂度 
                line >?= i;
                while (line < n && d[line][j] == '.' && d[line][j + s] == '.') {
                    if (con[line][j]) bottom = line;
                    line ++;
                }
                
                // 当前行能够容下矩形 
                if (con[i][j]) {
                    // 分别为上下两个矩形的面积 
                    int top_rec = max(0, i - top - 1) * (s - 1);
                    int bottom_rec = max(0, bottom - i - 1) * (s - 1);
                    
                    // 更新f 
                    f[i][j] = max(g[i][j], g[i][j + 1]);
                    f[i][j] >?= top_rec;
                    
                    // 更新答案 
                    res >?= (long long) f[i][j] * bottom_rec;
                }
            }
        }
    }
    
    printf("%I64d\n", res == 0LL ? -1LL : res);
    
    return 0;
}