#654694 (Rust) No.1526 Sum of Mex 2

提出ソース

結果

問題	No.1526 Sum of Mex 2
ユーザー	Rheo Tommy
提出日時	2021-05-06 23:54:00
言語	Rust (1.77.0 + proconio)
結果	AC
実行時間	285 ms / 3,000 ms
コード長	11,972 bytes
コンパイル時間	12,723 ms
コンパイル使用メモリ	379,408 KB
実行使用メモリ	24,024 KB
最終ジャッジ日時	2024-11-08 21:31:35
合計ジャッジ時間	19,257 ms
ジャッジサーバーID （参考情報）	judge4 / judge3

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
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,248 KB
testcase_04	AC	1 ms 5,248 KB
testcase_05	AC	1 ms 5,248 KB
testcase_06	AC	1 ms 5,248 KB
testcase_07	AC	1 ms 5,248 KB
testcase_08	AC	1 ms 5,248 KB
testcase_09	AC	1 ms 5,248 KB
testcase_10	AC	2 ms 5,248 KB
testcase_11	AC	1 ms 5,248 KB
testcase_12	AC	1 ms 5,248 KB
testcase_13	AC	33 ms 5,248 KB
testcase_14	AC	39 ms 6,400 KB
testcase_15	AC	51 ms 6,912 KB
testcase_16	AC	271 ms 23,360 KB
testcase_17	AC	185 ms 19,968 KB
testcase_18	AC	19 ms 5,248 KB
testcase_19	AC	16 ms 5,248 KB
testcase_20	AC	110 ms 12,032 KB
testcase_21	AC	258 ms 22,888 KB
testcase_22	AC	145 ms 13,568 KB
testcase_23	AC	271 ms 23,132 KB
testcase_24	AC	279 ms 23,552 KB
testcase_25	AC	256 ms 22,912 KB
testcase_26	AC	269 ms 23,424 KB
testcase_27	AC	265 ms 23,552 KB
testcase_28	AC	281 ms 23,808 KB
testcase_29	AC	278 ms 23,596 KB
testcase_30	AC	274 ms 24,024 KB
testcase_31	AC	283 ms 23,920 KB
testcase_32	AC	265 ms 23,296 KB
testcase_33	AC	285 ms 22,960 KB
evil_largest	AC	990 ms 74,880 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#![allow(unused_macros)]
#![allow(dead_code)]
#![allow(unused_imports)]

// # ファイル構成
// - use 宣言
// - lib モジュール
// - main 関数
// - basic モジュール
//
// 常に使うテンプレートライブラリは basic モジュール内にあります。
// 問題に応じて使うライブラリ lib モジュール内にコピペしています。
// ライブラリのコードはこちら → https://github.com/RheoTommy/at_coder
// Twitter はこちら → https://twitter.com/RheoTommy

use std::collections::*;
use std::io::{stdout, BufWriter, Write};

use crate::basic::*;
use crate::lib::*;

pub mod lib {
    pub trait LazyMonoid {
        type M: Clone;
        type L: Clone;
        fn e() -> Self::M;
        fn id() -> Self::L;
    }

    pub struct LazySegTree<M: LazyMonoid> {
        n: usize,
        height: usize,
        len: usize,
        data: Vec<(M::M, M::L)>,
        fold: Box<dyn FnMut(&M::M, &M::M) -> M::M>,
        eval: Box<dyn FnMut(&M::M, &M::L) -> M::M>,
        merge: Box<dyn FnMut(&M::L, &M::L) -> M::L>,
    }

    impl<M: LazyMonoid> std::fmt::Debug for LazySegTree<M>
        where
            <M as LazyMonoid>::M: std::fmt::Debug,
            <M as LazyMonoid>::L: std::fmt::Debug,
    {
        fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {
            for i in 0..=self.height {
                for j in 0..1 << i {
                    write!(f, "{:?}   ", &self.data[(1 << i) + j])?;
                }
                writeln!(f)?;
            }
            write!(f, "")
        }
    }

    impl<M: LazyMonoid> LazySegTree<M> {
        pub fn new(
            n: usize,
            fold: Box<dyn FnMut(&M::M, &M::M) -> M::M>,
            eval: Box<dyn FnMut(&M::M, &M::L) -> M::M>,
            merge: Box<dyn FnMut(&M::L, &M::L) -> M::L>,
        ) -> Self {
            let len = n;
            let n = n.next_power_of_two();
            let height = n.trailing_zeros() as usize;
            Self {
                n,
                height,
                len,
                fold,
                eval,
                merge,
                data: vec![(M::e(), M::id()); 2 * n],
            }
        }
        pub fn from_slice(
            a: &[M::M],
            mut fold: Box<dyn FnMut(&M::M, &M::M) -> M::M>,
            eval: Box<dyn FnMut(&M::M, &M::L) -> M::M>,
            merge: Box<dyn FnMut(&M::L, &M::L) -> M::L>,
        ) -> Self {
            let n = a.len().next_power_of_two();
            let height = n.trailing_zeros() as usize;
            let mut data = vec![(M::e(), M::id()); 2 * n];
            for i in 0..a.len() {
                data[i + n].0 = a[i].clone();
            }
            for i in (1..n).rev() {
                data[i].0 = fold(&data[i * 2].0, &data[i * 2 + 1].0);
            }
            Self {
                n,
                height,
                len: a.len(),
                data,
                fold,
                eval,
                merge,
            }
        }
        fn apply(&mut self, x: usize, op: &M::L) {
            let node = &mut self.data[x];
            node.0 = (self.eval)(&node.0, op);
            node.1 = (self.merge)(&node.1, op);
        }
        fn propagate_at(&mut self, x: usize) {
            let op = std::mem::replace(&mut self.data[x].1, M::id());
            self.apply(2 * x, &op);
            self.apply(2 * x + 1, &op);
        }
        fn save_at(&mut self, x: usize) {
            self.data[x].0 = (self.fold)(&self.data[x * 2].0, &self.data[x * 2 + 1].0);
        }
        fn propagate(&mut self, l: usize, r: usize) {
            let l = l + self.n;
            let r = r + self.n;
            for i in (1..=self.height).rev() {
                if (l >> i) << i != l {
                    self.propagate_at(l >> i);
                }
                if (r >> i) << i != r {
                    self.propagate_at((r - 1) >> i);
                }
            }
        }
        fn save(&mut self, l: usize, r: usize) {
            let l = l + self.n;
            let r = r + self.n;
            for i in 1..=self.height {
                if (l >> i) << i != l {
                    self.save_at(l >> i);
                }
                if (r >> i) << i != r {
                    self.save_at((r - 1) >> i);
                }
            }
        }
        pub fn set(&mut self, l: usize, r: usize, op: M::L) {
            assert!(l <= r && r <= self.n);
            if l == r {
                return;
            }
            self.propagate(l, r);
            let mut x = l + self.n;
            let mut y = r + self.n;
            while x < y {
                if x & 1 == 1 {
                    self.apply(x, &op);
                    x += 1;
                }
                if y & 1 == 1 {
                    y -= 1;
                    self.apply(y, &op);
                }
                x >>= 1;
                y >>= 1;
            }
            self.save(l, r)
        }
        pub fn fold(&mut self, l: usize, r: usize) -> M::M {
            assert!(l <= r && r <= self.n);
            if l == r {
                return M::e();
            }
            self.propagate(l, r);
            let mut x = l + self.n;
            let mut y = r + self.n;
            let mut p = M::e();
            let mut q = M::e();
            while x < y {
                if x & 1 == 1 {
                    p = (self.fold)(&p, &self.data[x].0);
                    x += 1;
                }
                if y & 1 == 1 {
                    y -= 1;
                    q = (self.fold)(&self.data[y].0, &q);
                }
                x >>= 1;
                y >>= 1;
            }
            (self.fold)(&p, &q)
        }

        pub fn find_rightest(&mut self, l: usize, mut p: impl FnMut(&M::M) -> bool) -> usize {
            assert!(l < self.n);
            self.propagate(l, self.n);
            if p(&self.fold(l, self.n)) {
                return self.len;
            }
            let mut now = M::e();
            let mut res = l + self.n;
            while res + 1 != res.next_power_of_two() {
                if res & 1 == 1 {
                    let tmp = (self.fold)(&now, &self.data[res].0);
                    if !p(&tmp) {
                        break;
                    }
                    now = tmp;
                    res += 1;
                }
                res >>= 1;
            }
            while res < self.n {
                self.propagate_at(res);
                res <<= 1;
                let tmp = (self.fold)(&now, &self.data[res].0);
                if p(&tmp) {
                    now = tmp;
                    res += 1;
                }
            }
            (res - self.n).min(self.len)
        }
        pub fn find_leftest(&mut self, r: usize, mut p: impl FnMut(&M::M) -> bool) -> usize {
            assert!(r <= self.n);
            self.propagate(0, r);
            if p(&self.fold(0, r)) {
                return 0;
            }
            let mut now = M::e();
            let mut res = r + self.n;
            while res > 1 {
                if res & 1 == 1 {
                    res -= 1;
                    let tmp = (self.fold)(&self.data[res].0, &now);
                    if !p(&tmp) {
                        break;
                    }
                    now = tmp;
                }
                res >>= 1;
            }
            while res < self.n {
                self.propagate_at(res);
                res <<= 1;
                res += 1;
                let tmp = (self.fold)(&self.data[res].0, &now);
                if p(&tmp) {
                    res -= 1;
                    now = tmp;
                }
            }
            (res - self.n + 1).min(self.len)
        }
    }
}

struct M;

impl LazyMonoid for M {
    // max, sum, len
    type M = (usize, usize, usize);
    type L = Option<usize>;

    fn e() -> Self::M {
        (0, 0, 1)
    }

    fn id() -> Self::L {
        None
    }
}

fn main() {
    let mut io = IO::new();
    let n = io.next_usize();
    let a = io.next_vec::<usize>(n);

    let mut ind = vec![vec![n]; n + 1];
    ind[0].push(0);
    for i in 0..n {
        ind[a[i]].push(i);
    }

    for i in ind.iter_mut() {
        i.sort();
        i.reverse();
    }

    let mut b = vec![0; n + 1];
    for i in 1..=n {
        b[i] = b[i - 1].max(*ind[i].last().unwrap());
    }

    let mut lst = LazySegTree::<M>::from_slice(
        &(0..=n).map(|i| (b[i], b[i], 1)).collect::<Vec<_>>(),
        Box::new(|&(a, b, c), &(d, e, f)| (a.max(d), b + e, c + f)),
        Box::new(|&(a, b, c), d| {
            if d.is_none() {
                (a, b, c)
            } else {
                let d = d.unwrap();
                (d, d * c, c)
            }
        }),
        Box::new(|&a, &b|
            match (a, b) {
                (_, Some(t)) => Some(t),
                _ => a,
            }),
    );

    let mut sum = 0;
    // eprintln!("{:?}", lst);

    for l in 0..n {
        sum += (n + 1) * n - lst.fold(0, n + 1).1;

        let al = a[l];
        ind[al].pop();
        let next_ind = *ind[al].last().unwrap();

        let left = lst.find_rightest(0, |&(c, _, _)| c <= next_ind);

        if al < left {
            lst.set(al, left, Some(next_ind));
        }
        lst.set(0, 1, Some(l + 1));
    }

    io.println(sum);
}

pub mod basic {
    pub const U_INF: u64 = (1 << 60) + (1 << 30);
    pub const I_INF: i64 = (1 << 60) + (1 << 30);

    pub struct IO {
        iter: std::str::SplitAsciiWhitespace<'static>,
        buf: std::io::BufWriter<std::io::StdoutLock<'static>>,
    }

    impl Default for IO {
        fn default() -> Self {
            Self::new()
        }
    }

    impl IO {
        pub fn new() -> Self {
            use std::io::*;
            let mut input = String::new();
            std::io::stdin().read_to_string(&mut input).unwrap();
            let input = Box::leak(input.into_boxed_str());
            let out = Box::new(stdout());
            IO {
                iter: input.split_ascii_whitespace(),
                buf: BufWriter::new(Box::leak(out).lock()),
            }
        }
        pub fn next_str(&mut self) -> &str {
            self.iter.next().unwrap()
        }
        pub fn read<T: std::str::FromStr>(&mut self) -> T
            where
                <T as std::str::FromStr>::Err: std::fmt::Debug,
        {
            self.iter.next().unwrap().parse().unwrap()
        }
        pub fn next_usize(&mut self) -> usize {
            self.read()
        }
        pub fn next_uint(&mut self) -> u64 {
            self.read()
        }
        pub fn next_int(&mut self) -> i64 {
            self.read()
        }
        pub fn next_float(&mut self) -> f64 {
            self.read()
        }
        pub fn next_chars(&mut self) -> std::str::Chars {
            self.next_str().chars()
        }
        pub fn next_vec<T: std::str::FromStr>(&mut self, n: usize) -> Vec<T>
            where
                <T as std::str::FromStr>::Err: std::fmt::Debug,
        {
            (0..n).map(|_| self.read()).collect::<Vec<_>>()
        }
        pub fn print<T: std::fmt::Display>(&mut self, t: T) {
            use std::io::Write;
            write!(self.buf, "{}", t).unwrap();
        }
        pub fn println<T: std::fmt::Display>(&mut self, t: T) {
            self.print(t);
            self.print("\n");
        }
        pub fn print_iter<T: std::fmt::Display, I: Iterator<Item=T>>(
            &mut self,
            mut iter: I,
            sep: &str,
        ) {
            if let Some(v) = iter.next() {
                self.print(v);
                for vi in iter {
                    self.print(sep);
                    self.print(vi);
                }
            }
            self.print("\n");
        }
        pub fn flush(&mut self) {
            use std::io::Write;
            self.buf.flush().unwrap();
        }
    }
}

yukicoder

結果

テストケース

ソースコード