This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include <bits/stdc++.h>
using namespace std;
#ifdef LOCAL
#include "debug.h"
#else
#define debug(...) 42
#endif
#define file ""
#define mp make_pair
#define fi first
#define se second
#define all(x) x.begin(), x.end()
#define bit(x) (1LL << (x))
#define getbit(x, i) (((x) >> (i)) & 1)
#define popcount __builtin_popcountll
mt19937_64 rd(chrono::steady_clock::now().time_since_epoch().count());
int rand(int l, int r) {
    return l + rd() % (r - l + 1);
}
const int N = 2e3 + 5;
const int mod = (int)1e9 + 7; // 998244353;
const int lg = 25; // lg + 1
const int oo = 1e9;
const long long ooo = 1e18;
template<class X, class Y> bool mini(X &a, Y b) {
    return a > b ? (a = b, true) : false;
}
template<class X, class Y> bool maxi(X &a, Y b) {
    return a < b ? (a = b, true) : false;
}
void add(int &a, int b) {
    a += b;
    if (a >= mod) a -= mod;
    if (a < 0) a += mod;
}
int m, n;
int a[N][N];
int h[N];
int l[N], r[N];
long long ans = 0;
void process() {
    cin >> m >> n;
    for (int i = 1; i <= m; ++i)
        for (int j = 1; j <= n; ++j) {
            char c;
            cin >> c;
            a[i][j] = c == '.';
        }
    for (int i = 1; i <= m; ++i) {
        for (int j = 1; j <= n; ++j) if (a[i][j] == 1) h[j] += 1; else h[j] = 0;
        stack<int> st;
        for (int j = 1; j <= n; ++j) {
            while (st.size() && h[st.top()] >= h[j]) st.pop();
            l[j] = st.size() ? st.top() : 0;
            st.push(j);
        }
        st = stack<int> {};
        for (int j = n; j >= 1; --j) {
            while (st.size() && h[st.top()] > h[j]) st.pop();
            r[j] = st.size() ? st.top() : n + 1;
            st.push(j);
        }
        auto sum = [&] (int x) {
            return x * (x + 1) / 2;
        };
        for (int j = 1; j <= n; ++j) {
            int lenl = j - l[j], lenr = r[j] - j;
            long long cur = 1LL * lenl * sum(lenr) + 1LL * lenr * sum(lenl) - 1LL * lenl * lenr;
            ans += 1LL * cur * sum(h[j]);
        }
    }
    cout << ans;
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    // freopen(file".inp", "r", stdin);
    // freopen(file".out", "w", stdout);
    int tc = 1;
    // cin >> tc;
    
    while (tc--) {
        process();
    }
    return 0;
}
/*
*/
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict | Execution time | Memory | Grader output | 
|---|
| Fetching results... |