Problem:
Score : 300 points

Problem Statement

Snuke has a fair N-sided die that shows the integers from 1 to N with equal probability and a fair coin. He will play the following game with them:

Throw the die. The current score is the result of the die.

As long as the score is between 1 and K-1 (inclusive), keep flipping the coin. The score is doubled each time the coin lands heads up, and the score becomes 0 if the coin lands tails up.

The game ends when the score becomes 0 or becomes K or above. Snuke wins if the score is K or above, and loses if the score is 0.

You are given N and K. Find the probability that Snuke wins the game.

Constraints

1 ≤ N ≤ 10^5

1 ≤ K ≤ 10^5

All values in input are integers.

Input

Input is given from Standard Input in the following format:

N K

Output

Print the probability that Snuke wins the game. The output is considered correct when the absolute or relative error is at most 10^{-9}.

Sample Input 1

3 10

Sample Output 1

0.145833333333

If the die shows 1, Snuke needs to get four consecutive heads from four coin flips to obtain a score of 10 or above. The probability of this happening is \frac{1}{3} \times (\frac{1}{2})^4 = \frac{1}{48}.

If the die shows 2, Snuke needs to get three consecutive heads from three coin flips to obtain a score of 10 or above. The probability of this happening is \frac{1}{3} \times (\frac{1}{2})^3 = \frac{1}{24}.

If the die shows 3, Snuke needs to get two consecutive heads from two coin flips to obtain a score of 10 or above. The probability of this happening is \frac{1}{3} \times (\frac{1}{2})^2 = \frac{1}{12}.

Thus, the probability that Snuke wins is \frac{1}{48} + \frac{1}{24} + \frac{1}{12} = \frac{7}{48} \simeq 0.1458333333.

Sample Input 2

100000 5

Sample Output 2

0.999973749998

Solution in Rust (1.15.1) (judge verdict: Wrong Answer)

#![allow(dead_code)]
#![allow(non_snake_case)]
// Original macro is in https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8
// Thanks @tanakh :bow:
macro_rules! input {
    ($($r:tt)*) => {
        let stdin = std::io::stdin();
        let mut bytes = std::io::Read::bytes(std::io::BufReader::new(stdin.lock()));
        let mut next = || -> String {
            bytes
                .by_ref()
                .map(|r|r.unwrap() as char)
                .skip_while(|c|c.is_whitespace())
                .take_while(|c|!c.is_whitespace())
                .collect()
        };
        input_inner!{next, $($r)*}
    };
}
 
macro_rules! input_inner {
    ($next:expr) => {};
    ($next:expr, ) => {};
 
    ($next:expr, $var:ident : $t:tt $($r:tt)*) => {
        let $var = read_value!($next, $t);
        input_inner!{$next $($r)*}
    };
    
    ($next:expr, mut $var:ident : $t:tt $($r:tt)*) => {
        let mut $var = read_value!($next, $t);
        input_inner!{$next $($r)*}
    };
}
 
macro_rules! read_value {
    ($next:expr, ( $($t:tt),* )) => {
        ( $(read_value!($next, $t)),* )
    };
 
    ($next:expr, [ $t:tt ; $len:expr ]) => {
        (0..$len).map(|_| read_value!($next, $t)).collect::<Vec<_>>()
    };

    ($next:expr, [$t:tt]) => {
        {
            let len = read_value!($next, usize);
            (0..len).map(|_| read_value!($next, $t)).collect::<Vec<_>>()
        }
    };
 
    ($next:expr, chars) => {
        read_value!($next, String).chars().collect::<Vec<char>>()
    };
 
    ($next:expr, usize1) => {
        read_value!($next, usize) - 1
    };
 
    ($next:expr, $t:ty) => {
        $next().parse::<$t>().expect("Parse error")
    };
}

// for adhoc read (ex. https://atcoder.jp/contests/abc112/tasks/abc112_a)
fn read_line() -> String {
    let mut s = String::new();
    std::io::stdin().read_line(&mut s).unwrap();
    s.trim().to_string()
}

#[derive(Debug)]
struct Matrix<T=f32> {
    row: usize,
    col: usize,
    matrix: Vec<Vec<T>>,
}

impl Matrix {
    fn rank(&mut self) -> usize {
        let mut rank = 0;

        for col in 0..self.col {
            let mut pivot = 1000000;
            let mut ma = 1e-10;
            for row in rank..self.row {
                if self.matrix[row][col].abs() > ma {
                    ma = self.matrix[row][col].abs();
                    pivot = row;
                }
            }
            if pivot == 1000000 {
                continue;
            }
            self.matrix.swap(pivot, rank);

            let div = self.matrix[rank][col];
            for col2 in 0..self.col { self.matrix[rank][col2] /= div; }
            for row in 0..self.row {
                if row != rank && self.matrix[row][col].abs() > 1e-10 {
                    let fac = self.matrix[row][col];
                    for col2 in 0..self.col {
                        self.matrix[row][col2] -= self.matrix[rank][col2] * fac;
                    }
                }
            }
            rank += 1;
            // println!("---------------------");
            // println!("col {}", col);
            // println!("pivot {}, rank {}", pivot, rank);
            // println!("  swapped matrix {:?}", self.matrix);
        }
        return rank;
    }
}
fn gcd(a: i64, b: i64) -> i64 {
    if a < b { return gcd(b, a); }
    // let q = a/b;
    let r = a%b;
    if r == 0 { return b; }
    return gcd(b, r);
}
fn factorial(n: u64) -> u64 {
    if n == 1 { return 1; }
    let tmp: u64 = n % 1000_000_007;
    return tmp*(factorial(n-1)%1000_000_007)%1000_000_007;
}
// https://practice.contest.atcoder.jp/tasks/practice_1
fn main() {
    input! {
        N: i32, K: f32
    }
    let mut ans : f32 = 0.0;
    for k in 1..N+1 {
        if k >= K as i32 {
            ans += 1 as f32/N as f32;
        } else {
            let t = K/k as f32;
            let a: f32 = t.log2();
            let r = a.ceil();
            ans += (1 as f32/N as f32) * 0.5f32.powf(r);
            // println!("k {}, t {}, a {}, r {}, ans {}", k, t, a, r, ans);
        }
    }
    println!("{}", ans);
}