結果

問題 No.1025 Modular Equation
コンテスト
ユーザー convexineq
提出日時 2020-04-10 22:46:17
言語 Python3
(3.14.3 + numpy 2.4.4 + scipy 1.17.1)
コンパイル:
python3 -mpy_compile _filename_
実行:
python3 _filename_
結果
TLE  
実行時間 -
コード長 912 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 446 ms
コンパイル使用メモリ 20,576 KB
実行使用メモリ 63,080 KB
最終ジャッジ日時 2026-04-05 02:42:00
合計ジャッジ時間 22,435 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge3_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 6 TLE * 1 -- * 25
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

# 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