#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
char s[2][2010][2010];
int n,m,up[2010][2010],L[2010],R[2010];
vector<int>l[2010],r[2010];
ll ans=0;
struct bit{
	int tr[2010],n;
	void init(int N){
		n=N;
		fill(tr,tr+n+1,0);
	}
	void add(int x,int c){
		for(int i=x;i<=n;i+=(i&-i)) tr[i]+=c;
	}
	int qry(int x){
		int res=0;
		for(int i=x;i;i-=(i&-i)) res+=tr[i];
		return res; 
	}
	int qry(int l,int r){
		return qry(r)-qry(l-1);
	}
}lt,rt;
void sol(int x,int y,int p,int q,int ty){
	if(y-x>r-l){
		ty^=1,swap(x,p),swap(y,q);
	}
	if(x==y) return;
	int mid=(p+q)>>1;
	sol(x,y,p,mid,ty),sol(x,y,mid+1,q,ty);
	for(int i=x;i<=y;i++) for(int j=p;j<=q;j++) up[i][j]=i>x&&s[ty][i-1][j]==s[ty][i][j]?up[i-1][j]:i;
	for(int i=x;i<=y;i++){
		for(int j=mid  ;j>=p&&s[ty][i][mid  ]==s[ty][i][j];j--) L[i]=j;
		for(int j=mid+1;j<=q&&s[ty][i][mid+1]==s[ty][i][j];j++) R[i]=j;
	}
	for(int i=x+1;i<=y;i++) if(s[ty][i][mid]==s[ty][i][mid+1]){
		for(int j=mid  ;j>=L[i];j--) if(up[i][j]<i) l[up[i][j]].push_back(j);
		for(int j=mid+1;j<=R[i];j++) if(up[i][j]<i) r[up[i][j]].push_back(j);
		for(int j=x;j<i;j++){
			for(int k:l[j]) lt.add(k,1);
			for(int k:r[j]) rt.add(k,1);
			if(s[ty][j][mid]!=s[ty][i][mid]||s[ty][j][mid+1]!=s[ty][i][mid]) continue;
			ans+=1ll*lt.qry(L[j],mid)*rt.qry(mid+1,R[j]);
		}
		for(int j=x;j<i;j++){
			for(int k:l[j]) lt.add(k,-1);
			for(int k:r[j]) rt.add(k,-1);
			l[j].clear(),r[j].clear();
		}
	}
}
int main(){
	scanf("%d%d",&n,&m),lt.init(max(n,m)),rt.init(max(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];
	sol(1,n,1,m,0);
	printf("%lld\n",ans);
}