結果
| 問題 |
No.854 公平なりんご分配
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-03-20 19:10:16 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 2,567 bytes |
| コンパイル時間 | 3,108 ms |
| コンパイル使用メモリ | 205,472 KB |
| 最終ジャッジ日時 | 2025-01-19 20:21:04 |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 WA * 1 |
| other | AC * 1 WA * 49 TLE * 30 MLE * 12 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
class factor {
vector<int> p;
public:
factor() : p(vector<int>(303)) {}
factor(vector<int> p): p(p) {}
int get(int i) const { return p[i]; }
};
const factor ZERO = factor();
factor operator+(const factor &a, const factor &b) {
vector<int> c(303);
for (int i = 0; i < 303; i++)
c[i] = a.get(i) + b.get(i);
return factor(c);
}
template<typename T>
class SegTree {
int n;
vector<T> seg;
T nil;
function<T(T,T)> binop;
void update(int k, int l, int r, int p, T x) {
if (p < l || r <= p) return;
if (r - l == 1) seg[k] = x;
else {
update(2*k+1, l, (l+r)/2, p, x);
update(2*k+2, (l+r)/2, r, p, x);
seg[k] = binop(seg[2*k+1], seg[2*k+2]);
}
}
T query(int k, int l, int r, int a, int b) {
if (b <= l || r <= a) return nil;
if (a <= l && r <= b) return seg[k];
return binop(query(2*k+1, l, (l+r)/2, a, b), query(2*k+2, (l+r)/2, r, a, b));
}
public:
SegTree(int n, T nil, function<T(T, T)> binop): n(n), seg(4*n), nil(nil), binop(binop) {}
void init() {
for (int i = 0; i < 4*n; i++)
seg[i] = nil;
}
// update p-th value to x
void update(int p, T x) {
update(0, 0, n, p, x);
}
// query for range [l, r)
T query(int l, int r) {
return query(0, 0, n, l, r);
}
};
bool isprime(long long n) {
if (n < 2) return false;
if (n == 2) return true;
if (n % 2 == 0) return false;
for (long long i = 3; i * i <= n; i += 2)
if (n % i == 0)
return false;
return true;
}
int main() {
ios_base::sync_with_stdio(0);
cin.tie(0);
vector<int> primes;
for (int i = 1; i <= 2000; i++)
if (isprime(i))
primes.push_back(i);
int n;
cin >> n;
function<factor(factor,factor)> plus = [] (factor a, factor b) { return a + b; };
SegTree tree(n, ZERO, plus);
for (int i = 0; i < n; i++) {
int a;
cin >> a;
vector<int> p(303);
for (int j = 0; j < 303; j++) {
while (a % primes[j] == 0) {
a /= primes[j];
++p[j];
}
}
tree.update(i, factor(p));
}
int q;
cin >> q;
for (int i = 0; i < q; i++) {
long long p;
int l, r;
cin >> p >> l >> r;
factor x = tree.query(l-1, r);
bool ck = true;
for (int j = 0; j < 303; j++) {
int c = 0;
while (p % primes[j] == 0) {
++c;
p /= primes[j];
}
if (c > x.get(j)) {
ck = false;
break;
}
}
if (ck)
cout << "Yes" << endl;
else
cout << "No" << endl;
}
return 0;
}