結果
問題 | No.2759 Take Pictures, Elements? |
ユーザー |
![]() |
提出日時 | 2024-05-17 22:44:45 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 37 ms / 2,000 ms |
コード長 | 4,972 bytes |
コンパイル時間 | 14,931 ms |
コンパイル使用メモリ | 406,740 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-20 14:46:36 |
合計ジャッジ時間 | 16,206 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
#![allow(non_snake_case, unused_imports)]use proconio::{fastout, input, marker::*};#[fastout]fn main() {input! {N: usize, Q: usize,A: [usize; N],B: [usize; Q],}let mut dp = vec![usize::MAX; N];dp[0] = 0;// margin// for i in 0..N {// fromleft.insert(i, N - 1 - i);// fromright.insert(i, i);// }let (z, zi) = coordinate_compression(&A);let mut pos = vec![vec![]; z.len()];for (i, a) in A.iter().enumerate() {let &p = zi.get(a).unwrap();pos[p].push(i);}for i in 0..Q {let &zb = zi.get(&B[i]).unwrap();let ps = &pos[zb];let mut fromleft = SegmentTree::new(N, |&x, &y| std::cmp::min(x, y), usize::MAX);let mut fromright = SegmentTree::new(N, |&x, &y| std::cmp::min(x, y), usize::MAX);if i == 0 {fromleft.insert(0, N - 1);fromright.insert(0, 0);} else {let &zbp = zi.get(&B[i - 1]).unwrap();let psp = &pos[zbp];for &p in psp {fromleft.insert(p, dp[p] + N - 1 - p);fromright.insert(p, dp[p] + p);}}dp = vec![usize::MAX; N];for &p in ps {let fl = fromleft.prod(..=p) - (N - 1 - p);let fr = fromright.prod(p..) - p;dp[p] = std::cmp::min(dp[p], fl);dp[p] = std::cmp::min(dp[p], fr);}}println!("{}", dp.iter().min().unwrap());}pub fn coordinate_compression<T: std::cmp::Ord + Clone + Copy>(values: &[T],) -> (Vec<T>, std::collections::BTreeMap<T, usize>) {let s: std::collections::BTreeSet<T> = values.iter().cloned().collect();let z: Vec<T> = s.iter().cloned().collect();let zinv: std::collections::BTreeMap<T, usize> =z.iter().enumerate().map(|(i, &v)| (v, i)).collect();(z, zinv)}pub struct SegmentTree<M> {size: usize,tree: Vec<M>,op: fn(&M, &M) -> M,id: M,}impl<M: Copy + PartialEq> SegmentTree<M> {/// self = [id; size], self.op = op, self.id = idpub fn new(size: usize, op: fn(&M, &M) -> M, id: M) -> Self {return Self {size: size,tree: vec![id; 2 * size],op: op,id: id,};}/// self = arr, self.op = op, self.id = idpub fn from(arr: &[M], op: fn(&M, &M) -> M, id: M) -> Self {let size = arr.len();let mut tree = vec![id; 2 * size];for i in 0..size {tree[i + size] = arr[i];assert!(op(&id, &arr[i]) == arr[i]);}for i in (1..size).rev() {tree[i] = op(&tree[i << 1], &tree[i << 1 | 1]);}return Self {size: size,tree: tree,op: op,id: id,};}/// self[pos] <- valuepub fn insert(&mut self, mut pos: usize, value: M) -> () {pos += self.size;self.tree[pos] = value;while pos > 1 {pos >>= 1;self.tree[pos] = (self.op)(&self.tree[pos << 1], &self.tree[pos << 1 | 1]);}}/// return self[pos]pub fn get_point(&self, pos: usize) -> M {return self[pos];}/// return Π_{i ∈ [left, right)} self[i]pub fn get(&self, left: usize, right: usize) -> M {let (mut l, mut r) = (left + self.size, right + self.size);let (mut vl, mut vr) = (self.id, self.id);while l < r {if l & 1 == 1 {vl = (self.op)(&vl, &self.tree[l]);l += 1;}if r & 1 == 1 {r -= 1;vr = (self.op)(&self.tree[r], &vr);}l >>= 1;r >>= 1;}return (self.op)(&vl, &vr);}/// return Π_{i ∈ range} self[i]pub fn prod<R: std::ops::RangeBounds<usize>>(&self, range: R) -> M {let left = match range.start_bound() {std::ops::Bound::Included(&l) => l,std::ops::Bound::Excluded(&l) => l + 1,std::ops::Bound::Unbounded => 0,};let right = match range.end_bound() {std::ops::Bound::Included(&r) => r + 1,std::ops::Bound::Excluded(&r) => r,std::ops::Bound::Unbounded => self.size,};return self.get(left, right);}}impl<M> std::ops::Index<usize> for SegmentTree<M> {type Output = M;fn index(&self, index: usize) -> &Self::Output {&self.tree[index + self.size]}}impl<M: std::fmt::Display> std::fmt::Display for SegmentTree<M> {fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {write!(f,"{}",self.tree[self.size..].iter().map(|x| x.to_string()).collect::<Vec<_>>().join(" "))}}