結果

問題 No.2207 pCr検査
ユーザー 👑 emthrm
提出日時 2023-02-04 14:12:34
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 2,551 bytes
コンパイル時間 2,519 ms
コンパイル使用メモリ 252,380 KB
実行使用メモリ 123,152 KB
最終ジャッジ日時 2024-07-03 11:53:20
合計ジャッジ時間 43,333 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 29 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
using namespace std;
#define FOR(i,m,n) for(int i=(m);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()
using ll = long long;
constexpr int INF = 0x3f3f3f3f;
constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;
constexpr double EPS = 1e-8;
constexpr int MOD = 998244353;
// constexpr int MOD = 1000000007;
constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1};
constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1};
constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1};
template <typename T, typename U>
inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; }
template <typename T, typename U>
inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; }
struct IOSetup {
IOSetup() {
std::cin.tie(nullptr);
std::ios_base::sync_with_stdio(false);
std::cout << fixed << setprecision(20);
}
} iosetup;
std::vector<int> prime_sieve(const int n, const bool get_only_prime) {
std::vector<int> smallest_prime_factor(n + 1), prime;
std::iota(smallest_prime_factor.begin(), smallest_prime_factor.end(), 0);
for (int i = 2; i <= n; ++i) {
if (smallest_prime_factor[i] == i) [[unlikely]] prime.emplace_back(i);
for (const int p : prime) {
if (i * p > n || p > smallest_prime_factor[i]) break;
smallest_prime_factor[i * p] = p;
}
}
return get_only_prime ? prime : smallest_prime_factor;
}
struct OsaK {
const std::vector<int> smallest_prime_factor;
explicit OsaK(const int n) : smallest_prime_factor(prime_sieve(n, false)) {}
std::vector<std::pair<int, int>> query(int n) const {
std::vector<std::pair<int, int>> res;
while (n > 1) {
const int prime = smallest_prime_factor[n];
int exponent = 0;
for (; smallest_prime_factor[n] == prime; n /= prime) {
++exponent;
}
res.emplace_back(prime, exponent);
}
return res;
}
};
int main() {
constexpr int P = 10000000;
int k; cin >> k;
int nums[P]{}, p_k = 0, nonzero = 0;
while (k--) {
int p, e; cin >> p >> e;
nums[p] = -e;
++nonzero;
chmax(p_k, p);
}
const OsaK osa_k(p_k);
FOR(r, 1, p_k) {
for (const auto& [p, e] : osa_k.query(p_k - r + 1)) {
if (nums[p] == 0) ++nonzero;
if ((nums[p] += e) == 0) --nonzero;
}
for (const auto& [p, e] : osa_k.query(r)) {
if (nums[p] == 0) ++nonzero;
if ((nums[p] -= e) == 0) --nonzero;
}
if (nonzero == 0) {
cout << p_k << ' ' << r << '\n';
return 0;
}
}
cout << "-1 -1\n";
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0