結果
| 問題 | No.2709 1975 Powers | 
| コンテスト | |
| ユーザー |  Cecil | 
| 提出日時 | 2024-03-31 16:23:00 | 
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) | 
| 結果 | 
                                TLE
                                 
                             | 
| 実行時間 | - | 
| コード長 | 1,091 bytes | 
| コンパイル時間 | 141 ms | 
| コンパイル使用メモリ | 12,800 KB | 
| 実行使用メモリ | 17,700 KB | 
| 最終ジャッジ日時 | 2024-09-30 21:04:55 | 
| 合計ジャッジ時間 | 3,656 ms | 
| ジャッジサーバーID (参考情報) | judge2 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| 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):
    for j in range(i, N):
        if i==j or A[i]>=A[j]:
            continue
        for k in range(j, N):
            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)
            
            
            
        