結果

問題 No.854 公平なりんご分配
コンテスト
ユーザー ahe100
提出日時 2019-07-26 22:41:25
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,472 bytes
コンパイル時間 1,646 ms
コンパイル使用メモリ 171,984 KB
実行使用メモリ 136,420 KB
最終ジャッジ日時 2024-07-02 08:35:38
合計ジャッジ時間 32,069 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 1
other AC * 1 WA * 79 TLE * 1 -- * 11
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define rep(i,n) for(ll i = 0;i<((ll)(n));i++)
#define reg(i,a,b) for(ll i = ((ll)(a));i<=((ll)(b));i++)
#define irep(i,n) for(ll i = ((ll)(n)-1);i>=0;i--)
#define ireg(i,a,b) for(ll i = ((ll)(b));i>=((ll)(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){
				t/=prime[j];
				fact[i][j]++;
			}
		}
	}
	zero[0]=0;
	rep(j,m)dp[0][j]=0;
	reg(i,1,n){
		zero[i]=zero[i-1];
		if(a[i]==0)zero[i]++;
		rep(j,m)dp[i][j] = dp[i-1][j] + fact[a[i]][j];
	}
	// reg(i,1,n){
	// 	rep(j,3)cerr<<dp[i][j]<<" ";
	// 	cerr<<endl;
	// }
}

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