#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<iomanip>
#include<math.h>
#include<complex>
#include<queue>
#include<deque>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<functional>
#include<assert.h>
#include<numeric>
using namespace std;
#define REP(i,m,n) for(int i=(int)(m) ; i < (int) (n) ; ++i )
#define rep(i,n) REP(i,0,n)
typedef long long ll;
typedef pair<int,int> pint;
typedef pair<ll,int> pli;
const int inf=1e9+7;
const ll longinf=1LL<<60 ;
const ll mod=1e9+7 ;


struct Xor128{
    unsigned x,y,z,w;
 
    Xor128(unsigned w_=88675123):x(123456789),y(362436069),z(521288629),w(w_){};
 
    inline unsigned xor128(){
        unsigned t;
        t = x^(x<<11);
        x = y;
        y = z;
        z = w;
        return w = (w^(w>>19))^(t^(t>>8));
    }
 
    int nextInt(int x,int y){
		if(x>y)swap(x,y);
        return xor128()%(y-x)+x;
    }
 
    double nextDouble(double a,double b){
        return (double)(xor128()&0xffff)/0xffff*(b-a)+a;
    }
 
};
int main(){
    string s;
    cin>>s;
    ll n=s.size();
    auto rnd=Xor128();
    set<ll> st;
    ll sum=0;
    ll cumsum[100001][26];
    rep(i, 100000) rep(j, 26) cumsum[i + 1][j] = cumsum[i][j] + (s[i] - 'a' == j ? 1 : 0);
    rep(i,100000){
        int l=rnd.nextInt(0, n);
        int r=rnd.nextInt(1, n+1);
        if(l==r)continue;
        if(l>r)swap(l,r);
        if(st.count(n*r+l))continue;
        st.insert(n*r+l);
        int cnt=0;
        rep(j, 26) if (cumsum[r][j] - cumsum[l][j] > 0) ++ cnt;
        sum+=cnt;
        if(st.size()==n*(n+1)/2)break;
    }
    cout<<1.0*sum/st.size()<<endl;

    return 0;
}