結果
問題 | No.1973 Divisor Sequence |
ユーザー | k1suxu |
提出日時 | 2022-06-11 09:22:02 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,712 ms / 2,000 ms |
コード長 | 2,661 bytes |
コンパイル時間 | 3,285 ms |
コンパイル使用メモリ | 254,072 KB |
実行使用メモリ | 178,948 KB |
最終ジャッジ日時 | 2023-10-21 16:33:26 |
合計ジャッジ時間 | 8,532 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge10 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,348 KB |
testcase_01 | AC | 2 ms
4,348 KB |
testcase_02 | AC | 160 ms
29,924 KB |
testcase_03 | AC | 10 ms
5,400 KB |
testcase_04 | AC | 94 ms
18,848 KB |
testcase_05 | AC | 29 ms
16,528 KB |
testcase_06 | AC | 99 ms
25,380 KB |
testcase_07 | AC | 21 ms
6,908 KB |
testcase_08 | AC | 54 ms
14,748 KB |
testcase_09 | AC | 83 ms
16,804 KB |
testcase_10 | AC | 105 ms
20,224 KB |
testcase_11 | AC | 14 ms
5,352 KB |
testcase_12 | AC | 81 ms
20,948 KB |
testcase_13 | AC | 25 ms
8,296 KB |
testcase_14 | AC | 67 ms
17,720 KB |
testcase_15 | AC | 141 ms
21,872 KB |
testcase_16 | AC | 137 ms
25,908 KB |
testcase_17 | AC | 87 ms
17,332 KB |
testcase_18 | AC | 7 ms
4,392 KB |
testcase_19 | AC | 123 ms
22,004 KB |
testcase_20 | AC | 21 ms
6,116 KB |
testcase_21 | AC | 123 ms
30,384 KB |
testcase_22 | AC | 114 ms
22,016 KB |
testcase_23 | AC | 1,712 ms
178,948 KB |
testcase_24 | AC | 314 ms
46,156 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i,n) for(int i = 0; i < (int)n; i++) #define FOR(n) for(int i = 0; i < (int)n; i++) #define repi(i,a,b) for(int i = (int)a; i < (int)b; i++) #define pb push_back #define m0(x) memset(x,0,sizeof(x)) #define fill(x,y) memset(x,y,sizeof(x)) #define bg begin() #define ed end() #define all(x) x.bg,x.ed //#define mp make_pair #define vi vector<int> #define vvi vector<vi> #define vll vector<ll> #define vvll vector<vll> #define vs vector<string> #define vvs vector<vs> #define vc vector<char> #define vvc vector<vc> #define pii pair<int,int> #define pllll pair<ll,ll> #define vpii vector<pair<int,int>> #define vpllll vector<pair<ll,ll>> #define vpis vector<pair<int,string>> #define vplls vector<pair<ll, string>> #define vpsi vector<pair<string, int>> #define vpsll vector<pair<string, ll>> template<typename T> void chmax(T &a, const T &b) {a = (a > b? a : b);} template<typename T> void chmin(T &a, const T &b) {a = (a < b? a : b);} using ll = long long; using ld = long double; using ull = unsigned long long; const ll INF = numeric_limits<long long>::max() / 2; const ld pi = 3.1415926535897932384626433832795028; const ll mod = 1e9 + 7; int dx[] = {-1, 0, 1, 0, -1, -1, 1, 1}; int dy[] = {0, -1, 0, 1, -1, 1, -1, 1}; #define int long long vector<pair<long long, long long> > prime_factorize(long long N) { vector<pair<long long, long long> > res; for (long long a = 2; a * a <= N; ++a) { if (N % a != 0) continue; long long ex = 0; // 指数 // 割れる限り割り続ける while (N % a == 0) { ++ex; N /= a; } // その結果を push res.push_back({a, ex}); } // 最後に残った数について if (N != 1) res.push_back({N, 1}); return res; } void solve() { int n, m; cin >> n >> m; vpii p = prime_factorize(m); int ans = 1; for(auto e : p) { int sz = e.second; vvi dp(n, vi(sz+1, 0)); rep(j, sz+1) dp[0][j] = 1; rep(i, n-1) { rep(j, sz+1) { rep(k, sz+1-j) { dp[i + 1][j] += dp[i][k]; dp[i + 1][j] %= mod; } } // int s = 0; // rep(j, sz+1) { // s += dp[i][sz-j]; // s %= mod; // dp[i+1][j] = s; // } } int sum = 0; rep(j, sz+1) (sum += dp[n-1][j]) %= mod; (ans *= sum) %= mod; } cout << ans << endl; } signed main() { cin.tie(nullptr); ios::sync_with_stdio(false); solve(); return 0; }