#include using namespace std; struct cww{cww(){ios::sync_with_stdio(false);cin.tie(0);}}star; #define all(c) (c).begin(), (c).end() #define rall(c) (c).rbegin(), (c).rend() #define vv(type, c, m, n, i) vector> c(m, vector(n, i)) #define mat(type, c, m, n) vector> c(m, vector(n));for(auto& r:c)for(auto& i:r)cin>>i; #define rep(i,a,n) for(int i=(a), i##_len=(n); ii##_len; --i) #define each(x,y) for(auto &(x):(y)) #define var(type, ...)type __VA_ARGS__;Scan(__VA_ARGS__); template void Scan(T& t) { cin >> t; } templatevoid Scan(First& first, Rest&...rest) { cin >> first; Scan(rest...); } #define vec(type, c, n) vector c(n);for(auto& i:c) cin>>i; #define len(x) ((int)(x).size()) #define mp make_pair #define eb emplace_back #define em emplace #define pb pop_back #define fi first #define se second #define get(a, i) get(a) #define bit_count __builtin_popcount typedef long long ll; typedef long double ld; typedef vector vi; typedef vector vd; typedef vector vll; typedef vector vs; typedef vector vb; typedef pair pii; #define mini min_element #define maxi max_element #define sum accumulate #define chmax(a,b) if(ab)a=b template inline void print(T t){cout << t << "\n";} template inline void print(H h, T... t){cout << h << " ";print(t...);} inline ll comb(int y, int x) { if (x == 0 || y == x) return 1; if (x < 0 || y < x) return 0; ll t = 1; ll maxt = ll(1e15) + 1; rep(i, 1, min(x, y - x) + 1) { t = t * (y + 1 - i) / i; if (t > maxt) return maxt; } return t; } #define MAX 70 ll sai[11][MAX]; ll ika[11][MAX]; ll A[MAX], B[MAX]; void f2(int n) { rep(i, 1, n) { rep(j, i * 4, i * 6 + 1) { rep(k, 4, 7) { ika[i + 1][j + k] += ika[i][j] * ika[1][k]; } } } } void f(int n) { rep(i, 1, n) { rep(j, i, i * 6 + 1) { rep(k, 1, 7) { sai[i + 1][j + k] += sai[i][j] * sai[1][k]; } } } } double solve() { var(int, n, k); sai[1][1] = sai[1][2] = sai[1][3] = sai[1][4] = sai[1][5] = sai[1][6] = 1; ika[1][4] = ika[1][5] = ika[1][6] = 2; f(n); if (k != 0) { f2(k); rep(i, n - k, (n - k) * 6 + 1) { rep(j, k * 4, k * 6 + 1) { A[i + j] += sai[n - k][i] * ika[k][j]; } } rrep(i, n * 6, n - 1) { B[i] = B[i + 1] + A[i]; } } else { rrep(i, n * 6, n - 1) { B[i] = B[i + 1] + sai[n][i]; } } ll cnt = 0; rep(i, n, n * 6 + 1) { cnt += sai[n][i] * B[i + 1]; } return cnt / pow(6, n * 2); } int main() { cout << setprecision(12) << fixed; print(solve()); return 0; }