結果

問題 No.2798 Multiple Chain
コンテスト
ユーザー bortik
提出日時 2024-06-28 22:59:35
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.89.0)
結果
TLE  
実行時間 -
コード長 3,173 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 12,453 ms
コンパイル使用メモリ 290,488 KB
最終ジャッジ日時 2025-02-22 01:19:51
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 TLE * 1 -- * 1
other -- * 51
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

// https://codeforces.com/blog/entry/96344
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt")

#include <bits/stdc++.h>

using namespace std;

using ll = long long;
using u32 = unsigned int;
using u64 = unsigned long long;
using i128 = __int128;
using u128 = unsigned __int128;
using f128 = __float128;

template<typename T>
bool chmax(T& a, const T& b) {
    bool res = a < b;
    a = max(a, b);
    return res;
}
template<typename T>
bool chmin(T& a, const T& b){
    bool res = a > b;
    a = min(a, b);
    return res;
}

typedef vector<long long> vl;
typedef pair<ll,ll> pll;
typedef vector<pair<ll, ll>> vll;
typedef vector<int> vi;
typedef vector<pair<int,int>> vii;
typedef pair<int,int> pii;

const int inf = 1000000009;
const ll linf = 4000000000000000009;


// https://trap.jp/post/1224/
template<class... T>
void input(T&... a){
    (cin >> ... >> a);
}

void print(){
    cout << '\n';
}
template<class T, class... Ts>
void print(const T& a, const Ts&... b){
    cout << a;
    (cout << ... << (cout << ' ', b));
    cout << '\n';
}

#define rep1(a)          for(int i = 0; i < a; i++)
#define rep2(i, a)       for(int i = 0; i < a; i++)
#define rep3(i, a, b)    for(int i = a; i < b; i++)
#define rep4(i, a, b, c) for(int i = a; i < b; i += c)
#define overload4(a, b, c, d, e, ...) e
#define rep(...) overload4(__VA_ARGS__, rep4, rep3, rep2, rep1)(__VA_ARGS__)


#define eb emplace_back
#define mp make_pair
#define mt make_tuple
#define fi first
#define se second
//---------------------------------


vector<ll> part(ll n, ll m){
    vector dp(m+1, vector<ll>(n+1, 0));
    dp[0][0] = 1;
    rep(i, 1, m+1){
        rep(j, n+1){
            if(j-i >= 0){
                dp[i][j] = dp[i-1][j] + dp[i][j-i];
            } else {
                dp[i][j] = dp[i-1][j];
            }
        }
    }

    vector<ll> res(m);
    rep(i, m){
        res[i] = dp[i][n];
    }

    return res;
}

template <typename T>
std::vector<std::pair<T, T>> prime_factor(T n) {
    std::vector<std::pair<T, T>> ret;
    for (T i = 2; i * i <= n; i++) {
        if (n % i != 0) continue;
        T tmp = 0;
        while (n % i == 0) {
            tmp++;
            n /= i;
        }
        ret.push_back(std::make_pair(i, tmp));
    }
    if (n != 1) ret.push_back(std::make_pair(n, 1));
    return ret;
}

void solve(){
    ll n; cin >> n;
    auto facts = prime_factor(n);
    vector<vector<ll>> parts;
    ll mx = 0;
    //cout << "HELLO" << endl;
    for(auto [p, e]: facts){
        if(e==1) continue;
        parts.push_back(part(e, 10));
        //print(p, e);
        chmax(mx, e);
    }
    ll ans = 1;
    int s = parts.size();
    rep(L, 2, mx+1){
        rep(j, s){
            ll cur = parts[j][L]-parts[j][L-1];
            if(cur < 1) continue;
            rep(i, j){
                cur *= parts[i][L-1];
            }
            rep(i, j+1, s){
                cur *= parts[i][L];
            }
            ans += cur;
        }
        //print(parts[0][L]-parts[0][L-1]);
    }
    print(ans);
}


int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    int t = 1;
    //cin >> t;
    rep(i,0,t) solve();
}
0