結果
| 問題 |
No.2709 1975 Powers
|
| コンテスト | |
| ユーザー |
Cecil
|
| 提出日時 | 2024-03-31 16:23:48 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,097 bytes |
| コンパイル時間 | 277 ms |
| コンパイル使用メモリ | 12,800 KB |
| 実行使用メモリ | 18,080 KB |
| 最終ジャッジ日時 | 2024-09-30 21:04:59 |
| 合計ジャッジ時間 | 3,688 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | TLE * 1 -- * 24 |
ソースコード
from collections import defaultdict as dd
import bisect
def pow(x, n, mod):
ans = 1
#n が 0 になるまで計算を続ける
while n:
if n % 2:
ans *= x
ans %= mod
x *= x
n >>= 1
return ans%mod
N,P,Q = map(int, input().split())
A = list(map(int, input().split()))
A = list(set(A))
A.sort()
N = len(A)
A_div = [ [0 for i in range(4)] for _ in range(N)]
num = [10,9,7,5]
for i in range(N):
for j in range(4):
A_div[i][j] = pow(num[j], A[i], P)
ans = 0
for i in range(N-3):
for j in range(i, N-2):
if i==j or A[i]>=A[j]:
continue
for k in range(j, N-1):
if j==k or A[j]>=A[k]:
continue
for l in range(k, N):
if l==k or A[k]>=A[l]:
continue
if A[i]<A[j]<A[k]<A[l]:
rem = A_div[i][0]+A_div[j][1]+A_div[k][2]+A_div[l][3]
rem %= P
if rem == Q:
ans += 1
else:
continue
print(ans)
Cecil