#Union-Findライブラリ(ac-library python) #https://github.com/not522/ac-library-python/blob/master/atcoder/dsu.py import typing class DSU: ''' Implement (union by size) + (path halving) Reference: Zvi Galil and Giuseppe F. Italiano, Data structures and algorithms for disjoint set union problems ''' def __init__(self, n: int = 0) -> None: self._n = n self.parent_or_size = [-1] * n def merge(self, a: int, b: int) -> int: assert 0 <= a < self._n assert 0 <= b < self._n x = self.leader(a) y = self.leader(b) if x == y: return x if -self.parent_or_size[x] < -self.parent_or_size[y]: x, y = y, x self.parent_or_size[x] += self.parent_or_size[y] self.parent_or_size[y] = x return x def same(self, a: int, b: int) -> bool: assert 0 <= a < self._n assert 0 <= b < self._n return self.leader(a) == self.leader(b) def leader(self, a: int) -> int: assert 0 <= a < self._n parent = self.parent_or_size[a] while parent >= 0: if self.parent_or_size[parent] < 0: return parent self.parent_or_size[a], a, parent = ( self.parent_or_size[parent], self.parent_or_size[parent], self.parent_or_size[self.parent_or_size[parent]] ) return a def size(self, a: int) -> int: assert 0 <= a < self._n return -self.parent_or_size[self.leader(a)] def groups(self) -> typing.List[typing.List[int]]: leader_buf = [self.leader(i) for i in range(self._n)] result: typing.List[typing.List[int]] = [[] for _ in range(self._n)] for i in range(self._n): result[leader_buf[i]].append(i) return list(filter(lambda r: r, result)) #ライブラリ終わり H, W = map(int, input().split()) d = DSU(H+W); for i in range(H): s = input() for j in range(W): if s[j] == '@': d.merge(i, j+H) res = d.groups() K = len(res) dpmin = [[10**9 for j in range(H+1)] for i in range(K+1)] dpmax = [[-10**9 for j in range(H+1)] for i in range(K+1)] dpmin[0][0] = 0 dpmax[0][0] = 0 for i in range(K): rcnt, ccnt = 0, 0 for v in res[i]: if v < H: rcnt+=1 else: ccnt+=1 for j in range(H): if dpmin[i][j] != 10**9: dpmin[i+1][j+rcnt] = min(dpmin[i+1][j+rcnt], dpmin[i][j] + ccnt) dpmax[i+1][j+rcnt] = max(dpmax[i+1][j+rcnt], dpmax[i][j] + ccnt) dpmin[i+1][j] = min(dpmin[i+1][j], dpmin[i][j]) dpmax[i+1][j] = max(dpmax[i+1][j], dpmax[i][j]) ans = 0 for i in range(H): if dpmin[K][i] != 10**9: ans = max(ans, i*(W-dpmin[K][i])+(H-i)*dpmin[K][i]) ans = max(ans, i*(W-dpmax[K][i])+(H-i)*dpmax[K][i]) print(ans)