#include <bits/stdc++.h>
using namespace std;

#define int long long

const int N=1e6+5;
const int Mod=1e9+7;
const int inv2=(Mod+1)/2;

char c[105], s[N];
int sum1[65][65], sum2[65][65], sum3[65][65], id[256], totcnt[256], cnt[256];

void add(int& x,int y)
{
	(x += y) %= Mod;
	if(x>=Mod)
		x -= Mod;
}

signed main()
{
	scanf("%s",s+1);
	int n=strlen(s+1);
	for(int i=1;i<=n;i++)
		totcnt[s[i]]++;

	for(int i=1;i<=26;i++)
		c[i]=i+'a'-1;

	for(int i=1;i<=26;i++)
		c[26+i]=i+'A'-1;

	for(int i=1;i<=10;i++)
		c[52+i]=i+'0'-1;

	int k=62;
	for(int i=1;i<=k;i++)
		id[c[i]]=i;

	int sum=0;
	for(int i=1;i<=k;i++)
		(sum += totcnt[c[i]]*totcnt[c[i]]) %= Mod;

	int ans=0;
	for(int i=n;i>=3;i--)
	{
		for(int j=1;j<=k;j++)
			add(sum3[id[s[i]]][j] , sum2[id[s[i]]][j]);

		for(int j=1;j<=k;j++)
			add(sum2[j][id[s[i]]] , sum1[j][id[s[i]]]);

		(sum -= totcnt[s[i]]*totcnt[s[i]]) %= Mod;
		cnt[s[i]]++, totcnt[s[i]]--;
		(sum += totcnt[s[i]]*totcnt[s[i]]) %= Mod;
		for(int j=1;j<=k;j++)
			if(j!=id[s[i]])
				add(sum1[id[s[i]]][j] , cnt[c[j]]);

		for(int j=1;j<=k;j++)
		{
			if(id[s[i]]==j)
				continue;

			int newsum=(sum-totcnt[s[i]]*totcnt[s[i]]-totcnt[c[j]]*totcnt[c[j]])%Mod;
			(ans += sum3[id[s[i]]][j]*(((i-1-totcnt[s[i]]-totcnt[c[j]])*(i-1-totcnt[s[i]]-totcnt[c[j]])-newsum)%Mod)) %= Mod;
		}
	}
	printf("%lld\n",(ans*inv2%Mod+Mod)%Mod);
	return 0;
}