結果

問題 No.206 数の積集合を求めるクエリ
ユーザー ayaoniayaoni
提出日時 2020-12-20 00:32:03
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 237 ms / 7,000 ms
コード長 2,072 bytes
コンパイル時間 678 ms
コンパイル使用メモリ 81,776 KB
実行使用メモリ 111,808 KB
最終ジャッジ日時 2023-10-21 09:53:37
合計ジャッジ時間 6,583 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
53,548 KB
testcase_01 AC 41 ms
53,548 KB
testcase_02 AC 54 ms
64,056 KB
testcase_03 AC 40 ms
53,548 KB
testcase_04 AC 41 ms
53,548 KB
testcase_05 AC 47 ms
61,408 KB
testcase_06 AC 64 ms
70,684 KB
testcase_07 AC 65 ms
70,700 KB
testcase_08 AC 65 ms
70,700 KB
testcase_09 AC 65 ms
70,700 KB
testcase_10 AC 40 ms
53,548 KB
testcase_11 AC 55 ms
64,056 KB
testcase_12 AC 79 ms
75,104 KB
testcase_13 AC 74 ms
74,540 KB
testcase_14 AC 74 ms
74,612 KB
testcase_15 AC 72 ms
74,068 KB
testcase_16 AC 73 ms
74,040 KB
testcase_17 AC 218 ms
111,532 KB
testcase_18 AC 199 ms
95,304 KB
testcase_19 AC 217 ms
108,080 KB
testcase_20 AC 197 ms
97,692 KB
testcase_21 AC 207 ms
99,780 KB
testcase_22 AC 208 ms
100,080 KB
testcase_23 AC 217 ms
109,804 KB
testcase_24 AC 237 ms
111,808 KB
testcase_25 AC 232 ms
108,124 KB
testcase_26 AC 217 ms
95,484 KB
testcase_27 AC 212 ms
92,324 KB
testcase_28 AC 225 ms
100,296 KB
testcase_29 AC 224 ms
100,052 KB
testcase_30 AC 214 ms
93,740 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

sys.setrecursionlimit(10**7)
def I(): return int(sys.stdin.readline().rstrip())
def MI(): return map(int,sys.stdin.readline().rstrip().split())
def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))
def LI2(): return list(map(int,sys.stdin.readline().rstrip()))
def S(): return sys.stdin.readline().rstrip()
def LS(): return list(sys.stdin.readline().rstrip().split())
def LS2(): return list(sys.stdin.readline().rstrip())


# 998244353 = 119*2**23+1

mod = 998244353
primitive_root = 3  # mod の原始根
roots = [pow(primitive_root,(mod-1) >> i,mod) for i in range(24)]
inv_roots = [pow(r,mod-2,mod) for r in roots]
# roots[i] = 1 の 2**i 乗根、inv_roots[i] = 1 の 2**i 乗根の逆元


# 順番は変わる

def ntt(A,n):
    for i in range(n):
        m = 1 << (n-i-1)
        for start in range(1 << i):
            w = 1
            start *= m*2
            for j in range(m):
                A[start+j],A[start+j+m] = (A[start+j]+A[start+j+m]) % mod,(A[start+j]-A[start+j+m])*w % mod
                w *= roots[n-i]
                w %= mod
    return A


def inv_ntt(A,n):
    for i in range(n):
        m = 1 << i
        for start in range(1 << (n-i-1)):
            w = 1
            start *= m*2
            for j in range(m):
                A[start+j],A[start+j+m] = (A[start+j]+A[start+j+m]*w) % mod,(A[start+j]-A[start+j+m]*w) % mod
                w *= inv_roots[i+1]
                w %= mod
    a = pow(2,n*(mod-2),mod)
    for i in range(1 << n):
        A[i] *= a
        A[i] %= mod
    return A


def convolution(A,B):
    a,b = len(A),len(B)
    deg = a+b-2
    n = deg.bit_length()
    N = 1 << n
    A += [0]*(N-a)  # A の次数を 2冪-1 にする
    B += [0]*(N-b)  # B の次数を 2冪-1 にする
    A = ntt(A,n)
    B = ntt(B,n)
    C = [(A[i]*B[i]) % mod for i in range(N)]
    C = inv_ntt(C,n)
    return C[:deg+1]


L,M,N = MI()
A = LI()
B = LI()

C = [0]*(N+1)
for a in A:
    C[a] = 1
D = [0]*(N+1)
for b in B:
    D[N-b] = 1

ANS = convolution(C,D)

Q = I()
for i in range(Q):
    print(ANS[N+i])
0