#include #define Mod 1000000007 long long div_mod(long long x, long long y, long long z) { if (x % y == 0) return x / y; else return (div_mod((1 + x / y) * y - x, (z % y), y) * z + x) / y; } int main() { char S[200001]; scanf("%s", S); int i, j, l, pos[26][200001] = {}, k[26] = {}; long long pow[200001], pow_inv[200001]; for (l = 0; S[l] != 0; l++) pos[S[l] - 'a'][k[S[l] - 'a']++] = l; for (i = 1, pow[0] = 1; i <= l; i++) pow[i] = pow[i-1] * 2 % Mod; for (i = l - 1, pow_inv[l] = div_mod(1, pow[l], Mod); i >= 0; i--) pow_inv[i] = pow_inv[i+1] * 2 % Mod; long long ans = l * pow[l-1] % Mod, tmp; for (i = 0; i < 26; i++) { for (j = 1, tmp = 0; j < k[i]; j++) { tmp = (tmp * pow_inv[pos[i][j] - pos[i][j-1]] % Mod + pow[l - (pos[i][j] - pos[i][j-1] + 1)]) % Mod; ans = (ans - tmp + Mod) % Mod; } } printf("%lld\n", ans); fflush(stdout); return 0; }