結果
| 問題 |
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 |
ソースコード
# 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)
convexineq