#include <iostream> // cin, cout, string
#include <algorithm> // sort, min/max, upper/lower_bound, transform, count/copy_if
#include <vector>
#include <utility> // pair, swap
#include <functional>
#include <iomanip> // setprecision
#include <map>
#include <set> // multiset
#include <unordered_map>
#include <unordered_set>
#include <tuple>
#include <stack>
#include <queue> // priority_queue
#include <deque>
#include <bitset>
#include <array>
#include <list>
#include <forward_list>
#include <complex>
#include <cassert>
#include <numeric> // accumulate, itoa
#include <random> // mt19937
#include <chrono>
#include <regex>
//#include <string>
//#include <iterator> // begin, end, next, size(C++17)
//#include <cmath>
//#include <cstdint>
//#include <limits>
//#include <type_traits>
//#include <optional> // C++17
#include <boost/range/irange.hpp>
#include <boost/optional.hpp>
//#include <boost/multiprecision/cpp_int.hpp>

using namespace std;

//#define int mnt
#define let const auto&
#define var auto&&

using mnt = int_fast64_t; // int_least64_t, intmax_t

using ll = long long;
using ld = long double;

template<typename T=mnt> using vec = vector<T>;
template<typename T=mnt> using mat = vec<vec<T>>;
template<typename T=mnt> using dic = map<string, T>;
template<typename T=mnt> using dict = unordered_map<string, T>;
template<typename T=mnt> using pqg = priority_queue<T, vector<T>, greater<T>>; // dijkstra
//template<typename T=mnt> using pq = priority_queue<T>;

using pii = pair<int, int>;
using vi = vector<int>;
using mi = vector<vector<int>>; // graph
using vp = vector<pii>;
using vb = vector<bool>;

//using cint = boost::multiprecision::cpp_int;
using boost::irange;
template<typename T=mnt> using opt = boost::optional<T>;

constexpr int IINF = 0x3f3f3f3f;
constexpr mnt LINF = 0x3f3f3f3f3f3f3f3fLL;
constexpr int INF = (int)LINF;

//constexpr mnt MNT_MIN = numeric_limits<mnt>::min();
//constexpr mnt MNT_MAX = numeric_limits<mnt>::max();

constexpr mnt MOD = 1e9+7;

#define _overload(_1,_2,_3,name,...) name
#define _rep1(i,a,b) for(int i=(a), _##i##_=(b);i<_##i##_;++i)
#define _rep2(i,n) _rep1(i,0,n)
#define _rep3(n) _rep2(i,n)
#define rep(...) _overload(__VA_ARGS__,_rep1,_rep2,_rep3)(__VA_ARGS__)
#define _rrep1(i,a,b) for(int i=(b)-1, _##i##_=(a);i>=_##i##_;--i)
#define _rrep2(i,n) _rrep1(i,0,n)
#define _rrep3(n) _rrep2(i,n)
#define rrep(...) _overload(__VA_ARGS__,_rrep1,_rrep2,_rrep3)(__VA_ARGS__)

#define eb emplace_back
#define pb push_back
#define pob pop_back
#define puf push_front
#define pof pop_front

#define np next_permutation
#define pp prev_permutation

#define mp make_pair
#define mt make_tuple
#define fat forward_as_tuple

#define ft first
#define sd second

#define all(a) begin(a), end(a)
#define rall(a) rbegin(a), rend(a)

#define elif else if
#define unless(a) if(!(a))
#define qt(a) #a

#define endl '\n'
#define decimal(n) fixed << setprecision(n)
#define zfl(n) setfill(‘0’) << setw(n)

template<typename T=int, class C> inline constexpr T sz(const C& a) { return (T)a.size(); };

template<class C> inline constexpr C& rev(C& a) { reverse(begin(a), end(a)); return a; }
template<class C> inline constexpr C& sort(C& a) { sort(begin(a), end(a)); return a; }

template<class C> inline constexpr auto& min(const C& a) { return *min_element(begin(a), end(a)); }
template<class C> inline constexpr auto& max(const C& a) { return *max_element(begin(a), end(a)); }

template<typename T> inline constexpr void rs(vector<T>& a, size_t n) { a.resize(n); }
// template<typename T, typename... Ts> inline constexpr void rs(vector<T>& a, size_t n, Ts&... m) { rs(a, n); for(auto&& i : a) rs(i, m...); }

inline constexpr char upper(char c) { return ('a' <= c && c <= 'z' ? c-'a'+'A' : c); }
inline constexpr char lower(char c) { return ('A' <= c && c <= 'Z' ? c-'A'+'a' : c); }

template<class C> inline constexpr auto sum(const C& c) { return accumulate(next(begin(c)), end(c), *begin(c)); }
template<class C, class F> inline constexpr auto foldl(const C& c, const F& f) { return accumulate(next(begin(c)), end(c), *begin(c), f); }

template<typename T=mnt, T l=0, T r=MOD>
auto rnd = bind(uniform_int_distribution<T>(l, r), mt19937(chrono::steady_clock::now().time_since_epoch().count()));

template<typename T> inline constexpr T gcd2(T a, T b) { if(a==0 || b==0) return 0; while(b) { a = a%b; swap(a, b); } return a; }
template<typename T> inline constexpr T extgcd2(T a, T b, T& x, T& y) { x = 1; y = 0; while (b) { T t = a / b; a -= t * b; swap(a, b); x -= t * y; swap(x, y); } return a; }
template<typename T> inline constexpr T lcm(T a, T b) { return a / gcd(a, b) * b; }

template<typename T> inline constexpr bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } return false; }
template<typename T> inline constexpr bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } return false; }

template<typename T, typename U> istream& operator>>(istream& is, pair<T, U>& p) { is >> p.first; is >> p.second; return is; }
template<typename T, typename U> ostream& operator<<(ostream& os, const pair<T, U>& p) { os << p.first << ' ' << p.second; return os; }

template<typename T> istream& operator>>(istream& is, vector<T>& v) { for(auto&& x : v) is >> x; return is; }
template<typename T> ostream& operator<<(ostream& os, const vector<T>& v) { for(size_t i=0;i<v.size();++i) { if(i>0) os << ' '; os << v[i]; } return os; }
template<typename T> ostream& operator<<(ostream& os, const vector<vector<T>>& v) { for(size_t i=0;i<v.size();++i) { if(i>0) os << '\n'; os << v[i]; } return os; }

template<typename T, typename U> inline constexpr auto& operator+=(pair<T, U>& p, const pair<T, U>& q) { p.first += q.first; p.second += q.second; return p; }
template<typename T, typename U> inline constexpr auto& operator-=(pair<T, U>& p, const pair<T, U>& q) { p.first -= q.first; p.second -= q.second; return p; }
template<typename T, typename U> inline constexpr auto operator+(const pair<T, U>& p, const pair<T, U>& q) { auto r = p; r += q; return r; }
template<typename T, typename U> inline constexpr auto operator-(const pair<T, U>& p, const pair<T, U>& q) { auto r = p; r -= q; return r; }

template<typename T> inline constexpr auto dot(const pair<T, T>& p, const pair<T, T>& q) { return p.first*q.first + p.second*q.second; }
template<typename T> inline constexpr long double norm(const pair<T, T>& p) { return sqrt(dot(p, p)); }

template <typename T> inline constexpr auto& operator+=(vector<T>& v, const vector<T>& w) { if(v.size() < w.size()) v.resize(w.size()); for(size_t i=0; i<w.size(); ++i) v[i] += w[i]; return v; }
template <typename T> inline constexpr auto& operator-=(vector<T>& v, const vector<T>& w) { if(v.size() < w.size()) v.resize(w.size()); for(size_t i=0; i<w.size(); ++i) v[i] -= w[i]; return v; }
template <typename T> inline constexpr auto operator+(const vector<T>& v, const vector<T>& w) { auto u = v; u += w; return u; }
template <typename T> inline constexpr auto operator-(const vector<T>& v, const vector<T>& w) { auto u = v; u -= w; return u; }

#define make_v(T,...) __make_v<T>(__VA_ARGS__)
template <typename T> auto __make_v(size_t a) { return vector<T>(a); }
template <typename T, typename... Ts> auto __make_v(size_t a, Ts... ts) { return vector<decltype(__make_v<T>(ts...))>(a, __make_v<T>(ts...)); }

template<typename T> struct has_itr {
  template<typename U> static auto check(U x) -> decltype(begin(x), end(x), true_type()); static false_type check(...);
  static bool constexpr value = decltype( check( declval<T>() ) )::value;
};
template <typename T, typename V, enable_if_t<!has_itr<T>::value, nullptr_t> = nullptr> void fill_v(T& t, const V& v) { t = v; }
template <typename T, typename V, enable_if_t< has_itr<T>::value, nullptr_t> = nullptr> void fill_v(T& t, const V& v) { for(auto&& e : t) fill_v(e, v); }

template <typename T, class F, enable_if_t<!has_itr<T>::value, nullptr_t> = nullptr> inline constexpr T& apply(T& t, const F& f) { return t = f(t); }
template <typename T, class F, enable_if_t< has_itr<T>::value, nullptr_t> = nullptr> inline constexpr auto& apply(T& t, const F& f) { for(auto&& e : t) apply(e, f); return t; }
template <typename T, class F, enable_if_t< has_itr<T>::value, nullptr_t> = nullptr> inline constexpr auto apply(T&& t, const F& f) { auto u = t; return apply(u, f); }

template <typename T> struct is_pair : false_type {}; template <typename T, typename U> struct is_pair<pair<T, U>> : true_type {};
template<typename T, enable_if_t< is_pair<T>::value, nullptr_t> = nullptr> inline constexpr T& dcr(T& t) { --t.ft; --t.sd; return t; }
template<typename T, enable_if_t<!is_pair<T>::value, nullptr_t> = nullptr> inline constexpr T& dcr(T& t) { return --t; }

struct InitIO {
  InitIO() {
    cin.tie(nullptr); cout.tie(nullptr);
    ios_base::sync_with_stdio(false);
    cout << fixed << setprecision(15);
    #ifndef LOCAL
      cerr.rdbuf(nullptr);
    #endif
  }

} initio;
#ifdef LOCAL
  #include <boost/type_index.hpp>
  #define trace(...) __trace(" " #__VA_ARGS__ ",", __VA_ARGS__)
  template <typename T> inline constexpr const char* __trace(const char* n, const T& t)
  { const char* c = strchr(n+1, ','); cerr.write(n, c-n) << " = " << t << "; // " << boost::typeindex::type_id_with_cvr<decltype(t)>().pretty_name() <<"\n"; return c; }
  template <typename T, typename... Ts> inline constexpr const char* __trace(const char* n, const T& t, const Ts&... ts)
  { const char* c = __trace(n, t); return __trace(c+1, ts...); }
#else
  #define trace(...)
#endif

#define in(T,...) __in<T>(__VA_ARGS__)
template<typename T> inline void sc(T& x) { cin >> x; }
template<typename T, typename... Ts> inline void sc(T& t, Ts&... ts) { sc(t); sc(ts...); }
template<typename T=mnt> inline T __in() { T x; sc(x); return x; }
template<typename T, typename... Ts> inline auto __in(Ts ...ts) { auto v = __make_v<T>(ts...); sc(v); return v; }
template<typename T> inline void out(const T& x) { cout << x; }
template<typename T, typename... Ts> inline void out(const T& t, const Ts&... ts) { out(t); out(' '); out(ts...); }
template<typename T> inline void err(const T& x) { cerr << x; }
template<typename T, typename... Ts> inline void err(const T& t, const Ts&... ts) { err(t); err(' '); err(ts...); }
template<typename... Ts> inline void print(const Ts&... ts) { out(ts...); out('\n'); }
template<typename... Ts> inline void debug(const Ts&... ts) { err(ts...); err('\n'); }
inline string input() { string str; cin >> ws; getline(cin, str); return str; }

inline constexpr auto yes(const bool& c) { return c ? "yes" : "no"; }
inline constexpr auto Yes(const bool& c) { return c ? "Yes" : "No"; }
inline constexpr auto YES(const bool& c) { return c ? "YES" : "NO"; }
inline constexpr auto Possible(const bool& c) { return c ? "Possible" : "Impossible"; }
inline constexpr auto POSSIBLE(const bool& c) { return c ? "POSSIBLE" : "IMPOSSIBLE"; }

constexpr int lim = 1e5+7;



signed main() {

  let N = in();
  var s = in(string, N);
  for(var t : s) sort(t);
  dict<> cnt;
  for(var t : s) ++cnt[t];
  mnt ans = 0;
  for(var c : cnt) ans += c.sd*(c.sd-1)/2;
  print(ans);

  return 0;
}


Measured memory_bytes: 9632