結果
問題 | No.2677 Minmax Independent Set |
ユーザー |
|
提出日時 | 2024-03-16 09:02:28 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 373 ms / 2,000 ms |
コード長 | 5,988 bytes |
コンパイル時間 | 13,960 ms |
コンパイル使用メモリ | 382,444 KB |
実行使用メモリ | 139,548 KB |
最終ジャッジ日時 | 2024-09-30 03:53:45 |
合計ジャッジ時間 | 24,630 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 61 |
コンパイルメッセージ
warning: unused import: `BufWriter` --> src/main.rs:5:22 | 5 | use std::io::{Write, BufWriter}; | ^^^^^^^^^ | = note: `#[warn(unused_imports)]` on by default warning: unused import: `Write` --> src/main.rs:5:15 | 5 | use std::io::{Write, BufWriter}; | ^^^^^
ソースコード
#[allow(unused_imports)]use std::cmp::*;#[allow(unused_imports)]use std::collections::*;use std::io::{Write, BufWriter};// https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8macro_rules! input {($($r:tt)*) => {let stdin = std::io::stdin();let mut bytes = std::io::Read::bytes(std::io::BufReader::new(stdin.lock()));let mut next = move || -> String{bytes.by_ref().map(|r|r.unwrap() as char).skip_while(|c|c.is_whitespace()).take_while(|c|!c.is_whitespace()).collect()};input_inner!{next, $($r)*}};}macro_rules! input_inner {($next:expr) => {};($next:expr,) => {};($next:expr, $var:ident : $t:tt $($r:tt)*) => {let $var = read_value!($next, $t);input_inner!{$next $($r)*}};}macro_rules! read_value {($next:expr, ( $($t:tt),* )) => { ($(read_value!($next, $t)),*) };($next:expr, [ $t:tt ; $len:expr ]) => {(0..$len).map(|_| read_value!($next, $t)).collect::<Vec<_>>()};($next:expr, chars) => {read_value!($next, String).chars().collect::<Vec<char>>()};($next:expr, usize1) => (read_value!($next, usize) - 1);($next:expr, [ $t:tt ]) => {{let len = read_value!($next, usize);read_value!($next, [$t; len])}};($next:expr, $t:ty) => ($next().parse::<$t>().expect("Parse error"));}trait Monoid {type V: Clone;type E: Clone;fn add(a: &Self::E, b: &Self::E) -> Self::E;fn e() -> Self::E;fn put_edge(a: &Self::V, eidx: EIdx) -> Self::E;fn put_vertex(a: &Self::E, vidx: usize) -> Self::V;}// Ported from https://trap.jp/post/1702/struct MonoidalReroot<M: Monoid> {#[allow(unused)]pub acc: Vec<Vec<M::E>>,#[allow(unused)]pub acc_rev: Vec<Vec<M::E>>,}type EIdx = usize;impl<M: Monoid> MonoidalReroot<M> {pub fn new(g: &[Vec<(usize, EIdx)>]) -> Self {let n = g.len();let mut acc = vec![vec![]; n];let mut acc_rev = vec![vec![]; n];let mut ch = vec![vec![]; n];for i in 0..n {acc[i] = vec![M::e(); g[i].len() + 1];acc_rev[i] = vec![M::e(); g[i].len() + 1];}Self::dfs1(0, n, &g, &mut ch);for i in 0..n {let l = ch[i].len();ch[i].sort_unstable_by_key(|&(eidx, _)| eidx);for j in 0..l {let val = M::add(&acc[i][j], &ch[i][j].1);acc[i][j + 1] = val;}for j in (0..l).rev() {let val = M::add(&acc_rev[i][j + 1], &ch[i][j].1);acc_rev[i][j] = val;}}Self::dfs2(0, n, &g, &mut ch, &acc, &acc_rev, M::e());for i in 0..n {let l = ch[i].len();ch[i].sort_unstable_by_key(|&(eidx, _)| eidx);for j in 0..l {let val = M::add(&acc[i][j], &ch[i][j].1);acc[i][j + 1] = val;}for j in (0..l).rev() {let val = M::add(&acc_rev[i][j + 1], &ch[i][j].1);acc_rev[i][j] = val;}}MonoidalReroot {acc: acc,acc_rev: acc_rev,}}fn dfs1(v: usize, par: usize, g: &[Vec<(usize, EIdx)>],ch: &mut [Vec<(EIdx, M::E)>],) -> M::V {let mut me = M::e();for i in 0..g[v].len() {let (w, eidx) = g[v][i];if w == par { continue; }let chval = Self::dfs1(w, v, g, ch);let chval = M::put_edge(&chval, eidx);ch[v].push((eidx, chval.clone()));me = M::add(&me, &chval);}M::put_vertex(&me, v)}fn dfs2(v: usize, par: usize, g: &[Vec<(usize, EIdx)>],ch: &mut [Vec<(EIdx, M::E)>],acc: &[Vec<M::E>],acc_rev: &[Vec<M::E>],passed: M::E,) {let mut parenteidx = None;for i in 0..g[v].len() {let (w, eidx) = g[v][i];if w == par {parenteidx = Some(eidx);continue;}let j = ch[v].binary_search_by_key(&eidx, |&(eidx, _)| eidx).unwrap();let inherited = M::add(&acc[v][j], &acc_rev[v][j + 1]);let inherited = if par >= g.len() {inherited} else {M::add(&inherited, &passed)};let inherited = M::put_vertex(&inherited, v);let inherited = M::put_edge(&inherited, eidx);Self::dfs2(w, v, g, ch, acc, acc_rev, inherited);}if let Some(eidx) = parenteidx {ch[v].push((eidx, passed.clone()));}}pub fn query(&self, v: usize) -> M::V {M::put_vertex(&self.acc_rev[v][0], v)}}enum IndepMonoid {}impl Monoid for IndepMonoid {type V = [i64; 2];type E = [i64; 2];fn add(a: &Self::E, b: &Self::E) -> Self::E {[a[0] + b[0], a[1] + b[1]]}fn e() -> Self::E {[0, 0]}fn put_edge(&a: &Self::V, _eidx: EIdx) -> Self::E {[a[0], std::cmp::max(a[0], a[1])]}fn put_vertex(a: &Self::E, _vidx: usize) -> Self::V {[a[1], a[0] + 1]}}fn main() {// In order to avoid potential stack overflow, spawn a new thread.let stack_size = 104_857_600; // 100 MBlet thd = std::thread::Builder::new().stack_size(stack_size);thd.spawn(|| solve()).unwrap().join().unwrap();}// Tags: rerootingfn solve() {input! {n: usize,uv: [(usize1, usize1); n - 1],}let mut g = vec![vec![]; n];for i in 0..n - 1 {let (u, v) = uv[i];g[u].push((v, i));g[v].push((u, i + n));}let rr = MonoidalReroot::<IndepMonoid>::new(&g);let mut ans = rr.query(0)[1];for i in 1..n {ans = std::cmp::min(ans, rr.query(i)[1]);}println!("{}", ans);}