1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
//! 削除可能ヒープ。

use std::collections::BinaryHeap;

/// 削除可能ヒープ。
///
/// # Idea
/// ヒープを二つ使う。
///
/// # Implementation notes
/// `peek()` を `&self` にするため、`remove()` と `pop()` の直後はとりあえず
/// `force()` している。
///
/// # Examples
/// ```
/// use nekolib::ds::RemovableHeap;
///
/// let mut heap = RemovableHeap::new();
/// heap.push(2);
/// heap.push(1);
/// heap.push(3);
/// heap.push(0);
/// assert_eq!(heap.peek(), Some(&3)); // {0, 1, 2, 3}
///
/// heap.remove(2);                    // {0, 1, 3}
/// assert_eq!(heap.pop(), Some(3));   // {0, 1}
/// assert_eq!(heap.peek(), Some(&1));
/// heap.remove(1);                    // {0}
/// assert_eq!(heap.peek(), Some(&0));
/// assert_eq!(heap.pop(), Some(0));   // {}
/// assert!(heap.is_empty());
/// ```
///
/// 削除する要素がヒープに入っていないとき、こまる。
/// ```
/// use nekolib::ds::RemovableHeap;
///
/// let mut heap = RemovableHeap::new();
/// heap.push(1);
/// heap.remove(2);
/// heap.push(2);
/// heap.push(3);                      // {1, 2, 3}?
/// assert_eq!(heap.pop(), Some(3));   // seems good?
/// assert_ne!(heap.peek(), Some(&2)); // but...
/// assert_eq!(heap.peek(), Some(&1));
/// ```
///
/// 未来の push 操作を打ち消すというわけでもない。
/// ```
/// use nekolib::ds::RemovableHeap;
///
/// let mut heap = RemovableHeap::new();
/// heap.push(1);
/// heap.push(2);
/// heap.remove(100);
/// heap.remove(3);
/// heap.push(101);
/// heap.push(3);                      // {1, 2, 101}?
/// assert_eq!(heap.pop(), Some(101)); // seems good, so far
/// assert_ne!(heap.peek(), Some(&2)); // but...
/// assert_eq!(heap.peek(), Some(&3));
/// ```
#[derive(Clone)]
pub struct RemovableHeap<T> {
    alive: BinaryHeap<T>,
    dead: BinaryHeap<T>,
    len: usize,
}

impl<T: Ord> RemovableHeap<T> {
    /// 空のヒープで初期化する。
    pub fn new() -> Self {
        Self {
            alive: BinaryHeap::new(),
            dead: BinaryHeap::new(),
            len: 0,
        }
    }

    /// 要素を追加する。
    pub fn push(&mut self, item: T) {
        self.len += 1;
        self.alive.push(item);
    }

    /// 要素を削除する。
    ///
    /// # Notes
    /// 削除処理を遅延させて行うため、`&T` ではなく `T` であることに注意。
    ///
    /// また、`item` はその時点でヒープに入っている要素のいずれかと等しい必要がある。
    /// 「未来に追加される要素を打ち消す」あるいは「単に無視する」という仕様は、
    /// どちらも操作順しだいでどちらもうまくいかなくなるはず。
    ///
    /// `BTreeSet` を使えば判定できるが、それなら最初から `BTreeSet` でやればよい。
    pub fn remove(&mut self, item: T) {
        self.len -= 1;
        self.dead.push(item);
        self.force();
    }

    /// 最大値を取り出す。
    pub fn pop(&mut self) -> Option<T> {
        self.force();
        let res = self.alive.pop();
        if res.is_some() {
            self.len -= 1;
            self.force();
        }
        res
    }

    /// 最大値を取得する。
    pub fn peek(&self) -> Option<&T> { self.alive.peek() }

    /// 空のとき `true` を返す。
    pub fn is_empty(&self) -> bool { self.len == 0 }

    /// 要素数を返す。
    pub fn len(&self) -> usize { self.len }

    fn force(&mut self) {
        while let (Some(x), Some(y)) = (self.alive.peek(), self.dead.peek()) {
            if x != y {
                break;
            }
            self.alive.pop();
            self.dead.pop();
        }
    }
}

#[test]
fn test() {
    let mut heap = RemovableHeap::new();
    assert_eq!(heap.pop(), None);
    heap.push(1); // {1}
    assert_eq!(heap.peek(), Some(&1));
    assert_eq!(heap.len(), 1);
    heap.push(2); // {1, 2}
    assert_eq!(heap.peek(), Some(&2));
    assert_eq!(heap.len(), 2);
    heap.remove(1); // {2}
    assert_eq!(heap.len(), 1);
    assert_eq!(heap.pop(), Some(2)); // {}
    assert_eq!(heap.len(), 0);
    assert_eq!(heap.pop(), None); // {}
    assert_eq!(heap.len(), 0);

    heap.push(2); // {2}
    heap.push(1); // {1, 2}
    assert_eq!(heap.len(), 2);
    heap.remove(2); // {1}
    assert_eq!(heap.len(), 1);
    assert_eq!(heap.pop(), Some(1)); // {}
    assert_eq!(heap.len(), 0);
}