結果

問題 No.2068 Restricted Permutation
ユーザー 👑 rin204rin204
提出日時 2022-09-02 22:02:18
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 174 ms / 2,000 ms
コード長 2,381 bytes
コンパイル時間 345 ms
コンパイル使用メモリ 87,296 KB
実行使用メモリ 127,824 KB
最終ジャッジ日時 2023-08-10 04:04:31
合計ジャッジ時間 6,067 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 123 ms
107,416 KB
testcase_01 AC 123 ms
107,244 KB
testcase_02 AC 122 ms
107,220 KB
testcase_03 AC 125 ms
107,400 KB
testcase_04 AC 124 ms
107,104 KB
testcase_05 AC 127 ms
107,224 KB
testcase_06 AC 127 ms
107,244 KB
testcase_07 AC 126 ms
107,400 KB
testcase_08 AC 125 ms
107,056 KB
testcase_09 AC 123 ms
107,232 KB
testcase_10 AC 125 ms
107,240 KB
testcase_11 AC 124 ms
107,000 KB
testcase_12 AC 122 ms
107,440 KB
testcase_13 AC 137 ms
107,552 KB
testcase_14 AC 130 ms
107,564 KB
testcase_15 AC 130 ms
107,468 KB
testcase_16 AC 163 ms
117,972 KB
testcase_17 AC 174 ms
127,824 KB
testcase_18 AC 162 ms
127,092 KB
testcase_19 AC 164 ms
127,428 KB
testcase_20 AC 163 ms
127,264 KB
testcase_21 AC 168 ms
127,424 KB
testcase_22 AC 167 ms
127,472 KB
testcase_23 AC 171 ms
127,604 KB
testcase_24 AC 147 ms
117,012 KB
testcase_25 AC 148 ms
116,964 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

# https://atcoder.jp/contests/code-festival-2016-qualc/submissions/33924036

from bisect import bisect


class Bit:
    def __init__(self, n):
        self.size = n
        self.n0 = 1 << (n.bit_length() - 1)
        self.tree = [0] * (n + 1)
    
    def range_sum(self, l, r):
        return self.sum(r - 1) - self.sum(l - 1)
        
    def sum(self, i):
        i += 1
        s = 0
        while i > 0:
            s += self.tree[i]
            i -= i & -i
        return s
        
    def get(self, i):
        return self.sum(i) - self.sum(i - 1)
 
    def add(self, i, x):
        i += 1
        while i <= self.size:
            self.tree[i] += x
            i += i & -i
         
    def lower_bound(self, x):
        pos = 0
        plus = self.n0
        while plus > 0:
            if pos + plus <= self.size and self.tree[pos + plus] < x:
                x -= self.tree[pos + plus]
                pos += plus
            plus //= 2
        return pos

MOD = 998244353
N = 2000000
fact = [0 for _ in range(N)]
invfact = [0 for _ in range(N)]
fact[0] = 1
for i in range(1, N):
    fact[i] = fact[i - 1] * i % MOD

invfact[N - 1] = pow(fact[N - 1], MOD - 2, MOD)

for i in range(N - 2, -1, -1):
    invfact[i] = invfact[i + 1] * (i + 1) % MOD

def nCk(n, k):
    if k < 0 or n < k:
        return 0
    else:
        return (fact[n] * invfact[k] % MOD) * invfact[n - k] % MOD

def nPk(n, k):
    if k < 0 or n < k:
        return 0
    else:
        return fact[n] * invfact[n - k] % MOD

def nHk(n, k):
    if n == k == 0:
        return 1
    return nCk(n + k - 1, k)

n, k, x = map(int, input().split())
P = [0] * n
P[k - 1] = x
P = P[::-1]
tf = [False] * (n + 1)
ind = []
for i, p in enumerate(P):
    tf[p] = True

lst = []
for i in range(1, n + 1):
    if not tf[i]:
        lst.append(i)

inv2 = pow(2, MOD - 2, MOD)
ans = 0
l = len(lst)
invl = pow(l, MOD - 2, MOD)
bit = Bit(n + 1)
x = 0
t = 0
for i, p in enumerate(P):
    if p == 0:
        ans += (x * fact[i]) * inv2 % MOD
        ans %= MOD
        x += 1

        ans += (t * fact[i] % MOD) * invl % MOD
        ans %= MOD
    else:
        c = bit.sum(p)
        ans += fact[i] * c % MOD
        ans %= MOD
        bit.add(p, 1)
        
        c = bisect(lst, p)
        ans += (c * invl % MOD) * (fact[i] * x % MOD) % MOD
        ans %= MOD
        t += l - c

ans *= fact[l]
ans %= MOD
print(ans)
0