結果

問題 No.2899 Taffy Permutation
ユーザー mattu34mattu34
提出日時 2024-09-24 15:03:31
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 720 ms / 2,000 ms
コード長 2,823 bytes
コンパイル時間 367 ms
コンパイル使用メモリ 82,916 KB
実行使用メモリ 124,192 KB
最終ジャッジ日時 2024-09-24 15:03:42
合計ジャッジ時間 10,170 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 147 ms
88,948 KB
testcase_01 AC 151 ms
89,060 KB
testcase_02 AC 147 ms
88,964 KB
testcase_03 AC 150 ms
89,092 KB
testcase_04 AC 149 ms
89,108 KB
testcase_05 AC 147 ms
88,904 KB
testcase_06 AC 147 ms
89,116 KB
testcase_07 AC 148 ms
89,140 KB
testcase_08 AC 145 ms
88,836 KB
testcase_09 AC 720 ms
124,016 KB
testcase_10 AC 666 ms
124,096 KB
testcase_11 AC 700 ms
124,192 KB
testcase_12 AC 653 ms
124,156 KB
testcase_13 AC 718 ms
124,148 KB
testcase_14 AC 681 ms
123,848 KB
testcase_15 AC 692 ms
123,784 KB
testcase_16 AC 686 ms
123,948 KB
testcase_17 AC 661 ms
123,912 KB
testcase_18 AC 146 ms
89,100 KB
testcase_19 AC 698 ms
124,116 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import *
import sys
import heapq
from heapq import heapify, heappop, heappush
import bisect
from bisect import bisect_left, bisect_right
import itertools
from functools import lru_cache
from types import GeneratorType
from fractions import Fraction
import math
import copy
import random

# import numpy as np

# sys.setrecursionlimit(int(1e7))
# @lru_cache(maxsize=None) # CPython特化
# @bootstrap # PyPy特化(こっちのほうが速い) yield dfs(), yield Noneを忘れずに


def bootstrap(f, stack=[]):  # yield
    def wrappedfunc(*args, **kwargs):
        if stack:
            return f(*args, **kwargs)
        else:
            to = f(*args, **kwargs)
            while True:
                if type(to) is GeneratorType:
                    stack.append(to)
                    to = next(to)
                else:
                    stack.pop()
                    if not stack:
                        break
                    to = stack[-1].send(to)
            return to

    return wrappedfunc


dxdy1 = ((0, 1), (0, -1), (1, 0), (-1, 0))  # 上下右左
dxdy2 = (
    (0, 1),
    (0, -1),
    (1, 0),
    (-1, 0),
    (1, 1),
    (-1, -1),
    (1, -1),
    (-1, 1),
)  # 8方向すべて
dxdy3 = ((0, 1), (1, 0))  # 右 or 下
dxdy4 = ((1, 1), (1, -1), (-1, 1), (-1, -1))  # 斜め
INF = float("inf")
_INF = 1 << 60
MOD = 998244353
mod = 998244353
MOD2 = 10**9 + 7
mod2 = 10**9 + 7
# memo : len([a,b,...,z])==26
# memo : 2^20 >= 10^6
# 小数の計算を避ける : x/y -> (x*big)//y  ex:big=10**9
# @:小さい文字, ~:大きい文字,None: 空の文字列
# ユークリッドの互除法：gcd(x,y)=gcd(x,y-x)
# memo : d 桁以下の p 進表記を用いると p^d-1 以下のすべての
#        非負整数を表現することができる
# memo : (X,Y) -> (X+Y,X−Y) <=> 点を原点を中心に45度回転し、√2倍に拡大
# memo : (x,y)のx正から見た偏角をラジアンで(-πからπ]：　math.atan2(y, x)
# memo : a < bのとき ⌊a⌋ ≦ ⌊b⌋

input = lambda: sys.stdin.readline().rstrip()
mi = lambda: map(int, input().split())
li = lambda: list(mi())
ii = lambda: int(input())
py = lambda: print("Yes")
pn = lambda: print("No")
pf = lambda: print("First")
ps = lambda: print("Second")


N = ii()
S = list(input())
dp = [[0] * (N + 1) for _ in range(N + 1)]
dp[0][0] = 1
acc_pre = [0] + [1] * N
for i in range(N):
    s = S[i]
    acc = [0]
    for j in range(N + 1):
        if s == "1":
            dp[i + 1][j] += dp[i][j] * (N - i - j)
        else:
            if j + 1 <= N:
                dp[i + 1][j] += acc_pre[j + 1]
            if j + 1 <= N:
                dp[i + 1][j] += dp[i][j + 1] * (j + 1)
        dp[i + 1][j] %= MOD
        acc.append((acc[-1] + dp[i + 1][j] % MOD))
    acc_pre = copy.deepcopy(acc)
# print(dp)
print(dp[-1][0])
0