Problem: There is a grid with H horizontal rows and W vertical columns, and there are obstacles on some of the squares.  Snuke is going to choose one of the squares not occupied by an obstacle and place a lamp on it. The lamp placed on the square will emit straight beams of light in four cardinal directions: up, down, left, and right. In each direction, the beam will continue traveling until it hits a square occupied by an obstacle or it hits the border of the grid. It will light all the squares on the way, including the square on which the lamp is placed, but not the square occupied by an obstacle.  Snuke wants to maximize the number of squares lighted by the lamp.  You are given H strings S_i (1 \leq i \leq H), each of length W. If the j-th character (1 \leq j \leq W) of S_i is `#`, there is an obstacle on the square at the i-th row from the top and the j-th column from the left; if that character is `.`, there is no obstacle on that square.  Find the maximum possible number of squares lighted by the lamp.

Submission attempts (in order):
import sys, math, heapq
stdin = sys.stdin
read_int = lambda : list(map(int,stdin.readline().split()))
read_str = lambda : stdin.readline().rstrip()

H, W = read_int()
grid = []
import numpy as np
L = np.zeros((H, W))
R = np.zeros((H, W))
U = np.zeros((H, W))
D = np.zeros((H, W))
res = np.zeros((H, W))
# L = [[0 for _ in range(W)] for _ in range(H)]
# R = [[0 for _ in range(W)] for _ in range(H)]
# U = [[0 for _ in range(W)] for _ in range(H)]
# D = [[0 for _ in range(W)] for _ in range(H)]
# res = [[0 for _ in range(W)] for _ in range(H)]

for _ in range(H):
    row = list(" ".join(input()).split())
    grid.append(row)

# r = []
# heapq.heapify(r)

#     heapq.heappush(r, 0)

def solve():
    ans = 0
    for i in range(H):
        for j in range(W):
            if grid[i][j] == "#":
                L[i][j] = U[i][j] = 0
            else:
                L[i][j] = 1 if j == 0 else L[i][j - 1] + 1
                U[i][j] = 1 if i == 0 else U[i - 1][j] + 1
                res[i][j] += L[i][j] + U[i][j]
                ans = max(res[i][j], ans)
                # r[i* W + j] -= L[i][j] + U[i][j]

            if grid[H - i - 1][W - j - 1] == "#":
                R[H - i - 1][W - j - 1] = D[H - i - 1][W - j - 1] = 0
            else:
                R[H - i - 1][W - j - 1] = 1 if j == 0 else R[H - i - 1][W - j] + 1
                D[H - i - 1][W - j - 1] = 1 if i == 0 else D[H - i][W - j - 1] + 1
                res[H - i - 1][W - j - 1] += R[H - i - 1][W - j - 1] + D[H - i - 1][W - j - 1]
                ans = max(res[H - i - 1][W - j - 1], ans)
                
    # print(res)
    return ans - 3

# if __name__ == "__main__":
#     print(solve())

if __name__ == "__main__":
    import timeit
    loop = 1
    n = 1000
    time = timeit.timeit(lambda: solve(), number=loop)
    print(time/loop)



import sys, math, heapq
stdin = sys.stdin
read_int = lambda : list(map(int,stdin.readline().split()))
H, W = read_int()


# nomal
# L = [[0 for _ in range(W)] for _ in range(H)]
# R = [[0 for _ in range(W)] for _ in range(H)]
# U = [[0 for _ in range(W)] for _ in range(H)]
# D = [[0 for _ in range(W)] for _ in range(H)]
res = [[0 for _ in range(W)] for _ in range(H)]
grid = []

# numpy
import numpy as np
L, R, U, D = [np.zeros((H, W), int) for _ in range(4)]
grid_np = np.zeros((H,W), int)


for i in range(H):
    row = list(" ".join(input()).split())
    grid.append(row)
    grid_np[i] = np.array(row) == '.'

def solve():
    ans = 0
    for i in range(H):
        for j in range(W):
            if grid[i][j] == ".":
                L[i][j] = 1 if j == 0 else L[i][j - 1] + 1
                U[i][j] = 1 if i == 0 else U[i - 1][j] + 1
                res[i][j] += L[i][j] + U[i][j]
                ans = max(res[i][j], ans)

            if grid[H - i - 1][W - j - 1] == ".":
                R[H - i - 1][W - j - 1] = 1 if j == 0 else R[H - i - 1][W - j] + 1
                D[H - i - 1][W - j - 1] = 1 if i == 0 else D[H - i][W - j - 1] + 1
                res[H - i - 1][W - j - 1] += R[H - i - 1][W - j - 1] + D[H - i - 1][W - j - 1]
                ans = max(res[H - i - 1][W - j - 1], ans)
    return ans - 3

def solve_np():
    ans = 0

    for j in range(W):
        
        if j == 0:
            L[:, j] = 1 * grid_np[:, j]
            R[:, -j - 1] = 1 * grid_np[:, -j - 1]
        else:
            L[:, j] = (L[:, j - 1] + 1) * grid_np[:, j]
            R[:, -j - 1] = (R[:, -j] + 1) * grid_np[:, -j - 1]

    for i in range(H):
        if i == 0:
            U[i] = 1 * grid_np[i]
            D[-i -1] = 1 * grid_np[-i -1]
        else:
            U[i] = (U[i - 1] + 1) * grid_np[i]
            D[-i - 1] = (D[-i] + 1) * grid_np[-i - 1]
    ans = np.max((L+R+U+D) -3)
    return ans

if __name__ == "__main__":
    # print(solve())
    print(solve_np())

# if __name__ == "__main__":

#     loop = 1
#     n = 1000

#     print(time/loop)



Verdicts: Wrong Answer

Accepted