#include<bits/stdc++.h>

using namespace std;

const int N = 305;
int n, m, cnt[3 * N][3 * N], ans;
char c[3 * N][3 * N], cc[3 * N][3 * N];

bool valid (int x, int y) { return x > 0 && y > 0 && x <= n && y <= m; }

void rotation () {
    for (int i = 0; i < 3 * N; ++i) {
        for (int j = 0; j < 3 * N; ++j) {
            cc[i][j] = '.';
        }
    }
    for (int i = 1; i <= m; ++i) {
        for (int j = 1; j <= n; ++j) {
            cc[i + m][j + n] = c[j + n][m + 1 - i + m];
        }
    }
    swap(n, m);
    for (int i = 0; i < 3 * N; ++i) {
        for (int j = 0; j < 3 * N; ++j) {
            c[i][j] = '.';
        }
    }
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            c[i + n][j + m] = cc[i + n][j + m];
        }
    }
}

void solve () {
    memset(cnt, 0, sizeof cnt);
    for (int i = 1; i <= 3 * n; ++i) {
        for (int j = 1; j <= 3 * m; ++j) {
            cnt[i][j] = cnt[i - 1][j + 1] + (c[i][j] == '#');
        }
    }
    for (int x = 1; x <= n; ++x) {
        for (int y = 1; y <= m; ++y) {
            for (int d = 1; d <= n; ++d) {
                int X = x - d, Y = y + d;
                if (!valid(X, Y) || c[x + n][y + m] == '.' || c[X + n][Y + m] == '.') continue ;
                ans += cnt[x - d + n][y - d + m] - cnt[x - 2 * d + n][y + m];
            }
        }
    }
}

int main () {
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    cin >> n >> m;
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            cin >> c[i + n][j + m];
        }
    }

    for (int _ = 0; _ < 4; ++_) {
        solve();
        rotation();
    }

    cout << ans;

    return 0;
}
