////////////////////////////// // Check before you submit. // #define _ATCODER_LIBRARY const long long MOD = 1e9+7; // const long long MOD = 998244353; ///////////////////////////// #include using namespace std; #include using namespace boost::multiprecision; #ifdef _ATCODER_LIBRARY #include using namespace atcoder; #endif // _ATCODER_LIBRARY const long long INF = 1LL << 60; const double PI = acos(-1); using ll = long long; using P = pair; #define FOR(i,a,b) for (ll i=(a);i<(ll)(b);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() #define SUM(v) accumulate(ALL(v),0ll) templateistream& operator>>(istream&i,vector&v){REP(j,v.size())i>>v[j];return i;} templatestring join(vector&v){stringstream s;REP(i,v.size())s<<' '<ostream& operator<<(ostream&o,vector&v){if(v.size())o<string join(vector>&vv){string s="\n";REP(i,vv.size())s+=join(vv[i])+"\n";return s;} templateostream& operator<<(ostream&o,vector>&vv){if(vv.size())o<istream& operator>>(istream&i,pair&v){return i>>v.first>>v.second;} templateostream& operator<<(ostream&o,pair&v){return o<T up(T a, T b){assert(b);return (a+b-1)/b;} templatebool eq(A const&... a){auto t={a...};assert(t.size());auto tar=*t.begin();for(const auto&e:t)if(tar!=e)return false;return true;} templatebool chmin(T &a, T b){if(a>b){a=b;return false;}return true;} templatebool chmax(T &a, T b){if(abool chmax(T &a, initializer_listl){return chmax(a,max(l));} templatebool chmin(T &a, initializer_listl){return chmin(a,min(l));} ////////////////////////////////////////////////////////////////// // My Library ////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////// // Contest Code ////////////////////////////////////////////////////////////////// vector normal{1, 2, 3, 4, 5, 6}; vector ikasama{4, 4, 5, 5, 6, 6}; double dp1[110][1000]{}; double dp2[110][1000]{}; int main(int argc, char **argv) { init_init_init(); ll N, K; cin >> N >> K; dp1[0][0] = dp2[0][0] = 1.f; REP(i, N) { REP(sum, 1000) { REP(j, 6) { if (sum-normal[j]>=0) dp1[i+1][sum] += dp1[i][sum-normal[j]] / 6; } REP(j, 6) { if (i < K) { if (sum-ikasama[j]>=0) dp2[i+1][sum] += dp2[i][sum-ikasama[j]] / 6; } else { if (sum-normal[j]>=0) dp2[i+1][sum] += dp2[i][sum-normal[j]] / 6; } } } } double res{0.f}; REP(i, 1000) for (ll j = i+1; j<1000; ++j) res += dp1[N][i] * dp2[N][j]; std::cout << res << std::endl; }