結果

問題 No.854 公平なりんご分配
ユーザー kwm_t
提出日時 2023-05-01 15:37:06
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 762 ms / 3,153 ms
コード長 2,052 bytes
コンパイル時間 1,959 ms
コンパイル使用メモリ 200,108 KB
最終ジャッジ日時 2025-02-12 16:21:52
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 92
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
//#include <atcoder/all>
using namespace std;
//using namespace atcoder;
//using mint = modint1000000007;
//const int mod = 1000000007;
//using mint = modint998244353;
//const int mod = 998244353;
//const int INF = 1e9;
//const long long LINF = 1e18;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i,l,r)for(int i=(l);i<(r);++i)
#define rrep(i, n) for (int i = (n-1); i >= 0; --i)
#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)
#define all(x) (x).begin(),(x).end()
#define allR(x) (x).rbegin(),(x).rend()
#define endl "\n"
#define P pair<int,int>
template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }
vector<long long> sieve(long long  n = 1000000) {
	vector<bool>p(n + 1, true);
	vector<long long>ret;
	p[0] = false;
	p[1] = false;
	for (int i = 2; i <= n; ++i) {
		if (!p[i]) continue;
		ret.push_back(i);
		for (int j = 2 * i; j <= n; j += i) p[j] = false;
	}
	return ret;
}
int x[100005][304];
int sx[100005][304];
int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int n; cin >> n;
	auto plist = sieve(2000);
	rep(i, n) {
		int a; cin >> a;
		if (0 == a) {
			x[i][303] = 1;
			continue;
		}
		rep(j, plist.size()) {
		int cnt = 0;
			while (0 == a % plist[j]) {
				cnt++;
				a /= plist[j];
			}
			x[i][j] = cnt;
		}
	}
	rep(i, n) rep(j, 304)sx[i + 1][j] = sx[i][j] + x[i][j];
	int q; cin >> q;
	while (q--) {
		int p, l, r; cin >> p >> l >> r;
		l--, r--;
		if (sx[r + 1][303] - sx[l][303]) {
			cout << "Yes" << endl;
			continue;
		}
		vector<int>cnt(304);
		rep(i, plist.size()) {
			while (0 == p % plist[i]) {
				cnt[i]++;
				p /= plist[i];
			}
		}
		if (1 != p) {
			cout << "NO" << endl;
			continue;
		}
		bool chk = true;
		rep(i, 303) {
			int x = sx[r + 1][i] - sx[l][i];
			if (x < cnt[i])chk = false;
		}
		if(chk)cout << "Yes" << endl;
		else cout << "NO" << endl;
	}
	return 0;
}
0