結果
問題 |
No.990 N×Mマス計算(Kの倍数)
|
ユーザー |
![]() |
提出日時 | 2020-02-14 23:47:53 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 845 ms / 2,000 ms |
コード長 | 1,031 bytes |
コンパイル時間 | 205 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 22,212 KB |
最終ジャッジ日時 | 2024-11-16 01:50:06 |
合計ジャッジ時間 | 5,781 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 19 |
ソースコード
import bisect N,M,K = map(int, input().split()) B = [i for i in input().split()] op = B.pop(0) A = [int(input()) for _ in range(N)] for i in range(M): B[i] = int(B[i]) def gcd(a, b): if a < b: a, b = b, a if a%b == 0: return b return gcd(b, a%b) ans = 0 if op == "+": for i in range(N): A[i] %= K for j in range(M): B[j] %= K B.sort() for a in A: l = bisect.bisect_left(B, (K-a)%K) r = bisect.bisect_right(B, (K-a)%K) ans += r-l if K == 1: ans = N*M else: dic = {} i = 1 while i < K**0.5: if K%i == 0: dic[i] = 0 dic[K//i] = 0 i += 1 if K%int(K**0.5) == 0: dic[int(K**0.5)] = 0 for b in B: b = gcd(b, K) dic[b] += 1 Key = sorted(dic) for i in range(len(Key)): for j in range(i+1, len(Key)): if Key[j]%Key[i] == 0: dic[Key[i]] += dic[Key[j]] for a in A: ans += dic[K//gcd(a, K)] print(ans)