Problem:


A Garden with Ponds

  Mr. Gardiner is a modern garden designer who is excellent at utilizing the terrain features.
  His design method is unique: he first decides the location of ponds and design them with the terrain features intact.

  According to his unique design procedure, all of his ponds are rectangular with simple aspect ratios.
  First, Mr. Gardiner draws a regular grid on the map of the garden site so that the land is divided into cells of unit square, and annotates every cell with its elevation.
  In his design method, a pond occupies a rectangular area consisting of a number of cells.
  Each of its outermost cells has to be higher than all of its inner cells.
  For instance, in the following grid map, in which numbers are elevations of cells, a pond can occupy the shaded area, where the outermost cells are shaded darker and the inner cells are shaded lighter.
  You can easily see that the elevations of the outermost cells are at least three and those of the inner ones are at most two.

  A rectangular area on which a pond is built must have at least one inner cell.
  Therefore, both its width and depth are at least three.

  When you pour water at an inner cell of a pond, the water can be kept in the pond until its level reaches that of the lowest outermost cells.
  If you continue pouring, the water inevitably spills over.  
  Mr. Gardiner considers the larger capacity the pond has, the better it is.
  Here, the capacity of a pond is the maximum amount of water it can keep.
  For instance, when a pond is built on the shaded area in the above map, its capacity is (3 &#8722; 1) + (3 &#8722; 0) + (3 &#8722; 2) = 6, where 3 is the lowest elevation of the outermost cells and 1, 0, 2 are the elevations of the inner cells.
  Your mission is to write a computer program that, given a grid map describing the elevation of each unit square cell, calculates the largest possible capacity of a pond built in the site.

  Note that neither of the following rectangular areas can be a pond.
  In the left one, the cell at the bottom right corner is not higher than the inner cell.
  In the right one, the central cell is as high as the outermost cells.

  

Input

The input consists of at most 100 datasets, each in the following format.

d w
e1, 1 ... e1, w
...
ed, 1 ... ed, w

The first line contains d and w, representing the depth and the width, respectively, of the garden site described in the map. 
They are positive integers between 3 and 10, inclusive.
Each of the following d lines contains w integers between 0 and 9, inclusive, separated by a space.
The x-th integer in the y-th line of the d lines is the elevation of the unit square cell with coordinates (x, y).

The end of the input is indicated by a line containing two zeros separated by a space.

Output

  For each dataset, output a single line containing the largest possible capacity of a pond that can be built in the garden site described in the dataset.
  If no ponds can be built, output a single line containing a zero.

Sample Input

3 3
2 3 2
2 1 2
2 3 1
3 5
3 3 4 3 3
3 1 0 2 3
3 3 4 3 2
7 7
1 1 1 1 1 0 0
1 0 0 0 1 0 0
1 0 1 1 1 1 1
1 0 1 0 1 0 1
1 1 1 1 1 0 1
0 0 1 0 0 0 1
0 0 1 1 1 1 1
6 6
1 1 1 1 2 2
1 0 0 2 0 2
1 0 0 2 0 2
3 3 3 9 9 9
3 0 0 9 0 9
3 3 3 9 9 9
0 0

Output for the Sample Input

0
3
1
9



Sample input:
3 3
2 3 2
2 1 2
2 3 1
3 5
3 3 4 3 3
3 1 0 2 3
3 3 4 3 2
7 7
1 1 1 1 1 0 0
1 0 0 0 1 0 0
1 0 1 1 1 1 1
1 0 1 0 1 0 1
1 1 1 1 1 0 1
0 0 1 0 0 0 1
0 0 1 1 1 1 1
6 6
1 1 1 1 2 2
1 0 0 2 0 2
1 0 0 2 0 2
3 3 3 9 9 9
3 0 0 9 0 9
3 3 3 9 9 9
0 0


Sample output:
0
3
1
9

Submission status: Accepted

Solution code:
import math,string,itertools,fractions,heapq,collections,re,array,bisect,sys,random,time,copy,functools

sys.setrecursionlimit(10**7)
inf = 10**20
eps = 1.0 / 10**13
mod = 10**9+7
dd = [(-1,0),(0,1),(1,0),(0,-1)]
ddn = [(-1,0),(-1,1),(0,1),(1,1),(1,0),(1,-1),(0,-1),(-1,-1)]

def LI(): return [int(x) for x in sys.stdin.readline().split()]
def LI_(): return [int(x)-1 for x in sys.stdin.readline().split()]
def LF(): return [float(x) for x in sys.stdin.readline().split()]
def LS(): return sys.stdin.readline().split()
def I(): return int(sys.stdin.readline())
def F(): return float(sys.stdin.readline())
def S(): return input()
def pf(s): return print(s, flush=True)


def main():
    rr = []

    def f(d,w):
        aa = [LI() for _ in range(d)]
        r = 0
        for i in range(d-2):
            for j in range(w-2):
                for k in range(i+2,d):
                    for l in range(j+2,w):
                        mi = inf
                        for m in range(i,k+1):
                            t = aa[m][j]
                            if mi > t:
                                mi = t
                            t = aa[m][l]
                            if mi > t:
                                mi = t
                        for m in range(j,l+1):
                            t = aa[i][m]
                            if mi > t:
                                mi = t
                            t = aa[k][m]
                            if mi > t:
                                mi = t
                        s = 0
                        for m in range(i+1,k):
                            for n in range(j+1,l):
                                if mi <= aa[m][n]:
                                    s = -inf
                                    break
                                s += mi - aa[m][n]
                        if r < s:
                            r = s
        return r


    while 1:
        n,m = LI()
        if n == 0:
            break
        rr.append(f(n,m))

    return '\n'.join(map(str, rr))


print(main())

Submission status: Runtime Error

Solution code:
while 1:
    d, w = map(int, input().split())
    if d == w == 0:
        break
    E = [[*map(int, input().split())] for i in range(d)]
    ans = 0
    for x2 in range(d):
        for y2 in range(w):
            for x1 in range(x2-1):
                for y1 in range(y2-1):
                    # [x1, x2] x [y1, y2]
                    co = 10**18
                    for x in range(x1, x2+1):
                        co = min(co, E[x][y1], E[x][y2])
                    for y in range(y1, y2+1):
                        co = min(co, E[x1][y], E[x2][y])
                    ci = 0
                    cnt = 0
                    for x in range(x1+1, x2):
                        for y in range(y1+1, y2):
                            ci = max(ci, E[x][y])
                            cnt += E[x][y]
                    if ci < co:
                        ans = max(ans, (x2-x1-1)*(y2-y1-1)*co - cnt)
    print(ans)

Submission status: Runtime Error

Solution code:
while 1:
    d, w = map(int, input().split())
    if d == w == 0:
        break
    E = [[*map(int, input().split())] for i in range(d)]
    ans = 0
    for x2 in range(d):
        for y2 in range(w):
            for x1 in range(x2-1):
                for y1 in range(y2-1):
                    co = 10**18
                    for x in range(x1, x2+1):
                        co = min(co, E[x][y1], E[x][y2])
                    for y in range(y1, y2+1):
                        co = min(co, E[x1][y], E[x2][y])
                    ci = 0
                    cnt = 0
                    for x in range(x1+1, x2):
                        for y in range(y1+1, y2):
                            ci = max(ci, E[x][y])
                            cnt += E[x][y]
                    if ci < co:
                        ans = max(ans, (x2-x1-1)*(y2-y1-1)*co - cnt)
    print(ans)

Submission status: Runtime Error

Solution code:
while 1:
    d, w = map(int, input().split())
    if d==w==0:
        break
    E=[list(map(int, input().split())) for i in range(d)]
    a=0
    for x2 in range(d):
        for y2 in range(w):
            for x1 in range(x2-1):
                for y1 in range(y2-1):
                    co=10**18
                    for x in range(x1, x2+1):
                        co=min(co, E[x][y1], E[x][y2])
                    for y in range(y1, y2+1):
                        co=min(co, E[x1][y], E[x2][y])
                    ci=cnt=0
                    for x in range(x1+1, x2):
                        for y in range(y1+1, y2):
                            ci=max(ci, E[x][y])
                            cnt+=E[x][y]
                    if ci < co:
                        a=max(ans, (x2-x1-1)*(y2-y1-1)*co - cnt)
    print(a)

Submission status: Runtime Error

Solution code:
while True:
    d, w = map(int, input().split())
    if d==w==0:
        break
    E=[list(map(int, input().split())) for i in range(d)]
    a=0
    for x2 in range(d):
        for y2 in range(w):
            for x1 in range(x2-1):
                for y1 in range(y2-1):
                    co=10**18
                    for x in range(x1, x2+1):
                        co=min(co, E[x][y1], E[x][y2])
                    for y in range(y1, y2+1):
                        co=min(co, E[x1][y], E[x2][y])
                    ci=cnt=0
                    for x in range(x1+1, x2):
                        for y in range(y1+1, y2):
                            ci=max(ci, E[x][y])
                            cnt+=E[x][y]
                    if ci < co:
                        a=max(ans, (x2-x1-1)*(y2-y1-1)*co - cnt)
    print(a)