結果
問題 | No.1973 Divisor Sequence |
ユーザー | HIcoder |
提出日時 | 2023-07-16 22:19:30 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,713 ms / 2,000 ms |
コード長 | 1,310 bytes |
コンパイル時間 | 1,083 ms |
コンパイル使用メモリ | 110,608 KB |
実行使用メモリ | 179,060 KB |
最終ジャッジ日時 | 2023-10-17 20:10:03 |
合計ジャッジ時間 | 6,064 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,348 KB |
testcase_01 | AC | 2 ms
4,348 KB |
testcase_02 | AC | 169 ms
30,024 KB |
testcase_03 | AC | 10 ms
5,500 KB |
testcase_04 | AC | 99 ms
18,948 KB |
testcase_05 | AC | 32 ms
16,640 KB |
testcase_06 | AC | 106 ms
25,480 KB |
testcase_07 | AC | 21 ms
7,008 KB |
testcase_08 | AC | 58 ms
14,860 KB |
testcase_09 | AC | 87 ms
16,904 KB |
testcase_10 | AC | 111 ms
20,336 KB |
testcase_11 | AC | 15 ms
5,452 KB |
testcase_12 | AC | 87 ms
20,900 KB |
testcase_13 | AC | 27 ms
8,396 KB |
testcase_14 | AC | 70 ms
17,832 KB |
testcase_15 | AC | 153 ms
21,972 KB |
testcase_16 | AC | 150 ms
26,008 KB |
testcase_17 | AC | 94 ms
17,432 KB |
testcase_18 | AC | 7 ms
4,432 KB |
testcase_19 | AC | 124 ms
22,032 KB |
testcase_20 | AC | 21 ms
6,216 KB |
testcase_21 | AC | 132 ms
30,496 KB |
testcase_22 | AC | 120 ms
22,116 KB |
testcase_23 | AC | 1,713 ms
179,060 KB |
testcase_24 | AC | 328 ms
46,268 KB |
ソースコード
#include<iostream> #include<set> #include<algorithm> #include<vector> #include<string> #include<set> #include<map> #include<numeric> #include<queue> #include<cmath> using namespace std; typedef long long ll; const ll INF=1LL<<60; typedef pair<ll,ll> P; typedef pair<int,P> PP; const ll MOD=1e9+7; vector<P> primediv(ll n){ vector<P> res; for(ll p=2;p*p<=n;p++){ if(n%p==0){ P tmp=make_pair(p,0); while(n%p==0){ tmp.second++; n/=p; } res.push_back(tmp); } } if(n>1){ res.emplace_back(n,1); n/=n; } return res; } int main(){ int N; ll M; cin>>N>>M; auto divM=primediv(M); ll ans=1; for(auto [v,num]:divM){ vector<vector<ll>> dp(N+1,vector<ll>(num+1,0)); for(int j=0;j<=num;j++){ dp[0][j]=1; } for(int i=1;i<N;i++){ for(int j=0;j<=num;j++){ for(int l=0;l<=num-j;l++){ dp[i][j]+=dp[i-1][l]; dp[i][j]%=MOD; } } } ll sum=0; for(int j=0;j<=num;j++){ sum+=dp[N-1][j]; sum%=MOD; } ans*=sum; ans%=MOD; } cout<<ans<<endl; }