結果

問題 No.1529 Constant Lcm
ユーザー jutamajutama
提出日時 2021-06-04 20:32:17
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 667 ms / 3,000 ms
コード長 6,244 bytes
コンパイル時間 3,090 ms
コンパイル使用メモリ 195,312 KB
実行使用メモリ 10,368 KB
最終ジャッジ日時 2024-04-30 00:00:28
合計ジャッジ時間 10,568 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 1 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 1 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 592 ms
9,728 KB
testcase_11 AC 216 ms
6,016 KB
testcase_12 AC 11 ms
5,376 KB
testcase_13 AC 511 ms
8,960 KB
testcase_14 AC 301 ms
7,040 KB
testcase_15 AC 339 ms
7,296 KB
testcase_16 AC 256 ms
6,400 KB
testcase_17 AC 134 ms
5,376 KB
testcase_18 AC 147 ms
5,376 KB
testcase_19 AC 315 ms
7,040 KB
testcase_20 AC 653 ms
10,368 KB
testcase_21 AC 656 ms
10,240 KB
testcase_22 AC 664 ms
10,240 KB
testcase_23 AC 667 ms
10,368 KB
testcase_24 AC 655 ms
10,240 KB
testcase_25 AC 641 ms
10,240 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <atcoder/all>
#include <bitset>
#include <fstream>
#include <functional>
#include <iostream>
#include <iomanip>
#include <limits>
#include <map>
#include <math.h>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <stdio.h>
#include <stdlib.h>
#include <unordered_map>
#include <unordered_set>
#include <vector>
using namespace std;
using namespace atcoder;

typedef long long ll;
typedef long double ld;
typedef std::pair<int, int> pii;
typedef std::pair<ll, ll> pll;
typedef std::pair<ld, ld> pid;
typedef std::pair<ll, std::string> pls;
typedef std::vector<bool> vb;
typedef std::vector<vb> vvb;
typedef std::vector<int> vi;
typedef std::vector<vi> vvi;
typedef std::vector<vvi> vvvi;
typedef std::vector<vvvi> vvvvi;
typedef std::vector<ll> vl;
typedef std::vector<vl> vvl;
typedef std::vector<vvl> vvvl;
typedef std::vector<vvvl> vvvvl;
typedef std::vector<ld> vd;
typedef std::vector<vd> vvd;
typedef std::vector<std::string> vs;
#define rep(i,n) for(auto i=0; i<n; ++i)
#define repm(i,s,n) for(auto i=s; i<n; ++i)
#define repd(i,n) for(auto i=n-1; i>=0; --i)
#define repdm(i,e,n) for(auto i=n-1; i>=e; --i)
#define all(a) (a).begin(), (a).end()
#define rall(a) (a).rbegin(), (a).rend()

template <class T> inline bool chmax(T& a, T b, int eq = 0) { if (a < b || (a == b && eq)) { a = b; return 1; } return 0; }
template <class T> inline bool chmin(T& a, T b, int eq = 0) { if (a > b || (a == b && eq)) { a = b; return 1; } return 0; }

template <class mint, internal::is_modint_t<mint>* = nullptr> constexpr istream& operator>>(istream& is, mint& x) noexcept {long long v = 0; std::cin >> v; x = v; return is;}
template <class mint, internal::is_modint_t<mint>* = nullptr> constexpr ostream& operator<<(ostream& os, const mint& x) noexcept {os << x.val(); return os;}
inline void _n() { std::cout << std::endl; }
template <class T> inline void _(const T a) { std::cout << a; }
template <class T> inline void _l(const T a) { _(a); _n(); }
template <class T> inline void _s(const T a) { _(a); _(' '); }
template <class T> inline void _v(const std::vector<T> v) { for(auto a : v) _(a); }
template <class T> inline void _vl(const std::vector<T> v) { for(auto a : v) _l(a); }
template <class T> inline void _vs(const std::vector<T> v) { for(auto a : v) _s(a); _n(); }
template <class T> inline void _vvl(const std::vector<std::vector<T>> v) { for(auto a : v) { _v(a); _n(); } }
template <class T> inline void _vvs(const std::vector<std::vector<T>> v) { for(auto a : v) { _vs(a); } }
inline void ynl(const bool b) {_l(b ? "yes" : "no");}
inline void yns(const bool b) {_l(b ? "Yes" : "No");}
inline void ynu(const bool b) {_l(b ? "YES" : "NO");}

constexpr int INF = numeric_limits<int>::max() >> 1;
constexpr long long INF_LL = numeric_limits<long long>::max() >> 1LL;
const long long MOD1 = 1000000007;
const long long MOD9 = 998244353;
using mint1 = atcoder::modint1000000007;
using mint9 = atcoder::modint998244353;

//* ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ *//


void input() {
    
}

struct sieve_E {

    public :
    sieve_E(long long n) : _n(n) {
        size = (_n/spp + 1) * srl;
        sieve = vector<bool>(size, true);
        sieving_prime = vector<unsigned long long>(size, 0);

        sieving();
    }

    bool is_prime(long long n) {
        assert(0 <= n && n <= _n);
        if (n == 0 || n == 1) return false;
        for (auto p : sp) {if(n == p) return true; if(n%p == 0) return false;}
        return sieve[ntoi(n)];
    }
    unsigned long long get_sieving_prime(long long n) {
        assert(0 <= n && n <= _n);
        if (n == 0 || n == 1) return 0;
        for (auto p : sp) {if(n%p == 0) return p;}
        if (is_prime(n)) return n;
        return sieving_prime[ntoi(n)];
    }
    unsigned long long get_sieve_size() {return _n;}

    private :
    long long _n;

    unsigned long long size;
    std::vector<bool> sieve;
    std::vector<unsigned long long> sieving_prime;

    const std::vector<unsigned int> sp = {2, 3, 5};                      // skip primes
    const std::vector<unsigned int> sr = {1, 7, 11, 13, 17, 19, 23, 29}; // stand reminders
    const unsigned int spp = 2 * 3 * 5;                                  // skip primes product
    const unsigned int srl = 8;                                          // stand reminders length

    void sieving() {
        unsigned long long max_i = ((unsigned long long)sqrt(_n)/spp + 1) * srl;
        for (unsigned long long i = 1; i <= max_i ; ++i) {
            if(!sieve[i]) continue;

            unsigned long long num_i = (unsigned long long)iton(i);
            for (unsigned long long j = num_i * num_i; j < (unsigned long long)_n; j += num_i) {
                long long idx_j = ntoi(j);
                if(idx_j > 0) {
                    sieve[idx_j] = false;
                    sieving_prime[idx_j] = num_i;
                }
            }
        }
    }

    unsigned long long iton(unsigned long long idx) {
        return (unsigned long long) ((idx/srl) * spp + sr[idx%srl]);
    }
    long long ntoi(unsigned long long num) {
        unsigned int sri = (unsigned int)(std::find(sr.begin(), sr.end(), num%spp) - sr.begin());
        return sri == srl ? -1 : (num/spp) * srl + sri;
    }
};



void solve() {
    
    
    
    int N; cin >> N;
    
    sieve_E sieve(N);
    
    map<ll, ll> mp;
    repm(i, 1, N) {
        ll tmp_s = i;
        ll tmp_b = N-i;
        if(tmp_s > tmp_b) break;
        
        map<ll, ll> tmp_mp;
        while(tmp_s > 1) {
            ll sp = sieve.get_sieving_prime(tmp_s);
            tmp_mp[sp]++;
            tmp_s /= sp;
        }
        while(tmp_b > 1) {
            ll sp = sieve.get_sieving_prime(tmp_b);
            tmp_mp[sp]++;
            tmp_b /= sp;
        }
        
        for(auto [k, v] : tmp_mp) {
            chmax(mp[k], v);
        }
    }
    
    mint9 ans = 1;
    for(auto [k, v] : mp) {
        ans *= pow_mod(k, v, MOD9);
    }
    _l(ans);
    
    
    
    
    
}

//* ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ *//

int main() {
    
    std::ifstream in("input.txt");
    std::cin.rdbuf(in.rdbuf());
    std::cin.tie(0);
    std::ios::sync_with_stdio(false);
    
    input();
    solve();
    
    return 0;
}
0