#include<bits/stdc++.h>
using namespace std;
#define V vector
#define pb push_back
#define fi first
#define se second
#define lb lower_bound
#define all1(x) (x).begin() + 1, (x).end()

typedef pair<int, int> pi;
typedef long long LL;
typedef pair<LL, LL> pl;

struct SA {
    int n;
    std::vector<int> sa, rk, ht;
    SA(const std::string &s) {
        n = s.size() - 1;
        sa.resize(n + 1);
        rk.resize(n + 1);
        ht.resize(n);
        std::iota(sa.begin() + 1, sa.end(), 1);
        std::sort(sa.begin() + 1, sa.end(), [&](int a, int b) {
            return s[a] < s[b];
        });

        rk[sa[1]] = 1;
        for (int i = 2; i <= n; ++i) {
            rk[sa[i]] = rk[sa[i - 1]] + (s[sa[i]] != s[sa[i - 1]]);
        }
        int k = 1;
        std::vector<int> tmp, cnt(n + 1);
        tmp.reserve(n + 1);
        while (rk[sa[n]] < n) {
            tmp.resize(1);
            for (int i = 1; i <= k; ++i) {
                tmp.push_back(n - k + i);
            }
            for (auto i : sa) {
                if (i > k) {
                    tmp.push_back(i - k);
                }
            }
            std::fill(cnt.begin() + 1, cnt.end(), 0);
            for (int i = 1; i <= n; ++i) {
                cnt[rk[i]] += 1;
            }
            for (int i = 2; i <= n; ++i) {
                cnt[i] += cnt[i - 1];
            }
            for (int i = n; i >= 1; --i) {
                sa[cnt[rk[tmp[i]]]--] = tmp[i];
            }
            std::swap(rk, tmp);
            rk[sa[1]] = 1;
            for (int i = 2; i <= n; ++i) {
                rk[sa[i]] =
                    rk[sa[i - 1]] + (tmp[sa[i - 1]] < tmp[sa[i]] ||
                                     sa[i - 1] + k == n + 1 ||
                                     tmp[sa[i - 1] + k] < tmp[sa[i] + k]);
            }
            k *= 2;
        }
        for (int i = 1, j = 0; i <= n; ++i) {
            if (rk[i] == 1) {
                j = 0;
            } else {
                for (j -= j > 0; i + j <= n && sa[rk[i] - 1] + j <= n &&
                                 s[i + j] == s[sa[rk[i] - 1] + j];)
                    ++j;
                ht[rk[i] - 1] = j;
            }
        }
    }
};


void solve() {
    int n; cin >> n;
    string s; cin >> s;
    s = " " + s;
    // string rs = s;
    // reverse(all1(rs));
    // SA sa(rs);
    V<int> cnt(26);
    for(int i = 1; i <= n; i++) {
        cnt[s[i] - 'a']++;
    }

    // V<LL> sub(n + 1);
    // map<int, int> mp;
    // LL cur = 1LL * (n + 1) * n / 2 - accumulate(all1(sa.ht), 0LL);
    // sub[1] = cur;
    // for(int i = 1; i <= n - 1; i++) {
    //     mp.insert({i, sa.ht[i]});
    // }
    // // for(int i = 1; i < n; i++) cout << sa.ht[i] << " \n"[i == n - 1];
    // for(int i = 2; i <= n; i++) {
    //     cout << i - 1 << endl;
    //     for(auto [x, y] : mp) cout << x << " " << y << endl;
    //     int x = sa.rk[i - 1];
    //     cur -= (n - i + 2);
    //     auto it = mp.lb(x);
    //     if(it == mp.end()) {
    //         cur += mp.rbegin()->se;
    //         mp.erase(prev(mp.end()));
    //         sub[i] = cur;
    //         continue;
    //     }

    //     cur += it->se;
    //     sub[i] = cur;
    //     if(it == mp.begin()) {
    //         mp.erase(it);
    //         continue;
    //     }

    //     int y = it->se;
    //     it--;
    //     cur += it->se;
    //     cur -= min(y, it->se);
    //     it->se = y;
    //     sub[i] = cur;
    //     ++it;
    //     mp.erase(it);
    // }

    // for(int i = n; i >= 1; i--) cout << sub[i] << " \n"[i == 1];

    // // cout << sub[n] << endl;
    
    // V<int> n_sub(n + 1);
    // sub.pb(0);
    // for(int i = 1; i <= n; i++) {
    //     n_sub[i] = sub[n + 1 - i] - sub[n + 1 - i + 1];
    // }

    LL ans = 0;
    for(int i = 0; i <= n; i++) {
        if(i) cnt[s[i] - 'a']--;
        LL add = 1;
        for(int j = 0; j < 26; j++) {
            if(i != n && j == s[i + 1] - 'a') continue;
            add += cnt[j];
        }
        ans += (add);
    }
    cout << ans - 1 << endl;

}
int main() {
    ios::sync_with_stdio(0), cin.tie(0);
    solve();
}