結果

問題 No.2207 pCr検査
ユーザー 👑 NachiaNachia
提出日時 2023-02-03 22:13:40
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,332 ms / 3,000 ms
コード長 4,209 bytes
コンパイル時間 1,147 ms
コンパイル使用メモリ 92,568 KB
実行使用メモリ 97,080 KB
最終ジャッジ日時 2024-07-02 20:15:23
合計ジャッジ時間 22,284 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 607 ms
91,880 KB
testcase_03 AC 359 ms
52,928 KB
testcase_04 AC 412 ms
68,168 KB
testcase_05 AC 686 ms
87,672 KB
testcase_06 AC 407 ms
82,228 KB
testcase_07 AC 291 ms
52,112 KB
testcase_08 AC 174 ms
30,512 KB
testcase_09 AC 500 ms
68,052 KB
testcase_10 AC 223 ms
51,040 KB
testcase_11 AC 476 ms
66,192 KB
testcase_12 AC 897 ms
76,564 KB
testcase_13 AC 506 ms
47,696 KB
testcase_14 AC 98 ms
14,468 KB
testcase_15 AC 1,052 ms
90,132 KB
testcase_16 AC 460 ms
43,240 KB
testcase_17 AC 794 ms
70,512 KB
testcase_18 AC 596 ms
56,552 KB
testcase_19 AC 285 ms
30,152 KB
testcase_20 AC 460 ms
45,712 KB
testcase_21 AC 1,101 ms
92,648 KB
testcase_22 AC 1,237 ms
96,952 KB
testcase_23 AC 1,332 ms
97,080 KB
testcase_24 AC 124 ms
18,088 KB
testcase_25 AC 459 ms
40,976 KB
testcase_26 AC 377 ms
37,980 KB
testcase_27 AC 789 ms
71,164 KB
testcase_28 AC 829 ms
95,292 KB
testcase_29 AC 841 ms
95,296 KB
testcase_30 AC 832 ms
95,296 KB
testcase_31 AC 846 ms
95,300 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#line 1 "Main.cpp"
#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <atcoder/modint>
#line 2 "nachia\\math\\prime-sieve-explicit.hpp"

#line 5 "nachia\\math\\prime-sieve-explicit.hpp"
#include <cassert>
#line 7 "nachia\\math\\prime-sieve-explicit.hpp"


namespace nachia{

namespace prime_sieve_explicit_internal{
    std::vector<bool> isprime = { false }; // a[x] := isprime(2x+1)

    void CalcIsPrime(int z){
        if((int)isprime.size() *2+1 < z+1){
            int new_z = isprime.size();
            while(new_z*2+1 < z+1) new_z *= 2;
            z = new_z-1;
            isprime.resize(z+1, true);
            for(int i=1; i*(i+1)*2<=z; i++) if(isprime[i]){
                for(int j=i*(i+1)*2; j<=z; j+=i*2+1) isprime[j] = false;
            }
        }
    }
    
    std::vector<int> prime_list = {2};
    int prime_list_max = 0;

    void CalcPrimeList(int z){
        while((int)prime_list.size() < z){
            if((int)isprime.size() <= prime_list_max + 1) CalcIsPrime(prime_list_max + 1);
            for(int p=prime_list_max+1; p<(int)isprime.size(); p++){
                if(isprime[p]) prime_list.push_back(p*2+1);
            }
            prime_list_max = isprime.size() - 1;
        }
    }

    void CalcPrimeListUntil(int z){
        if(prime_list_max < z){
            CalcIsPrime(z);
            for(int p=prime_list_max+1; p<(int)isprime.size(); p++){
                if(isprime[p]) prime_list.push_back(p*2+1);
            }
            prime_list_max = isprime.size() - 1;
        }
    }

}


bool IsprimeExplicit(int n){
    using namespace prime_sieve_explicit_internal;
    if(n == 2) return true;
    if(n % 2 == 0) return false;
    CalcIsPrime(n);
    return isprime[(n-1)/2];
}

int NthPrimeExplicit(int n){
    using namespace prime_sieve_explicit_internal;
    CalcPrimeList(n);
    return prime_list[n];
}

int PrimeCountingExplicit(int n){
    using namespace prime_sieve_explicit_internal;
    if(n < 2) return 0;
    CalcPrimeListUntil(n);
    auto res = std::upper_bound(prime_list.begin(), prime_list.end(), n) - prime_list.begin();
    return (int)res;
}

// [l, r)
std::vector<bool> SegmentedSieveExplicit(long long l, long long r){
    assert(0 <= l); assert(l <= r);
    long long d = r - l;
    if(d == 0) return {};
    std::vector<bool> res(d, true);
    for(long long p=2; p*p<=r; p++) if(IsprimeExplicit(p)){
        long long il = (l+p-1)/p, ir = (r+p-1)/p;
        if(il <= p) il = p;
        for(long long i=il; i<ir; i++) res[i*p-l] = false;
    }
    if(l < 2) for(long long p=l; p<2 && p<r; p++) res[l-p] = false;
    return res;
}


} // namespace nachia
#line 7 "Main.cpp"
using namespace std;
using i32 = int;
using u32 = unsigned int;
using i64 = long long;
using u64 = unsigned long long;
#define rep(i,n) for(int i=0; i<(int)(n); i++)
const i64 INF = 1001001001001001001;

using Modint = atcoder::static_modint<998244353>;

int main(){
    int k; cin >> k;
    vector<int> P(k), E(k);
    rep(i,k) cin >> P[i] >> E[i];
    int Pmax = P.back();
    int PCnt = nachia::PrimeCountingExplicit(Pmax);
    vector<int> Pex(PCnt);
    rep(i,k) Pex[nachia::PrimeCountingExplicit(P[i])-1] = E[i];
    vector<pair<int,int>> nx(Pmax+1);
    rep(i,PCnt){
        int p = nachia::NthPrimeExplicit(i);
        for(int q=1; q*p<=Pmax; q++) nx[p*q] = {i,q};
    }
    vector<int> Pnow(PCnt);
    int MatchCount = 0;
    rep(i,PCnt) if(Pnow[i] == Pex[i]) MatchCount++;
    for(int r=1; r<=Pmax; r++){
        for(int x=Pmax-r+1; x>1; x=nx[x].second){
            int p = nx[x].first;
            if(Pnow[p] == Pex[p]) MatchCount--;
            Pnow[p]++;
            if(Pnow[p] == Pex[p]) MatchCount++;
        }
        for(int x=r; x>1; x=nx[x].second){
            int p = nx[x].first;
            if(Pnow[p] == Pex[p]) MatchCount--;
            Pnow[p]--;
            if(Pnow[p] == Pex[p]) MatchCount++;
        }
        if(MatchCount == PCnt){
            cout << Pmax << ' ' << r << endl;
            return 0;
        }
    }
    cout << "-1 -1\n";
    return 0;
}



struct ios_do_not_sync{
    ios_do_not_sync(){
        ios::sync_with_stdio(false);
        cin.tie(nullptr);
    }
} ios_do_not_sync_instance;
0