結果

問題 No.1554 array_and_me
ユーザー tktk_snsntktk_snsn
提出日時 2021-06-18 23:06:19
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 371 ms / 2,000 ms
コード長 4,316 bytes
コンパイル時間 469 ms
コンパイル使用メモリ 87,188 KB
実行使用メモリ 102,128 KB
最終ジャッジ日時 2023-09-05 00:05:47
合計ジャッジ時間 11,475 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 80 ms
77,500 KB
testcase_01 AC 181 ms
80,532 KB
testcase_02 AC 178 ms
79,996 KB
testcase_03 AC 177 ms
80,656 KB
testcase_04 AC 179 ms
80,312 KB
testcase_05 AC 180 ms
80,420 KB
testcase_06 AC 371 ms
101,736 KB
testcase_07 AC 357 ms
102,044 KB
testcase_08 AC 357 ms
101,972 KB
testcase_09 AC 364 ms
101,876 KB
testcase_10 AC 357 ms
102,128 KB
testcase_11 AC 144 ms
99,700 KB
testcase_12 AC 140 ms
99,696 KB
testcase_13 AC 143 ms
99,660 KB
testcase_14 AC 141 ms
99,560 KB
testcase_15 AC 141 ms
99,684 KB
testcase_16 AC 141 ms
80,048 KB
testcase_17 AC 146 ms
79,564 KB
testcase_18 AC 143 ms
80,052 KB
testcase_19 AC 141 ms
79,632 KB
testcase_20 AC 143 ms
79,628 KB
testcase_21 AC 205 ms
81,068 KB
testcase_22 AC 204 ms
80,588 KB
testcase_23 AC 209 ms
80,696 KB
testcase_24 AC 207 ms
80,816 KB
testcase_25 AC 207 ms
80,760 KB
testcase_26 AC 272 ms
90,112 KB
testcase_27 AC 270 ms
90,516 KB
testcase_28 AC 277 ms
90,120 KB
testcase_29 AC 270 ms
89,940 KB
testcase_30 AC 271 ms
90,176 KB
testcase_31 AC 267 ms
90,276 KB
testcase_32 AC 264 ms
90,104 KB
testcase_33 AC 268 ms
90,664 KB
testcase_34 AC 269 ms
90,680 KB
testcase_35 AC 260 ms
90,204 KB
testcase_36 AC 246 ms
90,576 KB
testcase_37 AC 251 ms
90,140 KB
testcase_38 AC 248 ms
90,476 KB
testcase_39 AC 248 ms
90,536 KB
testcase_40 AC 251 ms
91,012 KB
testcase_41 AC 183 ms
80,052 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.buffer.readline
sys.setrecursionlimit(10 ** 7)
mod = 998244353
U = 10 ** 5 + 10
fact = [1] * (U+1)
inv = [1] * (U+1)
for i in range(1, U+1):
    fact[i] = fact[i-1] * i % mod
inv[U] = pow(fact[i], mod-2, mod)
for i in reversed(range(U)):
    inv[i] = inv[i+1] * (i + 1) % mod


class SegTree(object):
    def __init__(self, N, op_data, u_data):
        self._n = N
        self.log = (N-1).bit_length()
        self.size = 1 << self.log

        self.op = op_data
        self.e = u_data

        self.data = [u_data] * (2 * self.size)
        # self.len = [1] * (2 * self.size)

    def _update(self, i):
        self.data[i] = self.op(self.data[i << 1], self.data[i << 1 | 1])

    def initialize(self, arr=None):
        """ segtreeをarrで初期化する。len(arr) == Nにすること """
        if arr:
            for i, a in enumerate(arr, self.size):
                self.data[i] = a
        for i in reversed(range(1, self.size)):
            self._update(i)
            # self.len[i] = self.len[i << 1] + self.len[i << 1 | 1]

    def update(self, p, x):
        """ data[p] = x とする (0-indexed)"""
        p += self.size
        self.data[p] = x
        for i in range(1, self.log + 1):
            self._update(p >> i)

    def get(self, p):
        """ data[p]を返す """
        return self.data[p + self.size]

    def prod(self, l, r):
        """
        op_data(data[l], data[l+1], ..., data[r-1])を返す (0-indexed)
        """
        sml = self.e
        smr = self.e
        l += self.size
        r += self.size

        while l < r:
            if l & 1:
                sml = self.op(sml, self.data[l])
                l += 1
            if r & 1:
                r -= 1
                smr = self.op(self.data[r], smr)
            l >>= 1
            r >>= 1
        return self.op(sml, smr)

    def all_prod(self):
        """ op(data[0], data[1], ... data[N-1])を返す """
        return self.data[1]

    def max_right(self, l, func):
        """
        func(l, l+1, ..., r-1) = True,
        func(l, l+1, ..., r-1, r) = Falseとなる r を返す
        """
        if l == self._n:
            return self._n
        l += self.size
        sm = self.e
        while True:
            while l % 2 == 0:
                l >>= 1
            if not func(self.op(sm, self.data[l])):
                while l < self.size:
                    l <<= 1
                    if func(self.op(sm, self.data[l])):
                        sm = self.op(sm, self.data[l])
                        l += 1
                return l - self.size
            sm = self.op(sm, self.data[l])
            l += 1
            if (l & -l) == l:
                break
        return self._n

    def min_left(self, r, func):
        """
        func(     l, l+1, ..., r-1) = True,
        func(l-1, l, l+1, ..., r-1) = Falseとなる l を返す
        """
        if r == 0:
            return 0
        r += self.size
        sm = self.e
        while True:
            r -= 1
            while r > 1 and r & 1:
                r >>= 1
            if not func(self.op(self.data[r], sm)):
                while r < self.size:
                    r = r << 1 | 1
                    if func(self.op(self.data[r], sm)):
                        sm = self.op(self.data[r], sm)
                        r -= 1
                return r + 1 - self.size
            sm = self.op(self.data[r], sm)
            if (r & -r) == r:
                break
        return 0


def test():
    N, K = map(int, input().split())
    A = list(map(int, input().split()))
    S = sum(A) % mod
    S_inv = pow(S, mod-2, mod)

    u = (-1, -1, -1)

    def op(L, R):
        if L[0] == -1:
            return R
        if R[0] == -1:
            return L
        if L[0] * R[1] > L[1] * R[0]:
            return L
        return R

    count = [0] * N
    seg = SegTree(N, op, u)
    seg.initialize([(a, 1, i) for i, a in enumerate(A)])
    for _ in range(K):
        X, Y, idx = seg.all_prod()
        count[idx] += 1
        seg.update(idx, (X, Y+1, idx))
    res = 1
    for i in range(N):
        res = res * pow(A[i] * S_inv % mod, count[i], mod)
        res = res * inv[count[i]] % mod
    res *= fact[K]
    res %= mod
    return res


T = int(input())
for _ in range(T):
    print(test())
0