#include <bits/stdc++.h>
#define il inline
#define ll long long
#define int ll
#define Max 1000005
using namespace std;
il ll read()
{
	char c=getchar();
	ll x=0,f=1;
	while(c>'9'||c<'0')
	{
		if(c=='-') f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9')
	{
		x=x*10+c-'0';
		c=getchar();
	}
	return x*f;
}
#define Mod 998244353
struct Tree
{
	vector<int> t;
	int sz;
#define lowbit(x) x&-x
	Tree()
	{
		sz=0;
	}
	void resize(int n)
	{
		t.resize(n+1);
	}
	void insert(int x,int k)
	{
		sz=(sz+k)%Mod;
		for(int i=x;i<t.size();i+=lowbit(i)) t[i]=(t[i]+k)%Mod;
	}
	int qry(int x)
	{
		int res=0;
		for(int i=x;i;i-=lowbit(i))
		{
			res=(res+t[i])%Mod;
		}
		return (sz-res+Mod)%Mod;
	}
}t[3][Max*2];
int n,s1[3][Max],s2[3][Max];
char s[Max];
map<int,int> g[3];
int x[3][Max],y[3][Max];
vector<int> yy[3][Max*2];
char c[3]={'I','C','P'};
void precalc()
{
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<3;j++)
		{
			s1[j][i]=s1[j][i-1]+(s[i]==c[j]);
			s2[j][i]=s2[j][i-1]+(s[i]==c[(j+1)%3]);
			int w=s1[j][i]-s2[j][i];
			int k=i-3*s1[j][i];
			x[j][i]=w,y[j][i]=k;
			yy[j][w+Max].push_back(k);
		}
	}
	for(int j=0;j<3;j++)
		yy[j][Max].push_back(0);
//	puts("qwq1");
	for(int j=0;j<3;j++)
	{
		for(int i=0;i<2*Max;i++) 
			if(yy[j][i].size()) sort(yy[j][i].begin(),yy[j][i].end()),yy[j][i].erase(unique(yy[j][i].begin(),yy[j][i].end()),yy[j][i].end()),t[j][i].resize(yy[j][i].size());
	}
//	puts("qwq2");
}
int calc()
{
	for(int i=0;i<3;i++) 
	{
		int y=lower_bound(yy[i][Max].begin(),yy[i][Max].end(),0)-yy[i][Max].begin()+1;
		g[i][0]=1,t[i][0+Max].insert(y,1);
	}
	int ans=0;
	for(int i=1;i<=n;i++)
	{
		int res=0;
		for(int j=0;j<3;j++)
		{
			int w=s1[j][i]-s2[j][i];
			int k=i-3*s1[j][i];
			int x=w+Max,y=lower_bound(yy[j][x].begin(),yy[j][x].end(),k)-yy[j][x].begin()+1;
			int nw=(g[j][w]-t[j][x].qry(y-1)+Mod)%Mod;
			res=(res+nw)%Mod;
		//	cout<<c[j]<<' '<<c[(j+1)%3]<<' '<<nw<<" qwq\n";
		}
		ans=res;
		for(int j=0;j<3;j++)
		{
			int w=s1[j][i]-s2[j][i];
			int k=i-3*s1[j][i];
			int x=w+Max,y=lower_bound(yy[j][x].begin(),yy[j][x].end(),k)-yy[j][x].begin()+1;
			g[j][w]=(g[j][w]+res)%Mod;
			t[j][x].insert(y,res);
		}
	//	cout<<i<<' '<<res<<"\n";
	}
	return ans;
}
signed main()
{
	scanf("%s",s+1);
	n=strlen(s+1);
	precalc();
	cout<<calc()<<"\n";
}