S=input()
L=len(S)

M=[[float("inf")]*(L+1) for i in range(26)]

for i in range(L-1,-1,-1):
    for j in range(26):
        if ord(S[i])-97==j:
            M[j][i]=i
        else:
            M[j][i]=M[j][i+1]

ANS=0

for j in range(26):
    for i in range(L):
        k=M[j][i]-i
        if k==float("inf"):
            continue
        ANS+=(L-i)-k

print(ANS/(L*(L+1)/2))