結果

問題 No.854 公平なりんご分配
ユーザー LayCurseLayCurse
提出日時 2019-07-26 22:08:42
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 333 ms / 3,153 ms
コード長 4,128 bytes
コンパイル時間 2,028 ms
コンパイル使用メモリ 201,192 KB
実行使用メモリ 7,500 KB
最終ジャッジ日時 2023-09-15 01:40:56
合計ジャッジ時間 10,257 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge13
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,440 KB
testcase_01 AC 2 ms
5,472 KB
testcase_02 AC 2 ms
5,412 KB
testcase_03 AC 3 ms
5,432 KB
testcase_04 AC 2 ms
5,420 KB
testcase_05 AC 2 ms
5,416 KB
testcase_06 AC 2 ms
5,408 KB
testcase_07 AC 2 ms
5,436 KB
testcase_08 AC 2 ms
5,668 KB
testcase_09 AC 2 ms
5,428 KB
testcase_10 AC 2 ms
5,440 KB
testcase_11 AC 2 ms
5,440 KB
testcase_12 AC 2 ms
5,484 KB
testcase_13 AC 2 ms
5,468 KB
testcase_14 AC 2 ms
5,440 KB
testcase_15 AC 2 ms
5,488 KB
testcase_16 AC 2 ms
5,428 KB
testcase_17 AC 2 ms
5,532 KB
testcase_18 AC 2 ms
5,440 KB
testcase_19 AC 2 ms
5,444 KB
testcase_20 AC 2 ms
5,420 KB
testcase_21 AC 2 ms
5,428 KB
testcase_22 AC 4 ms
5,464 KB
testcase_23 AC 3 ms
5,444 KB
testcase_24 AC 5 ms
5,452 KB
testcase_25 AC 2 ms
5,440 KB
testcase_26 AC 5 ms
5,452 KB
testcase_27 AC 3 ms
5,516 KB
testcase_28 AC 3 ms
5,440 KB
testcase_29 AC 2 ms
5,412 KB
testcase_30 AC 3 ms
5,444 KB
testcase_31 AC 4 ms
5,428 KB
testcase_32 AC 39 ms
5,652 KB
testcase_33 AC 39 ms
5,696 KB
testcase_34 AC 72 ms
5,872 KB
testcase_35 AC 53 ms
5,756 KB
testcase_36 AC 42 ms
5,696 KB
testcase_37 AC 38 ms
5,644 KB
testcase_38 AC 29 ms
5,652 KB
testcase_39 AC 102 ms
6,184 KB
testcase_40 AC 44 ms
5,724 KB
testcase_41 AC 50 ms
5,832 KB
testcase_42 AC 66 ms
5,800 KB
testcase_43 AC 71 ms
6,080 KB
testcase_44 AC 78 ms
6,048 KB
testcase_45 AC 75 ms
6,200 KB
testcase_46 AC 96 ms
6,300 KB
testcase_47 AC 39 ms
5,656 KB
testcase_48 AC 62 ms
5,796 KB
testcase_49 AC 59 ms
5,928 KB
testcase_50 AC 32 ms
5,596 KB
testcase_51 AC 96 ms
6,084 KB
testcase_52 AC 132 ms
6,224 KB
testcase_53 AC 58 ms
5,824 KB
testcase_54 AC 129 ms
5,964 KB
testcase_55 AC 28 ms
5,792 KB
testcase_56 AC 27 ms
5,680 KB
testcase_57 AC 59 ms
5,848 KB
testcase_58 AC 100 ms
5,908 KB
testcase_59 AC 26 ms
5,664 KB
testcase_60 AC 88 ms
5,784 KB
testcase_61 AC 47 ms
5,680 KB
testcase_62 AC 95 ms
6,100 KB
testcase_63 AC 64 ms
5,860 KB
testcase_64 AC 28 ms
5,532 KB
testcase_65 AC 60 ms
5,716 KB
testcase_66 AC 62 ms
5,680 KB
testcase_67 AC 92 ms
5,908 KB
testcase_68 AC 105 ms
5,912 KB
testcase_69 AC 25 ms
5,736 KB
testcase_70 AC 25 ms
5,596 KB
testcase_71 AC 32 ms
5,636 KB
testcase_72 AC 67 ms
5,860 KB
testcase_73 AC 34 ms
5,720 KB
testcase_74 AC 111 ms
5,964 KB
testcase_75 AC 65 ms
5,748 KB
testcase_76 AC 74 ms
5,932 KB
testcase_77 AC 77 ms
6,012 KB
testcase_78 AC 124 ms
6,036 KB
testcase_79 AC 118 ms
5,972 KB
testcase_80 AC 111 ms
6,160 KB
testcase_81 AC 63 ms
5,828 KB
testcase_82 AC 27 ms
6,744 KB
testcase_83 AC 133 ms
7,500 KB
testcase_84 AC 119 ms
7,444 KB
testcase_85 AC 155 ms
7,480 KB
testcase_86 AC 24 ms
6,700 KB
testcase_87 AC 27 ms
6,836 KB
testcase_88 AC 28 ms
6,828 KB
testcase_89 AC 27 ms
6,704 KB
testcase_90 AC 27 ms
6,696 KB
testcase_91 AC 27 ms
6,764 KB
testcase_92 AC 332 ms
7,448 KB
testcase_93 AC 333 ms
7,428 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