結果

問題 No.854 公平なりんご分配
ユーザー chocopuu
提出日時 2019-07-26 22:46:06
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 1,858 bytes
コンパイル時間 1,863 ms
コンパイル使用メモリ 183,872 KB
実行使用メモリ 245,088 KB
最終ジャッジ日時 2024-07-02 08:42:28
合計ジャッジ時間 23,765 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 24 WA * 17 RE * 51
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define FOR(i, s, n) for (int i = (s); i < (n); i++)
#define RFOR(i, s, n) for (int i = (n) - 1; i >= (s); i--)
#define REP(i, n) FOR(i, 0, n)
#define RREP(i, n) RFOR(i, 0, n)
#define ALL(a) a.begin(), a.end()
const long long MOD = 1e9+7, INF = 1e18;
template<class T>inline bool CHMAX(T&a,T b){if(a<b){a=b;return true;}return false;}
template<class T>inline bool CHMIN(T&a,T b){if(a>b){a=b;return true;}return false;}

vector<int> P;
void make_primes(int n)
{
	vector<bool> primes;
	n++;
	primes.resize(n, true);
	primes[0] = primes[1] = false;
	FOR(i, 2, sqrt(n)) if (primes[i]) for (int j = 0; i * (j + 2) < n; j++)
		primes[i * (j + 2)] = false;

	FOR(i, 2, n) if (primes[i]) P.push_back(i);
}

vector<pair<long long, long long> > prime_factorize(long long n) {
    //計算量は√N!!!
    vector<pair<long long, long long> > res;
    for (long long p = 2; p * p <= n; ++p) {
        if (n % p != 0) continue;
        int num = 0;
        while (n % p == 0) { ++num; n /= p; }
        res.push_back(make_pair(p, num));
    }
    if (n != 1) res.push_back(make_pair(n, 1));
    return res;
}

signed main(){
	int N;
	cin >> N;
	make_primes(2000);
	vector<int>inv(2000,0);
	REP(i,P.size())inv[P[i]] = i;
	vector<int> a(N); for(auto& e:a)cin >> e;
	vector<vector<int>>cnt(N+1,vector<int>(P.size(),0));
	REP(i,N){
		auto tmp = prime_factorize(a[i]);
		for(auto e:tmp){
			cnt[i+1][inv[e.first]] += e.second;
		}
	}
	REP(i,N)REP(j,P.size())cnt[i+1][j] += cnt[i][j];
	int Q;
	cin >> Q;
	vector<int>ans;
	REP(i,Q){
		int p,l,r;
		cin >> p >> l >> r;
		l--;r--;
		auto tmp = prime_factorize(p);
		int flg = 0;
		for(auto e:tmp)if(e.second>cnt[r+1][inv[e.first]]-cnt[l][inv[e.first]])flg = 1;
		ans.push_back(flg);
	}
	for(auto e:ans){
		if(e)cout<<"NO"<<endl;
		else cout<<"Yes"<<endl;
	}
}
0