結果

問題 No.752 mod数列
ユーザー rsk0315rsk0315
提出日時 2020-04-23 18:06:49
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 59 ms / 2,000 ms
コード長 983 bytes
コンパイル時間 840 ms
コンパイル使用メモリ 63,872 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-04-21 19:52:22
合計ジャッジ時間 4,376 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 1 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 1 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 5 ms
5,376 KB
testcase_11 AC 3 ms
5,376 KB
testcase_12 AC 5 ms
5,376 KB
testcase_13 AC 5 ms
5,376 KB
testcase_14 AC 6 ms
5,376 KB
testcase_15 AC 42 ms
5,376 KB
testcase_16 AC 59 ms
5,376 KB
testcase_17 AC 55 ms
5,376 KB
testcase_18 AC 40 ms
5,376 KB
testcase_19 AC 47 ms
5,376 KB
testcase_20 AC 42 ms
5,376 KB
testcase_21 AC 42 ms
5,376 KB
testcase_22 AC 40 ms
5,376 KB
testcase_23 AC 41 ms
5,376 KB
testcase_24 AC 41 ms
5,376 KB
testcase_25 AC 17 ms
5,376 KB
testcase_26 AC 42 ms
5,376 KB
testcase_27 AC 31 ms
5,376 KB
testcase_28 AC 40 ms
5,376 KB
testcase_29 AC 47 ms
5,376 KB
testcase_30 AC 44 ms
5,376 KB
testcase_31 AC 2 ms
5,376 KB
testcase_32 AC 2 ms
5,376 KB
testcase_33 AC 5 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdint>
#include <cstdio>
#include <algorithm>
#include <vector>

int main() {
  intmax_t p;
  int q;
  scanf("%jd %d", &p, &q);

  std::vector<intmax_t> d{0};
  for (intmax_t i = 1; i*i <= p; ++i) {
    d.push_back(i);
    if (i*i < p) d.push_back(p/i);
  }

  std::sort(d.begin(), d.end());

  std::vector<intmax_t> s{0};
  for (size_t i = 1; i < d.size(); ++i) {
    intmax_t dl = d[i-1]+1;
    intmax_t dr = d[i];
    intmax_t sum = ((p % dl) + (p % dr)) * (dr-dl+1) / 2;
    s.push_back(sum);
  }

  s.insert(s.begin(), 0);
  for (size_t i = 1; i < s.size(); ++i) s[i] += s[i-1];

  auto f = [&](intmax_t r) -> intmax_t {
    if (r == 0) return 0;
    auto it = std::upper_bound(d.begin(), d.end(), r);
    size_t j = it - d.begin();
    intmax_t dl = it[-1]+1;
    intmax_t dr = r;
    return s[j] + ((p % dl) + (p % dr)) * (dr-dl+1) / 2;
  };

  for (int i = 0; i < q; ++i) {
    intmax_t l, r;
    scanf("%jd %jd", &l, &r);
    printf("%jd\n", f(r)-f(l-1));
  }
}
0