結果

問題 No.1025 Modular Equation
ユーザー convexineq
提出日時 2020-04-10 23:43:10
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 589 bytes
コンパイル時間 99 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 21,704 KB
最終ジャッジ日時 2024-09-16 03:11:14
合計ジャッジ時間 8,230 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
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


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 = mm = ((1<<32)-1)<<64
for _ in range(n):
    m <<= 96
    m ^= mm
pa = (1 << 64) - ((1 << 64) % MOD)

for ai in a:
    y = 0
    for i,c in enumerate(lst):
        y += c<<(i*ai%p*96)
    x = x*y
    x = (x>>NUM) + (x&ALL)
    x -= ((x & m) >> 64) * pa

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