#include <bits/stdc++.h>
#define f first
#define s second
using namespace std;
#define int long long
const int maxn = 1e6 + 69;
int a[maxn], b[maxn], c[maxn];
int ab[maxn], ac[maxn], bc[maxn], dp[maxn];
signed main()
{
   ios_base::sync_with_stdio(false);
   cin.tie(0);
   cout.tie(0);
   int tt = 1;
   //cin >> tt;
   while(tt--) {
		string s;
		cin >> s;
		int n = s.size();
		s = "#" + s;
		const int mod = 998244353;
		auto add = [&](int x, int y) {
			x += y;
			if (x >= mod)
				x -= mod;
			return x;
		};
		dp[0] = 1;
		for (int i = 1;i <= n;i++) {
			a[i] = a[i - 1] + (s[i] == 'I');
			b[i] = b[i - 1] + (s[i] == 'P');
			c[i] = c[i - 1] + (s[i] == 'C');
			ab[i] = a[i] - b[i];
			ac[i] = a[i] - c[i];
			bc[i] = b[i] - c[i];
		}
		map<int, vector<int>> z1, z2, z3;
		map<int, int> c1, c2, c3;
		for (int i = 0;i <= n;i++) {
			z1[ab[i]].push_back(3 * c[i] - i);
			z2[ac[i]].push_back(3 * b[i] - i);
			z3[bc[i]].push_back(3 * a[i] - i);
		}
		int p = 0;
		vector<vector<int>> f1, f2, f3;
		vector<vector<int>> r1, r2, r3;
		for (auto &x : z1) {
			sort(x.s.begin(), x.s.end());
			x.s.erase(unique(x.s.begin(), x.s.end()), x.s.end());
			f1.push_back(vector<int>(x.s.size() + 2, 0));
			r1.push_back(x.s);
			c1[x.f] = p++;
		}
		p = 0;
		for (auto &x : z2) {
			sort(x.s.begin(), x.s.end());
			x.s.erase(unique(x.s.begin(), x.s.end()), x.s.end());
			f2.push_back(vector<int>(x.s.size() + 2, 0));
			r2.push_back(x.s);
			c2[x.f] = p++;
		}
		p = 0;
		for (auto &x : z3) {
			sort(x.s.begin(), x.s.end());
			x.s.erase(unique(x.s.begin(), x.s.end()), x.s.end());
			f3.push_back(vector<int>(x.s.size() + 2, 0));
			r3.push_back(x.s);
			c3[x.f] = p++;
		}
		for (int i = 0;i <= n;i++) {
			int k1 = c1[ab[i]];
			auto u1 = lower_bound(r1[k1].begin(), r1[k1].end(), 3 * c[i] - i) - r1[k1].begin() + 1;
			for (int v = u1 - 1;v >= 1;v -= (v & -v))
				dp[i] = add(dp[i], f1[k1][v]);
				
			int k2 = c2[ac[i]];
			auto u2 = lower_bound(r2[k2].begin(), r2[k2].end(), 3 * b[i] - i) - r2[k2].begin() + 1;
			for (int v = u2 - 1;v >= 1;v -= (v & -v))
				dp[i] = add(dp[i], f2[k2][v]);
				
			int k3 = c3[bc[i]];
			auto u3 = lower_bound(r3[k3].begin(), r3[k3].end(), 3 * a[i] - i) - r3[k3].begin() + 1;
			for (int v = u3 - 1;v >= 1;v -= (v & -v))
				dp[i] = add(dp[i], f3[k3][v]);
			
			for (int v = u1;v < (int)f1[k1].size();v += (v & -v))
				f1[k1][v] = add(f1[k1][v], dp[i]);
			for (int v = u2;v < (int)f2[k2].size();v += (v & -v))
				f2[k2][v] = add(f2[k2][v], dp[i]);
			for (int v = u3;v < (int)f3[k3].size();v += (v & -v))
				f3[k3][v] = add(f3[k3][v], dp[i]);
		}
		cout << dp[n];
   }
   return 0;
}