結果

問題 No.752 mod数列
ユーザー rsk0315rsk0315
提出日時 2020-04-23 19:10:46
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 58 ms / 2,000 ms
コード長 1,047 bytes
コンパイル時間 654 ms
コンパイル使用メモリ 59,648 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-13 19:41:39
合計ジャッジ時間 4,412 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 1 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 1 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 1 ms
5,248 KB
testcase_09 AC 2 ms
5,248 KB
testcase_10 AC 4 ms
5,248 KB
testcase_11 AC 3 ms
5,248 KB
testcase_12 AC 4 ms
5,248 KB
testcase_13 AC 4 ms
5,248 KB
testcase_14 AC 4 ms
5,248 KB
testcase_15 AC 38 ms
5,248 KB
testcase_16 AC 58 ms
5,248 KB
testcase_17 AC 58 ms
5,248 KB
testcase_18 AC 41 ms
5,248 KB
testcase_19 AC 46 ms
5,248 KB
testcase_20 AC 41 ms
5,248 KB
testcase_21 AC 41 ms
5,248 KB
testcase_22 AC 41 ms
5,248 KB
testcase_23 AC 40 ms
5,248 KB
testcase_24 AC 43 ms
5,248 KB
testcase_25 AC 18 ms
5,248 KB
testcase_26 AC 41 ms
5,248 KB
testcase_27 AC 30 ms
5,248 KB
testcase_28 AC 38 ms
5,248 KB
testcase_29 AC 45 ms
5,248 KB
testcase_30 AC 44 ms
5,248 KB
testcase_31 AC 2 ms
5,248 KB
testcase_32 AC 2 ms
5,248 KB
testcase_33 AC 3 ms
5,248 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};
  {
    std::vector<intmax_t> tmp;
    for (intmax_t i = 1; i*i <= p; ++i) {
      d.push_back(i);
      if (i*i < p) tmp.push_back(p/i);
    }
    d.insert(d.end(), tmp.rbegin(), tmp.rend());
  }

  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