#include <iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iomanip>
#include<algorithm>
#include<deque>
#define ll long long
#define mod 998244353
#define sc(A) scanf("%d",&A)
#define rep(A,B,C) for(int C=A;C<=B;++C)
#define fep(A,B,C) for(int C=A;C>=B;--C)
#define put(x) printf("%d\n",x)
#define db double
using namespace std;
const int MAXN=1000010,maxn=77;
int n,x,maxx=62;
int f[MAXN][maxn];
int g[MAXN][maxn];
int c1[maxn][maxn],c2[maxn][maxn];
char s[MAXN];
int cnt[maxn];
int a[MAXN];
int v[MAXN];


int mux(int x,int y){return x-y<0?x-y+mod:x-y;}
int modify(char a)
{
	if(a>='a'&&a<='z')return a-'a'+1;
	if(a>='A'&&a<='Z')return a-'A'+27;
	return a-'0'+53;
}

int add(int x,int y){return x+y>=mod?x+y-mod:x+y;}

signed main() 
{
	//freopen("1.in","r",stdin);
	scanf("%s",s+1);
	n=strlen(s+1);
	fep(n,1,i)
	{
		int x=modify(s[i]);
		a[i]=x;
		rep(1,maxx,j)
		{
			if(j==x)continue;
			f[x][j]=cnt[j];
			g[i][j]=c2[x][j];
			c1[x][j]=add(c1[x][j],f[x][j]);
			c2[j][x]=add(c2[j][x],c1[j][x]);
		}
		++cnt[x];
	}

	memset(cnt,0,sizeof(cnt));
	int ans=0;
	int cnt1=0;
	rep(1,n,i)
	{
		rep(1,maxx,j)
		{
			if(j==x)continue;
			int res=i-1-cnt[j];
			int ww=(ll)res*(res-1)/2%mod;
			ww=mux(ww,cnt1);
			ww=add(ww,v[a[i]]);
			ww=add(ww,v[j]);
			ans=add(ans,(ll)ww*g[i][j]%mod);
		}
		++cnt[a[i]];
		cnt1=mux(cnt1,v[a[i]]);
		v[a[i]]=(ll)cnt[a[i]]*(cnt[a[i]]-1)/2%mod;
		cnt1=add(cnt1,v[a[i]]);
	}

	put(ans);
	return 0;
}