結果
問題 |
No.854 公平なりんご分配
|
ユーザー |
![]() |
提出日時 | 2020-12-21 23:24:06 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,336 bytes |
コンパイル時間 | 256 ms |
コンパイル使用メモリ | 82,512 KB |
実行使用メモリ | 114,148 KB |
最終ジャッジ日時 | 2024-09-21 13:18:03 |
合計ジャッジ時間 | 13,332 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 20 WA * 30 TLE * 1 -- * 41 |
ソースコード
from itertools import chain, accumulate import sys input = sys.stdin.buffer.readline sys.setrecursionlimit(10 ** 7) U = 2000 U = 100 def prime_set(N): """ Nまでの素数のsetを返す """ if N < 4: return ({}, {}, {2}, {2, 3})[N] Nsq = int(N ** 0.5 + 0.5) + 1 primes = {2, 3} | set(chain(range(5, N + 1, 6), range(7, N + 1, 6))) for i in range(5, Nsq, 2): if i in primes: primes -= set(range(i * i, N + 1, i * 2)) return primes prime = sorted(prime_set(U)) N = int(input()) A = list(map(int, input().split())) memo = [[0] * (N + 1) for _ in range(len(prime) + 1)] for i, a in enumerate(A, 1): if a == 0: memo[0][i] = 1 memo[0] = tuple(accumulate(memo[0])) for i, p in enumerate(prime, 1): for j, a in enumerate(A, 1): while a % p == 0: memo[i][j] += 1 a //= p memo[i] = tuple(accumulate(memo[i])) Q = int(input()) ans = [] for _ in range(Q): X, S, T = map(int, input().split()) S -= 1 ok = (memo[0][T] - memo[0][S] == 0) for i, p in enumerate(prime, 1): cnt = 0 while X % p == 0: cnt += 1 X //= p ok &= memo[i][T] - memo[i][S] >= cnt if ok and X == 1: ans.append("Yes") else: ans.append("NO") if ans: print("\n".join(ans))