結果
| 問題 |
No.1068 #いろいろな色 / Red and Blue and more various colors (Hard)
|
| コンテスト | |
| ユーザー |
とりゐ
|
| 提出日時 | 2021-10-18 02:23:35 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 2,695 ms / 3,500 ms |
| コード長 | 1,374 bytes |
| コンパイル時間 | 248 ms |
| コンパイル使用メモリ | 82,052 KB |
| 実行使用メモリ | 200,380 KB |
| 最終ジャッジ日時 | 2024-09-17 21:40:43 |
| 合計ジャッジ時間 | 51,834 ms |
|
ジャッジサーバーID (参考情報) |
judge6 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 29 |
ソースコード
p, g = 998244353, 3
invg = pow(g, p-2, p)
W = [pow(g, (p - 1) >> i, p) for i in range(24)]
iW = [pow(invg, (p - 1) >> i, p) for i in range(24)]
def fft(k, f):
for l in range(k)[::-1]:
d = 1 << l
u = 1
for i in range(d):
for j in range(i, 1 << k, 2*d):
f[j], f[j+d] = (f[j] + f[j+d]) % p, u * (f[j] - f[j+d]) % p
u = u * W[l+1] % p
def ifft(k, f):
for l in range(k):
d = 1 << l
u = 1
for i in range(d):
for j in range(i, 1 << k, 2*d):
f[j+d] *= u
f[j], f[j+d] = (f[j] + f[j+d]) % p, (f[j] - f[j+d]) % p
u = u * iW[l+1] % p
def convolve(a, b):
n0, n1 = len(a), len(b)
k = (max(n0, n1) - 1).bit_length() + 1
n = 1 << k
a = a + [0] * (n-n0)
b = b + [0] * (n-n1)
fft(k, a), fft(k, b)
for i in range(n):
a[i] = a[i] * b[i] % p
ifft(k, a)
invn = pow(n, p - 2, p)
return [a[i] * invn % p for i in range(n0 + n1 - 1)]
from collections import deque
def conv_all(l):
# 分割統治
q = deque(l)
while len(q)>=2:
a = q.popleft()
b = q.popleft()
q.append(convolve(a,b))
return q[0]
n,q=map(int,input().split())
a=list(map(int,input().split()))
FFT=conv_all([[(i-1)%p,1] for i in a])
for i in list(map(int,input().split())):
print(FFT[i])
とりゐ