from math import gcd, log10, sqrt
from sympy import isprime
def expected_value(n, k):
S = 0
for x in range(2, int(sqrt(n)) + 1):
S += (n - k) % x
S -= (n + k) % x
return 1 + S/(2*k)
while True:
n = int(input("n="))
k = int(input("k="))
E = expected_value(n, k)
print("EV=", E)