#include <bits/stdc++.h>
#include <iostream>

#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>

#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define rep_r(i, n) for (int i = n - 1; i >= 0; i--)
#define rep1(i, n) for (int i = 1; i <= (int)(n); i++)
#define REP(i, n) for (int i = 0, i##_len = (n); i < i##_len; ++i)
#define all(x) (x).begin(), (x).end()
#define SZ(x) ((ll)(x).size())
#define bit(n) (1LL << (n))
#define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end());

#define INF bit(60)
#define pb push_back
#define mod 1000000007

using namespace std;
using uif = uint_fast64_t;
using ll = long long int;
using tTree = __gnu_pbds::tree<ll, __gnu_pbds::null_type, less<ll>, __gnu_pbds::rb_tree_tag, __gnu_pbds::tree_order_statistics_node_update>;

ll dx[4] = {1, 0, -1, 0};
ll dy[4] = {0, 1, 0, -1};

#define FACSIZE 200002

ll invfac[FACSIZE];
ll fac[FACSIZE];

template <class T>
bool chmax(T &, const T &);
template <class T>
bool chmin(T &, const T &);
ll gcd(ll, ll);
ll powLL(ll x, ll y);

ll mod_pow(ll, ll);
ll mod_add(ll, ll);
ll mod_mul(ll, ll);
ll mod_div(ll, ll);
ll comb(ll, ll);
void make_fact(ll);
void make_invfact(ll);
void fact_init();

vector<ll> dp;
ll dir[4][2] = {
    {1, 0},
    {0, 1},
    {-1, 0},
    {0, -1}};

void init(ll);
ll root(ll);
void unite(ll, ll);
bool same(ll, ll);

bool range(ll y, ll x, ll h, ll w)
{
  return x >= 0 && y >= 0 && x < w && y < h;
}

ll enc(ll y, ll x, ll w)
{
  return w * y + x;
}

pair<ll, ll> dec(ll e, ll w)
{

  auto y = e / w;
  auto x = e % w;

  return make_pair(y, x);
}

int main(void)
{
  ll h, w;
  cin >> h >> w;
  vector<string> vec;
  rep(i, h)
  {
    string s;
    cin >> s;
    vec.push_back(s);
  }
  /*
  for (auto s : vec)
    cout << s << endl;
    */

  init(h * w);

  for (ll i = 0; i < h; i++)
  {
    for (ll j = 0; j < w; j++)
    {
      char c = vec[i][j];
      if (c == '#')
      {
        for (ll k = 0; k < 4; k++)
        {
          ll nextY = i + dir[k][0];
          ll nextX = j + dir[k][1];
          if (range(nextY, nextX, h, w))
          {
            if (vec[nextY][nextX] == '.')
              unite(enc(nextY, nextX, w), enc(i, j, w));
          }
        }
      }
      else
      {
        for (ll k = 0; k < 4; k++)
        {
          ll nextY = i + dir[k][0];
          ll nextX = j + dir[k][1];
          if (range(nextY, nextX, h, w))
          {
            if (vec[nextY][nextX] == '#')
              unite(enc(nextY, nextX, w), enc(i, j, w));
          }
        }
      }
    }
  }
  /*
  for (ll i = 0; i < h * w; i++)
  {
    cout << root(dp[i]) << endl;
  }
  */

  unordered_map<ll, ll> umB;
  unordered_map<ll, ll> umW;

  for (ll i = 0; i < h; i++)
  {
    for (ll j = 0; j < w; j++)
    {
      auto r = root(enc(i, j, w));

      //cout << vec[i][j] << endl;
      if (vec[i][j] == '#')
      {
        umB[r]++;
      }
      else
      {
        umW[r]++;
      }
    }
  }
  /*
  cout << "umB:" << endl;
  for (auto p : umB)
  {
    cout << p.first << "," << p.second << endl;
  }
  cout << endl;
  cout << "umW:" << endl;
  for (auto p : umW)
  {
    cout << p.first << "," << p.second << endl;
  }
  */

  ll ans = 0;
  for (auto p : umB)
  {

    auto key = p.first;
    if (umW.find(key) != end(umW))
    {
      ans += p.second * umW[key];
    }
  }

  cout << ans << endl;

  return 0;
}

void init(ll n)
{
  dp.resize(n);
  rep(i, n) dp[i] = i;
}
ll root(ll x)
{
  if (dp[x] == x)
    return x;
  else
    return dp[x] = root(dp[x]);
}
void unite(ll x, ll y)
{
  auto rx = root(x);
  auto ry = root(y);
  dp[ry] = dp[rx];
  root(y);
}
bool same(ll x, ll y)
{
  return root(x) == root(y);
}

ll mod_pow(ll x, ll r)
{
  if (r == 0)
    return 1;
  else if (r == 1)
    return x % mod;
  else if (r % 2 == 0)
  {
    ll t = mod_pow(x, r / 2) % mod;
    return mod_mul(t, t);
  }
  else
  {
    ll t = mod_pow(x, r / 2) % mod;
    ll k = mod_mul(t, t);
    return (k % mod) * (x % mod);
  }
}
ll mod_add(ll a, ll b)
{
  return ((a % mod) + (b % mod)) % mod;
}

ll mod_mul(ll a, ll b)
{
  return ((a % mod) * (b % mod)) % mod;
}

ll mod_div(ll a, ll b)
{
  return mod_mul(a, mod_pow(b, mod - 2));
}

void fact_init()
{
  make_fact(FACSIZE - 1);
  make_invfact(FACSIZE);
}

void make_fact(ll n)
{
  fac[0] = 1;
  rep(i, n)
  {
    fac[i + 1] = mod_mul(fac[i], i + 1);
  }
}

void make_invfact(ll n)
{
  invfac[n] = mod_pow(fac[n], mod - 2);
  for (int i = n - 1; i >= 0; i--)
  {
    invfac[i] = mod_mul(invfac[i + 1], i + 1);
  }
}

ll comb(ll n, ll r)
{

  return mod_mul(mod_mul(fac[n], invfac[r]), invfac[n - r]);
}
template <class T>
bool chmax(T &a, const T &b)
{
  if (a < b)
  {
    a = b;
    return 1;
  }
  return 0;
}
template <class T>
bool chmin(T &a, const T &b)
{
  if (b < a)
  {
    a = b;
    return 1;
  }
  return 0;
}

ll qp(ll a, ll b)
{
  ll ans = 1LL;
  do
  {
    if (b & 1LL)
      ans = 1LL * mod_mul(ans, a) % mod;
    a = 1LL * mod_mul(a, a) % mod;
  } while (b >>= 1LL);
  return ans;
}
ll qp(ll a, ll b, ll mo)
{
  ll ans = 1LL;
  do
  {
    if (b & 1LL)
      ans = 1LL * (ans % mo) * (a % mo);
    a = 1LL * (a % mo) * (a % mo);
  } while (b >>= 1LL);
  return ans;
}
ll gcd(ll a, ll b)
{
  return b ? gcd(b, a % b) : a;
}
ll powLL(ll x, ll y)
{
  ll ans = 1LL;
  for (ll i = 0LL; i < y; i++)
    ans *= x;
  return ans;
}

Measured memory_bytes: 8484