結果

問題 No.854 公平なりんご分配
コンテスト
ユーザー ahe100
提出日時 2019-07-26 23:21:14
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,864 bytes
コンパイル時間 1,668 ms
コンパイル使用メモリ 173,952 KB
実行使用メモリ 129,508 KB
最終ジャッジ日時 2024-07-02 09:47:03
合計ジャッジ時間 15,233 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 56 WA * 36
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define rep(i,n) for(int i = 0;i<((int)(n));i++)
#define reg(i,a,b) for(int i = ((int)(a));i<=((int)(b));i++)
#define irep(i,n) for(int i = ((int)(n)-1);i>=0;i--)
#define ireg(i,a,b) for(int i = ((int)(b));i>=((int)(a));i--)

/*
何故かエラー出る
	解決。配列範囲には気をつけよう
倍数かどうかだけ見れば良い。
*/

int n,m,a[100010],q,p[100010];
int l[100010],r[100010],fact[2010][305],dp[100010][305],zero[100010];
vector<ll> prime,pri2;

void prime_list(ll n,vector<ll>& prime_){
	ll visited[300010]={};
	reg(i,2,n){
		if(visited[i]==0){
			prime_.push_back(i);
			reg(j,2,n){
				if(j*i>n)break;
				visited[i*j]=1;
			}
		}
	}
}
void factoring(ll n,vector<pair<ll,ll>>& v,vector<ll>& p){
	rep(i,p.size()){
		int t=0;
		while(n%p[i]==0){
			n/=p[i];
			t++;
		}
		if(t>0)v.push_back({i,t});
	}
	if(n>1)v.push_back({10000,1});
}
ll modpow(__int128_t a, ll n,ll mo) {
	__int128_t r=1;
	while(n) r=r*((n%2)?a:1)%mo,a=a*a%mo,n>>=1;
	return r;
}

void init(){
	cin>>n;
	reg(i,1,n)cin>>a[i];
	cin>>q;
	rep(i,q)cin>>p[i]>>l[i]>>r[i];
	prime_list(2000,prime);
	m = prime.size();
	rep(j,m)fact[0][j]=0;
	reg(i,1,2000){
		ll t = i;
		rep(j,m){
			fact[i][j]=0;
			while(t%prime[j]==0){
				// if(prime[j]>=45)
				t/=prime[j];
				fact[i][j]++;
			}
		}
	}
	// 12ms
	zero[0]=0;
	rep(j,m)dp[0][j]=0;
	for(int i=1;i<=n;++i){
		zero[i]=zero[i-1];
		if(a[i]==0)zero[i]++;
		for(int j=0;j<m;++j)dp[i][j] = dp[i-1][j] + fact[a[i]][j];
	}
}

int main(void){
	init();
	for(int i=0;i<q;++i){
		bool ok=true;
		vector<pair<ll,ll>> v;
		factoring(p[i],v,prime);
		irep(j,v.size()){
			if(v[j].first>m){
				ok=false;
				break;
			}
			if(v[j].second > dp[r[i]][v[j].first] - dp[l[i]-1][v[j].first]){
				ok=false;
				break;
			}
		}
		cout<<(ok?"Yes":"NO")<<"\n";
	}
	return 0;
}
0