結果
問題 | No.2570 最大最大公約数 |
ユーザー |
![]() |
提出日時 | 2023-12-02 16:11:03 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 364 ms / 2,000 ms |
コード長 | 2,525 bytes |
コンパイル時間 | 4,389 ms |
コンパイル使用メモリ | 261,136 KB |
最終ジャッジ日時 | 2025-02-18 05:15:44 |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 28 |
ソースコード
#include<bits/stdc++.h>#include<atcoder/all>using namespace std;using int64 = long long;constexpr int mod = 998244353;constexpr int64 infll = (1LL << 62) - 1;constexpr int inf = (1 << 30) - 1;struct IoSetup {IoSetup() {cin.tie(nullptr);ios::sync_with_stdio(false);cout << fixed << setprecision(10);cerr << fixed << setprecision(10);}} iosetup;template<typename T1, typename T2>ostream &operator<<(ostream &os, const pair<T1, T2> &p) {os << p.first << " " << p.second;return os;}template<typename T1, typename T2>istream &operator>>(istream &is, pair<T1, T2> &p) {is >> p.first >> p.second;return is;}template<typename T>ostream &operator<<(ostream &os, const vector<T> &v) {for (int i = 0; i < (int) v.size(); i++) {os << v[i] << (i + 1 != v.size() ? " " : "");}return os;}template<typename T>istream &operator>>(istream &is, vector<T> &v) {for (T &in : v) is >> in;return is;}template<typename T1, typename T2>inline bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); }template<typename T1, typename T2>inline bool chmin(T1 &a, T2 b) { return a > b && (a = b, true); }template<typename T = int64>vector<T> make_v(size_t a) {return vector<T>(a);}template<typename T, typename... Ts>auto make_v(size_t a, Ts... ts) {return vector<decltype(make_v<T>(ts...))>(a, make_v<T>(ts...));}template<typename T, typename V>typename enable_if<is_class<T>::value == 0>::type fill_v(T &t, const V &v) {t = v;}template<typename T, typename V>typename enable_if<is_class<T>::value != 0>::type fill_v(T &t, const V &v) {for (auto &e : t) fill_v(e, v);}template<typename F>struct FixPoint : F {explicit FixPoint(F &&f) : F(forward<F>(f)) {}template<typename... Args>decltype(auto) operator()(Args &&... args) const {return F::operator()(*this, forward<Args>(args)...);}};template<typename F>inline decltype(auto) MFP(F &&f) {return FixPoint<F>{forward<F>(f)};}int main() {int N, K;cin >> N >> K;vector< int64 > A(N);cin >> A;vector< set< int64 > > dp(K + 1);dp[K].emplace(0);for(auto& a : A) {vector< set< int64 > > dp2(K + 1);for(int i = K; i >= 0; i--) {for(int64 j = max(0ll, a - i); j <= a + i; j++) {for(auto& p : dp[i]) {auto nxt = __gcd(p, j);dp2[i - abs(a - j)].emplace(__gcd(p, j));}}}dp = dp2;}int64 ret = 0;for(auto& p : dp) chmax(ret, *p.rbegin());cout << ret << endl;}