結果
問題 | No.469 区間加算と一致検索の問題 |
ユーザー |
|
提出日時 | 2017-12-12 12:46:43 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 476 ms / 5,000 ms |
コード長 | 1,355 bytes |
コンパイル時間 | 2,204 ms |
コンパイル使用メモリ | 167,476 KB |
実行使用メモリ | 44,596 KB |
最終ジャッジ日時 | 2024-06-12 23:01:35 |
合計ジャッジ時間 | 24,880 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 49 |
ソースコード
import std.algorithm, std.conv, std.range, std.stdio, std.string;import std.random;void main(){auto rd = readln.split.to!(int[]), n = rd[0], q = rd[1];auto p = primes!int(5*10^^7)[10000..$];auto b = p.randomSample(n).array;auto bc = new long[](n+1);foreach (i; 0..n) bc[i+1] = bc[i] + b[i];int[long] h;h[0] = 0;auto curr = 0L;foreach (i; 0..q) {auto line = readln.chomp;if (line[0] == '?') {writeln(h[curr]);} else {auto rd2 = line[2..$].splitter;auto l = rd2.front.to!int; rd2.popFront();auto r = rd2.front.to!int; rd2.popFront();auto k = rd2.front.to!int;curr += (bc[r] - bc[l]) * k;if (curr !in h)h[curr] = i+1;}}}pure T[] primes(T)(T n){import std.algorithm, std.bitmanip, std.conv, std.range;auto sieve = BitArray();sieve.length((n + 1) / 2);sieve = ~sieve;foreach (p; 1..((nsqrt(n) - 1) / 2 + 1))if (sieve[p])for (auto q = p * 3 + 1; q < (n + 1) / 2; q += p * 2 + 1)sieve[q] = false;auto r = sieve.bitsSet.map!(to!T).map!("a * 2 + 1").array;r[0] = 2;return r;}pure T nsqrt(T)(T n){import std.algorithm, std.conv, std.range, core.bitop;if (n <= 1) return n;T m = 1 << (n.bsr / 2 + 1);return iota(1, m).map!"a * a".assumeSorted!"a <= b".lowerBound(n).length.to!T;}