let rec f n p x b r s=if mod_float n p=0.then f(n/.p)p(x*.p)x r s
else if n>p*.p then f n(p+.1.+.s)1. 0.(r*.(x-.b))1.else Printf.printf"%0.f
"(r*.(x-.b)*.if n=1.then 1.else n-.1.);;while 1=1do
f(read_float())2. 1. 0. 1. 0.done