#include <algorithm>
#include <iostream>
constexpr int N(3010);
int main()
{
	static char s[N][N];
	std::ios::sync_with_stdio(false);
	std::cin.tie(nullptr);
	int n,m,c(0);std::cin>>n>>m;
	for(int i(1);i<=n;++i)
	{
		std::cin>>s[i]+1;
		for(int j(1);j<=m;++j) c+=(s[i][j]=='.');
	}
	static bool clu[N][N],crd[N][N];
	if(s[1][1]=='.') clu[1][1]=1;if(s[n][m]=='.') crd[n][m]=true;
	for(int i(1);i<=n;++i) for(int j(1);j<=m;++j) if(i!=1 || j!=1) {if(s[i][j]=='*') clu[i][j]=false;else clu[i][j]=clu[i-1][j]|clu[i][j-1];}
	for(int i(n);i;--i) for(int j(m);j;--j) if(i!=n || j!=m) {if(s[i][j]=='*') crd[i][j]=false;else crd[i][j]=crd[i+1][j]|crd[i][j+1];}
	if(!clu[n][m])
	{
		std::cout<<1ll*c*(c-1)/2<<std::endl;
		return 0;
	}
	static bool cap[N][N];
	for(int i(1);i<=n;++i) for(int j(1);j<=m;++j) cap[i][j]=clu[i][j]&crd[i][j];
	static int idx1[N][N],idx2[N][N];idx1[1][1]=idx2[1][1]=1;
	static std::pair<int,int> p1[N<<1],p2[N<<1];p1[1]=p2[1]={1,1};
	for(int i(1),x(1),y(1);i<=n+m-2;++i)
	{
		if(cap[x+1][y]) ++x;
		else ++y;
		idx1[x][y]=i+1;
		p1[i+1]={x,y};
	}
	for(int i(1),x(1),y(1);i<=n+m-2;++i)
	{
		if(cap[x][y+1]) ++y;
		else ++x;
		idx2[x][y]=i+1;
		p2[i+1]={x,y};
	}
	int cnt(0);
	for(int i(1);i<=n;++i) for(int j(1);j<=m;++j) cnt+=(idx1[i][j] && idx2[i][j]);
	int ans(cnt*(cnt-1)/2+cnt*(c-cnt));
	static int g[N][N],h[N][N];
	for(int i(1);i<=n;++i) for(int j(1);j<=m;++j) if(cap[i][j])
	{
		if(idx2[i][j]) g[i][j]=idx2[i][j];
		else
		{
			g[i][j]=1e9;
			if(cap[i-1][j]) g[i][j]=std::min(g[i][j],g[i-1][j]);
			if(cap[i][j-1]) g[i][j]=std::min(g[i][j],g[i][j-1]);
		}
	}
	for(int i(n);i;--i) for(int j(m);j;--j) if(cap[i][j])
	{
		if(idx2[i][j]) h[i][j]=idx2[i][j];
		else
		{
			h[i][j]=0;
			if(cap[i+1][j]) h[i][j]=std::max(h[i][j],h[i+1][j]);
			if(cap[i][j+1]) h[i][j]=std::max(h[i][j],h[i][j+1]);
		}
	}
	static int pre[N<<1],suf[N<<1];
	for(int i(1);i<=n+m-1;++i)
	{
		auto [x,y]=p2[i];
		if(idx1[x][y]) pre[i]=0;
		else
		{
			int r(1e9);
			if(cap[x-1][y]) r=std::min(r,g[x-1][y]);
			if(cap[x][y-1]) r=std::min(r,g[x][y-1]);
			pre[i]=pre[r]+1;
		}
	}
	for(int i(n+m-1);i;--i)
	{
		auto [x,y]=p2[i];
		if(idx1[x][y]) suf[i]=0;
		else
		{
			int r(0);
			if(cap[x+1][y]) r=std::max(r,h[x+1][y]);
			if(cap[x][y+1]) r=std::max(r,h[x][y+1]);
			suf[i]=suf[r]+1;
		}
	}
	for(int i(1),j(0);i<=n+m-1;++i)
	{
		auto [x,y]=p1[i];
		if(idx2[x][y]){j=i;continue;}
		int r(1e9),l(0);
		for(int t(y+1);t<=m;++t) if(cap[x-1][t])
		{
			r=std::min(r,g[x-1][t]),l=std::max(l,h[x-1][t]);
			if(idx2[x-1][t]) break;
		}
		// std::cerr<<x<<" "<<y<<" "<<l<<" "<<r<<" "<<pre[r]<<" "<<suf[l]<<" "<<pre[r]+suf[l]-(l==r)<<std::endl;
		ans+=pre[r]+suf[l]-(l==r);
	}
	std::cout<<ans<<std::endl;
	return 0;
}