結果

問題 No.854 公平なりんご分配
ユーザー ahe100
提出日時 2019-07-26 23:02:47
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,562 bytes
コンパイル時間 1,574 ms
コンパイル使用メモリ 170,448 KB
実行使用メモリ 133,668 KB
最終ジャッジ日時 2024-07-02 09:02:32
合計ジャッジ時間 22,370 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 81 TLE * 1 -- * 10
権限があれば一括ダウンロードができます

ソースコード

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;

void prime_list(ll n,vector<ll>& prime_){
	ll visited[2010]={};
	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;
			}
		}
	}
}
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=m-1;j>=0;--j)dp[i][j] = dp[i-1][j] + fact[a[i]][j];
	}
}

int main(void){
	init();
	for(int i=0;i<q;++i){
		ll tmp=1;
		if(zero[r[i]]-zero[l[i]-1]>0)tmp=0;
		for(int j=m-1;j>=0;--j){
			if(tmp==0)break;
			tmp *= modpow(prime[j], dp[r[i]][j] - dp[l[i]-1][j], p[i]);
			tmp %= p[i];
		}
		cout<<(tmp==0?"Yes":"NO")<<"\n";
	}
	return 0;
}
0