結果
問題 | No.752 mod数列 |
ユーザー | tsutaj |
提出日時 | 2019-03-29 22:19:45 |
言語 | C++14 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 864 ms / 2,000 ms |
コード長 | 2,363 bytes |
コンパイル時間 | 799 ms |
コンパイル使用メモリ | 101,184 KB |
実行使用メモリ | 81,920 KB |
最終ジャッジ日時 | 2024-04-24 16:46:24 |
合計ジャッジ時間 | 12,202 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 147 ms
81,664 KB |
testcase_01 | AC | 145 ms
81,920 KB |
testcase_02 | AC | 140 ms
81,920 KB |
testcase_03 | AC | 144 ms
81,792 KB |
testcase_04 | AC | 150 ms
81,792 KB |
testcase_05 | AC | 154 ms
81,920 KB |
testcase_06 | AC | 155 ms
81,664 KB |
testcase_07 | AC | 150 ms
81,792 KB |
testcase_08 | AC | 150 ms
81,920 KB |
testcase_09 | AC | 147 ms
81,664 KB |
testcase_10 | AC | 156 ms
81,920 KB |
testcase_11 | AC | 154 ms
81,664 KB |
testcase_12 | AC | 155 ms
81,920 KB |
testcase_13 | AC | 152 ms
81,792 KB |
testcase_14 | AC | 156 ms
81,792 KB |
testcase_15 | AC | 313 ms
81,664 KB |
testcase_16 | AC | 327 ms
81,664 KB |
testcase_17 | AC | 864 ms
81,664 KB |
testcase_18 | AC | 329 ms
81,920 KB |
testcase_19 | AC | 333 ms
81,664 KB |
testcase_20 | AC | 363 ms
81,664 KB |
testcase_21 | AC | 373 ms
81,664 KB |
testcase_22 | AC | 379 ms
81,792 KB |
testcase_23 | AC | 355 ms
81,664 KB |
testcase_24 | AC | 378 ms
81,920 KB |
testcase_25 | AC | 212 ms
81,664 KB |
testcase_26 | AC | 564 ms
81,792 KB |
testcase_27 | AC | 393 ms
81,920 KB |
testcase_28 | AC | 518 ms
81,792 KB |
testcase_29 | AC | 510 ms
81,920 KB |
testcase_30 | AC | 365 ms
81,664 KB |
testcase_31 | AC | 151 ms
81,920 KB |
testcase_32 | AC | 148 ms
81,920 KB |
testcase_33 | AC | 150 ms
81,792 KB |
ソースコード
// #define _GLIBCXX_DEBUG // for STL debug (optional) #include <iostream> #include <iomanip> #include <cstdio> #include <string> #include <cstring> #include <deque> #include <list> #include <queue> #include <stack> #include <vector> #include <utility> #include <algorithm> #include <map> #include <set> #include <complex> #include <cmath> #include <limits> #include <cfloat> #include <climits> #include <ctime> #include <cassert> #include <numeric> #include <fstream> #include <functional> #include <bitset> using namespace std; #define debug(...) fprintf(stderr, __VA_ARGS__) #define int long long int template<typename T> void chmax(T &a, T b) {a = max(a, b);} template<typename T> void chmin(T &a, T b) {a = min(a, b);} template<typename T> void chadd(T &a, T b) {a = a + b;} typedef pair<int, int> pii; typedef long long ll; int dx[] = {0, 0, 1, -1}; int dy[] = {1, -1, 0, 0}; const ll INF = 1001001001001001LL; const ll MOD = 1000000007LL; const int K = 10000000; int sum[K+10]; signed main() { int P, Q; cin >> P >> Q; for(int i=1; i<=K; i++) { int mod = P % i; sum[i] = sum[i-1] + mod; } while(Q--) { if(Q % 10000 == 0) { fprintf(stderr, "Q = %lld\n", Q); } int l, r; cin >> l >> r; int ans = 0; if(r <= K) { ans = sum[r] - sum[l-1]; cout << ans << endl; continue; } else if(l <= K and K <= r) { ans = sum[K] - sum[l-1]; l = K+1; } int vl = P / l, vr = P / r, cur = l; int cnt = 0; for(int div=vl; div>=vr; div--) { cnt++; if(div == 0) { int ll = max(l, P + 1); int rr = r; int add = 2 * P * (rr - ll + 1) / 2; ans += add; } else { int ll = max(P / (div + 1) + 1, l); int rr = min(P / div, r); int lll = P % ll, rrr = P % rr; int add = (lll + rrr) * (rr - ll + 1) / 2; ans += add; if(P / cur != P / (cur + 1)) { cur++; div = P / cur + 1; } } } // fprintf(stderr, "div cnt = %lld\n", cnt); cout << ans << endl; } return 0; }