#362690 (C++14) No.854 公平なりんご分配

提出ソース

結果

問題	No.854 公平なりんご分配
ユーザー	やむなく
提出日時	2019-07-26 22:26:22
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	662 ms / 3,153 ms
コード長	3,362 bytes
コンパイル時間	1,825 ms
コンパイル使用メモリ	184,532 KB
実行使用メモリ	279,812 KB
最終ジャッジ日時	2024-07-02 07:52:44
合計ジャッジ時間	15,451 ms
ジャッジサーバーID （参考情報）	judge3 / judge5

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,816 KB
testcase_01	AC	2 ms 6,816 KB
testcase_02	AC	2 ms 6,816 KB
testcase_03	AC	1 ms 6,816 KB
testcase_04	AC	2 ms 6,940 KB
testcase_05	AC	2 ms 6,944 KB
testcase_06	AC	1 ms 6,944 KB
testcase_07	AC	2 ms 6,944 KB
testcase_08	AC	2 ms 6,940 KB
testcase_09	AC	2 ms 6,944 KB
testcase_10	AC	2 ms 6,940 KB
testcase_11	AC	2 ms 6,940 KB
testcase_12	AC	2 ms 6,940 KB
testcase_13	AC	2 ms 6,940 KB
testcase_14	AC	2 ms 6,940 KB
testcase_15	AC	2 ms 6,944 KB
testcase_16	AC	2 ms 6,944 KB
testcase_17	AC	1 ms 6,944 KB
testcase_18	AC	1 ms 6,944 KB
testcase_19	AC	2 ms 6,940 KB
testcase_20	AC	2 ms 6,944 KB
testcase_21	AC	2 ms 6,940 KB
testcase_22	AC	4 ms 6,940 KB
testcase_23	AC	4 ms 6,940 KB
testcase_24	AC	5 ms 6,940 KB
testcase_25	AC	3 ms 6,940 KB
testcase_26	AC	6 ms 6,940 KB
testcase_27	AC	5 ms 6,940 KB
testcase_28	AC	4 ms 6,944 KB
testcase_29	AC	3 ms 6,940 KB
testcase_30	AC	4 ms 6,944 KB
testcase_31	AC	6 ms 6,944 KB
testcase_32	AC	84 ms 68,608 KB
testcase_33	AC	66 ms 36,224 KB
testcase_34	AC	125 ms 78,336 KB
testcase_35	AC	101 ms 70,912 KB
testcase_36	AC	45 ms 7,296 KB
testcase_37	AC	82 ms 67,200 KB
testcase_38	AC	47 ms 41,984 KB
testcase_39	AC	174 ms 80,384 KB
testcase_40	AC	76 ms 34,304 KB
testcase_41	AC	102 ms 65,408 KB
testcase_42	AC	118 ms 77,824 KB
testcase_43	AC	136 ms 66,432 KB
testcase_44	AC	144 ms 75,136 KB
testcase_45	AC	109 ms 19,072 KB
testcase_46	AC	168 ms 74,624 KB
testcase_47	AC	61 ms 40,064 KB
testcase_48	AC	109 ms 77,184 KB
testcase_49	AC	104 ms 77,952 KB
testcase_50	AC	79 ms 66,944 KB
testcase_51	AC	167 ms 76,288 KB
testcase_52	AC	101 ms 36,992 KB
testcase_53	AC	65 ms 36,736 KB
testcase_54	AC	105 ms 39,680 KB
testcase_55	AC	48 ms 36,864 KB
testcase_56	AC	43 ms 37,632 KB
testcase_57	AC	62 ms 34,688 KB
testcase_58	AC	101 ms 12,928 KB
testcase_59	AC	43 ms 33,152 KB
testcase_60	AC	82 ms 34,688 KB
testcase_61	AC	32 ms 7,936 KB
testcase_62	AC	89 ms 32,768 KB
testcase_63	AC	77 ms 33,408 KB
testcase_64	AC	27 ms 12,032 KB
testcase_65	AC	99 ms 73,472 KB
testcase_66	AC	50 ms 20,480 KB
testcase_67	AC	86 ms 41,984 KB
testcase_68	AC	82 ms 21,120 KB
testcase_69	AC	91 ms 80,512 KB
testcase_70	AC	44 ms 34,432 KB
testcase_71	AC	49 ms 35,584 KB
testcase_72	AC	57 ms 7,424 KB
testcase_73	AC	82 ms 67,456 KB
testcase_74	AC	130 ms 72,960 KB
testcase_75	AC	71 ms 35,712 KB
testcase_76	AC	115 ms 66,304 KB
testcase_77	AC	115 ms 76,928 KB
testcase_78	AC	99 ms 22,016 KB
testcase_79	AC	129 ms 63,232 KB
testcase_80	AC	122 ms 64,384 KB
testcase_81	AC	110 ms 65,408 KB
testcase_82	AC	430 ms 279,808 KB
testcase_83	AC	648 ms 279,680 KB
testcase_84	AC	626 ms 279,808 KB
testcase_85	AC	456 ms 279,812 KB
testcase_86	AC	394 ms 279,804 KB
testcase_87	AC	452 ms 279,804 KB
testcase_88	AC	455 ms 279,804 KB
testcase_89	AC	446 ms 279,808 KB
testcase_90	AC	477 ms 279,804 KB
testcase_91	AC	445 ms 279,808 KB
testcase_92	AC	662 ms 279,808 KB
testcase_93	AC	640 ms 279,808 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

//
// Created by yamunaku on 2019/07/26.
//

#include <bits/stdc++.h>

using namespace std;

#define rep(i, n) for(int i = 0; i < (n); i++)
#define repl(i, l, r) for(int i = (l); i < (r); i++)
#define per(i, n) for(int i = ((n)-1); i >= 0; i--)
#define perl(i, l, r) for(int i = ((r)-1); i >= (l); i--)
#define all(x) (x).begin(),(x).end()
#define MOD9 998244353
#define MOD1 1000000007
#define IINF 1000000000
#define LINF 1000000000000000000
#define SP <<" "<<
#define CYES cout<<"Yes"<<endl
#define CNO cout<<"NO"<<endl
#define CFS cin.tie(0);ios::sync_with_stdio(false)
#define CST(x) cout<<fixed<<setprecision(x)

typedef long long ll;
typedef long double ld;
typedef vector<int> vi;
typedef vector<vector<int>> mti;
typedef vector<ll> vl;
typedef vector<vector<ll>> mtl;

struct SegmentTree {
private:
    int n;
    vector<short> node;

public:
    void init(vi &v) {
        int sz=v.size();
        n = 1; while(n < sz) n *= 2;
        node.resize(2*n-1, 0);
        for(int i=0; i<sz; i++) node[i+n-1] = v[i];
        for(int i=n-2; i>=0; i--) node[i] = node[2*i+1]+node[2*i+2];
    }

    int getsum(int a, int b, int k=0, int l=0, int r=-1) {
        if(r < 0) r = n;
        if(r <= a || b <= l) return 0;
        if(a <= l && r <= b) return node[k];
        int vl = getsum(a, b, 2*k+1, l, (l+r)/2);
        int vr = getsum(a, b, 2*k+2, (l+r)/2, r);
        return vl+vr;
    }
};

struct SegmentTree2{
private:
    int n;
    vector<int> node;

public:
    void init(vi &v) {
        int sz=v.size();
        n = 1; while(n < sz) n *= 2;
        node.resize(2*n-1, IINF);
        for(int i=0; i<sz; i++) node[i+n-1] = v[i];
        for(int i=n-2; i>=0; i--) node[i] = min(node[2*i+1],node[2*i+2]);
    }

    int getmin(int a, int b, int k=0, int l=0, int r=-1) {
        if(r < 0) r = n;
        if(r <= a || b <= l) return IINF;
        if(a <= l && r <= b) return node[k];
        int vl = getmin(a, b, 2*k+1, l, (l+r)/2);
        int vr = getmin(a, b, 2*k+2, (l+r)/2, r);
        return min(vl,vr);
    }
};


int main(){
    CFS;
    vi isp(2001,1);
    isp[0]=isp[1]=0;
    rep(i,2001){
        if(isp[i]){
            for(int j=i+i;j<2001;j+=i){
                isp[j]=0;
            }
        }
    }
    vi prime;
    rep(i,2001) if(isp[i]) prime.push_back(i);
    int pn=prime.size();
    vector<SegmentTree> seg(pn);
    int n;
    cin >> n;
    vi a(n);
    rep(i,n) cin >> a[i];
    mti paa(pn,vi(n,0));
    rep(i,n){
        if(a[i]==0) continue;
        rep(j,pn){
            while(a[i]%prime[j]==0){
                a[i]/=prime[j];
                paa[j][i]++;
            }
            if(a[i]==1) break;
        }
    }
    rep(i,pn){
        seg[i].init(paa[i]);
    }
    SegmentTree2 seg2;
    seg2.init(a);
    int q;
    cin >> q;
    rep(_,q){
        int p,l,r;
        cin >> p >> l >> r;
        if(seg2.getmin(l-1,r)==0){
            CYES;
            continue;
        }
        rep(i,pn){
            if(p%prime[i]==0){
                int cnt=0;
                while(p%prime[i]==0){
                    p/=prime[i];
                    cnt++;
                }
                if(cnt>seg[i].getsum(l-1,r)){
                    CNO;
                    goto next;
                }
                if(p==1) break;
            }
        }
        if(p==1) CYES;
        else CNO;
        next:;
    }
    return 0;
}

yukicoder

結果

テストケース

ソースコード