#include<cstdio>
#include<algorithm>
#include<cstring>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
using namespace std;

const int N=610;
struct point{ int x,y; }p[90010];
int h,w,m=300,cnt,a[N][N],r[N][N],c[N][N];
long long ans;

int main(){
	scanf("%d%d",&h,&w);
	rep(i,1,h){
		char str[310]; scanf("%s",str+1);
		rep(j,1,w){
			if (str[j]=='#'){
				int x=i+j-1,y=i-j+m;
				p[cnt++]=(point){x,y}; a[x][y]=1;
			}
		}
	}
	m*=2;
	rep(i,1,m) rep(j,1,m) r[i][j]=r[i][j-1]+a[i][j],c[i][j]=c[i-1][j]+a[i][j];
	for (int i=0;i<cnt;i++){
		int x=p[i].x,y=p[i].y;
		for (int d=2;d<=m;d+=2){
			if (y+d<=m && a[x][y+d]){
				if (x+d<=m) ans+=r[x+d][y+d]-r[x+d][y-1];
				if (x-d>0) ans+=r[x-d][y+d]-r[x-d][y-1];
			}
			if (x+d<=m && a[x+d][y]){
				if (y+d<=m) ans+=c[x+d-1][y+d]-c[x][y+d];
				if (y-d>0) ans+=c[x+d-1][y-d]-c[x][y-d];
			}
		}
	}
	printf("%lld\n",ans);
	return 0;
}