結果
問題 | No.752 mod数列 |
ユーザー |
|
提出日時 | 2018-11-10 05:00:06 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,552 ms / 2,000 ms |
コード長 | 1,064 bytes |
コンパイル時間 | 2,071 ms |
コンパイル使用メモリ | 195,812 KB |
最終ジャッジ日時 | 2025-01-06 16:19:32 |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
#include <bits/stdc++.h> #define REP(i, n) for (int i = 0; (i) < (int)(n); ++ (i)) #define REP3(i, m, n) for (int i = (m); (i) < (int)(n); ++ (i)) using ll = long long; using namespace std; int main() { int p, q; cin >> p >> q; constexpr int K = 1000000; vector<ll> small(K); small[0] = 0; small[1] = 0; REP3 (i, 1, K - 1) { small[i + 1] = small[i] + (p % i); } auto solve1 = [&](int r) { ll large = 0; if (p + 1 < r) { large += p * (r - p - 1ll); r = p + 1; } for (int k = 1; r >= small.size(); ++ k) { int l = p / k + 1; if (r <= l) continue; ll a1 = p % l; ll a2 = p % (l + 1); ll d = a2 - a1; ll an = a1 + (r - l - 1) * d; large += (r - l) * (a1 + an) / 2; r = l; } return small[r] + large; }; REP (i, q) { int l, r; cin >> l >> r; ll answer = solve1(r + 1) - solve1(l); cout << answer << endl; } return 0; }