#include <iostream>
using namespace std;
typedef long long ll;

const int M=620,M2=310;
int h,w,a[M][M],b[M][M];

ll f(int B){
	ll res=0;
	for(int i=0;i<M;i++) for(int j=1;j<M;j++) b[i][j]=a[i][j]+b[i][j-1];
	for(int i=0;i<M;i++) for(int j=0;j<M;j++) if(a[i][j]){
		for(int k=j+1;k<M;k++) if(a[i][k]){
			int d=k-j;
			if(i+d<M) res+=b[i+d][k-B]-b[i+d][j-1+B];
			if(d<=i) res+=b[i-d][k-B]-b[i-d][j-1+B];
		}
	}
	for(int i=0;i<M;i++) for(int j=0;j<M;j++) b[M-j-1][i]=a[i][j];
	for(int i=0;i<M;i++) for(int j=0;j<M;j++) a[i][j]=b[i][j],b[i][j]=0;
	return res;
}

int main(){
	cin>>h>>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][i-j+M2+1]++;
	}
	cout<<f(0)+f(1)<<endl;
}