#include <bits/stdc++.h>
#define ll __int128
using namespace std;
const int Nmax=2005;
ll dp[Nmax][Nmax][5][3];
ll s[Nmax][3];
bitset<Nmax> b[Nmax];
int lin=0,col=1;
///lin->0
///col->0
void print(__int128 x) {
if (x < 0) {
putchar('-');
x = -x;
}
if (x > 9) print(x / 10);
putchar(x % 10 + '0');
}
signed main()
{
int n,m;
cin>>n>>m;
for(int i=1;i<=n;i++)
{
string str;
cin>>str;
for(int j=1;j<=m;j++)
{
if(str[j-1]=='1')
{
dp[i][j][1][0]=1;
dp[i][j][1][1]=1;
s[i][col]++;
s[j][lin]++;
b[i][j]=1;
}
else b[i][j]=0;
}
}
ll ans=0;
for(int k=2;k<=5;k++)
{
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(b[i][j])
{
dp[i][j][k][lin]=(s[i][col]-dp[i][j][k-1][col]);
dp[i][j][k][col]=(s[j][lin]-dp[i][j][k-1][lin]);
if(k==5)
{
ans=(ans+dp[i][j][k][lin]+dp[i][j][k][col]);
}
/// cout<<i<<' '<<j<<' '<<k<<' '<<dp[i][j][k][lin]<<' '<<dp[i][j][k][col]<<'\n';
}
}
}
if(k==5) break;
for(int i=1;i<=n;i++) s[i][col]=0;
for(int i=1;i<=m;i++) s[i][lin]=0;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
s[i][col]=(s[i][col]+dp[i][j][k][col]);
s[j][lin]=(s[j][lin]+dp[i][j][k][lin]);
}
}
}
///cout<<ans<<'\n';
for(int i=1;i<=n;i++)
{
for(int j=i+1;j<=n;j++)
{
bitset<Nmax> comb=(b[i]&b[j]);
ll cont=comb.count();
cont=(cont*(cont-1)/2);
//for(int k=1;k<=m;k++) cout<<comb[k]<<' ';
///cout<<'\n';
///cout<<i<<' '<<j<<' '<<cont<<'\n';
ans=(ans-8*cont);
}
}
print(ans);
return 0;
}
