import io
import sys
import bisect
import math
from itertools import permutations, combinations
from heapq import heappush, heappop
from collections import deque
from collections import defaultdict as dd
sys.setrecursionlimit(10**7+10)

# mod = 998244353
mod = 10**9+7

_INPUT = """\
aaaaa


"""



def main():
    T = input()
    N = len(T)
    dp = [-1 for _ in range(N+1)] # 区間l,rのときのケース数
    def dfs(l):
        if dp[l] != -1:
            return dp[l]
        else:
            r = N-l
            ret = 1
            for i in range(r-l):
                if l+i>=r-1-i:break
                if T[l:l+i+1] == T[r-1-i:r]:
                    ret += dfs(l+i+1)
                    ret %= mod
        dp[l] = ret
        return ret
    dfs(0)
    print(dp[0]%mod)

if __name__ == "__main__":
    # sys.stdin = io.StringIO(_INPUT)
    main()