結果

問題 No.1712 Read and Pile
ユーザー chineristACchineristAC
提出日時 2021-05-06 22:33:50
言語 PyPy3
(7.3.15)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 2,690 bytes
コンパイル時間 159 ms
コンパイル使用メモリ 81,984 KB
実行使用メモリ 132,008 KB
最終ジャッジ日時 2024-09-17 16:54:55
合計ジャッジ時間 21,083 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 2
other AC * 6 WA * 34
権限があれば一括ダウンロードができます

ソースコード

diff # 

class BIT():
    def __init__(self,n,mod=None):
        self.BIT=[0]*(n+1)
        self.num=n
        self.mod = mod

    def query(self,idx):
        res_sum = 0
        while idx > 0:
            res_sum += self.BIT[idx]
            if self.mod:
                res_sum %= self.mod
            idx -= idx&(-idx)
        return res_sum

    #Ai += x O(logN)
    def update(self,idx,x):
        while idx <= self.num:
            self.BIT[idx] += x
            if self.mod:
                self.BIT[idx] %= self.mod
            idx += idx&(-idx)
        return

import sys,random,bisect
from collections import deque,defaultdict
from heapq import heapify,heappop,heappush
from itertools import permutations
from math import log,gcd

input = lambda :sys.stdin.readline().rstrip()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

mod = 998244353
i2 = pow(2,mod-2,mod)

def solve(N,M,A):
    if N==1:
        return M
    if N==2:
        res = M
        pre = 1
        Z = A.count(0)
        for i in range(M):
            if A[i]==0:
                res += i2
                res %= mod
            else:
                if pre!=0 and A[i]!=pre:
                    res += 1
                    res %= mod
                elif pre==0:
                    res += i2
                    res %= mod
            pre = A[i]
        return res

    p = ((N-2)*pow(N,mod-2,mod)) % mod
    ip = pow(p,mod-2,mod)
    pow_p = [1 for i in range(N+M+1)]
    pow_ip = [1 for i in range(N+M+1)]
    for i in range(1,N+M+1):
        pow_p[i] = (pow_p[i-1] * p) % mod
        pow_ip[i] = (pow_ip[i-1] * ip) % mod

    A = [N-i for i in range(N)] + A
    A = [-1] + A

    zero = [A[i]==0 for i in range(N+M+1)]
    for i in range(1,N+M+1):
        zero[i] += zero[i-1]

    lastappear = [N-i+1 for i in range(N+1)]
    res = 0
    bit_p = BIT(N+M,mod=998244353)
    bit_cnt = BIT(N+M)
    for i in range(1,N+1):
        bit_cnt.update(i,1)
        bit_p.update(i,pow_ip[zero[i]])

    for i in range(N+1,N+M+1):
        if A[i]==0:
            res += (N-1)
            res %= mod
        else:
            pre = lastappear[A[i]]
            k = bit_cnt.query(pre-1)

            res = (res + N-1) % mod
            res += (bit_p.query(i)-bit_p.query(pre)) * pow_p[zero[i]] % mod
            res %= mod

            res -= k * pow_p[zero[i]-zero[pre]] % mod
            res %= mod

            bit_cnt.update(pre,-1)
            bit_p.update(pre,-pow_ip[zero[pre]])

            lastappear[A[i]] = i
            bit_cnt.update(i,1)
            bit_p.update(i,pow_ip[zero[i]])

    res *= i2
    res %= mod
    res +=  M
    res %= mod

    return res

N,M = mi()
A = li()

print(solve(N,M,A))
0