結果

問題 No.2164 Equal Balls
ユーザー tassei903tassei903
提出日時 2022-12-15 10:21:46
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 2,480 ms / 5,000 ms
コード長 4,959 bytes
コンパイル時間 151 ms
コンパイル使用メモリ 82,592 KB
実行使用メモリ 139,948 KB
最終ジャッジ日時 2024-11-13 23:50:32
合計ジャッジ時間 42,906 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 66 ms
79,088 KB
testcase_01 AC 73 ms
77,632 KB
testcase_02 AC 84 ms
85,664 KB
testcase_03 AC 66 ms
78,976 KB
testcase_04 AC 68 ms
78,572 KB
testcase_05 AC 68 ms
79,072 KB
testcase_06 AC 75 ms
77,880 KB
testcase_07 AC 69 ms
77,584 KB
testcase_08 AC 201 ms
94,812 KB
testcase_09 AC 211 ms
96,100 KB
testcase_10 AC 142 ms
93,172 KB
testcase_11 AC 316 ms
102,248 KB
testcase_12 AC 188 ms
95,200 KB
testcase_13 AC 158 ms
93,828 KB
testcase_14 AC 183 ms
94,668 KB
testcase_15 AC 309 ms
101,828 KB
testcase_16 AC 195 ms
94,932 KB
testcase_17 AC 124 ms
92,836 KB
testcase_18 AC 228 ms
96,680 KB
testcase_19 AC 352 ms
102,412 KB
testcase_20 AC 204 ms
95,756 KB
testcase_21 AC 112 ms
92,500 KB
testcase_22 AC 101 ms
92,668 KB
testcase_23 AC 170 ms
118,820 KB
testcase_24 AC 158 ms
106,624 KB
testcase_25 AC 142 ms
118,196 KB
testcase_26 AC 142 ms
93,936 KB
testcase_27 AC 166 ms
111,316 KB
testcase_28 AC 170 ms
111,172 KB
testcase_29 AC 129 ms
100,308 KB
testcase_30 AC 162 ms
111,376 KB
testcase_31 AC 175 ms
119,268 KB
testcase_32 AC 172 ms
114,380 KB
testcase_33 AC 135 ms
111,644 KB
testcase_34 AC 165 ms
120,828 KB
testcase_35 AC 96 ms
94,032 KB
testcase_36 AC 217 ms
125,076 KB
testcase_37 AC 184 ms
127,452 KB
testcase_38 AC 2,411 ms
139,184 KB
testcase_39 AC 2,412 ms
139,328 KB
testcase_40 AC 2,413 ms
139,252 KB
testcase_41 AC 2,428 ms
139,040 KB
testcase_42 AC 2,380 ms
139,180 KB
testcase_43 AC 2,399 ms
139,060 KB
testcase_44 AC 2,407 ms
139,308 KB
testcase_45 AC 2,415 ms
139,232 KB
testcase_46 AC 2,414 ms
139,172 KB
testcase_47 AC 2,404 ms
139,532 KB
testcase_48 AC 2,480 ms
139,948 KB
testcase_49 AC 1,693 ms
123,812 KB
testcase_50 AC 1,688 ms
123,808 KB
testcase_51 AC 1,640 ms
124,400 KB
testcase_52 AC 1,641 ms
124,468 KB
testcase_53 AC 1,627 ms
124,080 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = lambda :sys.stdin.readline()[:-1]
ni = lambda :int(input())
na = lambda :list(map(int,input().split()))
yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")
no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")
#######################################################################
def Primitive_Root(p):
    """Z/pZ上の原始根を見つける
 
    p:素数
    """
    if p==2:
        return 1
    if p==998244353:
        return 3
    if p==10**9+7:
        return 5
 
    fac=[]
    q=2
    v=p-1
 
    while v>=q*q:
        e=0
        while v%q==0:
            e+=1
            v//=q
 
        if e>0:
            fac.append(q)
        q+=1
 
    if v>1:
        fac.append(v)
 
    g=2
    while g<p:
        if pow(g,p-1,p)!=1:
            return None
 
        flag=True
        for q in fac:
            if pow(g,(p-1)//q,p)==1:
                flag=False
                break
 
        if flag:
            return g
 
        g+=1
 
#参考元 https://atcoder.jp/contests/practice2/submissions/16789717
def NTT(A):
    """AをMod を法とする数論変換を施す
 
    ※Modはグローバル変数から指定
    """
    primitive=Primitive_Root(Mod)
 
    N=len(A)
    H=(N-1).bit_length()
 
    if Mod==998_244_353:
        m=998_244_352
        u=119
        e=23
        S=[1,998244352,911660635,372528824,929031873,
           452798380,922799308,781712469,476477967,166035806,
           258648936,584193783,63912897,350007156,666702199,
           968855178,629671588,24514907,996173970,363395222,
           565042129,733596141,267099868,15311432]
    else:
        m=Mod-1
        e=((m&-m)-1).bit_length()
        u=m>>e
        S=[pow(primitive,(Mod-1)>>i,Mod) for i in range(e+1)]
 
    for l in range(H, 0, -1):
        d = 1 << l - 1
        U = [1]*(d+1)
        u = 1
        for i in range(d):
            u=u*S[l]%Mod
            U[i+1]=u
 
        for i in range(1 <<H - l):
            s=2*i*d
            for j in range(d):
                A[s],A[s+d]=(A[s]+A[s+d])%Mod, U[j]*(A[s]-A[s+d])%Mod
                s+=1
 
#参考元 https://atcoder.jp/contests/practice2/submissions/16789717
def Inverse_NTT(A):
    """AをMod を法とする逆数論変換を施す
 
    ※Modはグローバル変数から指定
    """
    primitive=Primitive_Root(Mod)
 
    N=len(A)
    H=(N-1).bit_length()
 
    if Mod==998244353:
        m=998_244_352
        e=23
        u=119
        S=[1,998244352,86583718,509520358,337190230,
           87557064,609441965,135236158,304459705,685443576,
           381598368,335559352,129292727,358024708,814576206,
           708402881,283043518,3707709,121392023,704923114,950391366,
           428961804,382752275,469870224]
    else:
        m=Mod-1
        e=(m&-m).bit_length()-1
        u=m>>e
 
        inv_primitive=pow(primitive,Mod-2,Mod)
        S=[pow(inv_primitive,(Mod-1)>>i,Mod) for i in range(e+1)]
 
    for l in range(1, H + 1):
        d = 1 << l - 1
        for i in range(1 << H - l):
            u = 1
            for j in range(i * 2 * d, (i * 2 + 1) * d):
                A[j+d] *= u
                A[j], A[j+d] = (A[j] + A[j+d]) % Mod, (A[j] - A[j+d]) % Mod
                u = u * S[l] % Mod
 
    N_inv=pow(N,Mod-2,Mod)
    for i in range(N):
        A[i]=A[i]*N_inv%Mod
 
#参考元 https://atcoder.jp/contests/practice2/submissions/16789717
def Convolution_Mod(A,B):
    """A,BをMod を法とする畳み込みを求める.
 
    ※Modはグローバル変数から指定
    """
    L=len(A)+len(B)-1
    H=L.bit_length()
    N=1<<H
 
    A=A+[0]*(N-len(A))
    B=B+[0]*(N-len(B))
 
    NTT(A)
    NTT(B)
 
    for i in range(N):
        A[i]=A[i]*B[i]%Mod
 
    Inverse_NTT(A)
 
    del A[L:]
    return A
#=================================================
Mod=998244353

mod = 998244353
nn = 10**6
fact = [1] * nn
for i in range(nn - 1):
    fact[i + 1] = fact[i] * (i + 1) % mod
invfact = [1] * nn
invfact[nn - 1] = pow(fact[nn - 1], mod - 2, mod)
for i in range(nn - 1)[::-1]:
    invfact[i] = invfact[i + 1] * (i + 1) % mod
 
def binom(x, y):
    if x < 0 or y < 0 or x - y < 0:
        return 0
    return fact[x] * invfact[y] % mod * invfact[x - y] % mod

n, m = na()

a = na()
b = na()
from collections import deque
dq = deque()


LA = [0] * m
LB = [0] * m
for i in range(m):
    A = []
    B = []
    la = float("inf")
    lb = float("inf")
    for j in range(i, n, m):
        la = min(la, a[j])
        lb = min(lb, b[j])
    LB[i] = lb
    LA[i] = la
S = sum(LB)
for i in range(m):
    la,lb = LA[i],LB[i]
    op = [1] * (la+lb+1)
    for j in range(i, n, m):
        for k in range(-lb, la+1):
            op[k+lb] *= binom(a[j]+b[j],k+b[j])
            op[k+lb] %= mod
    dq.append(op)

while len(dq)>1:
    l = dq.popleft()
    r = dq.popleft()
    z = Convolution_Mod(l, r)[:S+1]
    dq.append(z)
z = dq.popleft()
if len(z) <= S:
    print(0)
else:
    print(z[S])
0