結果

問題 No.162 8020運動
ユーザー ああいい
提出日時 2022-03-02 23:11:34
言語 PyPy3
(7.3.15)
結果
MLE  
実行時間 -
コード長 1,993 bytes
コンパイル時間 651 ms
コンパイル使用メモリ 81,972 KB
実行使用メモリ 735,484 KB
最終ジャッジ日時 2024-07-16 13:38:00
合計ジャッジ時間 16,655 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample TLE * 1 MLE * 1 -- * 1
other -- * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

A = int(input())
P = list(map(int,input().split()))
P = [p / 100 for p in P]
N = 80 - A
#dp = [[0] * (1 << 14) for _ in range(N+1)]
#dp[0][-1] = 1
"""
dp = [0] * (1 << 14)
dp[-1] = 1
for j in range(N):
    nx = [0] * (1 << 14)
    for bit in range(1 << 14):
        T = bit
        while True:
            tmp = 1
            for i in range(14):
                mask = 1 << i
                if bit & mask == 0:
                    continue
                count = 0
                if bit & (mask >> 1):
                    count += 1
                if bit & (mask << 1):
                    count += 1
                if T & mask:
                    tmp *= 1 - P[count]
                else:
                    tmp *= P[count]
            nx[T] += dp[bit] * tmp
            if T == 0:break
            T -= 1
            T &= bit
    dp = nx
"""
n = 1 << 14
from collections import defaultdict
l = defaultdict(lambda :defaultdict(int))
for bit in range(n):
    T = bit
    while True:
        tmp = 1
        for i in range(14):
            mask = 1 << i
            if bit & mask == 0:
                continue
            count = 0
            if bit & (mask >> 1):
                count += 1
            if bit & (mask << 1):
                count += 1
            if T & mask:
                tmp *= 1 - P[count]
            else:
                tmp *= P[count]
        l[T][bit] += tmp
        if T == 0:break
        T -= 1
        T &= bit
dp = [0] * n
dp[-1] = 1
for _ in range(N):
    nx = [0] * n
    for bit in range(n):
        T = bit
        while True:
            nx[T] += dp[bit] * l[T][bit]
            if T == 0:break
            T -= 1
            T &= bit
    dp = nx
"""
def seki(x,y):
    l = [[0] * n for _ in range(n)]
    for i in range(n):
        for j in range(n):
            for k in range(n):
                l[i][j] += x[i][k] * y[k][j]
    return l
"""


ans = 0
for bit in range(n):
    num = bin(bit).count("1")
    ans +=num * dp[bit]
print(ans * 2)
0