#include <map>
#include <set>
#include <array>
#include <cmath>
#include <deque>
#include <bitset>
#include <random>
#include <string>
#include <vector>
#include <cassert>
#include <complex>
#include <iomanip>
#include <iostream>
#include <algorithm>
#include <unordered_map>
#include <unordered_set>

using namespace std;

vector<int> get_z(string s) {
  int n = (int)s.size();
  int j = 0;
  vector<int> z(n);
  for (int i = 1; i < n; ++i) {
    if (j + z[j] > i)
      z[i] = min(z[i - j], j + z[j] - i);
    while (s[z[i]] == s[i + z[i]])
      ++z[i];
    if (i + z[i] > j + z[j])
      j = i;
  }
  return z;
}

const int MAXN = 3'000'100;
const int ALPHA = 3;

int go[MAXN][ALPHA];
int suff[MAXN];

int depth[MAXN];
int real_go[MAXN][ALPHA];
int non3term[MAXN];
int term[MAXN];
int nodes_cnt = 0;

int node() {
  fill(real_go[nodes_cnt], real_go[nodes_cnt] + ALPHA, -1);
  depth[nodes_cnt] = 0;
  term[nodes_cnt] = 0;
  ++nodes_cnt;
  return nodes_cnt - 1;
}

void dfs_cnt(int v) {
  for (int i = 0; i < ALPHA; ++i) {
    if (real_go[v][i] == -1)
      continue;
    dfs_cnt(real_go[v][i]);
    term[v] += term[real_go[v][i]];
  }
}

long long ans = 0;

vector<int> stack;

void dfs_karasik(int v) {
  stack.emplace_back(v);
  int szs[ALPHA];
  int periods[ALPHA];
  
  
  for (int i = 0; i < ALPHA; ++i) {
    periods[i] = depth[v] + 1 - depth[v == 0 ? 0 : go[suff[v]][i]];  
    if (real_go[v][i] == -1) {
      szs[i] = 0;
      continue;
    }
    szs[i] = term[real_go[v][i]];


    fill(go[v], go[v] + ALPHA, -1);
    go[v][i] = real_go[v][i];
    for (int j = 0; j < ALPHA; ++j) {
      if (go[v][j] == -1) {
        go[v][j] = (v == 0 ? 0 : go[suff[v]][j]);
      } else {
        suff[go[v][j]] = (v == 0 ? 0 : go[suff[v]][j]);
      }
    }
    dfs_karasik(go[v][i]);
  }
  ans += szs[0] * 1ll * szs[1] * 1ll * szs[2];

  //cerr << v << ' ' << szs[0] << ' ' << szs[1] << ' ' << szs[2] << '\n';

  for (int i = 0; i < ALPHA; ++i) {
    if (periods[i] * 3 <= depth[v])
      ans += non3term[stack[periods[i]]] * 1ll * szs[(i + 1) % 3] * 1ll * szs[(i + 2) % 3];
  }

  stack.pop_back();
}

long long solve(vector<string> lines) {
  nodes_cnt = 0;
  ans = 0;
  node();

  int n = lines.size();
  for (int i = 0; i < n; ++i) {
    string s = lines[i];

    int len = s.size();
    auto z = get_z(s);
    int per = 1;
    while (per < len && (len % per != 0 || z[per] != len - per))
      ++per;
    s.resize(per);

    vector<int> cur;
    for (char c : s) {
      if (c == 'P')
        cur.emplace_back(0);
      else if (c == 'N')
        cur.emplace_back(1);
      else
        cur.emplace_back(2);
    }
    int cur_node = 0;
    for (int x : cur) {
      if (real_go[cur_node][x] == -1) {
        real_go[cur_node][x] = node();
        depth[real_go[cur_node][x]] = depth[cur_node] + 1;
      }
      cur_node = real_go[cur_node][x];
    }
    ++non3term[cur_node];
    for (int j = 0; j < 2; ++j) {
      for (int x : cur) {
        if (real_go[cur_node][x] == -1) {
          real_go[cur_node][x] = node();
          depth[real_go[cur_node][x]] = depth[cur_node] + 1;
        }
        cur_node = real_go[cur_node][x];
      }
    }
    ++term[cur_node];
  }
  dfs_cnt(0);
  dfs_karasik(0);
  return ans;
}

long long dumb(vector<string> lines) {
  long long ans = 0;
  for (int a = 0; a < lines.size(); ++a) {
    for (int b = a + 1; b < lines.size(); ++b) {
      for (int c = b + 1; c < lines.size(); ++c) {
        string s1 = lines[a], s2 = lines[b], s3 = lines[c];
        int lol = s1.size() + s2.size() + s3.size();
        for (int i = 0; i < lol; ++i) {
          char c1 = s1[i % s1.size()];
          char c2 = s2[i % s2.size()];
          char c3 = s3[i % s3.size()];
          int kek = (c1 != c2) + (c2 != c3) + (c3 != c1);
          if (kek) {

            if (kek == 3) {
              //cerr << s1 << ' ' << s2 << ' ' << s3 << '\n';
            }

            ans += kek == 3;
            break;
          }
        }
      }
    }
  }
  return ans;
}

//#define STRESS

int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);

#ifdef STRESS
  char symbols[3] = {'P', 'N', 'K'};
  for (int _ = 0; _ < 100100; ++_) {
    int n = rand() % 7 + 1;
    vector<string> lines(n);
    for (int i = 0; i < n; ++i) {
      int len = rand() % 5 + 1;
      for (int j = 0; j < len; ++j)
        lines[i] += symbols[rand() % 3];
    }
    auto ans1 = solve(lines);
    auto ans2 = dumb(lines);
    if (ans1 != ans2) {
      cerr << "BUG!!!\n";
      cerr << n << '\n';
      for (auto line : lines)
        cerr << line << '\n';
      cerr << "wrong: " << ans1 << " right: " << ans2 << "\n\n\n";
    }
  }

#else
  int n;
  cin >> n;
  vector<string> lines(n);
  for (int i = 0; i < n; ++i)
    cin >> lines[i];
  cout << solve(lines) << '\n';
#endif
  return 0;
}