結果
問題 | No.878 Range High-Element Query |
ユーザー | akakimidori |
提出日時 | 2019-09-06 23:12:41 |
言語 | Rust (1.77.0 + proconio) |
結果 |
AC
|
実行時間 | 44 ms / 2,000 ms |
コード長 | 5,739 bytes |
コンパイル時間 | 12,111 ms |
コンパイル使用メモリ | 408,216 KB |
実行使用メモリ | 18,984 KB |
最終ジャッジ日時 | 2024-06-24 21:08:29 |
合計ジャッジ時間 | 13,649 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
5,248 KB |
testcase_01 | AC | 1 ms
5,248 KB |
testcase_02 | AC | 1 ms
5,248 KB |
testcase_03 | AC | 1 ms
5,376 KB |
testcase_04 | AC | 1 ms
5,376 KB |
testcase_05 | AC | 1 ms
5,376 KB |
testcase_06 | AC | 1 ms
5,376 KB |
testcase_07 | AC | 1 ms
5,376 KB |
testcase_08 | AC | 1 ms
5,376 KB |
testcase_09 | AC | 1 ms
5,376 KB |
testcase_10 | AC | 1 ms
5,376 KB |
testcase_11 | AC | 40 ms
17,536 KB |
testcase_12 | AC | 31 ms
15,360 KB |
testcase_13 | AC | 32 ms
13,952 KB |
testcase_14 | AC | 26 ms
11,520 KB |
testcase_15 | AC | 32 ms
15,232 KB |
testcase_16 | AC | 41 ms
18,176 KB |
testcase_17 | AC | 44 ms
18,984 KB |
testcase_18 | AC | 43 ms
18,816 KB |
ソースコード
// 1-indexedなBIT // 座標xへの加算、[1,x]の和, ([1,x]の和)>=s となる最小のxの探索 // ---------- begin fenwick tree ---------- #[allow(dead_code)] mod fenwick { type T = i32; const ZERO: T = 0; pub struct Tree { a: Vec<T>, } impl Tree { pub fn new(n: usize) -> Tree { Tree { a: vec![ZERO; n + 1], } } pub fn add(&mut self, mut x: usize, v: T) { while x < self.a.len() { self.a[x] += v; x += x & (!x + 1); } } pub fn sum(&self, mut x: usize) -> T { debug_assert!(x < self.a.len()); let mut res = ZERO; while x > 0 { res += self.a[x]; x -= x & (!x + 1); } res } pub fn search(&self, s: T) -> usize { debug_assert!(self.sum(self.a.len() - 1) >= s); let mut k = 1; while 2 * k < self.a.len() { k *= 2; } let mut x = 0; let mut w = ZERO; while k > 0 { if k + x < self.a.len() && w + self.a[x + k] < s { w += self.a[x + k]; x += k; } k >>= 1; } x + 1 } } } // ---------- end fenwick tree ---------- //---------- begin SegmentTree Point update Range query ---------- mod segment_tree { pub struct PURQ<T: Clone, F: Fn(T, T) -> T> { n: usize, a: Vec<T>, id: T, op: F, } #[allow(dead_code)] impl<T: Clone, F: Fn(T, T) -> T> PURQ<T, F> { pub fn new(n: usize, id: T, op: F) -> PURQ<T, F> { let mut k = 1; while k < n { k *= 2; } PURQ { n: k, a: vec![id.clone(); 2 * k], id: id, op: op, } } pub fn update(&mut self, x: usize, v: T) { let mut k = self.n + x; let a = &mut self.a; a[k] = v; k >>= 1; while k > 0 { a[k] = (self.op)(a[2 * k].clone(), a[2 * k + 1].clone()); k >>= 1; } } pub fn update_tmp(&mut self, x: usize, v: T) { self.a[x + self.n] = v; } pub fn update_all(&mut self) { for k in (1..(self.n)).rev() { self.a[k] = (self.op)(self.a[2 * k].clone(), self.a[2 * k + 1].clone()); } } pub fn find(&self, mut l: usize, mut r: usize) -> T { let mut p = self.id.clone(); let mut q = self.id.clone(); l += self.n; r += self.n; while l < r { if (l & 1) == 1 { p = (self.op)(p, self.a[l].clone()); l += 1; } if (r & 1) == 1 { r -= 1; q = (self.op)(self.a[r].clone(), q); } l >>= 1; r >>= 1; } (self.op)(p, q) } } } //---------- end SegmentTree Point update Range query ---------- //https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8 より macro_rules! input { (source = $s:expr, $($r:tt)*) => { let mut iter = $s.split_whitespace(); input_inner!{iter, $($r)*} }; ($($r:tt)*) => { let s = { use std::io::Read; let mut s = String::new(); std::io::stdin().read_to_string(&mut s).unwrap(); s }; let mut iter = s.split_whitespace(); input_inner!{iter, $($r)*} }; } macro_rules! input_inner { ($iter:expr) => {}; ($iter:expr, ) => {}; ($iter:expr, $var:ident : $t:tt $($r:tt)*) => { let $var = read_value!($iter, $t); input_inner!{$iter $($r)*} }; } macro_rules! read_value { ($iter:expr, ( $($t:tt),* )) => { ( $(read_value!($iter, $t)),* ) }; ($iter:expr, [ $t:tt ; $len:expr ]) => { (0..$len).map(|_| read_value!($iter, $t)).collect::<Vec<_>>() }; ($iter:expr, chars) => { read_value!($iter, String).chars().collect::<Vec<char>>() }; ($iter:expr, usize1) => { read_value!($iter, usize) - 1 }; ($iter:expr, $t:ty) => { $iter.next().unwrap().parse::<$t>().expect("Parse error") }; } // ここまで use std::io::Write; fn run() { let out = std::io::stdout(); let mut out = std::io::BufWriter::new(out.lock()); input! { n: usize, q: usize, a: [usize; n], p: [(u8, usize, usize); q], } type T = (usize, usize); let f = |a: T, b: T| if a.0 > b.0 {a} else {b}; let mut s = segment_tree::PURQ::new(n, (0, 0), f); for i in 0..n { s.update_tmp(i, (a[i], i)); } s.update_all(); let mut range = vec![(0, 0); n]; let mut stack = vec![]; stack.push((0, n)); while let Some((l, r)) = stack.pop() { if l >= r { continue; } let (_, k) = s.find(l, r); range[k] = (l, r); stack.push((l, k)); stack.push((k + 1, r)); } let mut g = vec![vec![]; n]; for (i, &(_, l, r)) in p.iter().enumerate() { g[r - 1].push((l, i)); } let mut ans = vec![0; q]; let mut bit = fenwick::Tree::new(n); for i in 0..n { bit.add(range[i].0 + 1, 1); bit.add(i + 2, -1); for &(l, k) in g[i].iter() { ans[k] = bit.sum(l); } } for v in ans { writeln!(out, "{}", v).unwrap(); } } fn main() { run(); }