from fractions import Fraction
import sys

def chmax(a, b):
    if a < b:
        a = b
        return True
    return False

def chmin(a, b):
    if b < a:
        a = b
        return True
    return False

INTINF = 1000001000
INTMAX = 2147483647
LLMAX = 9223372036854775807
LLINF = 1000000000000000000

def main():
    N = int(input().strip())
    memo = {}

    memo[1] = Fraction(0)

    def solve(n):
        if n in memo:
            return memo[n]

        divsum = Fraction(0)
        d = 2
        while d * d <= n:
            if n % d == 0:
                divsum += solve(d)

                if d != n // d:
                    divsum += solve(n // d)
            d += 1

        ret = Fraction(n, n - 1) + divsum / Fraction(n - 1)
        memo[n] = ret
        return ret

    result = solve(N)
    print(f"{float(result):.10f}")

if __name__ == "__main__":
    main()