結果

問題 No.2770 Coupon Optimization
コンテスト
ユーザー srjywrdnprkt
提出日時 2023-11-30 11:01:33
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 2,636 ms / 3,000 ms
コード長 1,878 bytes
コンパイル時間 478 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 329,224 KB
最終ジャッジ日時 2024-12-14 20:23:40
合計ジャッジ時間 33,271 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

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

N, M = map(int,input().split())
A = list(map(int,input().split()))
B = list(map(int,input().split()))
for i in range(N):
  A[i] //= 100
for i in range(M):
  B[i] = 100-B[i]
  
for i in range(max(0, N-M)):
  B.append(100)
A.sort()
B.sort()

C = convolution(A,B)
for i in range(N):
  print(C[i])
0