Problem: Nick likes strings very much, he likes to rotate them, sort them, rearrange characters within a string... Once he wrote a random string of characters a, b, c on a piece of paper and began to perform the following operations: 

  * to take two adjacent characters and replace the second character with the first one, 
  * to take two adjacent characters and replace the first character with the second one 



To understand these actions better, let's take a look at a string «abc». All of the following strings can be obtained by performing one of the described operations on «abc»: «bbc», «abb», «acc». Let's denote the frequency of a character for each of the characters a, b and c as the number of occurrences of this character in the string. For example, for string «abc»: |a| = 1, |b| = 1, |c| = 1, and for string «bbc»: |a| = 0, |b| = 2, |c| = 1. 

While performing the described operations, Nick sometimes got balanced strings. Let's say that a string is balanced, if the frequencies of each character differ by at most 1. That is  - 1 ≤ |a| - |b| ≤ 1,  - 1 ≤ |a| - |c| ≤ 1 и  - 1 ≤ |b| - |c| ≤ 1. 

Would you help Nick find the number of different balanced strings that can be obtained by performing the operations described above, perhaps multiple times, on the given string s. This number should be calculated modulo 51123987.

Input

The first line contains integer n (1 ≤ n ≤ 150) — the length of the given string s. Next line contains the given string s. The initial string can be balanced as well, in this case it should be counted too. The given string s consists only of characters a, b and c.

Output

Output the only number — the number of different balanced strings that can be obtained by performing the described operations, perhaps multiple times, on the given string s, modulo 51123987.

Examples

Input

4
abca


Output

7


Input

4
abbc


Output

3


Input

2
ab


Output

1

Note

In the first sample it is possible to get 51 different strings through the described operations, but only 7 of them are balanced: «abca», «bbca», «bcca», «bcaa», «abcc», «abbc», «aabc». In the second sample: «abbc», «aabc», «abcc». In the third sample there is only one balanced string — «ab» itself.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int INF = 0x3f3f3f3f;
int n;
string s;
const int MOD = 51123987;
long long C[200][200], dp[2][66][66][66][8];
void MMOD(long long &tmp) { tmp %= MOD; }
bool check(int A, int B) { return abs(A - B) <= 1; }
int main() {
  ios::sync_with_stdio(false);
  cin >> n >> s;
  s = '$' + s;
  for (int i = 0; i < 200; ++i) C[i][0] = 1;
  for (int i = 1; i <= 199; ++i)
    for (int j = 1; j <= 199; ++j)
      C[i][j] = C[i - 1][j - 1] + C[i - 1][j], C[i][j] %= MOD;
  dp[0][0][0][0][0] = 1;
  int now = 1, prev = 0;
  for (int i = 1; i <= n; ++i) {
    int Key = s[i] - 'a';
    for (int a = 0; a <= 58; ++a)
      for (int b = 0; b <= 58; ++b) {
        for (int c = 0; c <= 58; ++c) {
          for (int mask = 0; mask < 8; ++mask) {
            if (dp[prev][a][b][c][mask] == 0) continue;
            if (mask & (1 << Key)) {
              dp[now][a][b][c][mask] += dp[prev][a][b][c][mask];
              MMOD(dp[now][a][b][c][mask]);
            } else {
              if (Key == 0) {
                dp[now][a + 1][b][c][1 << Key] += dp[prev][a][b][c][mask];
                MMOD(dp[now][a + 1][b][c][1 << Key]);
              } else {
                if (Key == 1)
                  dp[now][a][b + 1][c][1 << Key] += dp[prev][a][b][c][mask],
                      MMOD(dp[now][a][b + 1][c][1 << Key]);
                else
                  dp[now][a][b][c + 1][1 << Key] += dp[prev][a][b][c][mask],
                      MMOD(dp[now][a][b][c + 1][1 << Key]);
              }
              dp[now][a][b][c][mask + (1 << Key)] += dp[prev][a][b][c][mask],
                  MMOD(dp[now][a][b][c][mask + (1 << Key)]);
            }
          }
        }
      }
    for (int a = 0; a < 60; ++a)
      for (int b = 0; b < 60; ++b)
        for (int c = 0; c < 60; ++c)
          for (int mask = 0; mask < 8; ++mask) dp[prev][a][b][c][mask] = 0;
    swap(now, prev);
  }
  vector<pair<pair<int, int>, int> > pre;
  for (int i = 0; i <= 60; ++i)
    for (int j = 0; j <= 60; ++j)
      for (int k = 0; k <= 60; ++k) {
        if (i + j + k == n && check(i, j) && check(k, j) && check(i, k)) {
          pre.push_back(make_pair(make_pair(i, j), k));
        }
      }
  long long res = 0;
  for (auto it : pre) {
    for (int a = 0; a < 60; ++a)
      for (int b = 0; b < 60; ++b)
        for (int c = 0; c < 60; ++c)
          for (int mask = 0; mask < 8; ++mask) {
            int fa, fb, fc;
            fa = it.first.first;
            fb = it.first.second;
            fc = it.second;
            if (fa > 0 && a == 0) continue;
            if (fb > 0 && b == 0) continue;
            if (fc > 0 && c == 0) continue;
            long long tmp = 1;
            if (fa > 0) {
              tmp *= C[fa - 1][a - 1];
              MMOD(tmp);
            }
            if (fb > 0) {
              tmp *= C[fb - 1][b - 1];
              MMOD(tmp);
            }
            if (fc > 0) {
              tmp *= C[fc - 1][c - 1];
              MMOD(tmp);
            }
            res += tmp * dp[prev][a][b][c][mask] % MOD;
            MMOD(res);
          }
  }
  cout << res << endl;
  return 0;
}