結果
問題 |
No.1581 Multiple Sequence
|
ユーザー |
|
提出日時 | 2021-07-06 16:29:28 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 123 ms / 2,000 ms |
コード長 | 1,192 bytes |
コンパイル時間 | 1,786 ms |
コンパイル使用メモリ | 173,432 KB |
実行使用メモリ | 11,212 KB |
最終ジャッジ日時 | 2024-07-01 12:09:35 |
合計ジャッジ時間 | 4,017 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define REP(i,n) for (int i = 0; i < (n); ++i) #define DREP(i,s,n) for(int i = (s); i < (n); i++) template<class T> inline bool chmin(T& a, T b) {if (a > b) {a = b;return true;}return false;} template<class T> inline bool chmax(T& a, T b) {if (a < b) {a = b;return true;}return false;} using ll = long long; using P = pair<int,int>; using Pl = pair<long long,long long>; using veci = vector<int>; using vecl = vector<long long>; using vecveci = vector<vector<int>>; using vecvecl = vector<vector<long long>>; const int MOD = 1000000007; const double pi = acos(-1); ll gcd(ll a, ll b) {if(b == 0) return a; else return gcd(b,a%b);} ll lcm(ll a, ll b) {return a*b/gcd(a,b);} veci div(int N) { veci res; for(int i = 1; i*i <= N; i++) { if(N%i == 0) { res.push_back(i); if(i != N/i) res.push_back(N/i); } } sort(res.begin(),res.end()); return res; } int main() { int M; cin >> M; vecl dp(1000010); dp[0] = 1; for(int i = 1; i <= M; i++) { for(auto j : div(i)) { dp[i] += dp[i/j-1]; dp[i] %= MOD; } } cout << dp[M] << endl; }