結果

問題 No.206 数の積集合を求めるクエリ
ユーザー Navier_BoltzmannNavier_Boltzmann
提出日時 2023-11-24 16:08:12
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 6,257 ms / 7,000 ms
コード長 1,834 bytes
コンパイル時間 277 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 71,088 KB
最終ジャッジ日時 2024-09-26 08:39:27
合計ジャッジ時間 89,487 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 29 ms
11,136 KB
testcase_01 AC 29 ms
11,136 KB
testcase_02 AC 34 ms
11,264 KB
testcase_03 AC 31 ms
11,264 KB
testcase_04 AC 29 ms
11,264 KB
testcase_05 AC 30 ms
11,264 KB
testcase_06 AC 153 ms
12,928 KB
testcase_07 AC 152 ms
12,928 KB
testcase_08 AC 157 ms
12,800 KB
testcase_09 AC 154 ms
12,800 KB
testcase_10 AC 30 ms
11,264 KB
testcase_11 AC 35 ms
11,264 KB
testcase_12 AC 158 ms
13,056 KB
testcase_13 AC 154 ms
12,928 KB
testcase_14 AC 164 ms
12,800 KB
testcase_15 AC 158 ms
12,672 KB
testcase_16 AC 158 ms
12,800 KB
testcase_17 AC 6,110 ms
71,088 KB
testcase_18 AC 6,257 ms
67,408 KB
testcase_19 AC 6,199 ms
70,548 KB
testcase_20 AC 6,049 ms
61,484 KB
testcase_21 AC 6,085 ms
66,876 KB
testcase_22 AC 6,087 ms
65,812 KB
testcase_23 AC 6,046 ms
70,568 KB
testcase_24 AC 6,159 ms
70,972 KB
testcase_25 AC 6,236 ms
69,756 KB
testcase_26 AC 6,147 ms
67,280 KB
testcase_27 AC 6,025 ms
63,776 KB
testcase_28 AC 6,218 ms
66,460 KB
testcase_29 AC 6,164 ms
67,012 KB
testcase_30 AC 6,052 ms
63,448 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import cmath
def convolution(f,g):

    n = len(bin(len(f)+len(g)-1)) - 2
    fft_length = 1<<n

    f = f + [0]*(fft_length - len(f))
    g = g + [0]*(fft_length - len(g))

    def fft(a):

        for i in range(fft_length):
            j = 0
            for k in range(n):
                j |= ((i>>k)&1) << (n - 1 - k)

            if i<j:
                a[i],a[j] = a[j],a[i]
        
        for nn in range(n):
            b = 1<<nn
            for j in range(b):

                w = cmath.rect(1,(-cmath.pi)*j/b)

                for k in range(0,fft_length,2*b):

                    s = a[j+k]
                    t = a[j+k+b]*w
                    a[j+k] = s+t
                    a[j+k+b] = s-t
        
        return a
    
    def ifft(a):
        for i in range(fft_length):
            j = 0
            for k in range(n):
                j |= ((i>>k)&1) << (n - 1 - k)

            if i<j:
                a[i],a[j] = a[j],a[i]
        
        for nn in range(n):
            b = 1<<nn
            for j in range(b):

                w = cmath.rect(1,(cmath.pi)*j/b)

                for k in range(0,fft_length,2*b):

                    s = a[j+k]
                    t = a[j+k+b]*w
                    a[j+k] = s+t
                    a[j+k+b] = s-t
        
        return [i/fft_length for i in a]
            
    F = fft(f)
    G = fft(g)
    H = [i*j for i,j in zip(F,G)]
    return [round(i.real) for i in ifft(H)]


from collections import *
from itertools import *
from functools import *
from heapq import *
import sys,math
input = sys.stdin.readline

L,M,N = map(int,input().split())
A = list(map(int,input().split()))
B = list(map(int,input().split()))
X = [0]*(N+1)
Y = [0]*(N+1)
for a in A:
    X[a] += 1
for b in B:
    Y[N-b] += 1

Q = int(input())
Z = convolution(X,Y)
print(*Z[N:N+Q],sep='\n')
0