# include <bits/stdc++.h>
using namespace std;
typedef long long ll;
char s[2][2010][2010],(*t)[2010];
int L[2010],R[2010],dwn[2010][2010];
vector<int> v1[2010],v2[2010];
ll ans;
struct BIT
{
	int n,sum[2010];
	void init(int _n){n=_n;for(int i=1;i<=n;i++) sum[i]=0;}
	void update(int x,int y){for(int i=x;i<=n;i+=i&(-i)) sum[i]+=y;}
	int query(int x)
	{
		int ans=0;
		for(int i=x;i;i-=i&(-i)) ans+=sum[i];
		return ans;
	}
}B1,B2;
void solve(int l1,int r1,int l2,int r2,int op)
{
	if(l1==r1 || l2==r2) return;
	if(r1-l1>r2-l2) swap(l1,l2),swap(r1,r2),op^=1;
	int mid=(l2+r2)/2;
	t=s[op];
	for(int i=l1;i<=r1;i++)
	{
		int p=mid;
		for(;p>l2 && t[i][p-1]==t[i][p];p--);
		L[i]=mid-p+1;
		p=mid+1;
		for(;p<r2 && t[i][p+1]==t[i][p];p++);
		R[i]=p-mid;
	}
	for(int i=r1;i>=l1;i--)
		for(int j=l2;j<=r2;j++)
			if(i==r1 || t[i+1][j]!=t[i][j]) dwn[i][j]=1;
			else dwn[i][j]=dwn[i+1][j]+1;
	for(int i=l1;i<=r1;i++)
	{
		if(t[i][mid]!=t[i][mid+1]) continue;
		B1.init(mid-l2+1);B2.init(r2-mid);
		for(int j=mid;j>=l2 && mid-j+1<=L[i];j--)
			B1.update(mid-j+1,1),
			v1[i+dwn[i][j]].push_back(mid-j+1);
		for(int j=mid+1;j<=r2 && j-mid<=R[i];j++)
			B2.update(j-mid,1),
			v2[i+dwn[i][j]].push_back(j-mid);
		for(int j=i+1;j<=r1;j++)
		{
			for(int k:v1[j]) B1.update(k,-1);
			for(int k:v2[j]) B2.update(k,-1);
			ans+=(ll)B1.query(L[j])*B2.query(R[j]);
		}
		for(int j=i+1;j<=r1+1;j++) v1[j].clear(),v2[j].clear();
	}
	solve(l1,r1,l2,mid,op);solve(l1,r1,mid+1,r2,op);
}
int main()
{
	int n,m;
	cin>>n>>m;
	for(int i=1;i<=n;i++) scanf("%s",s[0][i]+1);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			s[1][j][i]=s[0][i][j];
	solve(1,n,1,m,0);
	cout<<ans<<endl;
	return 0;
}