#include <bits/stdc++.h>
using namespace std;
int n, i, j, x, sol, ok;
char s[100005];
int sum[100005], prv[100005][3], d[100005];
int main() {
  cin >> s + 1;
  n = strlen(s + 1);
  prv[0][0] = prv[0][1] = prv[0][2] = -1;
  for (i = 1; i <= n; i++) {
    sum[i] = (sum[i - 1] + s[i] - 'a' + 1) % 3;
    if (s[i] == s[i - 1]) {
      ok = 1;
    }
    prv[i][0] = prv[i - 1][0];
    prv[i][1] = prv[i - 1][1];
    prv[i][2] = prv[i - 1][2];
    prv[i][sum[i - 1]] = i - 1;
  }
  if (ok == 0) {
    cout << 1;
    return 0;
  }
  d[0] = 1;
  for (i = 1; i <= n; i++) {
    for (j = 1; j <= 2; j++) {
      x = (sum[i] - j + 3) % 3;
      if (prv[i][x] != -1) {
        d[i] = (d[i] + d[prv[i][x]]) % 1000000007;
      }
    }
    if (sum[i] == sum[n]) {
      d[i] = (d[i] + 1) % 1000000007;
    }
  }
  cout << (d[n] - 1 + 1000000007) % 1000000007;
  ;
}
