Advertisement
bibaboba12345

д1

Aug 24th, 2021
113
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.09 KB | None | 0 0
  1. #pragma GCC target ("avx2")
  2. #pragma GCC optimization ("O3")
  3. #pragma GCC optimization ("unroll-loops")
  4. #include <iostream>
  5. #include <cstdio>
  6. #include <string>
  7. #include <algorithm>
  8. #include <cmath>
  9. #include <vector>
  10. #include <queue>
  11. #include <stack>
  12. #include <deque>
  13. #include <set>
  14. #include <map>
  15. #include <climits>
  16. #include <cstdlib>
  17. #include<time.h>
  18. #include<iomanip>
  19. #include <ctime>
  20. using namespace std;
  21. const int N = 2e5 + 7;
  22. long long dp[N], mod, dn_s, n, suffixsum[N], pows[30];
  23. int main() {
  24. std::ios::sync_with_stdio(false);
  25. cin.tie(0);
  26. cout.tie(0);
  27. cin >> n >> mod;
  28. dp[n] = 1;
  29. dn_s = 0;
  30. pows[0] = 1;
  31. for (int i = 1; i < 30; i++) {
  32. pows[i] = pows[i - 1] * 2;
  33. pows[i] %= mod;
  34. }
  35. for (int i = n; i >= 1; i--) {
  36. dp[i] += dn_s;
  37. dp[i] %= mod;
  38. dn_s += dp[i];
  39. dn_s %= mod;
  40. for (int j = 2; j <= i; j++) {
  41. for (int pw = 29; pw >= 0; pw--) {
  42. if (i / j == (i / (j + pows[pw]))) {
  43. dp[i / j] += pows[pw] * dp[i];
  44. dp[i / j] %= mod;
  45. j += pows[pw];
  46. }
  47. }
  48. dp[i / j] += dp[i];
  49. dp[i / j] %= mod;
  50. }
  51. }
  52. cout << dp[1];
  53.  
  54. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement