#include<bits/stdc++.h>
#define int long long
using namespace std;
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-f;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int buf[30];
inline void print(int x,char ch=' '){
	if(x<0) putchar('-'),x=-x;
	int tot=0;
	do{
		buf[++tot]=x%10;
		x/=10;
	}while(x);
	for(int i=tot;i;i--) putchar(buf[i]+'0');
	putchar(ch);
}
const int MAXN=2e3+5;
int n,m,ans;
char s[MAXN][MAXN];
int down[MAXN][MAXN],up[MAXN][MAXN];
int pre[MAXN][MAXN],suf[MAXN][MAXN];
void init(){
	for(int i=1;i<=m;i++) down[n][i]=1;
	for(int i=n-1;i>=1;i--)
		for(int j=1;j<=m;j++)
			down[i][j]=1+down[i+1][j]*(s[i][j]==s[i+1][j]);
	
	for(int i=1;i<=m;i++) up[1][i]=1;
	for(int i=2;i<=n;i++)
		for(int j=1;j<=m;j++)
			up[i][j]=1+up[i-1][j]*(s[i][j]==s[i-1][j]);
	
	for(int i=1;i<=n;i++) pre[i][1]=1;
	for(int j=2;j<=m;j++)
		for(int i=1;i<=n;i++)
			pre[i][j]=1+pre[i][j-1]*(s[i][j]==s[i][j-1]);
	
	for(int i=1;i<=n;i++) suf[i][m]=1;
	for(int j=m-1;j>=1;j--)
		for(int i=1;i<=n;i++)
			suf[i][j]=1+suf[i][j+1]*(s[i][j]==s[i][j+1]);
}

int t1[MAXN][MAXN],t2[MAXN][MAXN],t3[MAXN][MAXN],t4[MAXN][MAXN];
void dfs(int x1,int y1,int x2,int y2){
	if(x1>=x2||y1>=y2)
		return ;
	int n=x2-x1+1,m=y2-y1+1;
	if(n<=m){
		int mid=(y1+y2)>>1;
		for(int i=x1;i<=x2;i++){
			for(int j=y1;j<=mid;j++){
				if(mid-j+1<=pre[i][mid]){
					t1[i][min(n,down[i][j])]++;
					t3[i][min(n,up[i][j])]++;
				}
			}
			for(int j=mid;j<=y2;j++){
				if(j-mid+1<=suf[i][mid]){
					t2[i][min(n,down[i][j])]++;
					t4[i][min(n,up[i][j])]++;
				}
			}
			for(int j=n-1;j>=1;j--){
				t1[i][j]+=t1[i][j+1];
				t2[i][j]+=t2[i][j+1];
				t3[i][j]+=t3[i][j+1];
				t4[i][j]+=t4[i][j+1];
			}
		}
		
		for(int i=x1;i<=x2;i++){
			for(int j=i+1;j<=x2;j++){
				int A,B;
				if(pre[i][mid]>=pre[j][mid]) A=t3[j][j-i+1];
				else A=t1[i][j-i+1];
				if(suf[i][mid]>=suf[j][mid]) B=t4[j][j-i+1];
				else B=t2[i][j-i+1];
				ans+=A*B;
				if(down[i][mid]>=j-i+1) ans--;
			}
		}
		
		for(int i=x1;i<=x2;i++)
			for(int j=1;j<=n;j++)
				t1[i][j]=t2[i][j]=t3[i][j]=t4[i][j]=0;
			
		dfs(x1,y1,x2,mid-1);
		dfs(x1,mid+1,x2,y2);
	}
	else{
		int mid=(x1+x2)>>1;
		for(int j=y1;j<=y2;j++){
			for(int i=x1;i<=mid;i++){
				if(mid-i+1<=up[mid][j]){
					t1[j][min(m,pre[i][j])]++;
					t3[j][min(m,suf[i][j])]++;
				}
			}
			for(int i=mid;i<=x2;i++){
				if(i-mid+1<=down[mid][j]){
					t2[j][min(m,pre[i][j])]++;
					t4[j][min(m,suf[i][j])]++;
				}
			}
			for(int i=m-1;i>=1;i--){
				t1[j][i]+=t1[j][i+1];
				t2[j][i]+=t2[j][i+1];
				t3[j][i]+=t3[j][i+1];
				t4[j][i]+=t4[j][i+1];
			}
		}
		
		for(int i=y1;i<=y2;i++){
			for(int j=i+1;j<=y2;j++){
				int A,B;
				if(up[mid][i]>=up[mid][j]) A=t1[j][j-i+1];
				else A=t3[i][j-i+1];
				if(down[mid][i]>=down[mid][j]) B=t2[j][j-i+1];
				else B=t4[i][j-i+1];
				ans+=A*B;
				if(suf[mid][i]>=j-i+1) ans--;
			}
		}
		
		for(int j=y1;j<=y2;j++)
			for(int i=1;i<=m;i++)
				t1[j][i]=t2[j][i]=t3[j][i]=t4[j][i]=0;
				
		dfs(x1,y1,mid-1,y2);
		dfs(mid+1,y1,x2,y2);
	}
}

signed main(){
	n=read(),m=read();
	for(int i=1;i<=n;i++) scanf("%s",s[i]+1);
	init();
	dfs(1,1,n,m);
	print(ans);
	return 0;
}