結果

問題 No.1621 Sequence Inversions
ユーザー kyaneko999kyaneko999
提出日時 2021-07-23 07:51:45
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 2,023 ms / 3,000 ms
コード長 3,629 bytes
コンパイル時間 357 ms
コンパイル使用メモリ 87,168 KB
実行使用メモリ 281,152 KB
最終ジャッジ日時 2023-09-25 05:27:57
合計ジャッジ時間 31,158 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 134 ms
77,852 KB
testcase_01 AC 139 ms
77,920 KB
testcase_02 AC 157 ms
79,408 KB
testcase_03 AC 164 ms
79,136 KB
testcase_04 AC 1,187 ms
234,324 KB
testcase_05 AC 138 ms
77,728 KB
testcase_06 AC 155 ms
79,540 KB
testcase_07 AC 175 ms
79,272 KB
testcase_08 AC 1,827 ms
280,552 KB
testcase_09 AC 1,670 ms
280,784 KB
testcase_10 AC 1,664 ms
280,708 KB
testcase_11 AC 1,814 ms
280,612 KB
testcase_12 AC 1,910 ms
280,844 KB
testcase_13 AC 1,958 ms
280,996 KB
testcase_14 AC 1,987 ms
280,792 KB
testcase_15 AC 2,005 ms
280,932 KB
testcase_16 AC 1,479 ms
223,784 KB
testcase_17 AC 1,932 ms
272,680 KB
testcase_18 AC 1,423 ms
212,072 KB
testcase_19 AC 155 ms
79,340 KB
testcase_20 AC 530 ms
119,988 KB
testcase_21 AC 2,004 ms
280,748 KB
testcase_22 AC 2,023 ms
280,856 KB
testcase_23 AC 2,012 ms
281,152 KB
testcase_24 AC 136 ms
78,048 KB
testcase_25 AC 141 ms
78,836 KB
testcase_26 AC 136 ms
77,808 KB
testcase_27 AC 138 ms
77,720 KB
testcase_28 AC 138 ms
78,040 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from sys import exit, stdin, setrecursionlimit
from collections import deque, defaultdict, Counter
from copy import deepcopy
from bisect import bisect_left, bisect_right, insort_left, insort_right
from heapq import heapify, heappop, heappush
from itertools import product, permutations, combinations, combinations_with_replacement
from functools import reduce
from math import gcd, sin, cos, tan, asin, acos, atan, atan2, degrees, radians, ceil, floor, sqrt, factorial
from math import pi as PI
from random import randint
# from decimal import Decimal

setrecursionlimit(500000)
INF = (1<<61)-1
EPS = 1e-10
MOD = 10**9+7
MOD = 998244353

def input():
    return stdin.readline().strip('\n')
def intput():
    return int(input())
def minput():
    return input().split()
def linput():
    return input().split()
def mint():
    return map(int,input().split())
def lint():
    return list(map(int,input().split()))
def ilint():
    return intput(),lint()
def lcm(x,y):
    return x*y//gcd(x,y)
def lgcd(l):
    return reduce(gcd,l)
def llcm(l):
    return reduce(lcm,l)
def powmod(n,i,mod=MOD):
    return pow(n,mod-1+i,mod) if i<0 else pow(n,i,mod)
def div2(x):
    return x.bit_length()
def div10(x):
    return len(str(x))-(x==0)
def popcount(x):
    return bin(x).count('1')
def digit(x,i,max_len=None):
    s = str(x)
    if max_len:
        i -= max_len-len(s)
    return int(s[i-1]) if i>0 else 0
def digitsum(x):
    ans = 0
    for i in range(div10(x)):
        ans += digit(x,i+1)
    return ans
def pf(x,mode='counter'):
    C = Counter()
    p = 2
    while x>1:
        k = 0
        while x%p==0:
            x //= p
            k += 1
        if k>0:
            C[p] += k
        p = p+2-(p==2) if p*p<x else x
    if mode=='counter':
        return C
    S = set([1])
    for k in C:
        T = set()
        for x in S:
            for i in range(C[k]+1):
                T.add(x*(k**i))
        S = T
    if mode=='set':
        return S
    if mode=='list':
        return sorted(S)
def isprime(x):
    if x<2:
        return False
    return len(pf(x,'set'))==2
def matmul(A, B):
    # import numpy
    A1, A2 = A >> 15, A & (1 << 15) - 1
    B1, B2 = B >> 15, B & (1 << 15) - 1
    X = np.dot(A1, B1) % MOD
    Y = np.dot(A2, B2)
    Z = np.dot(A1 + A2, B1 + B2) - X - Y
    return ((X << 30) + (Z << 15) + Y) % MOD
def matpow(A, N):
    P = np.eye(A.shape[0], dtype=np.int64)
    while N:
        if N & 1:
            P = matmul(P, A)
        A = matmul(A, A)
        N >>= 1
    return P
def zash(S):
    lis = sorted(S)
    dic = {}
    for i,x in enumerate(lis):
        dic[x] = i
    return lis, dic
def pr(*x):
    print(*x, sep='', end='') if len(x) else print()
def lprint(l):
    for x in l: print(x)
def ston(c, c0='a'):
    return ord(c)-ord(c0)
def ntos(x, c0='a'):
    return chr(x+ord(c0))
def judge(x, l=['Yes', 'No']):
    print(l[0] if x else l[1])
def debug(*x, flag=1):
    if flag: print(*x)

######################################################

N,K=mint()
M=N*(N-1)//2
A=lint()
C=Counter(A)
dp=[0]*(M+1)
dp[0]=1
cnt=0
M0=N*N//4
dpp=[[[0]*(M0+1) for _ in range(N+1)] for _ in range(N+1)]
for c in range(N+1):
    dpp[0][c][0]=1
for j in range(1,N+1):
    for c in range(N+1):
        for x in range(M0+1):
            dpp[j][c][x]=dpp[j-1][c][x]
            if x-j>=0:
                dpp[j][c][x]+=dpp[j][c-1][x-j]
                dpp[j][c][x]%=MOD
for key in sorted(C.keys()):
    v=C[key]
    Mp=cnt*v
    nxt=[0]*(M+1)
    for x in range(Mp+1):
        for k in range(M-x+1):
            nxt[k+x]+=dp[k]*dpp[v][cnt][x]
            nxt[k+x]%=MOD
    dp=nxt
    cnt+=v
print(dp[K])
0