#!/usr/bin/env python3
c,t,r = [int(i) for i in input().split()]
p = float(input())

DP = [0]
for x in range(1,c+1):
    DP.append(min([t + DP[x-i] + (r+1/p)*((1-p)**-i - 1) for i in range(1,x+1)]))
#print(*[round(i,2) for i in DP])
print(DP[c])