結果

問題 No.854 公平なりんご分配
ユーザー LayCurseLayCurse
提出日時 2019-08-21 05:39:32
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 293 ms / 3,153 ms
コード長 4,543 bytes
コンパイル時間 2,865 ms
コンパイル使用メモリ 216,228 KB
実行使用メモリ 8,276 KB
最終ジャッジ日時 2024-10-07 14:57:51
合計ジャッジ時間 11,063 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
6,816 KB
testcase_01 AC 2 ms
6,820 KB
testcase_02 AC 2 ms
6,816 KB
testcase_03 AC 2 ms
6,824 KB
testcase_04 AC 2 ms
6,824 KB
testcase_05 AC 2 ms
7,628 KB
testcase_06 AC 2 ms
7,508 KB
testcase_07 AC 2 ms
6,824 KB
testcase_08 AC 2 ms
6,820 KB
testcase_09 AC 2 ms
6,820 KB
testcase_10 AC 2 ms
7,632 KB
testcase_11 AC 2 ms
6,816 KB
testcase_12 AC 3 ms
6,816 KB
testcase_13 AC 2 ms
6,816 KB
testcase_14 AC 2 ms
7,508 KB
testcase_15 AC 2 ms
6,820 KB
testcase_16 AC 3 ms
7,632 KB
testcase_17 AC 2 ms
6,820 KB
testcase_18 AC 2 ms
6,816 KB
testcase_19 AC 3 ms
7,628 KB
testcase_20 AC 2 ms
7,632 KB
testcase_21 AC 2 ms
6,824 KB
testcase_22 AC 4 ms
7,624 KB
testcase_23 AC 3 ms
6,816 KB
testcase_24 AC 4 ms
6,816 KB
testcase_25 AC 3 ms
6,816 KB
testcase_26 AC 4 ms
6,820 KB
testcase_27 AC 4 ms
7,508 KB
testcase_28 AC 3 ms
6,816 KB
testcase_29 AC 3 ms
7,756 KB
testcase_30 AC 4 ms
6,820 KB
testcase_31 AC 4 ms
6,820 KB
testcase_32 AC 33 ms
6,816 KB
testcase_33 AC 34 ms
6,820 KB
testcase_34 AC 63 ms
6,820 KB
testcase_35 AC 47 ms
6,816 KB
testcase_36 AC 40 ms
6,816 KB
testcase_37 AC 32 ms
6,816 KB
testcase_38 AC 24 ms
6,824 KB
testcase_39 AC 88 ms
6,820 KB
testcase_40 AC 40 ms
7,764 KB
testcase_41 AC 44 ms
6,816 KB
testcase_42 AC 60 ms
7,764 KB
testcase_43 AC 66 ms
6,820 KB
testcase_44 AC 69 ms
6,816 KB
testcase_45 AC 69 ms
6,816 KB
testcase_46 AC 86 ms
6,820 KB
testcase_47 AC 35 ms
6,816 KB
testcase_48 AC 55 ms
6,820 KB
testcase_49 AC 52 ms
6,820 KB
testcase_50 AC 28 ms
6,820 KB
testcase_51 AC 84 ms
7,756 KB
testcase_52 AC 121 ms
6,816 KB
testcase_53 AC 53 ms
6,820 KB
testcase_54 AC 118 ms
6,820 KB
testcase_55 AC 25 ms
7,640 KB
testcase_56 AC 24 ms
6,816 KB
testcase_57 AC 54 ms
6,820 KB
testcase_58 AC 90 ms
6,816 KB
testcase_59 AC 21 ms
6,820 KB
testcase_60 AC 79 ms
6,820 KB
testcase_61 AC 43 ms
7,624 KB
testcase_62 AC 85 ms
6,820 KB
testcase_63 AC 56 ms
6,816 KB
testcase_64 AC 24 ms
7,636 KB
testcase_65 AC 52 ms
7,632 KB
testcase_66 AC 56 ms
7,508 KB
testcase_67 AC 80 ms
6,816 KB
testcase_68 AC 97 ms
6,816 KB
testcase_69 AC 21 ms
6,820 KB
testcase_70 AC 22 ms
7,504 KB
testcase_71 AC 27 ms
6,820 KB
testcase_72 AC 61 ms
6,820 KB
testcase_73 AC 29 ms
7,620 KB
testcase_74 AC 98 ms
6,816 KB
testcase_75 AC 57 ms
6,816 KB
testcase_76 AC 59 ms
6,816 KB
testcase_77 AC 66 ms
7,628 KB
testcase_78 AC 111 ms
7,628 KB
testcase_79 AC 105 ms
6,816 KB
testcase_80 AC 98 ms
6,820 KB
testcase_81 AC 54 ms
6,820 KB
testcase_82 AC 25 ms
6,820 KB
testcase_83 AC 117 ms
6,860 KB
testcase_84 AC 106 ms
7,244 KB
testcase_85 AC 140 ms
7,752 KB
testcase_86 AC 22 ms
7,636 KB
testcase_87 AC 25 ms
6,824 KB
testcase_88 AC 28 ms
7,636 KB
testcase_89 AC 27 ms
7,636 KB
testcase_90 AC 23 ms
7,764 KB
testcase_91 AC 24 ms
6,816 KB
testcase_92 AC 293 ms
6,864 KB
testcase_93 AC 292 ms
8,276 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize ("Ofast")
#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;
}
template<class T> inline void walloc1d(T **arr, int x, void **mem = &wmem){
  static int skip[16]={0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1};
  (*mem) = (void*)( ((char*)(*mem)) + skip[((unsigned long long)(*mem)) & 15] );
  (*arr)=(T*)(*mem);
  (*mem)=((*arr)+x);
}
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;
  walloc1d(&isprime, r, &mem);
  walloc1d(&sf, r, &mem);
  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 20190820-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