#include <algorithm>
#include <cassert>
#include <cctype>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <functional>
#include <iostream>
#include <list>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <vector>

using namespace std;

#define pb push_back
#define mp make_pair
#define rep(i,n) for(int i = 0; i < (n); i++)
#define repr(i,b,e) for(int i = (b); i <= (e); i++)

#define pr(x) do{cerr << (#x) << " = " << (x) << endl;}while(0)
#define pri(x,i) do{cerr << (#x) << "[" << i << "] = " << (x[i]) << endl;}while(0)

namespace std {
	template<class S, class T> pair<S, T> operator + (const pair<S, T> &a, const pair<S, T> &b) {return mp<S, T>(a.first + b.first, a.second + b.second);}
	template<class S, class T> pair<S, T> operator - (const pair<S, T> &a, const pair<S, T> &b) {return mp<S, T>(a.first - b.first, a.second - b.second);}
	template<class S, class T> pair<S, T>& operator += (pair<S, T> &a, const pair<S, T> &other) {return a = a + other;}
	template<class S, class T> pair<S, T>& operator -= (pair<S, T> &a, const pair<S, T> &other) {return a = a - other;}
	template<class S, class T> ostream& operator << (ostream& ost, pair<S, T> p) {ost << "(" << p.first << ", " << p.second << ")";return ost;}
	template<class T> ostream& print_iterable_container(ostream& ost, T& container) {
		for (auto&& iter = begin(container); iter != end(container); ++iter) ost << (iter == begin(container) ? " [ " : ", ") << (*iter);
		return ost << " ]";
	}
#	define APPLY_PRINT_ITERABLE_CONTAINER(container_type)\
	template<class T> ostream& operator << (ostream& ost, container_type<T>& container) {\
		return print_iterable_container(ost, container);\
	}
	APPLY_PRINT_ITERABLE_CONTAINER(vector);
	APPLY_PRINT_ITERABLE_CONTAINER(list);
	APPLY_PRINT_ITERABLE_CONTAINER(deque);
}

template<class T> inline bool CHMAX(T &l,const T &r){return l<r?l=r,1:0;}
template<class T> inline bool CHMIN(T &l,const T &r){return l>r?l=r,1:0;}

using ll = long long;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
using vi = vector<int>;
using vll = vector<ll>;
using vvi = vector<vi>;
using vvll = vector<vll>;
using vpii = vector<pii>;
using vpll = vector<pll>;

using ld = long double;

template<class T> using VV = vector<vector<T> >;
template<class T> using maxPQ = priority_queue<T>;
template<class T> using minPQ = priority_queue<T, vector<T>, greater<T> >;
template<class T, class U> using uomap = unordered_map<T, U>;
template<class T> using uoset = unordered_set<T>;

const ll INF = 1001001001;
const ll INFll = (ll)INF * INF;
const ll MOD = 1000000007;
const ll MOD_FFT = 998244353;
const ld EPS = 1e-15;

//function<void(int,int)> tensi_no_3p=[&](int v,int d){
//};

int H, W, T;
vvi B;

int main() {
	cin >> H >> W;
	T = H + W - 1;
	B.resize(T);
	rep(i, T) B[i].resize(T);
	rep(y, H) {
		string s;
		cin >> s;
		rep(x, W) {
			B[y - x + W - 1][y + x] = s[x] == '#' ? 1 : 0;
		}
	}

	vvi Yc = B, Xc = B;
	repr(yy, 1, T - 1) rep(xx, T) Yc[yy][xx] += Yc[yy - 1][xx];
	rep(yy, T) repr(xx, 1, T - 1) Xc[yy][xx] += Xc[yy][xx - 1];

	ll ans = 0;
	rep(yy, T) rep(xx, T) {
		for (int ww = 2; yy + ww < T && xx + ww < T; ww += 2) {
			if (B[yy][xx] && B[yy][xx + ww]) ans += Xc[yy + ww][xx + ww - 1] - Xc[yy + ww][xx];
			if (B[yy][xx] && B[yy + ww][xx]) ans += Yc[yy + ww - 1][xx + ww] - Yc[yy][xx + ww];
			if (B[yy + ww][xx] && B[yy + ww][xx + ww]) ans += Xc[yy][xx + ww - 1] - Xc[yy][xx];
			if (B[yy][xx + ww] && B[yy + ww][xx + ww]) ans += Yc[yy + ww - 1][xx] - Yc[yy][xx];
			if (B[yy][xx] && B[yy][xx + ww] && B[yy + ww][xx]) ans++;
			if (B[yy][xx] && B[yy][xx + ww] && B[yy + ww][xx + ww]) ans++;
			if (B[yy][xx] && B[yy + ww][xx] && B[yy + ww][xx + ww]) ans++;
			if (B[yy][xx + ww] && B[yy + ww][xx] && B[yy + ww][xx + ww]) ans++;
		}
	}
	cout << ans << endl;
	return 0;
}


Measured memory_bytes: 4608