結果

問題 No.473 和と積の和
ユーザー 🍡yurahuna
提出日時 2016-12-19 22:50:05
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 405 ms / 3,000 ms
コード長 2,278 bytes
コンパイル時間 1,852 ms
コンパイル使用メモリ 181,536 KB
実行使用メモリ 19,328 KB
最終ジャッジ日時 2024-12-16 03:41:17
合計ジャッジ時間 5,058 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 43
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
#define int long long   // <-----!!!!!!!!!!!!!!!!!!!

#define rep(i,n) for (int i=0;i<(n);i++)
#define rep2(i,a,b) for (int i=(a);i<(b);i++)
#define rrep(i,n) for (int i=(n)-1;i>=0;i--)
#define rrep2(i,a,b) for (int i=(a)-1;i>=b;i--)
#define chmin(a,b) (a)=min((a),(b));
#define chmax(a,b) (a)=max((a),(b));
#define all(a) (a).begin(),(a).end()
#define rall(a) (a).rbegin(),(a).rend()
#define printV(v) cout<<(#v)<<":";for(auto(x):(v)){cout<<" "<<(x);}cout<<endl;
#define printVS(vs) cout<<(#vs)<<":"<<endl;for(auto(s):(vs)){cout<<(s)<< endl;}
#define printVV(vv) cout<<(#vv)<<":"<<endl;for(auto(v):(vv)){for(auto(x):(v)){cout<<" "<<(x);}cout<<endl;}
#define printP(p) cout<<(#p)<<(p).first<<" "<<(p).second<<endl;
#define printVP(vp) cout<<(#vp)<<":"<<endl;for(auto(p):(vp)){cout<<(p).first<<" "<<(p).second<<endl;}

inline void output(){ cout << endl; }
template<typename First, typename... Rest>
inline void output(const First& first, const Rest&... rest) {
    cout << first << " "; output(rest...);
}

using ll = long long;
using Pii = pair<int, int>;
using vi = vector<int>;
using vvi = vector<vi>;
using vvvi = vector<vvi>;
const int inf = 1ll << 60;
const int mod = 1e9 + 7;
using Graph = vector<vector<int>>;

using State = tuple<int, int, int>;
map<State, int> dp;
vi divisors;

// nの2以上sqrt(n)以下の約数を列挙
void setDivisor(int n) {
    for (int i = 2; i * i <= n; i++) {
        if (n % i == 0) {
            divisors.emplace_back(i);
            // if (i != n / i) res.emplace_back(n / i);
        }
    }
    sort(all(divisors));
}

// y以上の整数n個からなる集合で、すべての要素の積がxとなるものの個数
int dfs(int n, int x, int y) {
    // output(n, x, y);
    if (x < y) return 0;
    if (n == 1) return 1;
    State s(n, x, y);
    if (dp.count(s)) return dp[s];
    int &ret = dp[s];
    for (int i = 0; i < divisors.size(); ++i) {
        int p = divisors[i];
        if (p * p > x) break;
        if (p < y || x % p != 0) continue;
        ret += dfs(n - 1, x / p, p);
    }
    return ret;
}

signed main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(0);

    int n, x;
    cin >> n >> x;
    setDivisor(x + 1);
    cout << dfs(n, x + 1, 2) << endl;
}
0