結果

問題 No.1973 Divisor Sequence
ユーザー magsta
提出日時 2022-06-06 22:07:18
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,183 bytes
コンパイル時間 1,876 ms
コンパイル使用メモリ 176,256 KB
実行使用メモリ 178,940 KB
最終ジャッジ日時 2024-09-21 05:50:16
合計ジャッジ時間 6,385 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2 WA * 1
other WA * 22
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
using namespace std;
#define MOD 1000000007
vector<pair<long long, long long>> fact(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;
}
res.push_back({ a, ex });
}
if (N != 1) res.push_back({ N, 1 });
return res;
}
int main() {
long long N, M; cin >> N >> M;
auto factor = fact(M);
long long N_ = factor.size();
long long ans = 1;
for (int q = 0; q < N_; q++) {
vector<vector<long long>> dp(N + 1, vector<long long>(factor[q].second + 1, 0));
dp[0][0] = 1;
for (int i = 1; i <= N; i++) {
for (int j = factor[q].second; j >= 0; j--) {
for (int k = 0; k <= j; k++) {
dp[i][j] = (dp[i][j] + dp[i - 1][k]) % MOD;
}
}
}
long long sum = 0;
for (int i = 0; i <= factor[q].second; i++) {
sum = (sum + dp[N][i]) % MOD;
}
ans = ans * sum % MOD;
}
cout << ans << endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0