#include <bits/stdc++.h>
#define cs const
#define pb push_back
using namespace std; 
cs int N = 5e4 + 5; 
cs int mod = 998244353; 

int add(int a, int b) {
    return a + b >= mod ? a + b - mod : a + b;
}
int dec(int a, int b) {
    return a - b < 0 ? a - b + mod : a - b;
}
int mul(int a, int b) {
    return 1ll * a * b % mod;
}
void Add(int &a, int b){
    a = add(a, b);
}
void Dec(int &a, int b) {
    a = dec(a, b);
}
void Mul(int &a, int b) {
    a = mul(a, b);
}
int ksm(int a, int b) {
    int c = 1; 
    for(; b; b >>= 1, Mul(a, a))
    if(b & 1) Mul(c, a);
    return c; 
}

int F(char c) {
    if(c == 'Y' || c == 'R') return 1; 
    if(c == 'P' || c == 'B') return 2; 
    return 3; 
}

int n; 
char S[N * 2];
int dp[2][N / 2][2][2];

int main() {
    #ifdef zqj
    freopen("1.in","r",stdin);
    #endif
    cin >> n; 
    scanf("%s", S);
    int now = 0, nxt = 1; 
    dp[0][0][0][0] = dp[0][0][0][1] = 1; 
    S[n * 2] = S[0];
    for(int i = 1; i <= n; i++) {
        int a = F(S[i * 2 - 1]), b = F(S[i * 2]), c; 
        for(int j = min(i - 1, n / 2); j >= 0; j--) 
        for(int o : {0, 1}) for(int e : {0, 1}) dp[nxt][j][o][e] = 0; 
        for(int j = min(i - 1, n / 2); j >= 0; j--)  {
            for(int o : {0, 1})
            for(int e : {0, 1}) 
            if(c = dp[now][j][o][e]) {
                if(b & (e + 1)) Add(dp[nxt][j][o][e], c + (a == 3 ? c : 0));
                if(j == 0) {
                    assert(o == 0);
                    if(b & ((e ^ 1) + 1)) {
                        if(a & 1) Add(dp[nxt][j + 1][0][e ^ 1], c);
                        if(a & 2) Add(dp[nxt][j + 1][1][e ^ 1], c);
                    }
                } else {
                    if(b & ((e ^ 1) + 1)) {
                        if(a & 1) {
                            if(o == 0) Add(dp[nxt][j - 1][j == 1 ? 0 : 1][e ^ 1], c);
                            else Add(dp[nxt][j + 1][0][e ^ 1], c);
                        }
                        if(a & 2) {
                            if(o == 0) Add(dp[nxt][j + 1][1][e ^ 1], c); 
                            else Add(dp[nxt][j - 1][0][e ^ 1], c);
                        }
                    }
                }
            }
        } swap(now, nxt);
        // cout<<a<<' '<<b<<endl;
        // for(int j=0; j<=i;j++)
        // for(int k:{0,1})
        // for(int l:{0,1}){
        //     if(dp[i][j][k][l])cerr<<i<<' '<<j<<' '<<k<<' '<<l<<' '<<dp[i][j][k][l]<<endl;
        // }
    } int sum = 0;
    for(int i = 0; i < 2; i++)
    for(int j = 0; j < 2; j++)
    Add(sum, dp[now][0][i][j]);
    cout << sum << '\n';
    return 0; 
}