結果

問題 No.2389 Cheating Code Golf
ユーザー cleantted
提出日時 2023-07-19 20:40:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 173 ms / 2,000 ms
コード長 2,051 bytes
コンパイル時間 155 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 91,008 KB
最終ジャッジ日時 2024-09-19 19:06:40
合計ジャッジ時間 10,224 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 50
権限があれば一括ダウンロードができます

ソースコード

diff # 

# import bisect
import copy
import heapq
import itertools
import math
import operator
import random
import sys
from bisect import bisect, bisect_left, bisect_right, insort
from collections import Counter, deque
from fractions import Fraction
from functools import cmp_to_key, lru_cache, partial
from inspect import currentframe
from math import ceil, gcd, log10, pi, sqrt
from typing import Iterable, Iterator, List, Tuple, TypeVar, Union

# import pypyjit
# pypyjit.set_param('max_unroll_recursion=-1')
# import string
# import networkx as nx
input = sys.stdin.readline
sys.setrecursionlimit(10000000)
# mod = 10 ** 9 + 7
mod = 998244353
# mod = 1 << 128
# mod = 10 ** 30 + 1
INF = 1 << 61
DIFF = 10 ** -9
DX = [1, 0, -1, 0, 1, 1, -1, -1]
DY = [0, 1, 0, -1, 1, -1, 1, -1]

def read_values(): return tuple(map(int, input().split()))
def read_index(): return tuple(map(lambda x: int(x) - 1, input().split()))
def read_list(): return list(read_values())
def read_lists(N): return [read_list() for _ in range(N)]
def dprint(*values): print(*values, file=sys.stderr)
def dprint2(*values):
    names = {id(v): k for k, v in currentframe().f_back.f_locals.items()}
    dprint(", ".join(f"{names.get(id(value), '???')}={repr(value)}" for value in values))


def cmp(l, r):
    a1, b1, p1 = l
    a2, b2, p2 = r

    t1 = (a1 - b1) * a2 * b2 * (p2 - 1)
    t2 = (a2 - b2) * a1 * b1 * (p1 - 1)
    if t1 > t2:
        return -1
    elif t1 < t2:
        return 1
    return 0


def main():
    N, M = read_values()

    L = read_lists(N)
    L.sort(key=cmp_to_key(cmp))
    for i in range(N):
        L[i][1] = min(L[i][0], L[i][1])

    dp = [[0] * (M + 1) for _ in range(N + 1)] 

    for i in reversed(range(N)):
        a, b, p = L[i]
        for m in range(M + 1):
            r = 0
            t = 1
            for k in range(m):
                r += 1 / p * t * (1 / b + dp[i + 1][m - k])
                t *= (p - 1) / p
            r += t * (1 / a + dp[i + 1][0])
            dp[i][m] = r

    print(dp[0][M])


if __name__ == "__main__":
    main()
0