結果

問題 No.1025 Modular Equation
ユーザー convexineqconvexineq
提出日時 2020-04-10 22:46:17
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 912 bytes
コンパイル時間 193 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 95,340 KB
最終ジャッジ日時 2024-09-15 22:55:57
合計ジャッジ時間 21,604 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge6
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,040 ms
58,960 KB
testcase_01 AC 1,039 ms
59,212 KB
testcase_02 AC 1,137 ms
59,800 KB
testcase_03 AC 1,144 ms
59,084 KB
testcase_04 AC 1,140 ms
58,576 KB
testcase_05 AC 1,106 ms
58,960 KB
testcase_06 AC 1,109 ms
58,952 KB
testcase_07 AC 1,127 ms
58,704 KB
testcase_08 AC 1,130 ms
59,088 KB
testcase_09 TLE -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

# coding: utf-8
# Your code here!
import sys
sys.setrecursionlimit(10**6)
readline = sys.stdin.readline
read = sys.stdin.read

import numpy
from scipy.fftpack import fft, ifft
def FastConvolution(A,B): # A,B: 2^c 成分を想定
    return ifft(fft(A)*fft(B))

p,n,k,b,*a = [int(i) for i in read().split()]

M = 1
while M <= p: M <<= 1
M <<= 1
V = (-M)%p
q = (M+V)//p

lst = [0]*p
for i in range(p): lst[pow(i,k,p)] += 1

x = [1]+[0]*(M-1)

MOD = 10**9+7
#print(lst,"lst")
#print(p,q)
for ai in a:
    y = [0]*M
    for i,c in enumerate(lst):
        y[i*ai%p] += c

    #print(x,y)    
    x = FastConvolution(x,y)
    x = numpy.mod(numpy.round(numpy.real(x)).astype(int),MOD)
    z = numpy.reshape(numpy.append(x,[0]*V),[q,p])
    x = numpy.append(numpy.mod(numpy.sum(z, axis=0),MOD),[0]*(M-p))
    #x = numpy.sum(numpy.reshape(numpy.append(x,[0]*V),[p,q]))
    #print(x)
    
    #print(x)


print(x[b])
0