結果

問題 No.1307 Rotate and Accumulate
ユーザー WizistWizist
提出日時 2020-12-05 05:04:40
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,130 bytes
コンパイル時間 509 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 305,932 KB
最終ジャッジ日時 2024-09-15 10:25:35
合計ジャッジ時間 16,091 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

#!/usr/bin/env python3
#
# No.1307 Rotate and Accumulate
#
import sys, os, math

def fft(a, inverse=False):
	n = 1
	while n < len(a): n <<= 1
	if len(a) < n: a += [0] * (n - len(a))
	j, k = 0, 1
	for i in range(1, n):
		bit = n >> 1
		while j & bit: j ^= bit; bit >>= 1
		j ^= bit
		if i < j: a[i], a[j] = a[j], a[i]
	while k < n:
		wk = math.cos(math.pi / k) + math.sin(-math.pi / k if inverse else math.pi / k) * 1j
		for i in range(0, n, 2 * k):
			w = 1 + 0j
			for j in range(k):
				u, v = a[i + j], a[i + j + k] * w
				a[i + j] = u + v
				a[i + j + k] = u - v
				w *= wk
		k *= 2
	if inverse:
		for i in range(n): a[i] /= n

def multiply(a, b):
	n = 1
	while n < len(a) + len(b): n <<= 1
	f, g, r = [0] * n, [0] * n, [0] * n
	for i, v in enumerate(a): f[i] = v
	for i, v in enumerate(b): g[i] = v
	fft(f); fft(g)
	for i in range(n): f[i] *= g[i]
	fft(f, True)
	for i, v in enumerate(f): r[i] = round(v.real)
	return r

def read_ints(): return list(map(int, input().split()))

n, q = read_ints()
a, r = read_ints(), read_ints()

f = a + a
g = [0] * (n + 1)
for x in r: g[n - x] += 1

c = multiply(f, g)
print(c[n:n + n])
0