結果
問題 | No.980 Fibonacci Convolution Hard |
ユーザー | Navier_Boltzmann |
提出日時 | 2023-11-17 22:35:45 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,571 bytes |
コンパイル時間 | 216 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 431,528 KB |
最終ジャッジ日時 | 2024-09-26 05:33:29 |
合計ジャッジ時間 | 8,769 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | TLE | - |
testcase_01 | -- | - |
testcase_02 | -- | - |
testcase_03 | -- | - |
testcase_04 | -- | - |
testcase_05 | -- | - |
testcase_06 | -- | - |
testcase_07 | -- | - |
testcase_08 | -- | - |
testcase_09 | -- | - |
testcase_10 | -- | - |
testcase_11 | -- | - |
testcase_12 | -- | - |
testcase_13 | -- | - |
testcase_14 | -- | - |
testcase_15 | -- | - |
testcase_16 | -- | - |
ソースコード
from collections import * from itertools import * from functools import * from heapq import * import sys,math input = sys.stdin.readline mod = 10**9 + 7 mem = [-1]*(2*10**6 + 1) mem[0]=0 mem[1]=0 mem[2]=1 import numpy as np p = int(input()) Q = int(input()) q = [int(input()) for _ in range(Q)] def convolve(f, g): tf = np.array(f, np.int64) tg = np.array(g, np.int64) fft_len = 1 while 2 * fft_len < len(tf) + len(tg) - 1: fft_len *= 2 fft_len *= 2 # フーリエ変換 Ff = np.fft.rfft(tf, fft_len) Fg = np.fft.rfft(tg, fft_len) # 各点積 Fh = Ff * Fg # フーリエ逆変換 h = np.fft.irfft(Fh, fft_len) # 小数になっているので、整数にまるめる h = np.rint(h).astype(np.int64) return h[:len(f) + len(g) - 1] def convolve2(f,g,p): f1,f2 = np.divmod(f,1<<15) g1,g2 = np.divmod(g,1<<15) a = convolve(f1,g1)%p c = convolve(f2,g2)%p b = (convolve(f1+f2,g1+g2) - a - c)%p h = (a<<30) + (b<<15) + c return h%p def convolve_pow(f,n,p): nbit = list(str(bin(n))[2:]) nbit = [int(i) for i in nbit] N = len(f) C = [1] + [0]*(N-1) B = f for i in range(len(nbit)): if nbit[-1-i] == 1: C = convolve2(C,B,p) B = convolve2(B,B,p) return C def f(x): if mem[x]!=-1: return mem[x] mem[x] = (p*f(x-1) + f(x-2))%mod return mem[x] X = [f(i) for i in range(2*10**6 + 1)] Z = convolve2(X,X,mod) # print(X) for qq in q: print(Z[qq])