結果

問題 No.1621 Sequence Inversions
ユーザー kyaneko999kyaneko999
提出日時 2021-07-22 22:36:39
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 4,053 bytes
コンパイル時間 171 ms
コンパイル使用メモリ 82,568 KB
実行使用メモリ 115,960 KB
最終ジャッジ日時 2024-07-17 18:52:31
合計ジャッジ時間 8,400 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 57 ms
72,476 KB
testcase_01 AC 55 ms
65,120 KB
testcase_02 AC 76 ms
74,656 KB
testcase_03 AC 57 ms
64,948 KB
testcase_04 AC 60 ms
67,752 KB
testcase_05 AC 57 ms
66,532 KB
testcase_06 AC 72 ms
72,932 KB
testcase_07 AC 113 ms
78,156 KB
testcase_08 AC 2,994 ms
92,552 KB
testcase_09 TLE -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
権限があれば一括ダウンロードができます

ソースコード

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=0):
    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
for key in sorted(C.keys()):
    v=C[key]
    Mp=cnt*v
    dpp=[[0]*(Mp+1) for k in range(v+1)]
    for k in range(v+1):
        dpp[k][0]=1
    for c in range(1,cnt+1):
        nxt=[[0]*(Mp+1) for k in range(v+1)]
        debug(key,c,dpp)
        for k in range(Mp//c+1):
            for cc in range(Mp-c*k+1):
                for kk in range(v-k+1):
                    debug(kk,k,cc,c)
                    nxt[kk+k][cc+c*k]+=dpp[kk][cc]
                    nxt[kk+k][cc+c*k]%=MOD
        dpp=nxt
#     for _ in range(v):
#         nxt=[0]*(Mp+1)
#         su=0
#         for x in range(Mp+1):
#             if x>cnt:
#                 su-=dpp[x-cnt]
#             su+=dpp[x]
#             su%=MOD
#             nxt[x]=su
#         dpp=nxt
    nxt=[0]*(M+1)
    for c in range(Mp+1):
        for k in range(M-c+1):
            nxt[k+c]+=dp[k]*dpp[v][c]
            nxt[k+c]%=MOD
    dp=nxt
    cnt+=v
    #debug(dp)
# C=Counter(A)
# for v in C.values():
#     K+=v*(v-1)//2
print(dp[K])
0