結果

問題 No.1529 Constant Lcm
ユーザー yakki
提出日時 2021-06-04 21:08:33
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,942 bytes
コンパイル時間 1,404 ms
コンパイル使用メモリ 130,812 KB
最終ジャッジ日時 2025-01-21 23:09:23
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 16 TLE * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<string>
#include<vector>
#include<algorithm>
#include<bitset>
#include<set>
#include<map>
#include<stack>
#include<queue>
#include<deque>
#include<list>
#include<iomanip>
#include<cmath>
#include<cstring>
#include<functional>
#include<cstdio>
#include<cstdlib>
#include<numeric>
//#include<atcoder/all>
using namespace std;
//using namespace atcoder;

#define repr(i, a, b) for (int i = (int)(a); i < (int)(b); i++)
#define rep(i, n) repr(i, 0, n)
#define INF 2e9
//#define MOD 1000000007
#define MOD 998244353
#define LINF (long long)4e18
#define jck 3.141592
#define PI acos(-1.0)

const double EPS = 1e-10;

using ll = long long;
using Pi = pair<int,int>;
using Pl = pair<ll,ll>;
//using mint = modint998244353;

int dh[] = {-1,0,1,0};
int dw[] = {0,1,0,-1};

vector<pair<ll,int>> factorize(ll n){
    vector<pair<ll,int>> div;
    ll d = n;
    for(ll i = 2; i*i <= n; i++){
        pair<ll,int> p;
        if(d % i == 0){
            p = make_pair(i,0);
            while(d % i == 0){
                d /= i;
                p.second++;
            }
            div.emplace_back(p);
        }
    }
    if(d != 1) div.emplace_back(d,1);
    return div;
}

ll modpow(ll a, ll n, ll m){
    if(n == 0) return 1;
    ll half = modpow(a,n/2,m);
    ll res = half * half % m;
    if(n & 1) res = res * (a%m) % m;
    return res;
}

int main(){
    int n; cin >> n;
    vector<int> cnt(n+1);
    ll ans = 1;
    for(int i = 1; i <= n-i; i++){
        auto div1 = factorize(i);
        auto div2 = factorize(n-i);
        map<int,int> cnt2;
        for(auto [u,v] : div1) cnt2[u] += v;
        for(auto [u,v] : div2) cnt2[u] += v;
        for(auto [u,v] : cnt2){
            if(cnt[u] < cnt2[u]){
                ans *= modpow(u,cnt2[u]-cnt[u],MOD);
                ans %= MOD;
                cnt[u] = cnt2[u];
            }
        }
        //cout << i << " " << ans << endl;
    }

    cout << ans << endl;
}
0