Advertisement
GerONSo

F спбгу need upd

Jan 10th, 2019
121
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.97 KB | None | 0 0
  1. /*
  2. ┓┏┓┏┓┃
  3. ┛┗┛┗┛┃
  4. ┓┏┓┏┓┃
  5. ┛┗┛┗┛┃
  6. ┓┏┓┏┓┃\○/
  7. ┛┗┛┗┛┃ / /
  8. ┓┏┓┏┓┃ノ
  9. ┛┗┛┗┛┃
  10. ┓┏┓┏┓┃
  11. ┛┗┛┗┛┃
  12. ┓┏┓┏┓┃
  13. ┛┗┛┗┛┃
  14. ┓┏┓┏┓┃
  15. ┛┗┛┗┛┃
  16. ┓┏┓┏┓┃┓
  17. ┛┗┛┗┛┃┃
  18. MIPTCLASSICMIPTCLASSICMIPTCLASSICMIPTCLASSICMIPTCLASSICMIPTCLASSIC
  19. */
  20.  
  21. // #define pragma
  22.  
  23. #ifdef pragma
  24. #pragma GCC optimize("Ofast")
  25. #pragma GCC optimize("no-stack-protector")
  26. #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
  27. #pragma GCC optimize("unroll-loops")
  28. #pragma GCC diagnostic ignored "-Wunused-result"
  29. #endif // pragma
  30.  
  31. #include<bits/stdc++.h>
  32. #include <ext/pb_ds/assoc_container.hpp>
  33. #include <ext/pb_ds/tree_policy.hpp>
  34.  
  35. #define ll long long
  36. #define all(x) begin(x), end(x)
  37. #define pb push_back
  38. #define x first
  39. #define y second
  40. #define int long long
  41. #define zero(two) memset(two, 0, sizeof(two))
  42.  
  43. using namespace std;
  44. using namespace __gnu_pbds;
  45.  
  46.  
  47. typedef vector<int> vi;
  48. typedef vector<bool> vb;
  49. typedef pair<int, int> pii;
  50. typedef long double ld;
  51. typedef vector<vi> matrix;
  52. template<typename T>
  53. using kawaii_set = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
  54.  
  55. const ld PI = atan2(0, -1);
  56.  
  57. void seriy() {
  58. ios::sync_with_stdio(0);
  59. cin.tie(0);
  60. cout.tie(0);
  61. cout << fixed << setprecision(10);
  62. #if 1
  63. freopen("input", "r", stdin);
  64. freopen("output", "w", stdout);
  65. #endif
  66. }
  67.  
  68. const int MAXN = 1e6 + 1;
  69. const int INF = 1e9 + 7;
  70.  
  71. vi pr(MAXN, 1);
  72. vi p;
  73. vector<deque<int>> vm;
  74. matrix prec(21);
  75.  
  76. int mul(deque<int> &st) {
  77. int ans = 1;
  78. for(auto i : st) {
  79. ans *= i;
  80. }
  81. return ans;
  82. }
  83.  
  84. int cnt(int x, int i) {
  85. deque<int> cp = vm[i];
  86. int ans = 0;
  87. int ml = mul(cp);
  88. while(ans < x) {
  89. ans++;
  90. // cerr << cp.front() << '\n';
  91. ml = ml / cp.front() * (cp.front() - 1);
  92. int nw = cp.front();
  93. cp.pop_front();
  94. cp.push_front(nw - 1);
  95. if(nw - 1 == 1) {
  96. cp.pop_front();
  97. }
  98. }
  99. return ml;
  100. }
  101.  
  102. signed main() {
  103. seriy();
  104. for(int i = 2; i < MAXN; i++) {
  105. if(pr[i]) {
  106. p.pb(i);
  107. if(i * 1ll * i < MAXN) {
  108. for(int j = i * i; j < MAXN; j += i) {
  109. pr[j] = 0;
  110. }
  111. }
  112. }
  113. }
  114. pr[0] = 0;
  115. pr[1] = 0;
  116. int n;
  117. cin >> n;
  118. int a;
  119. vm.resize(n);
  120. vi ap;
  121. map<int, char> c;
  122. for(int i = 0; i < n; i++) {
  123. cin >> a;
  124. int k = a;
  125. c.clear();
  126. if(pr[k]) {
  127. vm[i].push_back(2);
  128. continue;
  129. }
  130. for(int j = 0; j < (int)p.size(); j++) {
  131. if(p[j] > k) {
  132. break;
  133. }
  134. while(k % p[j] == 0) {
  135. if(c[p[j]] == 0) {
  136. c[p[j]] = 2;
  137. }
  138. else {
  139. c[p[j]]++;
  140. }
  141. k /= p[j];
  142. }
  143. }
  144. for(auto &[k, v] : c) {
  145. ap.pb(v);
  146. }
  147. sort(all(ap));
  148. for(auto j : ap) {
  149. if(j > 1) vm[i].push_back(j);
  150. }
  151. ap.clear();
  152. }
  153.  
  154. for(int i = 0; i < 21; i++) {
  155. t[i].resize(4 * MAXN);
  156. prec[i].resize(n);
  157. for(int j = 0; j < n; j++) {
  158. prec[i][j] = cnt(i, j);
  159. }
  160. build(i, 0, 0, n - 1);
  161. }
  162.  
  163. int q;
  164. cin >> q;
  165. vi l(q), r(q), x(q);
  166. for(int i = 0; i < q; i++) {
  167. cin >> l[i] >> r[i] >> x[i];
  168. l[i]--;
  169. r[i]--;
  170. x[i] = min(240, x[i]);
  171. int mn = 0;
  172. for(int j = 0; j < 21; j++) {
  173. int cur = get(j, 0, 0, n - 1, l[i], r[i]);
  174. if(cur <= x) {
  175. mn = cur;
  176. break;
  177. }
  178. }
  179.  
  180. }
  181. return 0;
  182. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement