結果

問題 No.1025 Modular Equation
ユーザー convexineq
提出日時 2020-04-10 23:18:03
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 694 bytes
コンパイル時間 257 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 22,308 KB
最終ジャッジ日時 2024-09-16 00:32:02
合計ジャッジ時間 8,040 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other WA * 6 TLE * 1 -- * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

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()]


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

MOD = 10**9+7


x = 1
M = 1<<96
NUM = 96*p
ALL = (1<<NUM)-1
m = int(("1" * 2 + "0" * 30) * p, 2)
pa = (1 << 30) - ((1 << 30) % MOD)

for ai in a:
    y = 0
    for i,c in enumerate(lst):
        if c!=0:
            y += c<<(i*ai%p*96)
    #print(bin(y)[2:])
    x = x*y
    x = (x>>(NUM)) + (x&ALL)
    x -= ((x & m) >> 96) * pa
    

ans = (x>>(b*96))&((1<<32)-1)
print(ans%MOD)
0