結果

問題 No.854 公平なりんご分配
ユーザー LayCurseLayCurse
提出日時 2019-07-26 22:08:42
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 321 ms / 3,153 ms
コード長 4,128 bytes
コンパイル時間 2,294 ms
コンパイル使用メモリ 203,520 KB
実行使用メモリ 6,784 KB
最終ジャッジ日時 2024-07-02 07:27:25
合計ジャッジ時間 9,939 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 3 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 2 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 2 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
testcase_18 AC 2 ms
5,376 KB
testcase_19 AC 2 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
testcase_21 AC 2 ms
5,376 KB
testcase_22 AC 3 ms
5,376 KB
testcase_23 AC 3 ms
5,376 KB
testcase_24 AC 4 ms
5,376 KB
testcase_25 AC 3 ms
5,376 KB
testcase_26 AC 5 ms
5,376 KB
testcase_27 AC 3 ms
5,376 KB
testcase_28 AC 3 ms
5,376 KB
testcase_29 AC 2 ms
5,376 KB
testcase_30 AC 3 ms
5,376 KB
testcase_31 AC 4 ms
5,376 KB
testcase_32 AC 36 ms
5,376 KB
testcase_33 AC 38 ms
5,376 KB
testcase_34 AC 72 ms
5,376 KB
testcase_35 AC 51 ms
5,376 KB
testcase_36 AC 42 ms
5,376 KB
testcase_37 AC 35 ms
5,376 KB
testcase_38 AC 26 ms
5,376 KB
testcase_39 AC 99 ms
5,376 KB
testcase_40 AC 43 ms
5,376 KB
testcase_41 AC 48 ms
5,376 KB
testcase_42 AC 63 ms
5,376 KB
testcase_43 AC 68 ms
5,376 KB
testcase_44 AC 74 ms
5,376 KB
testcase_45 AC 74 ms
5,376 KB
testcase_46 AC 93 ms
5,376 KB
testcase_47 AC 38 ms
5,376 KB
testcase_48 AC 58 ms
5,376 KB
testcase_49 AC 56 ms
5,376 KB
testcase_50 AC 30 ms
5,376 KB
testcase_51 AC 93 ms
5,376 KB
testcase_52 AC 131 ms
5,376 KB
testcase_53 AC 56 ms
5,376 KB
testcase_54 AC 125 ms
5,376 KB
testcase_55 AC 27 ms
5,376 KB
testcase_56 AC 25 ms
5,376 KB
testcase_57 AC 58 ms
5,376 KB
testcase_58 AC 102 ms
5,376 KB
testcase_59 AC 24 ms
5,376 KB
testcase_60 AC 85 ms
5,376 KB
testcase_61 AC 45 ms
5,376 KB
testcase_62 AC 96 ms
5,376 KB
testcase_63 AC 65 ms
5,376 KB
testcase_64 AC 27 ms
5,376 KB
testcase_65 AC 57 ms
5,376 KB
testcase_66 AC 59 ms
5,376 KB
testcase_67 AC 88 ms
5,376 KB
testcase_68 AC 104 ms
5,376 KB
testcase_69 AC 23 ms
5,376 KB
testcase_70 AC 24 ms
5,376 KB
testcase_71 AC 30 ms
5,376 KB
testcase_72 AC 68 ms
5,376 KB
testcase_73 AC 32 ms
5,376 KB
testcase_74 AC 107 ms
5,376 KB
testcase_75 AC 63 ms
5,376 KB
testcase_76 AC 68 ms
5,376 KB
testcase_77 AC 74 ms
5,376 KB
testcase_78 AC 126 ms
5,376 KB
testcase_79 AC 118 ms
5,376 KB
testcase_80 AC 109 ms
5,376 KB
testcase_81 AC 61 ms
5,376 KB
testcase_82 AC 27 ms
5,888 KB
testcase_83 AC 130 ms
6,656 KB
testcase_84 AC 118 ms
6,528 KB
testcase_85 AC 157 ms
6,656 KB
testcase_86 AC 23 ms
6,016 KB
testcase_87 AC 27 ms
5,888 KB
testcase_88 AC 27 ms
5,888 KB
testcase_89 AC 28 ms
5,888 KB
testcase_90 AC 26 ms
5,888 KB
testcase_91 AC 27 ms
6,016 KB
testcase_92 AC 316 ms
6,784 KB
testcase_93 AC 321 ms
6,784 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
void *wmem;
template<class S, class T> inline S min_L(S a,T b){
  return a<=b?a:b;
}
inline void rd(int &x){
  int k, m=0;
  x=0;
  for(;;){
    k = getchar_unlocked();
    if(k=='-'){
      m=1;
      break;
    }
    if('0'<=k&&k<='9'){
      x=k-'0';
      break;
    }
  }
  for(;;){
    k = getchar_unlocked();
    if(k<'0'||k>'9'){
      break;
    }
    x=x*10+k-'0';
  }
  if(m){
    x=-x;
  }
}
inline void wt_L(char a){
  putchar_unlocked(a);
}
inline void wt_L(const char c[]){
  int i=0;
  for(i=0;c[i]!='\0';i++){
    putchar_unlocked(c[i]);
  }
}
int Prime_L(int N, int res[], void *mem=wmem){
  bool *isprime;
  const int r=23000;
  int a, b, i, *sf, ss=1, sz=1;
  isprime = (bool*)mem;
  sf = (int*)(isprime + r);
  N /= 2;
  res[0] = 2;
  b =min_L(r, N);
  for(i=1;i<b;i++){
    isprime[i] = 1;
  }
  for(i=1;i<b;i++){
    if(isprime[i]){
      res[sz++] = 2*i+1;
      sf[ss] = 2*i*(i+1);
      if(sf[ss] < N){
        while(sf[ss] < r){
          isprime[sf[ss]] = 0;
          sf[ss] += res[ss];
        }
        ss++;
      }
    }
  }
  for(a=r; a<N; a+=r){
    b =min_L(a + r, N);
    isprime -= r;
    for(i=a;i<b;i++){
      isprime[i] = 1;
    }
    for(i=1;i<ss;i++){
      while(sf[i] < b){
        isprime[sf[i]] = 0;
        sf[i] += res[i];
      }
    }
    for(i=a;i<b;i++){
      if(isprime[i]){
        res[sz++] = 2*i+1;
      }
    }
  }
  return sz;
}
char memarr[96000000];
int N;
int A[100000];
int Q;
int P[100000];
int L[100000];
int R[100000];
int res[100000];
int ok[100001];
int fc[100001];
int ps;
int prm[1000];
int cnt[100000];
int main(){
  int i, j, k;
  wmem = memarr;
  ps =Prime_L(2000, prm);
  rd(N);
  {
    int Lj4PdHRW;
    for(Lj4PdHRW=0;Lj4PdHRW<N;Lj4PdHRW++){
      rd(A[Lj4PdHRW]);
    }
  }
  rd(Q);
  {
    int KL2GvlyY;
    for(KL2GvlyY=0;KL2GvlyY<Q;KL2GvlyY++){
      rd(P[KL2GvlyY]);
      rd(L[KL2GvlyY]);L[KL2GvlyY] += (-1);
      rd(R[KL2GvlyY]);
    }
  }
  for(i=0;i<Q;i++){
    res[i] = 1;
  }
  for(i=0;i<N;i++){
    if(A[i]==0){
      ok[i+1] = 1;
    }
  }
  for(i=0;i<N;i++){
    ok[i+1] += ok[i];
  }
  for(i=0;i<Q;i++){
    if(ok[R[i]]-ok[L[i]]){
      res[i] = 2;
      continue;
    }
  }
  for(k=0;k<ps;k++){
    j = 0;
    for(i=0;i<Q;i++){
      if(res[i]==1){
        cnt[i] = 0;
        while(P[i]%prm[k]==0){
          P[i] /= prm[k];
          cnt[i]++;
          j++;
        }
      }
    }
    if(j==0){
      continue;
    }
    for(i=0;i<N+1;i++){
      fc[i] = 0;
    }
    for(i=0;i<N;i++){
      if(ok[i+1]-ok[i]){
        continue;
      }
      while(A[i]%prm[k]==0){
        A[i] /= prm[k];
        fc[i+1]++;
      }
    }
    for(i=0;i<N;i++){
      fc[i+1] += fc[i];
    }
    for(i=0;i<Q;i++){
      if(res[i]==1){
        if(fc[R[i]]-fc[L[i]] < cnt[i]){
          res[i] = 0;
        }
      }
    }
  }
  for(i=0;i<Q;i++){
    if(res[i]==1 && P[i] > 1){
      res[i] = 0;
    }
  }
  for(i=0;i<Q;i++){
    if(res[i]){
      wt_L("Yes");
      wt_L('\n');
    }
    else{
      wt_L("NO");
      wt_L('\n');
    }
  }
  return 0;
}
// cLay varsion 20190721-1

// --- original code ---
// int N, A[1d5];
// 
// int Q, P[1d5], L[1d5], R[1d5];
// int res[1d5];
// 
// int ok[100001], fc[100001];
// 
// int ps, prm[1000];
// int cnt[1d5];
// 
// {
//   int i, j, k;
//   
//   ps = Prime(2000, prm);
// 
//   rd(N,A(N),Q,(P,L--,R)(Q));
//   rep(i,Q) res[i] = 1;
// 
//   rep(i,N) if(A[i]==0) ok[i+1] = 1;
// 
//   rep(i,N) ok[i+1] += ok[i];
//   rep(i,Q) if(ok[R[i]]-ok[L[i]]){ res[i] = 2; continue; }
// 
//   rep(k,ps){
//     j = 0;
//     rep(i,Q) if(res[i]==1){
//       cnt[i] = 0;
//       while(P[i]%prm[k]==0) P[i] /= prm[k], cnt[i]++, j++;
//     }
//     if(j==0) continue;
// 
//     rep(i,N+1) fc[i] = 0;
//     rep(i,N){
//       if(ok[i+1]-ok[i]) continue;
//       while(A[i]%prm[k]==0) A[i] /= prm[k], fc[i+1]++;
//     }
//     rep(i,N) fc[i+1] += fc[i];
// 
//     rep(i,Q) if(res[i]==1){
//       if(fc[R[i]]-fc[L[i]] < cnt[i]) res[i] = 0;
//     }
//   }
// 
//   rep(i,Q) if(res[i]==1 && P[i] > 1) res[i] = 0;
// 
//   rep(i,Q) wt(if[res[i],"Yes","NO"]);
// }
0