#include<iostream>
#include<cstring>
#include<cstdio>
#define next dskjf

using namespace std;

const int mo=1e9+7;

int f[100001],next[100001][3],sum[100001],ans=0;
char s[100001];

int main()
{
	scanf("%s",s+1);int len=strlen(s+1);
	bool bo=false;
	sum[0]=0;sum[1]=(s[1]=='a')?1:2;
	for(int i=2;i<=len;i++) 
	{
		sum[i]=(sum[i-1]+((s[i]=='a')?1:2))%3;
		if(s[i]==s[i-1]) bo=true;
	}
	if(!bo) {cout<<1<<endl;return 0;}
	memset(next,-1,sizeof(next));
	for(int i=len;i>=0;i--)
	{
		for(int j=0;j<=2;j++) next[i][j]=next[i+1][j];
		if(i<len) next[i][sum[i+1]]=i+1;
	}
	f[0]=1;
	for(int i=0;i<=len;i++)
	{
		for(int j=1;j<=2;j++)
		{
			int yu=next[i][(sum[i]+j)%3]; 
			if(yu!=-1) f[yu]=(f[yu]+f[i])%mo;
		}
		if((i>0)&&(sum[i]==sum[len])) ans=(ans+f[i])%mo;
	}
	cout<<ans<<endl;
	return 0;
}