Problem:
長方形探索
縦に
H
行、横に
W
列並べられた、合計
W × H
のマス目があります。いくつかのマス目には印がついています。各マス目の印の状態を読み込み、印のついていないマス目だけからなる最大の長方形の面積を出力するプログラムを作成してください。
入力データは 1 行
W
文字から構成され、
H
行が与えられます。たとえば以下のようなデータが与えられます。
..*....**.
..........
**....***.
....*.....
..*.......
...**.....
.*.*......
..........
..**......
.*..*.....
入力データの一行が、一行のマス目を表現します。入力データの文字列のうち、. (ピリオド) は印のついていないマス目、* (アスタリスク) は印のついているマス目を示しています。入力データの文字列には、ピリオド、アスタリスク、改行以外の文字は含まれません。
上記の例では、下図の 0 で示される長方形が最大となります。
..*....**.
..........
**....***.
....*00000
..*..00000
...**00000
.*.*.00000
.....00000
..**.00000
.*..*00000
よって、35 と出力すれば正解になります。なお、すべてのマス目に印がついている場合には、0 を出力してください。
Input
複数のデータセットが与えられます。各データセットはスペースで区切られた
H
と
W
からなる行から始まり、つづいて
H × W
の長方形が与えられます。
H, W
はともに 500　以下とします。
入力は２つの 0 を含む行で終わります。データセットの数は 20 を超えません。
Output
各データセットごとに、最大の長方形の面積を１行に出力してください。
Sample Input
10 10
...*....**
..........
**....**..
........*.
..*.......
**........
.*........
..........
....*..***
.*....*...
10 10
..*....*..
.*.*...*..
*****..*..
*...*..*..
*...*..*..
..........
****.*...*
..*..*...*
.*...*...*
****..***.
2 3
...
...
0 0
Output for the Sample Input
28
12
6

Submitted solutions:

code_snippet: #pragma region Macros

#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,fma,abm,mmx,avx,avx2")

#include <bits/extc++.h>
// #include <atcoder/all>
// using namespace atcoder;
using namespace std;
using namespace __gnu_pbds;

// #include <boost/multiprecision/cpp_dec_float.hpp>
// #include <boost/multiprecision/cpp_int.hpp>
// namespace mp = boost::multiprecision;
// using Bint = mp::cpp_int;
// using Bdouble = mp::number<mp::cpp_dec_float<256>>;
// Bdouble Beps = 0.00000000000000000000000000000001; // 1e-32
// const bool equals(Bdouble a, Bdouble b) { return mp::fabs(a - b) < Beps; }

#define pb emplace_back
#define int ll
#define endl '\n'
// #define unordered_map<int, int> gp_hash_table<int, int, custom_hash> 

#define sqrt __builtin_sqrtl
#define cbrt __builtin_cbrtl
#define hypot __builtin_hypotl

#define next asdnext
#define prev asdprev

using ll = long long;
using ld = long double;
const ld PI = acosl(-1);
const int INF = 1 << 30;
const ll INFL = 1LL << 61;
// const int MOD = 998244353;
const int MOD = 1000000007;

const ld EPS = 1e-10;
const bool equals(ld a, ld b) { return fabs((a) - (b)) < EPS; }

const vector<int> dx = {0, 1, 0, -1, 1, 1, -1, -1}; // → ↓ ← ↑ ↘ ↙ ↖ ↗
const vector<int> dy = {1, 0, -1, 0, 1, -1, -1, 1};

struct Edge {
    int from, to, cost;
    Edge() : from(-1), to(-1), cost(-1) {}
    Edge(int to, int cost) : to(to), cost(cost) {}
    Edge(int from, int to, int cost) : from(from), to(to), cost(cost) {}
    bool operator ==(const Edge& e) {
        return this->from == e.from && this->to == e.to && this->cost == e.cost;
    }
    bool operator !=(const Edge& e) {
        return this->from != e.from or this->to != e.to or this->cost != e.cost;
    }
};

chrono::system_clock::time_point start;
__attribute__((constructor))
void constructor() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout << fixed << setprecision(10);
    start = chrono::system_clock::now();
}

random_device seed_gen;
mt19937_64 rng(seed_gen());
uniform_int_distribution<int> dist_x(0, 1e9);
struct RNG {
    unsigned Int(unsigned l, unsigned r) {
        return dist_x(rng) % (r - l + 1) + l;
    }
    ld Double() {
        return ld(dist_x(rng)) / 1e9;
    }
} rnd;

using i64 = ll;
// using i64 = uint64_t;
// bit演算, x==0の場合は例外処理した方がよさそう. 区間は [l, r)
i64 lrmask(i64 l, i64 r) { return (1LL << r) - (1LL << l); }
i64 sub_bit(i64 x, i64 l, i64 r) { i64 b = x & ((1LL << r) - (1LL << l)); return b >> l; } // r溢れ可
i64 bit_width(i64 x) { return 64 - __builtin_clzll(x) + (x == 0); }

i64 popcount(i64 x) { return __builtin_popcountll(x); }
i64 popcount(i64 x, i64 l, i64 r) { return __builtin_popcountll(sub_bit(x, l, r)); }
i64 unpopcount(i64 x) { return bit_width(x) - __builtin_popcountll(x); }
i64 unpopcount(i64 x, i64 l, i64 r) { return r - l - __builtin_popcountll(sub_bit(x, l, r)); }
bool is_pow2(i64 x) { return __builtin_popcountll(x) == 1; }

i64 top_bit(i64 x) { return 63 - __builtin_clzll(x);} // 2^kの位 (x > 0)
i64 bot_bit(i64 x) { return __builtin_ctz(x);} // 2^kの位 (x > 0)
// i64 next_bit(i64 x, i64 k) { return 0; }
// i64 prev_bit(i64 x, i64 k) { return 0; }
// i64 kth_bit(i64 x, i64 k) { return 0; }
i64 MSB(i64 x) { if (x == 0) return 0; return 1LL << (63 - __builtin_clzll(x)); } // mask
i64 LSB(i64 x) { return (x & -x); } // mask

i64 countl_zero(i64 x) { return __builtin_clzll(x); }
i64 countl_one(i64 x) {
    i64 ret = 0, k = 63 - __builtin_clzll(x);
    while (k != -1 && (x & (1LL << k))) { k--; ret++; }
    return ret;
}
i64 countr_zero(i64 x) { return __builtin_ctzll(x); } // x==0のとき64が返ることに注意
i64 countr_one(i64 x) { i64 ret = 0; while (x & 1) { x >>= 1; ret++; } return ret; }

i64 floor_log2(i64 x) { if (x == 0) return 0; return 63 - __builtin_clzll(x); }
i64 bit_floor(i64 x) { if (x == 0) return 0; return 1LL << (63 - __builtin_clzll(x)); } // MSBと同じ
i64 ceil_log2(i64 x) { if (x == 0) return 0; return 64 - __builtin_clzll(x - 1); }
i64 bit_ceil(i64 x) { if (x == 0) return 0; return 1LL << (64 - __builtin_clzll(x - 1)); }

i64 rotl(i64 x, i64 k) { // 有効bit内でrotate. オーバーフロー注意
    i64 w = bit_width(x); k %= w;
    return ((x << k) | (x >> (w - k))) & ((1LL << w) - 1);
}
// i64 rotl(i64 x, i64 l, i64 m, i64 r) {}
i64 rotr(i64 x, i64 k) {
    i64 w = bit_width(x); k %= w;
    return ((x >> k) | (x << (w - k))) & ((1LL << w) - 1);
}
// i64 rotr(i64 x, i64 l, i64 m, i64 r) {}
i64 bit_reverse(i64 x) { // 有効bit内で左右反転
    i64 r = 0, w = bit_width(x);
    for (i64 i = 0; i < w; i++) r |= ((x >> i) & 1) << (w - i - 1);
    return r;
}
// i64 bit_reverse(i64 x, i64 l, i64 r) { return 0; }

bool is_palindrome(i64 x) { return x == bit_reverse(x); }
bool is_palindrome(i64 x, i64 l, i64 r) { i64 b = sub_bit(x, l, r); return b == bit_reverse(b); }
i64 concat(i64 a, i64 b) { return (a << bit_width(b)) | b; } // オーバーフロー注意
i64 erase(i64 x, i64 l, i64 r) { return x>>r<<l | x&(1LL<<l - 1); } // [l, r) をカット

i64 hamming(i64 a, i64 b) { return __builtin_popcountll(a ^ b); }
i64 hamming(i64 a, i64 b, i64 l, i64 r) { return __builtin_popcountll(sub_bit(a, l, r) ^ sub_bit(b, l, r)); }
i64 compcount(i64 x) { return (__builtin_popcountll(x ^ (x >> 1)) + (x & 1)) / 2; }
i64 compcount2(i64 x) { return compcount(x & (x >> 1)); } // 長さ2以上の連結成分の個数
i64 adjacount(i64 x) { return __builtin_popcountll(x & (x >> 1)); } // 隣接する1のペアの個数

i64 next_combination(i64 x) {
    i64 t = x | (x - 1); return (t + 1) | (((~t & -~t) - 1) >> (__builtin_ctz(x) + 1));
}


__int128_t POW(__int128_t x, int n) {
    __int128_t ret = 1;
    assert(n >= 0);
    if (x == 1 or n == 0) ret = 1;
    else if (x == -1 && n % 2 == 0) ret = 1; 
    else if (x == -1) ret = -1; 
    else if (n % 2 == 0) {
        assert(x < INFL);
        ret = POW(x * x, n / 2);
    } else {
        assert(x < INFL);
        ret = x * POW(x, n - 1);
    }
    return ret;
}
int per(int x, int y) { // x = qy + r (0 <= r < y) を満たすq
    assert(y != 0);
    if (x >= 0 && y > 0) return x / y;
    if (x >= 0 && y < 0) return x / y - (x % y < 0);
    if (x < 0 && y < 0) return x / y + (x % y < 0);
    return x / y - (x % y < 0); //  (x < 0 && y > 0) 
}
int mod(int x, int y) { // x = qy + r (0 <= r < y) を満たすr
    assert(y != 0);
    if (x >= 0) return x % y;
    __int128_t ret = x % y; // (x < 0)
    ret += (__int128_t)abs(y) * INFL;
    ret %= abs(y);
    return ret;
}
int floor(int x, int y) { // (ld)x / y 以下の最大の整数
    assert(y != 0);
    if (y < 0) x = -x, y = -y;
    return x >= 0 ? x / y : (x + 1) / y - 1;
}
int ceil(int x, int y) { // (ld)x / y 以上の最小の整数
    assert(y != 0);
    if (y < 0) x = -x, y = -y;
    return x > 0 ? (x - 1) / y + 1 : x / y;
}
int round(int x, int y) {
    assert(y != 0);
    return (x * 2 + y) / (y * 2);
}
int round(int x, int y, int k) { // (ld)(x/y)を10^kの位に関して四捨五入
    assert(y != 0); // TODO
    return INF;
}
int round2(int x, int y) { // 五捨五超入 // 未verify
    assert(y != 0);
    if (y < 0) y = -y, x = -x;
    int z = z / y;
    if ((z * 2 + 1) * y <= y * 2) z++;
    return z;
}
// int round(ld x, int k) { // xを10^kの位に関して四捨五入
// }
// int floor(ld x, int k) { // xを10^kの位に関してflooring
// }
// int ceil(ld x, int k) { // xを10^kの位に関してceiling
// }
// int kth(int x, int y, int k) { // x / yの10^kの位の桁
// }
int floor(ld x, ld y) { // 誤差対策TODO
    assert(!equals(y, 0));
    return floor(x / y);
    // floor(x) = ceil(x - 1) という話も
}
int ceil(ld x, ld y) { // 誤差対策TODO // ceil(p/q) = -floor(-(p/q))らしい
    assert(!equals(y, 0));
    return ceil(x / y);
    // ceil(x) = floor(x + 1)
}
int perl(ld x, ld y) { // x = qy + r (0 <= r < y, qは整数) を満たす q
    // 未verify. 誤差対策TODO. EPS外してもいいかも。
    assert(!equals(y, 0));
    if (x >= 0 && y > 0) return floor(x / y)+EPS;
    if (x >= 0 && y < 0) return -floor(x / fabs(y));
    if (x < 0 && y < 0) return floor(x / y) + (x - floor(x/y)*y < -EPS);
    return floor(x / y) - (x - floor(x/y)*y < -EPS); //  (x < 0 && y > 0) 
}
ld modl(ld x, ld y) { // x = qy + r (0 <= r < y, qは整数) を満たす r
    // 未verify. 誤差対策TODO. -0.0が返りうる。
    assert(!equals(y, 0));
    if (x >= 0) return x - fabs(y)*fabs(per(x, y));
    return x - fabs(y)*floor(x, fabs(y));
}
int seisuu(ld x) { return (int)x; } // 整数部分. 誤差対策TODO
int modf(ld x) {
	if (x < 0) return ceill(x);
	else return floorl(x);
}
// 正なら+EPS, 負なら-EPSしてから、文字列に直して小数点以下を捨てる?
int seisuu(int x, int y) {
    assert(y != 0);
    return x / y;
}
int seisuu(ld x, ld y) { // 誤差対策TODO
    assert(!equals(y, 0));
    return (int)(x / y);
}

template <class T> pair<T, T> max(const pair<T, T> &a, const pair<T, T> &b) {
    if (a.first > b.first or a.first == b.first && a.second > b.second) return a;
    return b;
}
template <class T> pair<T, T> min(const pair<T, T> &a, const pair<T, T> &b) {
    if (a.first < b.first or a.first == b.first && a.second < b.second) return a;
    return b;
}

template <class T> bool chmax(T &a, const T& b) {
    if (a < b) { a = b; return true; } return false;
}
template <class T> bool chmin(T &a, const T& b) {
    if (a > b) { a = b; return true; } return false;
}
template <class T> T mid(T a, T b, T c) { // 誤差対策TODO
    return a + b + c - max({a, b, c}) - min({a, b, c});
}
template <class T> void Sort(T &a, T &b, bool rev = false) { 
    if (rev == false) {  // TODO テンプレート引数
        if (a > b) swap(a, b);
    } else {
        if (b > a) swap(b, a);
    }
}
template <class T> void sort(T &a, T &b, T &c, bool rev = false) {
    if (rev == false) { 
        if (a > b) swap(a, b); if (a > c) swap(a, c); if (b > c) swap(b, c);
    } else {
        if (c > b) swap(c, b); if (c > a) swap(c, a); if (b > a) swap(b, a);
    }
}
template <class T> void sort(T &a, T &b, T &c, T &d, bool rev = false) {
    if (rev == false) { 
        if (a > b) swap(a, b); if (a > c) swap(a, c);  if (a > d) swap(a, d);
        if (b > c) swap(b, c); if (b > d) swap(b, d); if (c > d) swap(c, d);
    } else {
        if (d > c) swap(d, c); if (d > b) swap(d, b); if (d > a) swap(d, a);
        if (c > b) swap(c, b); if (c > a) swap(c, a); if (b > a) swap(b, a);
    }
}

struct custom_hash {
    static uint64_t splitmix64(uint64_t x) {
        x += 0x9e3779b97f4a7c15;
        x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9;
        x = (x ^ (x >> 27)) * 0x94d049bb133111eb;
        return x ^ (x >> 31);
    }

    size_t operator()(uint64_t x) const {
        static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count();
        return splitmix64(x + FIXED_RANDOM);
    }
};

class Compress {
public:
    int sz = 0;
    // gp_hash_table<int, int, custom_hash> Z, UZ;
    unordered_map<int, int> Z;    // 元の値 -> 圧縮した値
    unordered_map<int, int> UZ;   // 圧縮した値 -> 元の値

    Compress(const vector<int> &V, int base = 0) {
        this->sz = base;
        set<int> s(V.begin(), V.end());

        for (int x : s) {
            this->Z[x] = this->sz;
            this->UZ[this->sz] = x;
            this->sz++;
        }
    }
    
    Compress(const vector<int> &V1, const vector<int> &V2, int base = 0) {
        this->sz = base;
        vector<int> V3 = V2;
        V3.insert(V3.end(), V1.begin(), V1.end());
        set<int> s(V3.begin(), V3.end());

        for (int x : s) {
            this->Z[x] = this->sz;
            this->UZ[this->sz] = x;
            this->sz++;
        }
    }

    Compress(const vector<int> &V1, const vector<int> &V2, const vector<int> &V3, int base = 0) {
        this->sz = base;
        vector<int> V4 = V1;
        V4.insert(V4.end(), V2.begin(), V2.end());
        V4.insert(V4.end(), V3.begin(), V3.end());
        set<int> s(V4.begin(), V4.end());

        for (int x : s) {
            this->Z[x] = this->sz;
            this->UZ[this->sz] = x;
            this->sz++;
        }
    }

    Compress(const vector<int> &V1, const vector<int> &V2,
            const vector<int> &V3, const vector<int> &V4, int base = 0) {
        this->sz = base;
        vector<int> V5 = V1;
        V5.insert(V5.end(), V2.begin(), V2.end());
        V5.insert(V5.end(), V3.begin(), V3.end());
        V5.insert(V5.end(), V4.begin(), V4.end());
        set<int> s(V5.begin(), V5.end());

        for (int x : s) {
            this->Z[x] = this->sz;
            this->UZ[this->sz] = x;
            this->sz++;
        }
    }

    vector<int> zip(const vector<int> &V) {
        vector<int> ret = V;
        for (int i = 0; i < (int)V.size(); i++) {
            ret[i] = Z[ret[i]];
        }
        return ret;
    }

    vector<int> unzip(const vector<int> &V) {
        vector<int> ret = V;
        for (int i = 0; i < (int)V.size(); i++) {
            ret[i] = UZ[ret[i]];
        }
        return ret;
    }

    int size() { return sz; }

    int encode(int x) { return Z[x]; }
    int decode(int x) {
        if (UZ.find(x) == UZ.end()) return -1; // xが元の配列に存在しないとき
        return UZ[x];
    }
};

class UnionFind {
public:
	UnionFind() = default;
    UnionFind(int N) : par(N), sz(N, 1) {
        iota(par.begin(), par.end(), 0);
    }
	int root(int x) {
		if (par[x] == x) return x;
		return (par[x] = root(par[x]));
	}
	bool unite(int x, int y) {
		int rx = root(x);
		int ry = root(y);
        if (rx == ry) return false;
		if (sz[rx] < sz[ry]) swap(rx, ry);
		sz[rx] += sz[ry];
		par[ry] = rx;
        return true;
	}
	bool issame(int x, int y) { return (root(x) == root(y)); }
	int size(int x) { return sz[root(x)]; }
    vector<vector<int>> groups(int N) {
        vector<vector<int>> G(N);
        for (int x = 0; x < N; x++) {
            G[root(x)].push_back(x);
        }
		G.erase( remove_if(G.begin(), G.end(),
            [&](const vector<int>& V) { return V.empty(); }), G.end());
        return G;
    }
private:
	vector<int> par, sz;
};

template<typename T>
struct BIT {
    int N;             // 要素数
    vector<T> bit[2];  // データの格納先
    BIT(int N_, int x = 0) {
        N = N_ + 1;
        bit[0].assign(N, 0); bit[1].assign(N, 0);
        if (x != 0) {
            for (int i = 0; i < N; i++) add(i, x);
        }
    }
    BIT(const vector<int> &A) {
        N = A.size() + 1;
        bit[0].assign(N, 0); bit[1].assign(N, 0);
        for (int i = 0; i < (int)A.size(); i++) add(i, A[i]);
    }
    void add_sub(int p, int i, T x) {
        while (i < N) { bit[p][i] += x; i += (i & -i); }
    }
    void add(int l, int r, T x) {
        add_sub(0, l + 1, -x * l); add_sub(0, r + 1, x * r);
        add_sub(1, l + 1, x); add_sub(1, r + 1, -x);
    }
    void add(int i, T x) { add(i, i + 1, x); }
    T sum_sub(int p, int i) {
        T ret = 0;
        while (i > 0) { ret += bit[p][i]; i -= (i & -i); }
        return ret;
    }
    T sum(int i) { return sum_sub(0, i) + sum_sub(1, i) * i; }
    T sum(int l, int r) { return sum(r) - sum(l); }
    T get(int i) { return sum(i, i + 1); }
    void set(int i, T x) { T s = get(i); add(i, -s + x); }
};

template<int mod> class Modint {
public:
    int val = 0;
    Modint(int x = 0) { while (x < 0) x += mod; val = x % mod; }
    Modint(const Modint &r) { val = r.val; }

    Modint operator -() { return Modint(-val); } // 単項
    Modint operator +(const Modint &r) { return Modint(*this) += r; }
    Modint operator +(const int &q) { Modint r(q); return Modint(*this) += r; }
    Modint operator -(const Modint &r) { return Modint(*this) -= r; }
    Modint operator -(const int &q) { Modint r(q); return Modint(*this) -= r; }
    Modint operator *(const Modint &r) { return Modint(*this) *= r; }
    Modint operator *(const int &q) { Modint r(q); return Modint(*this) *= r; }
    Modint operator /(const Modint &r) { return Modint(*this) /= r; }
    Modint operator /(const int &q) { Modint r(q); return Modint(*this) /= r; }
    
    Modint& operator ++() { val++; if (val >= mod) val -= mod; return *this; } // 前置
    Modint operator ++(signed) { ++*this; return *this; } // 後置
    Modint& operator --() { val--; if (val < 0) val += mod; return *this; }
    Modint operator --(signed) { --*this; return *this; }
    Modint &operator +=(const Modint &r) { val += r.val; if (val >= mod) val -= mod; return *this; }
    Modint &operator +=(const int &q) { Modint r(q); val += r.val; if (val >= mod) val -= mod; return *this; }
    Modint &operator -=(const Modint &r) { if (val < r.val) val += mod; val -= r.val; return *this; }
    Modint &operator -=(const int &q) { Modint r(q);  if (val < r.val) val += mod; val -= r.val; return *this; }
    Modint &operator *=(const Modint &r) { val = val * r.val % mod; return *this; }
    Modint &operator *=(const int &q) { Modint r(q); val = val * r.val % mod; return *this; }
    Modint &operator /=(const Modint &r) {
        int a = r.val, b = mod, u = 1, v = 0;
        while (b) {int t = a / b; a -= t * b; swap(a, b); u -= t * v; swap(u, v);}
        val = val * u % mod; if (val < 0) val += mod;
        return *this;
    }
    Modint &operator /=(const int &q) {
        Modint r(q); int a = r.val, b = mod, u = 1, v = 0;
        while (b) {int t = a / b; a -= t * b; swap(a, b); u -= t * v; swap(u, v);}
        val = val * u % mod; if (val < 0) val += mod;
        return *this;
    }
    bool operator ==(const Modint& r) { return this -> val == r.val; }
    bool operator <(const Modint& r) { return this -> val < r.val; }
    bool operator >(const Modint& r) { return this -> val > r.val; }
    bool operator !=(const Modint& r) { return this -> val != r.val; }
};

using mint = Modint<MOD>;

istream &operator >>(istream &is, mint& x) {
    int t; is >> t; x = t; return (is);
}
ostream &operator <<(ostream &os, const mint& x) {
    return os << x.val;
}
mint modpow(const mint &x, int n) {
    if (n < 0) return (mint)1 / modpow(x, -n); // 未verify
    assert(n >= 0);
    if (n == 0) return 1;
    mint t = modpow(x, n / 2);
    t = t * t;
    if (n & 1) t = t * x;
    return t;
}

int modpow(__int128_t x, int n, int mod) {
    assert(n >= 0 && mod > 0); // TODO: n <= -1
    __int128_t ret = 1;
    while (n > 0) {
        if (n % 2 == 1) ret = ret * x % mod;
        x = x * x % mod;
        n /= 2;
    }
    return ret;
}
int modinv(__int128_t x, int mod) {
    assert(mod > 0);
    // assert(x > 0);
    if (x == 1 or x == 0) return 1;
    return mod - modinv(mod % x, mod) * (mod / x) % mod;
}

istream &operator >>(istream &is, __int128_t& x) {
    string S; is >> S;
    __int128_t ret = 0;
    int f = 1;
    if (S[0] == '-') f = -1; 
    for (int i = 0; i < S.length(); i++)
        if ('0' <= S[i] && S[i] <= '9')
            ret = ret * 10 + S[i] - '0';
    x = ret * f;
    return (is);
}
ostream &operator <<(ostream &os, __int128_t x) {
    ostream::sentry s(os);
    if (s) {
        __uint128_t tmp = x < 0 ? -x : x;
        char buffer[128]; char *d = end(buffer);
        do {
            --d; *d = "0123456789"[tmp % 10]; tmp /= 10;
        } while (tmp != 0);
        if (x < 0) { --d; *d = '-'; }
        int len = end(buffer) - d;
        if (os.rdbuf()->sputn(d, len) != len) os.setstate(ios_base::badbit);
    }
    return os;
}

__int128_t stoll(string &S) {
    __int128_t ret = 0; int f = 1;
    if (S[0] == '-') f = -1; 
    for (int i = 0; i < S.length(); i++)
        if ('0' <= S[i] && S[i] <= '9') ret = ret * 10 + S[i] - '0';
    return ret * f;
}
__int128_t gcd(__int128_t a, __int128_t b) { return b ? gcd(b, a % b) : a; }
__int128_t lcm(__int128_t a, __int128_t b) {
    return a / gcd(a, b) * b;
    // lcmが__int128_tに収まる必要あり
}

string to_string(ld x, int k) { // xの小数第k位までをstring化する
    assert(k >= 0);
    stringstream ss;
    ss << setprecision(k + 2) << x;
    string s = ss.str();
    if (s.find('.') == string::npos) s += '.';
    int pos = s.find('.');
    for (int i = 0; k >= (int)s.size() - 1 - pos; i++) s += '0';
    s.pop_back();
    if (s.back() == '.') s.pop_back();
    return s;

    // stringstream ss; // 第k+1位を四捨五入して第k位まで返す
    // ss << setprecision(k + 1) << x;
    // string s = ss.str();
    // if (s.find('.') == string::npos) s += '.';
    // int pos = s.find('.');
    // for (int i = 0; k > (int)s.size() - 1 - pos; i++) s += '0';
    // if (s.back() == '.') s.pop_back();
    // return s;
}
string to_string(__int128_t x) {
    string ret = "";
    if (x < 0) { ret += "-"; x *= -1; }
    while (x) { ret += (char)('0' + x % 10); x /= 10; }
    reverse(ret.begin(), ret.end());
    return ret;
}
string to_string(char c) { string s = ""; s += c; return s; }

template<class T> size_t HashCombine(const size_t seed,const T &v) {
    return seed^(hash<T>()(v)+0x9e3779b9+(seed<<6)+(seed>>2));
}
template<class T,class S> struct hash<pair<T,S>>{
    size_t operator()(const pair<T,S> &keyval) const noexcept {
        return HashCombine(hash<T>()(keyval.first), keyval.second);
    }
};
template<class T> struct hash<vector<T>>{
    size_t operator()(const vector<T> &keyval) const noexcept {
        size_t s=0;
        for (auto&& v: keyval) s=HashCombine(s,v);
        return s;
    }
};
template<int N> struct HashTupleCore{
    template<class Tuple> size_t operator()(const Tuple &keyval) const noexcept{
        size_t s=HashTupleCore<N-1>()(keyval);
        return HashCombine(s,get<N-1>(keyval));
    }
};
template <> struct HashTupleCore<0>{
    template<class Tuple> size_t operator()(const Tuple &keyval) const noexcept{ return 0; }
};
template<class... Args> struct hash<tuple<Args...>>{
    size_t operator()(const tuple<Args...> &keyval) const noexcept {
        return HashTupleCore<tuple_size<tuple<Args...>>::value>()(keyval);
    }
};

vector<mint> _fac, _finv, _inv;
void COMinit(int N) {
    _fac.resize(N + 1); _finv.resize(N + 1);  _inv.resize(N + 1);
    _fac[0] = _fac[1] = 1; _finv[0] = _finv[1] = 1; _inv[1] = 1;
    for (int i = 2; i <= N; i++) {
        _fac[i] = _fac[i-1] * mint(i);
        _inv[i] = -_inv[MOD % i] * mint(MOD / i);
        _finv[i] = _finv[i - 1] * _inv[i];
    }
}

mint FAC(int N) {
    if (N < 0) return 0; return _fac[N];
}
mint COM(int N, int K) {
    if (N < K) return 0; if (N < 0 or K < 0) return 0;
    return _fac[N] * _finv[K] * _finv[N - K];
}
mint PERM(int N, int K) {
    if (N < K) return 0; if (N < 0 or K < 0) return 0;
    return _fac[N] *  _finv[N - K];
}
mint NHK(int N, int K) { // initのサイズに注意
    if (N == 0 && K == 0)  return 1;
    return COM(N + K - 1, K);
}

#pragma endregion

vector<string> split_str(string S, char c) {
    vector<string> ret;
    stringstream ss;
    for (int i = 0; i < S.size(); i++) {
        if (S[i] == c) {
            ret.pb(ss.str());
            ss.str("");
            ss.clear();
        } else {
            ss << S[i];
        }
    }
    ret.pb(ss.str());
    return ret;
}

signed main() {
    int H, W;
    while (cin >> H >> W, H) {
        vector<vector<int>> A(H, vector<int>(W));
        for (int i = 0; i < H; i++) {
            for (int j = 0; j < W; j++) {
                char c;
                cin >> c;
                if (c == '.') A[i][j] = 1;
            }
        }

        // A[i][j] := マス(i, j)を含め、そこから上に何マスの白マスが続いているか
        for (int j = 0; j < W; j++) {
            for (int i = 0; i + 1 < H; i++) {
                if (A[i + 1][j] == 1) A[i + 1][j] += A[i][j];
            }
        }

        int Ans = 0;
        for (int i = 0; i < H; i++) {
            vector<pair<int,int>> st; // stackに(A[i][j], i)を積む
            st.pb(0, -1); // 番兵
            vector<int> left(W); // 左方向に伸ばせる限界の位置
            for (int j = 0; j < W; j++) {
                while (st.size() && st.back().first >= A[i][j]) {
                    st.pop_back();
                }
                if (st.size()) left[j] = st.back().second;
                st.pb(A[i][j], j);
            }

            vector<pair<int,int>> st2;
            st2.pb(0, W);
            vector<int> right(W); // 右方向に伸ばせる限界の位置
            for (int j = W - 1; j >= 0; j--) {
                while (st2.size() && st2.back().first >= A[i][j]) {
                    st2.pop_back();
                }
                if (st2.size()) right[j] = st2.back().second;
                st2.pb(A[i][j], j);
            }

            int ans = 0;
            for (int j = 0; j < W; j++) {
                ans = max(ans, A[i][j] * (right[j] - left[j] - 1));
            }

            Ans = max(Ans, ans);
        }
        cout << Ans << endl;
    }
}
accuracy: 1
time_ms: 10
memory_kb: 5200
final_rank: 11

code_snippet: #pragma GCC optimize("Ofast")
#include <bits/stdc++.h>
using namespace std;
typedef long long int ll;
typedef unsigned long long int ull;

mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
ll myRand(ll B) { return (ull)rng() % B; }

template<typename T> T largest_rectangle_in_histogram(const vector<T> &a) {
    int n = a.size();
    using P = pair<T, int>;
    vector<int> l(n, 0), r(n, 0);
    vector<P> st_l, st_r;
    st_l.push_back({-1, -1});
    for (int i = 0; i < n; ++i) {
        while (a[i] <= st_l.back().first) {
            st_l.pop_back();
        }
        l[i] = st_l.back().second;
        st_l.push_back({a[i], i});
    }
    st_r.push_back({-1, n});
    for (int i = n - 1; i >= 0; --i) {
        while (a[i] <= st_r.back().first) {
            st_r.pop_back();
        }
        r[i] = st_r.back().second;
        st_r.push_back({a[i], i});
    }
    T res = 0;
    for (int i = 0; i < n; ++i) {
        res = max(res, a[i] * (r[i] - l[i] - 1));
    }
    return res;
}

template<typename T> T largest_square_in_histogram(const vector<T> &a) {
    int n = a.size();
    using P = pair<T, int>;
    vector<int> l(n, 0), r(n, 0);
    vector<P> st_l, st_r;
    st_l.push_back({-1, -1});
    for (int i = 0; i < n; ++i) {
        while (a[i] <= st_l.back().first) {
            st_l.pop_back();
        }
        l[i] = st_l.back().second;
        st_l.push_back({a[i], i});
    }
    st_r.push_back({-1, n});
    for (int i = n - 1; i >= 0; --i) {
        while (a[i] <= st_r.back().first) {
            st_r.pop_back();
        }
        r[i] = st_r.back().second;
        st_r.push_back({a[i], i});
    }
    T res = 0;
    for (int i = 0; i < n; ++i) {
        T len = min(a[i], T(r[i] - l[i] - 1));
        res = max(res, len * len);
    }
    return res;
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int h,w;
    while (cin >> h >> w, h+w) {
        vector<vector<char>> c(h, vector<char>(w));
        for (int i = 0; i < h; ++i) {
            for (int j = 0; j < w; ++j) {
                cin >> c[i][j];
            }
        }
        int res = 0;
        vector<int> a(w);
        for (int i = 0; i < h; ++i) {
            for (int j = 0; j < w; ++j) {
                if (a[j] == 0) {
                    int ii = i;
                    while (ii < h and c[ii][j] == '.') {
                        ii += 1, a[j] += 1;
                    }
                } else {
                    a[j] -= 1;
                }
            }
            res = max(res, largest_rectangle_in_histogram(a));
        }
        cout << res << endl;
    }
}
accuracy: 1
time_ms: 10
memory_kb: 3464
final_rank: 1

code_snippet: #line 1 "test/aoj_0116.test.cpp"
#define PROBLEM "https://onlinejudge.u-aizu.ac.jp/problems/0116"

#ifndef call_include
#define call_include
#include <bits/stdc++.h>
using namespace std;
#endif

#line 2 "dp/largest_histogram_rectangle.cpp"

#ifndef call_include
#define call_include
#include <bits/stdc++.h>
using namespace std;
#endif

template<typename T> long long largest_histogram_rectangle(vector<T> &v) {
	v.emplace_back(0);
	int n = v.size();
	long long mxrec = 0;
	stack<pair<T, int>> s;
	for(int i = 0; i < n; i++) {
		int ni = i;
		while(s.size() > 0 and s.top().first > v[i]) {
			long long h = s.top().first, w = i - s.top().second;
			long long rec = h * w;
			if(mxrec < rec) mxrec = rec;
			ni = s.top().second;
			s.pop();
		}
		s.emplace(v[i], ni);
	}
	return mxrec;
}
#line 2 "structure/2d_array.cpp"

#ifndef call_include
#define call_include
#include <bits/stdc++.h>
using namespace std;
#endif

template<typename T> struct v2d {
private:
	vector<vector<T>> m;

public:
	v2d() {}
	v2d(int h, int w) : m(h, vector<T>(w)) {}
	v2d(int h, int w, const T &init) : m(h, vector<T>(w, init)) {}
	v2d(const initializer_list<initializer_list<T>> m_init) : m(m_init.begin(), m_init.end()) {}

	void assign(int h, int w) {
		m.assign(h, vector<T>(w));
	}
	void assign(int h, int w, const T init) {
		m.assign(h, vector<T>(w, init));
	}

	inline int size() const {
		return m.size();
	}

	void in() {
		for(vector<T> &v : m)
			for(T &val : v) cin >> val;
	}
	void in(int h, int w) {
		m.resize(h, vector<T>(w));
		in();
	}

	void out() {
		int h = m.size();
		for(vector<T> &v : m) {
			int sz = v.size();
			for(int j = 0; j < sz; j++) {
				cout << v[j] << (j == sz - 1 ? '\n' : ' ');
			}
		}
		cout << flush;
	}

	inline vector<T> &operator[](int idx) {
		assert(0 <= idx && idx < m.size());
		return m[idx];
	}

	bool rangeout(int x, int y) {
		if(x < 0 || y < 0 || x >= size() || y >= m[x].size()) return true;
		return false;
	}
};
#line 11 "test/aoj_0116.test.cpp"

int main() {
	while(1) {
		int H, W;
		cin >> H >> W;
		if(H == 0 and W == 0) return 0;

		v2d<char> m;
		m.in(H, W);
		v2d<int> hist(H, W);

		for(int j = 0; j < W; j++) {
			hist[0][j] = (m[0][j] == '*' ? 0 : 1);
		}
		for(int i = 1; i < H; i++) {
			for(int j = 0; j < W; j++) {
				hist[i][j] = (m[i][j] == '*' ? 0 : hist[i - 1][j] + 1);
			}
		}

		int ans = 0;
		for(int i = 0; i < H; i++) ans = max(ans, (int)largest_histogram_rectangle(hist[i]));
		cout << ans << endl;
	}
}
accuracy: 1
time_ms: 10
memory_kb: 5440
final_rank: 12

code_snippet: #include<bits/stdc++.h>
using namespace std;
int main(){
    int h, w, sum, piri[550][550], l[550], r[550];
    string masu;
    while(cin >> h >> w, h!=0){
        memset(piri, 0, sizeof(piri));
        sum = 0;
        for(int i=0;i<h;i++){
            cin >> masu;
            piri[i][0] = 0;
            piri[i][w+1] = 0;
            for(int j=0;j<w;j++){
                if(masu[j]=='.'){
                    piri[i][j+1] = 1;
                    if(i > 0) piri[i][j+1] += piri[i-1][j+1];
                }
            }
            stack<int> st;
            for(int j=0;j<=w;j++){
                while(!st.empty() && piri[i][st.top()]>=piri[i][j]) st.pop();
                if(st.empty()) l[j] = j;
                else l[j] = st.top()+1;
                st.push(j);
            }
            st = stack<int>();
            for(int j=w+1;j>0;j--){
                while(!st.empty() && piri[i][st.top()]>=piri[i][j]) st.pop();
                if(st.empty()) r[j] = j;
                else r[j] = st.top()-1;
                st.push(j);
                sum = max(sum, piri[i][j]*(r[j]-l[j]+1));
            }
        }
        cout << sum << endl;
    }
    return(0);
}
accuracy: 1
time_ms: 20
memory_kb: 4260
final_rank: 5

code_snippet: #include <stdio.h>
#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000000


int main(){
	int h,w;
	while(cin>>h>>w){
		if(h==0)break;
		
		vector<string> s(h);
		rep(i,h)cin>>s[i];
		
		int ans = 0;
		
		vector<int> v(w,0);
		rep(i,h){
			rep(j,w){
				if(s[i][j]=='*')v[j]=0;
				else v[j]++;
			}
			stack<pair<int,int>> S;
			S.emplace(0,-1);
			rep(j,w+1){
				int x;
				if(j<w)x =v[j];
				else x = 0;
				if(S.top().first < x){
					S.emplace(x,j);
				}
				else{
					int last;
					while(S.size()>0 && S.top().first>=x){
						ans = max(ans,S.top().first * (j-S.top().second));
						last = S.top().second;
						S.pop();
					}
					S.emplace(x,last);
				}
				
				
			}
		}
		
		cout<<ans<<endl;		
	}
    return 0;
}
accuracy: 1
time_ms: 10
memory_kb: 3464
final_rank: 1

code_snippet: #include<bits/stdc++.h>
using namespace std;
using ll = long long;

int main(){
    int h, w;
    while(cin >> h >> w, h){
        vector<string> A(h);
        vector<vector<int>> dp(h, vector<int>(w));
        int ans = 0;
        stack<pair<int,int>> ST;
        ST.emplace(-1, -1);
        for(int y = 0; y < h; y++){
            cin >> A[y];
            for(int x = 0; x < w; x++){
                if(A[y][x] == '.')dp[y][x] = (y == 0 ? 0 : dp[y - 1][x]) + 1;
                int p = x;
                while(ST.top().first > dp[y][x]){
                    ans = max(ans, ST.top().first * (x - ST.top().second));
                    p = ST.top().second;
                    ST.pop();
                }
                ST.emplace(dp[y][x], p);
            }
            while(ST.size() > 1){
                ans = max(ans, ST.top().first * (w - ST.top().second));
                ST.pop();
            }
        }
        cout << ans << '\n';
    }
}
accuracy: 1
time_ms: 10
memory_kb: 4528
final_rank: 8

code_snippet: #define _USE_MATH_DEFINES
#include<iostream>
#include<string>
#include<vector>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<list>
#include<algorithm>
#include<numeric>
#include<cmath>
#include<iomanip>
#include<cctype>
#include<sstream>
#include<regex>
#include<bitset>
#include<random>
using namespace std;
#define int long long
#define uint unsigned long long
const int mod=1e9+7;
const int Mod=998244353;
const int inf=mod*mod;
const double eps=1e-9;
string s[500];
int dp[500][500],l[500],r[500];
signed main(){
	int h,w;
	while(cin>>h>>w,h){
		for(int i=0;i<h;i++)
			cin>>s[i];
		for(int i=0;i<h;i++){
			for(int j=0;j<w;j++){
				dp[i][j]=0;
				if(s[i][j]=='*')
					continue;
				if(i!=0)
					dp[i][j]=dp[i-1][j];
				dp[i][j]++;
			}
		}
		int ans=0;
		for(int i=0;i<h;i++){
			stack<int>st;
			for(int j=0;j<w;j++){
				while(!st.empty()&&dp[i][st.top()]>=dp[i][j])
					st.pop();
				l[j]=st.empty()?0:st.top()+1;
				st.push(j);
			}
			st={};
			for(int j=w-1;j>=0;j--){
				while(!st.empty()&&dp[i][st.top()]>=dp[i][j])
					st.pop();
				r[j]=st.empty()?w:st.top();
				st.push(j);
			}
			for(int j=0;j<w;j++)
				ans=max(ans,dp[i][j]*(r[j]-l[j]));
		}
		cout<<ans<<endl;
	}
}
accuracy: 1
time_ms: 10
memory_kb: 5660
final_rank: 13

code_snippet: #include<bits/stdc++.h>
#define rep(i,n) for(int i=0; i< (n); i++)
using namespace std;
typedef long long ll;
const int mod = 1000000007;
const int inf = (1<<30);
int w,h;
struct p{
	int hei;
	int wid;
};

int dp(int tebl[]){
	// rep(i,w+1) cout<<tebl[i]<<" ";
	int ans=0; 
	stack<p> s;
	p k;
	k.hei=tebl[0];
	k.wid=0;
	s.push(k);
	//cout<<s.top().hei<<" "<<s.top().wid<<" "<<s.size()<<endl;
	for(int i=1; i<w+1; i++){
		if(s.top().hei < tebl[i]){
			//cout<<i<<" s.top().hei < tebl[i]"<<endl;
			k.hei=tebl[i]; k.wid=i;
			s.push(k);
		}else if(s.top().hei>tebl[i]){
			int g=i;
			while(s.size() && s.top().hei>tebl[i] ){
				p r=s.top();
				s.pop();
				ans=max(ans,(r.hei*(i-r.wid)));
				g=min(g,r.wid);
			}
			k.hei=tebl[i]; k.wid=g;
			s.push(k);
		}
		//cout<<s.top().hei<<" "<<s.top().wid<<" "<<s.size()<<endl;
	}
	return ans;
}
			


int main(){
	while(true){
		cin>>h>>w;
		if(h==0 && w==0) break;
		int tebl[h][w+1]={};
		rep(i,h){
			rep(j,w){
				char k; cin>>k;
				if(k=='.'){
					tebl[i][j]=1;
				}
			}
		}
		rep(j,w){
			int a=0;
			rep(i,h){
				if(tebl[i][j]==1){
					tebl[i][j]=++a;
				}else{
					a=0;
					tebl[i][j]=a;
				}
			}
		}
		/*rep(i,h){
			rep(j,w){
				cout<<tebl[i][j]<<" ";
			}
			cout<<endl;
		}*/
		int ans=0;
		rep(i,h){
			ans=max(ans,dp(tebl[i]));
			//cout<<ans<<endl;
		}
		cout<<ans<<endl;
	}
}
accuracy: 1
time_ms: 20
memory_kb: 4072
final_rank: 3

code_snippet: #include <iostream>
#include <vector>
#include <algorithm>
#include <utility>
using namespace std;

int solve_sub(vector<int> &a){
    int n = a.size();
    vector<pair<int,int>> wait;
    int res = 0;
    for(int i=0; i<n; i++){
        int lm = -1;
        while(!wait.empty() and a[i]<wait.back().second){
            res = max(res, (i-wait.back().first)*wait.back().second);
            lm = wait.back().first;
            wait.pop_back();
        }
        if(lm != -1){
            wait.emplace_back(lm, a[i]);
        }
        if(wait.empty() or a[i]>wait.back().second){
            wait.emplace_back(i, a[i]);
        }
    }
    return res;
}
int solve(vector<vector<int>> a){
    int n = a.size();
    int m = a[0].size();
    for(int j=0; j<m; j++){
        for(int i=0; i<n-1; i++){
            if(a[i+1][j] == 0) continue;
            a[i+1][j] += a[i][j];
        }
    }
    int res = 0;
    for(int i=0; i<n; i++){
        a[i].push_back(0);
        res = max(res, solve_sub(a[i]));
    }
    return res;
}

int main(){
    while(1){
        int h,w;
        cin >> h >> w;
        if(h == 0) break;

        vector<vector<int>> a(h, vector<int>(w));
        for(int i=0; i<h; i++){
            for(int j=0; j<w; j++){
                char c;
                cin >> c;
                if(c == '.') a[i][j] = 1;
                else a[i][j] = 0;
            }
        }
        cout << solve(a) << endl;
    }
}
accuracy: 1
time_ms: 20
memory_kb: 5756
final_rank: 14

code_snippet: #include <bits/stdc++.h>
using namespace std;
#define MAX 1400

struct Rectangle{
    // indexの地点の高さ
    int height;
    // 始点のX
    int pos;
};

// buf = Tの一行の配列
// size = Tの一行のサイズ
int getLargestRectangle(int w,int buf[]){
    // stackとして長方形を保存
    stack<Rectangle> s;
    int maxV = 0;
    // 配列の最後を0にすることで、範囲外後に最大面積を計算するようにする
    buf[w] = 0;
    // stoperまで回す
    for(int i = 0;i <= w;++i){
        Rectangle rect;
        // 高さ代入
        rect.height = buf[i];
        // 横の位置は計算を始める時点の位置
        rect.pos = i;
        // stack空ならば追加
        if(s.empty()){
            s.push(rect);
        }
        else{
            // stackのトップにあるrectの高さより高い場合
            if(s.top().height < rect.height){
                s.push(rect);
            }
            // stackのトップにあるrectの高さよりも低い場合(保存してあった最大の高さが無効になるので、無効になる高さの面積を計算する必要がある)
            else if(s.top().height > rect.height){
                // 新しく積むrectの幅を出すために保存しておく。
                int targetX = i;
                // スタックに入っているrectが対象のrectの高さ以上の場合最大面積を計算する
                while(!s.empty() && s.top().height >= rect.height){
                    Rectangle pre = s.top();
                    s.pop();
                    // 面積を計算する
                    int area = pre.height * (i - pre.pos);
                    maxV = max(maxV,area);
                    // 始点Xを入れ替える(stackに積むrectの高さは取り出したrectの高さよりも低いので、取り出した地点が、長方形を計算する際のX座標になる)
                    targetX = pre.pos;
                }
                rect.pos = targetX;
                s.push(rect);
            }
        }
    }
    return maxV;
}

int H,W;
// タイル情報を2次元で保存0:綺麗なタイル,1:汚れたタイル
char tiles[MAX][MAX];
// 高さを格納する2次元配列
int T[MAX][MAX];

int getLargestRectangle(){
    // 各々の地点における高さを格納するプロセス
    for(int i = 0;i < H;++i){
        for(int j = 0;j < W;++j){
            // 汚れたタイル
            if(tiles[i][j] == '*'){
                T[i][j] = 0;
            }
            else{
                // 対象のindex地点の高さを格納する
                // 上の地点 + 1の高さになる
                T[i][j] = (i > 0) ? T[i - 1][j] + 1 : 1;    
            }
        }
    }
    
    // 最大面積保存用
    int maxV = 0;

    for(int i = 0;i < H;++i){
        maxV = max(maxV,getLargestRectangle(W,T[i]));
    }
    return maxV;
}

int main(){
    while(1){
        scanf("%d %d",&H,&W);
        if(H == 0 && W == 0) break;
        
        for(int i = 0;i < H;++i){
            for(int j = 0;j < W;++j){
                scanf(" %c",&tiles[i][j]);
            }
        }
        printf("%d\n",getLargestRectangle());
    }
}
accuracy: 1
time_ms: 20
memory_kb: 6600
final_rank: 16

code_snippet: #include <bits/stdc++.h>
using namespace std;
#define MAX 1400

struct Rectangle{
    // indexの地点の高さ
    int height;
    // 始点のX
    int pos;
};

// buf = Tの一行の配列
// size = Tの一行のサイズ
int getLargestRectangle(int w,int buf[]){
    // stackとして長方形を保存
    stack<Rectangle> s;
    int maxV = 0;
    // 配列の最後を0にすることで、範囲外後に最大面積を計算するようにする
    buf[w] = 0;
    // stoperまで回す
    for(int i = 0;i <= w;++i){
        Rectangle rect;
        // 高さ代入
        rect.height = buf[i];
        // 横の位置は計算を始める時点の位置
        rect.pos = i;
        // stack空ならば追加
        if(s.empty()){
            s.push(rect);
        }
        else{
            // stackのトップにあるrectの高さより高い場合
            if(s.top().height < rect.height){
                s.push(rect);
            }
            // stackのトップにあるrectの高さよりも低い場合(保存してあった最大の高さが無効になるので、無効になる高さの面積を計算する必要がある)
            else if(s.top().height > rect.height){
                // 新しく積むrectの幅を出すために保存しておく。
                int targetX = i;
                // スタックに入っているrectが対象のrectの高さ以上の場合最大面積を計算する
                while(!s.empty() && s.top().height >= rect.height){
                    Rectangle pre = s.top();
                    s.pop();
                    // 面積を計算する
                    int area = pre.height * (i - pre.pos);
                    maxV = max(maxV,area);
                    // 始点Xを入れ替える(stackに積むrectの高さは取り出したrectの高さよりも低いので、取り出した地点が、長方形を計算する際のX座標になる)
                    targetX = pre.pos;
                }
                rect.pos = targetX;
                s.push(rect);
            }
        }
    }
    return maxV;
}

int H,W;
// タイル情報を2次元で保存0:綺麗なタイル,1:汚れたタイル
char tiles[MAX][MAX];
// 高さを格納する2次元配列
int T[MAX][MAX];

int getLargestRectangle(){
    // 各々の地点における高さを格納するプロセス
    for(int i = 0;i < H;++i){
        for(int j = 0;j < W;++j){
            // 汚れたタイル
            if(tiles[i][j] == '*'){
                T[i][j] = 0;
            }
            else{
                // 対象のindex地点の高さを格納する
                // 上の地点 + 1の高さになる
                T[i][j] = (i > 0) ? T[i - 1][j] + 1 : 1;    
            }
        }
    }
    
    // for(int i = 0;i < H;++i){
    //     for(int j = 0;j < W;++j){
    //         printf("%d",T[i][j]);
    //     }
    //     printf("\n");
    // }    
    // 最大面積保存用
    int maxV = 0;

    for(int i = 0;i < H;++i){
        maxV = max(maxV,getLargestRectangle(W,T[i]));
    }
    return maxV;
}

int main(){
    while(1){
        scanf("%d %d",&H,&W);
        if(H == 0 && W == 0) break;
        
        for(int i = 0;i < H;++i){
            for(int j = 0;j < W;++j){
                scanf(" %c",&tiles[i][j]);
                //printf("%c",tiles[i][j]);
            }
            //printf("\n");
        }
        printf("%d\n",getLargestRectangle());
    }
}
accuracy: 1
time_ms: 10
memory_kb: 6596
final_rank: 15

code_snippet: #include <iostream>	

struct Rect {
	int width;
	int height;
	int max_width;
	int max_height;
	Rect() : width(0), height(0), max_width(0), max_height(0) {}
};

void show(Rect** DP, int H, int W)
{
	for (int i=0; i<H; i++) {
		for (int j=0; j<W; j++) std::cout << DP[i][j].width << "," << DP[i][j].height << " ";
		std::cout << std::endl;
	}
}

int main()
{
	int W, H;
	while (std::cin >> H >> W) {
		
		if (W==0) break;
	
		char** board = new char*[H];
		for (int i=0; i<H; i++) board[i] = new char[W];
		for (int i=0; i<H; i++) {
			for (int j=0; j<W; j++) std::cin >> board[i][j];
		}
		
		Rect** DP = new Rect*[H];
		for (int i=0; i<H; i++) DP[i] = new Rect[W];
		
		int max = 0;
		for (int i=0; i<H; i++) {
			for (int j=0; j<W; j++) {
				if (board[i][j] == '*') continue;
				if (i>0) {
					if (j>0) {
						DP[i][j].max_width = DP[i][j-1].max_width + 1;
						DP[i][j].max_height = DP[i-1][j].max_height + 1;
						int w = DP[i][j].max_width;
						int max_area = DP[i][j].max_height; // DP[i][j].max_height * 1
						int max_w = 1;
						int max_h = DP[i][j].max_height;
						for (int h=2; h<=DP[i][j].max_height; h++) {
							w = std::min(w, DP[i-h+1][j].max_width);
							if (max_area < h*w) {
								max_area = h*w;
								max_w = w;
								max_h = h;
							}
						}
						DP[i][j].width = max_w;
						DP[i][j].height = max_h;
					}
					else { // j==0
						DP[i][j].height = DP[i-1][j].height + 1;
						DP[i][j].max_height = DP[i-1][j].max_height + 1;
						DP[i][j].width = 1;
						DP[i][j].max_width = 1;
					}
				}
				else { // i==0
					DP[i][j].height = 1;
					DP[i][j].max_height = 1;
					if (j>0) {
						DP[i][j].width = DP[i][j-1].width + 1;
						DP[i][j].max_width = DP[i][j-1].max_width + 1;
					} else {
						DP[i][j].width = 1;
						DP[i][j].max_width = 1;
					}
				}
				
				if (max < DP[i][j].width * DP[i][j].height) {
					max = DP[i][j].width * DP[i][j].height;
				}
			}
		}
		//show(DP, H, W);
		std::cout << max << std::endl;
		
		for (int i=0; i<H; i++) delete[] board[i];
		delete[] board;
		for (int i=0; i<H; i++) delete[] DP[i];
		delete[] DP;
	}
	
	return 0;
}
accuracy: 1
time_ms: 50
memory_kb: 6952
final_rank: 17

code_snippet: #include <iostream>
#include <algorithm>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <numeric>
#include <bitset>
#include <cmath>

static const int MOD = 1000000007;
using ll = long long;
using u32 = unsigned;
using u64 = unsigned long long;
using namespace std;

template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208;

int main() {
    int h, w;
    while(cin >> h >> w, h){
        vector<vector<int>> v(h, vector<int>(w, 0));
        for (int i = 0; i < h; ++i) {
            string s;
            cin >> s;
            for (int j = 0; j < w; ++j) {
                v[i][j] = (s[j] == '.');
            }
            if(i){
                for (int j = 0; j < w; ++j) {
                    if(v[i][j]) v[i][j] += v[i-1][j];
                }
            }
        }

        int ans = 0;
        for (int i = 0; i < h; ++i) {
            stack<int> s;
            v[i].emplace_back(0);
            vector<int> l(w+1);
            for (int j = 0; j <= w; ++j) {
                while(!s.empty() && v[i][j] <= v[i][s.top()]){
                    int k = s.top(); s.pop();
                    ans = max(ans, (j-l[k]-1)*v[i][k]);
                }
                l[j] = s.empty() ? -1 : s.top();
                s.emplace(j);
            }
        }
        cout << ans << "\n";
    }
    return 0;
}
accuracy: 1
time_ms: 10
memory_kb: 4844
final_rank: 10

code_snippet: #include<iostream>
#include<algorithm>
#include <map>
#include <vector>
#include <stack>
#include<string>

using namespace std;

char cs[501][501];
int p[501][501][2];

int main()
{
	int H, W;
	while (true)
	{
		cin >> H >> W;
		if (H == 0 && W == 0)
			break;
		for (size_t i = 0; i <= H; i++)
		{
			p[i][0][0] = 0;
			p[i][0][1] = 0;
		}
		for (size_t j = 0; j <= W; j++)
		{
			p[0][j][0] = 0;
			p[0][j][1] = 0;
		}

		int max = 0;
		for (size_t i = 1; i <= H; i++)
		{
			for (size_t j = 1; j <= W; j++)
			{
				cin >> cs[i][j];
				if (cs[i][j] == '*')
				{
					p[i][j][0] = 0;
					p[i][j][1] = 0;
					continue;
				}
				else
				{
					p[i][j][0] = p[i - 1][j][0] + 1;
					p[i][j][1] = p[i][j - 1][1] + 1;
					if (p[i][j][0] * p[i][j][1] <= max)
						continue;
					int hei = 501;
					for (int c = j ; c >= 1 ; c--)
					{
						hei = min(hei, p[i][c][0]);
						int area = hei * (j - c + 1);
						if (area > max)
							max = area;
					}
					int wid = 501;
					for (int r = i; r >= 1; r--)
					{
						wid = min(wid, p[r][j][1]);
						int area = wid * (i - r + 1);
						if (area > max)
							max = area;
					}
				}
			}
		}
		cout << max << endl;
	}
	//system("pause");
}
accuracy: 1
time_ms: 230
memory_kb: 5252
final_rank: 18

code_snippet: #include<iostream>
#include<algorithm>
#include <map>
#include <vector>
#include <stack>
#include<string>

using namespace std;

char cs[501][501];
int p[501][501][2];

int main()
{
	int H, W;
	while (true)
	{
		cin >> H >> W;
		if (H == 0 && W == 0)
			break;
		for (size_t i = 0; i <= H; i++)
		{
			p[i][0][0] = 0;
			p[i][0][1] = 0;
		}
		for (size_t j = 0; j <= W; j++)
		{
			p[0][j][0] = 0;
			p[0][j][1] = 0;
		}

		int max = 0;
		for (size_t i = 1; i <= H; i++)
		{
			for (size_t j = 1; j <= W; j++)
			{
				cin >> cs[i][j];
				if (cs[i][j] == '*')
				{
					p[i][j][0] = 0;
					p[i][j][1] = 0;
					continue;
				}
				else
				{
					p[i][j][0] = p[i - 1][j][0] + 1;
					p[i][j][1] = p[i][j - 1][1] + 1;
					int hei = 501;
					for (int c = j ; c >= 1 ; c--)
					{
						hei = min(hei, p[i][c][0]);
						int area = hei * (j - c + 1);
						if (area > max)
							max = area;
					}
					int wid = 501;
					for (int r = i; r >= 1; r--)
					{
						wid = min(wid, p[r][j][1]);
						int area = wid * (i - r + 1);
						if (area > max)
							max = area;
					}
				}
			}
		}
		cout << max << endl;
	}
	//system("pause");
}
accuracy: 1
time_ms: 390
memory_kb: 5208
final_rank: 19

code_snippet: #include <iostream>
#include <vector>
#include <algorithm>
using namespace std;

int rectangularSearching(short** table,int W,int H);

struct rect{
  unsigned short w;
  unsigned short h;
  bool operator<(const rect &another) const{
    return w*20000+h<another.w*20000+another.h;
  };
  bool operator==(const rect &another) const{
    return w==another.w and h==another.h;
  };
};


int main(void){
  int W,H;
  string s;
  while(1){
    cin>>H>>W;
    if(H==0 and W==0)break;
    short** table;
    table=(short **)malloc(sizeof(short*)*(H+2));
    for(int i=0;i<H+2;i++){
      table[i]=(short*)malloc(sizeof(short)*(W+2));
    }
    for(int i=0;i<H+2;i++){
      table[i][0]=0;
      table[i][W+1]=0;
    }
    for(int j=0;j<W+2;j++){
      table[0][j]=0;
      table[H-1][j]=0;
    }

    for(int i=0;i<H;i++){
      cin>>s;
      for(int j=0;j<W;j++){
        if(s[j]=='*'){
          table[i+1][j+1]=0;
        }
        else{
          table[i+1][j+1]=1;
        }
      }
    }
    cout<<rectangularSearching(table,W,H)<<endl;
    for(int i=0;i<H+2;i++){
      free(table[i]);
    }
    free(table);
  }
  return 0;
}

int rectangularSearching(short** table,int W,int H){
  vector<rect> dp[H][W];
  int retval=0;
  for(int i=0;i<H;i++){
    for(int j=0;j<W;j++){
      if(table[i+1][j+1]==0){
        continue;
      }

      if(table[i+1][j]==0 and table[i][j+1]==0){
        rect r;
        r.h=1;
        r.w=1;
        dp[i][j].push_back(r);
      }
      else if(table[i+1][j]==0 and table[i][j+1]==1){
        unsigned short height=0;
        while(table[i+1-height][j+1]==1){
          height++;
        }
        rect r;
        r.h=height;
        r.w=1;
        dp[i][j].push_back(r);
      }
      else if(table[i+1][j]==1 and table[i][j+1]==0){
        unsigned short width=0;
        while(table[i+1][j+1-width]==1){
          width++;
        }
        rect r;
        r.h=1;
        r.w=width;
        dp[i][j].push_back(r);
      }
      else{
        unsigned short width=0;
        while(table[i+1][j-width+1]==1){
          width++;
        }
        unsigned short height=0;
        while(table[i-height+1][j+1]==1){
          height++;
        }

        vector<rect> tmp;
        tmp.clear();
        tmp.shrink_to_fit();
        for(unsigned short k=0;k<dp[i][j-1].size();k++){
          rect r;
          unsigned short h=min(dp[i][j-1][k].h,height);
          unsigned short w=dp[i][j-1][k].w+1;
          r.h=h;
          r.w=w;
          tmp.push_back(r);
        }
        for(int k=0;k<dp[i-1][j].size();k++){
          rect r;
          unsigned short h=dp[i-1][j][k].h+1;
          unsigned short w=min(dp[i-1][j][k].w,width);
          r.h=h;
          r.w=w;
          tmp.push_back(r);
        }
        /*
        rect r,r2;
        for(int k=0;k<tmp.size();k++){
          r=tmp[k];
          bool contain=true;
          for(int l=0;l<tmp.size();l++){
            if(k==l)continue;
            r2=tmp[l];
            if(r.w<r2.w and r.h<r2.h){
              contain=false;
              break;
            }
          }
          if(contain){
            dp[i][j].push_back(r);
          }
        }
        */
        sort(tmp.begin(),tmp.end());
        tmp.erase(unique(tmp.begin(),tmp.end()),tmp.end());

        dp[i][j]=tmp;
      }
      vector<rect> v=dp[i][j];
      for(int k=0;k<v.size();k++){
        retval=max((int)v[k].h*(int)v[k].w,retval);
      }
      if(i>=1){
        dp[i-1][j].clear();
        dp[i-1][j].shrink_to_fit();
      }
    }
  }

  return retval;
}
accuracy: 1
time_ms: 500
memory_kb: 9952
final_rank: 20

code_snippet: #include <bits/stdc++.h>

using namespace std;

int H, W;
int DpOfRow[500][500];

int calculation() {
	int len, ans;

	ans = 0;
	for (int i = 0; i < H; i++) {
		for (int j = W - 1; j > -1; j--) {
			if (DpOfRow[i][j] == 0) {
				continue;
			}
			len = INT_MAX;
			for (int k = i; k < H; k++) {
				len = min(len, DpOfRow[k][j]);
				if (len == 0) {
					break;
				}
				ans = max(ans, len*(k - i + 1));
			}
		}
	}
	
	return ans;
}


int main() {
	char c;

	while (cin >> H >> W, H) {
		memset(DpOfRow, 0, sizeof(DpOfRow));
		for (int i = 0; i < H; i++) {
			for (int j = 0; j < W; j++) {
				cin >> c;
				if (c == '.') {
					DpOfRow[i][j] = (j == 0 ? 1 : DpOfRow[i][j - 1] + 1);
				}
			}
		}
		cout << calculation() << endl;
	}

	return 0;
}
accuracy: 1
time_ms: 50
memory_kb: 3980
final_rank: 7

code_snippet: #include <bits/stdc++.h>

using namespace std;

int H, W;
int DpOfRow[500][500];
char MapData[11][11];

void show() {
	cout << endl;
	for (int i = 0; i < H; i++) {
		for (int j = 0; j < W; j++) {
			cout << MapData[i][j];
		}
		cout << endl;
	}
}

int calculation() {
	int len, ans;

	ans = 0;
	for (int i = 0; i < H; i++) {
		for (int j = W - 1; j > -1; j--) {
			if (DpOfRow[i][j] == 0) {
				continue;
			}
			len = INT_MAX;
			for (int k = i; k < H; k++) {
				len = min(len, DpOfRow[k][j]);
				if (len == 0) {
					break;
				}
				ans = max(ans, len*(k - i + 1));
			}
		}
	}
	
	return ans;
}


int main() {
	char c;

	while (cin >> H >> W, H) {
		memset(DpOfRow, 0, sizeof(DpOfRow));
		for (int i = 0; i < H; i++) {
			for (int j = 0; j < W; j++) {
				cin >> c;
				if (c == '.') {
					DpOfRow[i][j] = (j == 0 ? 1 : DpOfRow[i][j - 1] + 1);
				}
				//MapData[i][j] = c;
			}
		}
		//show();
		cout << calculation() << endl;
	}

	return 0;
}
accuracy: 1
time_ms: 50
memory_kb: 4064
final_rank: 9

code_snippet: #include <bits/stdc++.h>

using namespace std;

int H, W;
int DpOfRow[500][500];
char MapData[11][11];

void show() {
	cout << endl;
	for (int i = 0; i < H; i++) {
		for (int j = 0; j < W; j++) {
			cout << MapData[i][j];
		}
		cout << endl;
	}
}

int calculation() {
	int len, ans;

	ans = 0;
	for (int i = 0; i < H; i++) {
		for (int j = W - 1; j > -1; j--) {
			if (DpOfRow[i][j] == 0) {
				continue;
			}
			len = INT_MAX;
			for (int k = i; k < W; k++) {
				len = min(len, DpOfRow[k][j]);
				if (len == 0) {
					break;
				}
				ans = max(ans, len*(k - i + 1));
			}
		}
	}
	
	return ans;
}


int main() {
	char c;

	while (cin >> H >> W, H) {
		memset(DpOfRow, 0, sizeof(DpOfRow));
		for (int i = 0; i < H; i++) {
			for (int j = 0; j < W; j++) {
				cin >> c;
				if (c == '.') {
					DpOfRow[i][j] = (j == 0 ? 1 : DpOfRow[i][j - 1] + 1);
				}
				//MapData[i][j] = c;
			}
		}
		//show();
		cout << calculation() << endl;
	}

	return 0;
}
accuracy: 1
time_ms: 50
memory_kb: 3976
final_rank: 6

code_snippet: #define _USE_MATH_DEFINES
#include <bits/stdc++.h>
using namespace std;

#define rep(i, n) REP(i, 0, n)
#define ALL(v) v.begin(), v.end()
#define MSG(a) cout << #a << " " << a << endl;
#define REP(i, x, n) for (int i = x; i < n; i++)
#define OP(m) cout << m << endl
typedef long long ll;
typedef unsigned long long ull;

int main()
{
    int w, h;
    while (cin >> h >> w, h)
    {
        int max_area = 0;
        int m[h][w];

        rep(i, h)
        {
            rep(j, w)
            {
                char c;
                cin >> c;
                m[i][j] = (c == '.' ? 1 : 0);
            }
        }
        int a[w + 1] = {};
        a[w] = -1;
        rep(i, h)
        {
            rep(j, w) a[j] = (m[i][j] == 1 ? a[j] + 1 : 0);
            stack<int> st;
            rep(j, w + 1)
            {
                while (!st.empty() && a[st.top()] > a[j])
                {
                    int t = st.top();
                    st.pop();
                    max_area = max(max_area, a[t] * (st.empty() ? j : j - st.top() - 1));
                }
                st.push(j);
            }
        }
        OP(max_area);
    }
    return 0;
}
accuracy: 1
time_ms: 20
memory_kb: 4084
final_rank: 4