#include "bits/stdc++.h"
using namespace std;
typedef long long ll;
const ll INF = 1LL << 60;
const ll MOD = 998244353;
typedef pair<int, int> P;
bool exist[600][600];
int x[600][600];
int y[600][600];
int main() {
	int H, W;
	cin >> H >> W;
	vector<string> s(H);
	for (int i = 0; i < H; i++) cin >> s[i];
	vector<P> coin;
	for (int i = 0; i < H; i++) {
		for (int j = 0; j < W; j++) {
			if (s[i][j] == '#') {
				int nx = j - i + 300, ny = j + i;
				x[ny][nx]++;
				y[ny][nx]++;
				exist[ny][nx] = 1;
				coin.emplace_back(nx, ny);
			}
		}
	}
	for (int i = 0; i < 600; i++) {
		for (int j = 1; j < 600; j++) {
			x[i][j] += x[i][j - 1];
		}
	}
	for (int i = 0; i < 600; i++) {
		for (int j = 1; j < 600; j++) {
			y[j][i] += y[j - 1][i];
		}
	}
	ll ans = 0;
	ll t = 0;
	for (int d = 1; d < 600; d++) {
		for (P p : coin) {
			if (p.first + d < 600 && exist[p.second][p.first + d]) {
				if (p.second + d < 600) {
					ans += x[p.second + d][p.first + d];
					if (p.first > 0) ans -= x[p.second + d][p.first - 1];
				}
				if (p.second - d >= 0) {
					ans += x[p.second - d][p.first + d];
					if (p.first > 0) ans -= x[p.second - d][p.first - 1];
				}
			}
			if (p.second + d < 600 && exist[p.second + d][p.first]) {
				if (p.first + d < 600) {
					ans += y[p.second + d][p.first + d];
					if (p.second > 0) ans -= y[p.second - 1][p.first + d];
				}
				if (p.first - d >= 0) {
					ans += y[p.second + d][p.first - d];
					if (p.second > 0) ans -= y[p.second - 1][p.first - d];
				}
			}
			for (int i = -1; i <= 1; i += 2) {
				int x = p.first + d * i;
				if (x < 0 || x >= 600 || !exist[p.second][x]) continue;
				for (int j = -1; j <= 1; j += 2) {
					int y = p.second + d * j;
					if (y < 0 || y >= 600 || !exist[y][p.first]) continue;
					t++;
				}
			}
		}
	}
	ans -= t;
	cout << ans << endl;
 }