結果

問題 No.1529 Constant Lcm
ユーザー rokahikou1rokahikou1
提出日時 2021-06-04 20:58:04
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 706 ms / 3,000 ms
コード長 13,352 bytes
コンパイル時間 2,173 ms
コンパイル使用メモリ 195,288 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-04-30 00:45:35
合計ジャッジ時間 10,036 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 3 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 616 ms
10,240 KB
testcase_11 AC 217 ms
6,144 KB
testcase_12 AC 10 ms
5,376 KB
testcase_13 AC 511 ms
9,216 KB
testcase_14 AC 297 ms
7,168 KB
testcase_15 AC 341 ms
7,424 KB
testcase_16 AC 262 ms
6,528 KB
testcase_17 AC 133 ms
5,376 KB
testcase_18 AC 146 ms
5,376 KB
testcase_19 AC 322 ms
7,424 KB
testcase_20 AC 706 ms
10,880 KB
testcase_21 AC 692 ms
10,880 KB
testcase_22 AC 685 ms
11,008 KB
testcase_23 AC 703 ms
10,880 KB
testcase_24 AC 686 ms
10,880 KB
testcase_25 AC 669 ms
10,880 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'std::pair<long long int, long long int> math::crt2(const std::pair<long long int, long long int>&, const std::pair<long long int, long long int>&)':
main.cpp:336:10: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
  336 |     auto [r1, m1] = rm1;
      |          ^
main.cpp:337:10: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
  337 |     auto [r2, m2] = rm2;
      |          ^
main.cpp: In function 'long long int math::garner(const std::vector<std::pair<long long int, long long int> >&)':
main.cpp:359:14: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
  359 |         auto [r, m] = rm[i];
      |              ^
main.cpp: In function 'long long int math::garner_mod(std::vector<std::pair<long long int, long long int> >, long long int)':
main.cpp:380:14: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
  380 |         auto [r, m] = rm[i];
      |              ^
main.cpp:387:1: warning: no return statement in function returning non-void [-Wreturn-type]
  387 | }
      | ^
main.cpp: In function 'int main()':
main.cpp:454:18: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
  454 |         for(auto [x, c] : p) {
      |                  ^
main.cpp:457:18: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
  457 |         for(auto [x, c] : q) {
      |                  ^
main.cpp:464:14: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
  464 |     for(auto [x, c] : mp) {
      |              ^

ソースコード

diff #

#pragma region Macros
#include <bits/stdc++.h>
#define rep(i, n) for(int(i) = 0; (i) < (n); (i)++)
#define rrep(i, n) for(int(i) = (n)-1; (i) >= 0; (i)--)
#define FOR(i, m, n) for(int(i) = (m); (i) < (n); (i)++)
#define ROF(i, m, n) for(int(i) = (n)-1; (i) >= (m); (i)--)
#define ALL(v) (v).begin(), (v).end()
#define LLA(v) (v).rbegin(), (v).rend()
#define SZ(v) (int)(v).size()
#define INT(...)     \
    int __VA_ARGS__; \
    read(__VA_ARGS__)
#define LL(...)     \
    ll __VA_ARGS__; \
    read(__VA_ARGS__)
#define DOUBLE(...)     \
    double __VA_ARGS__; \
    read(__VA_ARGS__)
#define CHAR(...)     \
    char __VA_ARGS__; \
    read(__VA_ARGS__)
#define STRING(...)     \
    string __VA_ARGS__; \
    read(__VA_ARGS__)
#define VEC(type, name, size) \
    vector<type> name(size);  \
    read(name)
#define VEC2(type, name, height, width)                     \
    vector<vector<type>> name(height, vector<type>(width)); \
    read(name)
#define DVEC(type, name1, name2, size)     \
    vector<type> name1(size), name2(size); \
    read(name1, name2);
#define TVEC(type, name1, name2, name3, size)           \
    vector<type> name1(size), name2(size), name3(size); \
    read(name1, name2, name3);
using namespace std;
using ll = long long;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
const int INF = 1 << 30;
const ll LINF = 1LL << 60;
const int MOD = 1e9 + 7;
const char newl = '\n';
const int dx[] = {1, 0, -1, 0};
const int dy[] = {0, 1, 0, -1};
template <class T> inline bool between(T x, T l, T r) { return l <= x && x < r; }
template <class T> inline vector<T> make_vec(size_t a, T val) { return vector<T>(a, val); }
template <class... Ts> inline auto make_vec(size_t a, Ts... ts) { return vector<decltype(make_vec(ts...))>(a, make_vec(ts...)); }
void read() {}
template <class T> inline void read(T &a) { cin >> a; }
template <class T, class S> inline void read(pair<T, S> &p) { read(p.first), read(p.second); }
template <class T> inline void read(vector<T> &v) {
    for(auto &&a : v) read(a);
}
template <class T, class U> inline void read(vector<T> &a, vector<U> &b) {
    for(int i = 0; i < a.size(); i++) {
        read(a[i]);
        read(b[i]);
    }
}
template <class T, class U, class V> inline void read(vector<T> &a, vector<U> &b, vector<V> &c) {
    for(int i = 0; i < a.size(); i++) {
        read(a[i]);
        read(b[i]);
        read(c[i]);
    }
}
template <class Head, class... Tail> inline void read(Head &head, Tail &...tail) { read(head), read(tail...); }
template <class T> void write(const T &a) { cout << a << '\n'; }
template <class T> void write(const vector<T> &a) {
    for(int i = 0; i < a.size(); i++) cout << a[i] << (i + 1 == a.size() ? '\n' : ' ');
}
template <class Head, class... Tail> void write(const Head &head, const Tail &...tail) {
    cout << head << ' ';
    write(tail...);
}
template <class T> void writel(const T &a) { cout << a << '\n'; }
template <class T> void writel(const vector<T> &a) {
    for(int i = 0; i < a.size(); i++) cout << a[i] << '\n';
}
template <class Head, class... Tail> void writel(const Head &head, const Tail &...tail) {
    cout << head << '\n';
    write(tail...);
}
template <class T> auto sum(const vector<T> &a) { return accumulate(ALL(a), T(0)); }
template <class T> auto min(const vector<T> &a) { return *min_element(ALL(a)); }
template <class T> auto max(const vector<T> &a) { return *max_element(ALL(a)); }
template <class T, class U> void msort(vector<T> &a, vector<U> &b) {
    assert(a.size() == b.size());
    vector<pair<T, U>> ab(a.size());
    for(int i = 0; i < a.size(); i++) ab[i] = {a[i], b[i]};
    sort(ALL(ab));
    for(int i = 0; i < a.size(); i++) {
        a[i] = ab[i].first;
        b[i] = ab[i].second;
    }
}
template <class T, class U, class V> void msort(vector<T> &a, vector<U> &b, vector<V> &c) {
    assert(a.size() == b.size() && b.size() == c.size());
    vector<tuple<T, U, V>> abc(a.size());
    for(int i = 0; i < a.size(); i++) abc[i] = {a[i], b[i], c[i]};
    sort(ALL(abc));
    for(int i = 0; i < a.size(); i++) {
        a[i] = get<0>(abc[i]);
        b[i] = get<1>(abc[i]);
        c[i] = get<2>(abc[i]);
    }
}
template <class T, class U> inline bool chmax(T &a, U b) {
    if(a < b) {
        a = b;
        return 1;
    }
    return 0;
}
template <class T, class U> inline bool chmin(T &a, U b) {
    if(a > b) {
        a = b;
        return 1;
    }
    return 0;
}
inline int bsf(int v) { return __builtin_ctz(v); } // 最下位の1が下から何番目か
inline int bsf(ll v) { return __builtin_ctzll(v); }
inline int bsr(int v) { return 31 - __builtin_clz(v); } // 最上位の1が下から何番目か
inline int bsr(ll v) { return 63 - __builtin_clzll(v); }
inline int lsb(int v) { return v & -v; } // 最上位の1だけ残す
inline ll lsb(ll v) { return v & -v; }
inline int msb(int v) { return 1 << bsr(v); } // 最上位の1だけ残す
inline ll msb(ll v) { return 1LL << bsr(v); }
struct IO {
    IO() {
        ios::sync_with_stdio(false);
        cin.tie(nullptr);
        cout << fixed << setprecision(10);
    }
} io;
#pragma endregion

/* ModInt */
template <uint_fast64_t MOD> class ModInt {
    using u64 = uint_fast64_t;

  public:
    u64 val;

    ModInt(const u64 x = 0) : val((x + MOD) % MOD) {}
    constexpr u64 &value() { return val; }
    static const u64 get_mod() { return MOD; }
    constexpr ModInt operator-() { return val ? MOD - val : 0; }
    constexpr ModInt operator+(const ModInt &rhs) const {
        return ModInt(*this) += rhs;
    }
    constexpr ModInt operator-(const ModInt &rhs) const {
        return ModInt(*this) -= rhs;
    }
    constexpr ModInt operator*(const ModInt &rhs) const {
        return ModInt(*this) *= rhs;
    }
    constexpr ModInt operator/(const ModInt &rhs) const {
        return ModInt(*this) /= rhs;
    }
    constexpr ModInt &operator+=(const ModInt &rhs) {
        val += rhs.val;
        if(val >= MOD) {
            val -= MOD;
        }
        return *this;
    }
    constexpr ModInt &operator-=(const ModInt &rhs) {
        if(val < rhs.val) {
            val += MOD;
        }
        val -= rhs.val;
        return *this;
    }
    constexpr ModInt &operator*=(const ModInt &rhs) {
        val = val * rhs.val % MOD;
        return *this;
    }

    constexpr ModInt &operator/=(const ModInt &rhs) {
        *this *= rhs.inv();
        return *this;
    }

    constexpr bool operator==(const ModInt &rhs) {
        return this->val == rhs.val;
    }
    constexpr bool operator!=(const ModInt &rhs) {
        return this->val != rhs.val;
    }
    friend constexpr ostream &operator<<(ostream &os, const ModInt<MOD> &x) {
        return os << x.val;
    }
    friend constexpr istream &operator>>(istream &is, ModInt<MOD> &x) {
        return is >> x.val;
    }

    constexpr ModInt inv() const { return ModInt(*this).pow(MOD - 2); }

    constexpr ModInt pow(long long e) const {
        u64 x = 1, p = val;
        while(e > 0) {
            if(e % 2 == 0) {
                p = (p * p) % MOD;
                e /= 2;
            } else {
                x = (x * p) % MOD;
                e--;
            }
        }
        return ModInt(x);
    }
};

// using mint = ModInt<1000000007UL>;
using mint = ModInt<998244353UL>;

namespace math {

std::set<long long> divisor(long long n) {
    std::set<long long> ret;
    for(long long i = 1; i * i <= n; i++) {
        if(n % i == 0) {
            ret.insert(i);
            if(i * i != n)
                ret.insert(n / i);
        }
    }
    return ret;
}

// 素数篩+素因数分解
// 初期化O(N),素因数分解O(logN)
// たくさん素因数分解するときはこっち
struct Sieve {
    int N;
    std::vector<int> sieve;
    Sieve(int n) : N(n + 1), sieve(n + 1) { init(); }
    void init() {
        std::iota(sieve.begin(), sieve.end(), 0);
        for(int i = 2; i * i <= N; i++) {
            if(sieve[i] < i)
                continue;
            for(int j = i * i; j <= N; j += i) {
                if(sieve[j] == j)
                    sieve[j] = i;
            }
        }
    }

    bool is_prime(int x) {
        assert(x <= N);
        return sieve[x] == x;
    }

    std::map<long long, int> prime_factorize(long long n) {
        assert(n <= N);
        std::map<long long, int> ret;
        while(n > 1) {
            ret[sieve[n]]++;
            n = n / sieve[n];
        }
        return ret;
    }
};

// 素因数分解
// O(sqrt(N))
std::map<long long, int> prime_factor(long long n) {
    std::map<long long, int> ret;
    for(long long i = 2; i * i <= n; i++) {
        while(n % i == 0) {
            ret[i]++;
            n /= i;
        }
    }
    if(n != 1)
        ret[n] = 1;
    return ret;
}

long long mod_pow(long long x, long long n, long long mod) {
    if(n == 0)
        return 1;
    long long res = mod_pow(x * x % mod, n / 2, mod);
    if(n & 1)
        res = res * x % mod;
    return res;
}

long long euler_phi(long long n) {
    long long ret = n;
    for(long long i = 2; i * i <= n; i++) {
        if(n % i == 0) {
            ret -= ret / i;
            while(n % i == 0)
                n /= i;
        }
    }
    if(n > 1)
        ret -= ret / n;
    return ret;
}

long long extgcd(long long a, long long b, long long &x, long long &y) {
    if(b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    long long d = extgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

long long mod_inv(long long a, long long mod) {
    long long b = mod, u = 1, v = 0;
    while(b) {
        long long t = a / b;
        std::swap(a -= t * b, b);
        std::swap(u -= t * b, v);
    }
    u %= mod;
    if(u < 0)
        u += mod;
    return u;
}

std::pair<long long, long long>
crt2(const std::pair<long long, long long> &rm1,
     const std::pair<long long, long long> &rm2) {
    long long p, q;
    auto [r1, m1] = rm1;
    auto [r2, m2] = rm2;
    long long d = extgcd(r1, m2, p, q);
    if((r2 - r1) % d != 0)
        return {0, -1};
    long long m = m1 * (m2 / d);
    long long tmp = (r2 - r1) / d * p % (m2 / d);
    long long c = (r1 + m1 * tmp);
    long long r = (c % m + m) % m;
    return {r, m};
}

/**
 * @brief  Garnerのアルゴリズム. x % m[i] == b[i] % m[i]を満たす最小のxを求める.
 *
 * @param rm {あまり,法}の配列. !!法はすべて互いに素!!
 * @return long long x % m[i] == b[i] % m[i]を満たす最小のx.
 */
long long garner(const std::vector<std::pair<long long, long long>> &rm) {
    int n = rm.size();
    long long m_prod = 1;
    long long res = rm[0].first % rm[0].second;
    for(int i = 1; i < n; i++) {
        auto [r, m] = rm[i];
        m_prod *= rm[i - 1].second;
        long long t = (((r - res) * mod_inv(m_prod, m)) % m + m) % m;
        res += t * m_prod;
    }
    return res;
}

/**
 * @brief  Garnerのアルゴリズム. x % m[i] == b[i] % m[i]を満たす最小のxを求める.
 *
 * @param rm {あまり,法}の配列. !!法はすべて互いに素!!
 * @param mod 答えの法.
 * @return long long x % m[i] == b[i] % m[i]を満たす最小のx (mod M).
 */
long long garner_mod(std::vector<std::pair<long long, long long>> rm,
                     long long mod) {
    rm.emplace_back(0, mod);
    int n = rm.size();
    std::vector<long long> coffs(n, 1), constants(n, 0);
    for(int i = 0; i < n - 1; i++) {
        auto [r, m] = rm[i];
        long long t = ((r - constants[i]) * mod_inv(coffs[i], m) % m + m) & m;
        for(int j = i + 1; j < n; j++) {
            (constants[j] += coffs[j] * t) %= rm[j].second;
            (coffs[j] *= m) %= rm[j].second;
        }
    }
}

long long floor_div(long long a, long long b) {
    assert(b != 0);
    if(b < 0)
        a = -a, b = -b;
    return a >= 0 ? a / b : (a - (b - 1)) / b;
}

long long ceil_div(long long a, long long b) {
    assert(b != 0);
    if(b < 0)
        a = -a, b = -b;
    return a >= 0 ? (a + b - 1) / b : a / b;
}
// left_most_bit
long long lmb(long long x) {
    long long a = 1;
    while(a <= x) {
        a *= 2;
    }
    return a / 2;
}

int bit_length(long long x) {
    assert(x >= 0);
    int ret = 0;
    while(x > 0) {
        ret++;
        x /= 2;
    }
    return ret;
}

long long ll_sqrt(long long x) {
    long long ret = sqrt(x) - 1;
    while((ret + 1) * (ret + 1) <= x)
        ret++;
    return ret;
}

// multiple of d in [left,right]
// multiple of d in [left,right],left>0
long long cmul(long long left, long long right, long long d) {
    return right / d - (left - 1) / d;
}

// multiple of d in [left,right]
long long count_multiple(long long left, long long right, long long d) {
    if(right < 0)
        return cmul(-right, -left, d);
    if(left > 0)
        return cmul(left, right, d);
    return cmul(0, right, d) + cmul(0, -left, d) - 1;
}

}; // namespace math

using namespace math;

int main() {
    INT(n);
    Sieve sieve(n + 1);
    map<int, int> mp;
    for(int i = 1; i <= n / 2; i++) {
        auto p = sieve.prime_factorize(i);
        auto q = sieve.prime_factorize(n - i);
        for(auto [x, c] : p) {
            q[x] += c;
        }
        for(auto [x, c] : q) {
            if(!mp.count(x) || c > mp[x]) {
                mp[x] = c;
            }
        }
    }
    mint res = 1;
    for(auto [x, c] : mp) {
        res *= mint(x).pow(c);
    }
    write(res);
}
0