結果
問題 | No.206 数の積集合を求めるクエリ |
ユーザー | ああいい |
提出日時 | 2022-02-12 19:40:39 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,203 ms / 7,000 ms |
コード長 | 1,752 bytes |
コンパイル時間 | 291 ms |
コンパイル使用メモリ | 87,248 KB |
実行使用メモリ | 284,616 KB |
最終ジャッジ日時 | 2023-09-11 07:47:51 |
合計ジャッジ時間 | 17,399 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 77 ms
71,752 KB |
testcase_01 | AC | 84 ms
75,952 KB |
testcase_02 | AC | 87 ms
77,036 KB |
testcase_03 | AC | 77 ms
71,944 KB |
testcase_04 | AC | 79 ms
75,652 KB |
testcase_05 | AC | 88 ms
76,868 KB |
testcase_06 | AC | 124 ms
78,480 KB |
testcase_07 | AC | 114 ms
79,488 KB |
testcase_08 | AC | 114 ms
79,336 KB |
testcase_09 | AC | 114 ms
79,324 KB |
testcase_10 | AC | 79 ms
75,984 KB |
testcase_11 | AC | 91 ms
77,100 KB |
testcase_12 | AC | 127 ms
79,296 KB |
testcase_13 | AC | 114 ms
78,876 KB |
testcase_14 | AC | 119 ms
78,764 KB |
testcase_15 | AC | 118 ms
78,488 KB |
testcase_16 | AC | 110 ms
78,768 KB |
testcase_17 | AC | 1,173 ms
284,552 KB |
testcase_18 | AC | 1,166 ms
280,632 KB |
testcase_19 | AC | 1,150 ms
282,368 KB |
testcase_20 | AC | 437 ms
134,748 KB |
testcase_21 | AC | 1,148 ms
280,792 KB |
testcase_22 | AC | 459 ms
135,836 KB |
testcase_23 | AC | 1,154 ms
282,340 KB |
testcase_24 | AC | 1,189 ms
284,616 KB |
testcase_25 | AC | 1,186 ms
282,564 KB |
testcase_26 | AC | 1,189 ms
280,772 KB |
testcase_27 | AC | 463 ms
138,480 KB |
testcase_28 | AC | 1,189 ms
281,560 KB |
testcase_29 | AC | 1,203 ms
281,108 KB |
testcase_30 | AC | 463 ms
139,968 KB |
ソースコード
L,M,N = map(int,input().split()) A = list(map(int,input().split())) B = list(map(int,input().split())) Q = int(input()) A.sort() B.sort() h = [A[0] - B[i] for i in reversed(range(M))] delta = [A[i] - A[0] for i in range(L)] import sys if h[0] < 0: sabun = -h[0] else: sabun = 0 n = h[-1] + sabun x = [0] * (n + 1) for i in h: x[sabun + i] += 1 m = delta[-1] y = [0] * (m + 1) for i in delta: y[i] += 1 import cmath #FFT 配列のサイズは2べきのみ #逆変換はinverse = True #配列Aは壊れる def fft(A,inverse = False): N = len(A) h = 0 while 1 << h < N: h += 1 for i in range(N): j = 0 for k in range(h): j |= (i >> k & 1) << (h-1-k) if i < j: A[i],A[j] = A[j],A[i] #↑バタフライ演算用配列整備 b = 1 if inverse:c = -1 else:c = 1 while b < N: for j in range(b): w = cmath.rect(1.0,-cmath.pi * j / b * c) for k in range(0,N,2*b): s = A[j+k] t = A[j+k+b] * w A[j+k] = s + t A[j+k+b] = s - t b *= 2 if inverse: for i in range(N): A[i] /= N return A def convolve(a,b): s = len(a) + len(b) - 1 t = 1 while t < s: t <<= 1 A = a + [0] * (t - len(a)) B = b + [0] * (t - len(b)) A = fft(A) B = fft(B) for i in range(t): A[i] *= B[i] A = fft(A,True) return A[:s] #↓整数版、配列は壊れない def Zconvolve(a,b): A = convolve(a,b) for i in range(len(A)): A[i] = round(A[i].real) return A C = Zconvolve(x,y) for i in range(Q): if sabun + i >= len(C): print(0) else: print(C[sabun+i])