結果

問題 No.1621 Sequence Inversions
ユーザー kyaneko999kyaneko999
提出日時 2021-07-23 07:51:45
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,807 ms / 3,000 ms
コード長 3,629 bytes
コンパイル時間 181 ms
コンパイル使用メモリ 82,372 KB
実行使用メモリ 279,380 KB
最終ジャッジ日時 2024-07-18 04:39:00
合計ジャッジ時間 26,088 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 59 ms
65,972 KB
testcase_01 AC 55 ms
65,256 KB
testcase_02 AC 69 ms
73,720 KB
testcase_03 AC 76 ms
74,952 KB
testcase_04 AC 999 ms
224,828 KB
testcase_05 AC 52 ms
65,512 KB
testcase_06 AC 72 ms
73,924 KB
testcase_07 AC 100 ms
81,852 KB
testcase_08 AC 1,599 ms
278,968 KB
testcase_09 AC 1,548 ms
279,320 KB
testcase_10 AC 1,454 ms
279,152 KB
testcase_11 AC 1,609 ms
279,060 KB
testcase_12 AC 1,807 ms
278,964 KB
testcase_13 AC 1,735 ms
279,028 KB
testcase_14 AC 1,770 ms
278,824 KB
testcase_15 AC 1,771 ms
278,940 KB
testcase_16 AC 1,269 ms
221,932 KB
testcase_17 AC 1,725 ms
271,132 KB
testcase_18 AC 1,238 ms
215,948 KB
testcase_19 AC 77 ms
74,804 KB
testcase_20 AC 417 ms
110,632 KB
testcase_21 AC 1,745 ms
279,012 KB
testcase_22 AC 1,765 ms
279,380 KB
testcase_23 AC 1,766 ms
278,864 KB
testcase_24 AC 53 ms
66,172 KB
testcase_25 AC 59 ms
69,396 KB
testcase_26 AC 52 ms
66,400 KB
testcase_27 AC 51 ms
65,364 KB
testcase_28 AC 51 ms
65,116 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