#include using namespace std; using ll = long long; using VI = vector; using VV = vector; using VS = vector; // tourist set template string to_string(pair p); template string to_string(tuple p); template string to_string(tuple p); string to_string(const string& s) { return '"' + s + '"'; } string to_string(const char* s) { return to_string((string) s); } string to_string(bool b) { return (b ? "true" : "false"); } string to_string(vector v) { bool first = true; string res = "{"; for (int i = 0; i < static_cast(v.size()); i++) { if (!first) { res += ", "; } first = false; res += to_string(v[i]); } res += "}"; return res; } template string to_string(bitset v) { string res = ""; for (size_t i = 0; i < N; i++) { res += static_cast('0' + v[i]); } return res; } template string to_string(A v) { bool first = true; string res = "{"; for (const auto &x : v) { if (!first) { res += ", "; } first = false; res += to_string(x); } res += "}"; return res; } template string to_string(pair p) { return "(" + to_string(p.first) + ", " + to_string(p.second) + ")"; } template string to_string(tuple p) { return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " + to_string(get<2>(p)) + ")"; } template string to_string(tuple p) { return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " + to_string(get<2>(p)) + ", " + to_string(get<3>(p)) + ")"; } void debug_out() { cerr << '\n'; } template void debug_out(Head H, Tail... T) { cerr << " " << to_string(H); debug_out(T...); } #ifdef LOCAL #define debug(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__) #else #define debug(...) 42 #endif // tourist set end templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b()) void no(){p_no(); exit(0);} void yes(){p_yes(); exit(0);} const ll mod = 1e9 + 7; const ll inf = 1e18; const double PI = acos(-1); double dp[15][100]; // taro double dp2[15][100]; int main(){ cin.tie(0); ios::sync_with_stdio(false); // input ll N, K; cin >> N >> K; double p = (double)K/N; // イカサマのサイコロを使う確率 double q = 1-p; //普通のサイコロを使う確率 // taro dp[0][0] = 1; // rep(i, N){ // rep(pos, 70){ // // イカサマのサイコロ // FOR(v, 4, 7){ // dp[i+1][pos+v] += p * 1.0/3 * dp[i][pos]; // } // // 普通のサイコロ // FOR(v, 1, 7){ // dp[i+1][pos+v] += q * 1.0/6 * dp[i][pos]; // } // } // } FOR(i, 0, K){ rep(pos, 70){ // イカサマのサイコロ FOR(v, 4, 7){ dp[i+1][pos+v] += 1.0/3 * dp[i][pos]; } } } FOR(i, K, N){ rep(pos, 70){ // 普通のサイコロ FOR(v, 1, 7){ dp[i+1][pos+v] += 1.0/6 * dp[i][pos]; } } } for(int i=1; i<13; i++){ cout << dp[N][i] << endl; } // jiro dp2[0][0] = 1; rep(i, N){ rep(pos, 70){ // 普通のサイコロ FOR(v, 1, 7){ dp2[i+1][pos+v] += 1.0/6 * dp2[i][pos]; } } } debug(dp2[2][2], 1.0/36); double ans = 0; // i : jiro's value FOR(i, N, 61){ double prob = dp2[N][i]; FOR(j, i+1, 61){ ans += prob * dp[N][j]; } } cout << setprecision(20); cout << ans << endl; return 0; }