fork download
copy 
  1. /*
  2. * @Author: hungeazy
  3. * @Date:   2025-11-18 14:22:35
  4. * @Last Modified by:   hungeazy
  5. * @Last Modified time: 2025-11-18 16:39:00
  6. */
  7. #include <bits/stdc++.h>
  8. #include <ext/pb_ds/assoc_container.hpp> 
  9. #include <ext/pb_ds/tree_policy.hpp> 
  10. // #pragma GCC optimize("O3")  
  11. // #pragma GCC optimize("unroll-loops")  
  12. // #pragma GCC target("avx2,bmi,bmi2,popcnt,lzcnt")  
  13. using namespace std;
  14. using namespace __gnu_pbds; 
  15. bool M1;
  16. #define fast ios_base::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL);
  17. #define int long long
  18. #define ll long long 
  19. #define ull unsigned long long
  20. #define sz(x) x.size()
  21. #define sqr(x) (1LL * (x) * (x))
  22. #define all(x) x.begin(), x.end()
  23. #define fill(f,x) memset(f,x,sizeof(f))
  24. #define FOR(i,l,r) for(int i=l;i<=r;i++)
  25. #define FOD(i,r,l) for(int i=r;i>=l;i--)
  26. #define debug(x) cout << #x << " = " << x << '\n'
  27. #define ii pair<int,int>
  28. #define iii pair<int,ii>
  29. #define di pair<ii,ii>
  30. #define vi vector<int>
  31. #define vii vector<ii>
  32. #define mii map<int,int>
  33. #define fi first
  34. #define se second
  35. #define pb push_back
  36. #define MOD 1000000007
  37. #define __lcm(a,b) (1ll * ((a) / __gcd((a), (b))) * (b))
  38. #define YES cout << "YES\n"
  39. #define NO cout << "NO\n"
  40. #define MASK(i) (1LL << (i))
  41. #define c_bit(i) __builtin_popcountll(i)
  42. #define BIT(x,i) ((x) & MASK(i))
  43. #define SET_ON(x,i) ((x) | MASK(i))
  44. #define SET_OFF(x,i) ((x) & ~MASK(i))
  45. #define oo 1e18
  46. #define name ""
  47. #define endl '\n'
  48. #define memory() cerr << abs(&M2-&M1)/1024.0/1024 << " MB" << endl
  49. #define time() cerr << endl << "-------------Time:" << 1000.0 * clock() / CLOCKS_PER_SEC << "ms." << endl
  50. template<typename T> bool maximize(T &res, const T &val) { if (res < val){ res = val; return true; }; return false; }
  51. template<typename T> bool minimize(T &res, const T &val) { if (res > val){ res = val; return true; }; return false; }
  52. template <class T> using ordered_set = tree <T, null_type, less_equal <T>, rb_tree_tag,tree_order_statistics_node_update>;
  53. const int N = (int)1e5+10;
  54. int n,q,a[N];
  55.  
  56. struct Query {
  57. 	int l,r,M;
  58. } query[N];
  59.  
  60. namespace sub1 {
  61.  
  62. 	bool approved() {
  63. 		return n <= 3e3 and q <= 3e3;
  64. 	}
  65.  
  66. 	void solve(void)
  67. 	{
  68. 		FOR(i,1,q)
  69. 		{
  70. 			int l = query[i].l, r = query[i].r, M = query[i].M;
  71. 			int ans = 0;
  72. 			FOR(i,l,r) ans += a[i]%M;
  73. 			cout << ans << endl;
  74. 		}
  75. 	}
  76.  
  77. }
  78.  
  79. namespace sub2 {
  80.  
  81. 	const int S = 1e3;
  82. 	int pre[N],ans[N],upd[N*3],sum[S+1],mx;
  83. 	vector<array<int,3>> vec[N];
  84.  
  85. 	/*
  86. 		Tổng i từ l -> r với a[i]%M
  87. 		<-> a[i]-M*[a[i]/M] với phép [] làm tròn xuống
  88. 	*/
  89.  
  90. 	void add(int l, int r, int val)
  91. 	{
  92. 		if (l > r or r > mx) return;
  93. 		if (r > mx) r = mx;
  94. 		upd[l] += val;
  95. 		upd[r+1] -= val;
  96. 		sum[(l-1)/S] += val;
  97. 		sum[r/S] -= val;
  98. 	}
  99.  
  100. 	int get(int pos)
  101. 	{
  102. 		if (pos > mx) pos = mx;
  103. 		int block = (pos-1)/S, ans = 0;
  104. 		FOR(bl,0,block-1) ans += sum[bl];
  105. 		FOR(i,block*S+1,pos) ans += upd[i];
  106. 		return ans;
  107. 	}
  108.  
  109. 	void solve(void)
  110. 	{
  111. 		mx = *max_element(a+1,a+n+1);
  112. 		FOR(i,1,n) pre[i] = pre[i-1]+a[i];
  113. 		FOR(i,1,q)
  114. 		{
  115. 			int l = query[i].l, r = query[i].r, M = query[i].M;
  116. 			if (M <= (N-10)*3)
  117. 			{
  118. 				vec[r].pb({M,i,1});
  119. 				vec[l-1].pb({M,i,-1});
  120. 			}
  121. 			ans[i] = pre[r]-pre[l-1];
  122. 		}	
  123. 		FOR(i,1,n)
  124. 		{
  125. 			int m = 1;
  126. 			while (m <= a[i])
  127. 			{
  128. 				int v = a[i]/m, r = a[i]/v;
  129. 				add(m,r,v);
  130. 				m = r+1;
  131. 			}
  132. 			for (auto [m,id,delta] : vec[i])
  133. 				ans[id] -= delta*m*get(m);
  134. 		}
  135. 		FOR(i,1,q) cout << ans[i] << endl;
  136. 	}
  137.  
  138. }
  139.  
  140. bool M2;
  141. signed main()
  142. {
  143.     fast;
  144.     if (fopen(name".inp","r"))
  145.     {
  146.     	freopen(name".inp","r",stdin);
  147.     	freopen(name".out","w",stdout);
  148.     }
  149.     cin >> n >> q;
  150.     FOR(i,1,n) cin >> a[i];
  151.     FOR(i,1,q) 
  152.     {
  153.     	cin >> query[i].l >> query[i].r >> query[i].M;
  154.     	query[i].l++; query[i].r++;
  155.     }
  156.     if (sub1::approved()) return sub1::solve(), time(), memory(), 0;
  157.     sub2::solve();
  158.     time();
  159.     memory();
  160.     return 0;
  161. }
  162. // ██░ ██  █    ██  ███▄    █   ▄████
  163. //▓██░ ██▒ ██  ▓██▒ ██ ▀█   █  ██▒ ▀█▒
  164. //▒██▀▀██░▓██  ▒██░▓██  ▀█ ██▒▒██░▄▄▄░
  165. //░▓█ ░██ ▓▓█  ░██░▓██▒  ▐▌██▒░▓█  ██▓
  166. //░▓█▒░██▓▒▒█████▓ ▒██░   ▓██░░▒▓███▀▒
  167. // ▒ ░░▒░▒░▒▓▒ ▒ ▒ ░ ▒░   ▒ ▒  ░▒   ▒
  168. // ▒ ░▒░ ░░░▒░ ░ ░ ░ ░░   ░ ▒░  ░   ░
  169. // ░  ░░ ░ ░░░ ░ ░    ░   ░ ░ ░ ░   ░
  170. // ░  ░  ░   ░              ░       ░
Success #stdin #stdout #stderr 0.01s 8140KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
Standard output is empty
stderr
copy 
-------------Time:6.744ms.
9.16397 MB
https://ideone.com/E6ixF7
language:
C++14 (gcc 8.3)
created:
6 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!