結果

問題 No.752 mod数列
コンテスト
ユーザー te-sh
提出日時 2018-11-16 12:25:21
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 70 ms / 2,000 ms
コード長 893 bytes
コンパイル時間 750 ms
コンパイル使用メモリ 103,288 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-13 01:55:32
合計ジャッジ時間 3,569 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.container, std.conv, std.math, std.range, std.typecons, std.stdio, std.string;

auto rdsp(){return readln.splitter;}
void pick(R,T)(ref R r,ref T t){t=r.front.to!T;r.popFront;}
void readV(T...)(ref T t){auto r=rdsp;foreach(ref v;t)pick(r,v);}

void main()
{
  long p; int q; readV(p, q);

  auto maxM = min(p, 10L^^4), maxK = p/(maxM+1);

  auto b = new long[](maxK+1);
  foreach (k; 1..maxK+1)
    b[k] = b[k-1] + p%k;

  auto c = new long[](maxM+2); c[$-1] = b[$-1];
  foreach_reverse (m; 1..maxM+1) {
    auto u = p/(m+1)+1, v = p/m;
    c[m] = c[m+1] + p*(v-u+1) - m*(v*(v+1)-u*(u-1))/2;
  }

  auto calc(long x)
  {
    if (x <= maxK) {
      return b[x];
    } else {
      auto m = p/x, u = p/(m+1)+1;
      return c[m+1] + p*(x-u+1) - m*(x*(x+1)-u*(u-1))/2;
    }
  }

  foreach (_; 0..q) {
    long l, r; readV(l, r);
    writeln(calc(r) - calc(l-1));
  }
}
0