Submission #3412956

---

Source Code Expand

import bisect
import collections
import heapq
import itertools
import operator
import sys

import numpy as np


def main():
    n, d = map(int, sys.stdin.readline().split())
    ns = []
    while True:
        for i in [2, 3, 5]:
            div, mod = divmod(d, i)
            if mod == 0:
                ns.append(i)
                d = div
                break
        else:
            if d != 1:
                print("{:.8f}".format(0))
                return
            break

    dp = np.zeros((101, 61, 40, 30))
    dp[0][0][0][0] = 1
    for l in range(n):
        for i in range(59):
            for j in range(39):
                for k in range(29):
                    dp[l + 1][i][j][k] += dp[l][i][j][k]
                    dp[l + 1][i + 1][j][k] += dp[l][i][j][k]
                    dp[l + 1][i][j + 1][k] += dp[l][i][j][k]
                    dp[l + 1][i + 2][j][k] += dp[l][i][j][k]
                    dp[l + 1][i][j][k + 1] += dp[l][i][j][k]
                    dp[l + 1][i + 1][j + 1][k] += dp[l][i][j][k]

    ans = 0
    n2 = ns.count(2)
    n3 = ns.count(3)
    n5 = ns.count(5)
    for i in range(n2, 61):
        for j in range(n3, 40):
            for k in range(n5, 30):
                ans += dp[n][i][j][k]

    print("{:.8f}".format(float(ans / 6 ** n)))


if __name__ == "__main__":
    main()

Submission Info

Judge Result