diff --git a/src/data_structures/fenwick.md b/src/data_structures/fenwick.md
index 82755b452..56220ad61 100644
--- a/src/data_structures/fenwick.md
+++ b/src/data_structures/fenwick.md
@@ -300,6 +300,18 @@ $$h(i) = i + (i ~\&~ (-i)).$$
As you can see, the main benefit of this approach is that the binary operations complement each other very nicely.
+**Note**: for improved portability, the extraction of the last set bit can alternatively be computed as $i - (i ~\&~ (i-1))$. Accordingly:
+
+$$g(i) = i ~\&~ (i-1)$$
+
+and
+
+$$h(i) = 2i - (i ~\&~ (i-1)).$$
+
+This also allows the code for $g(i)$ to be succinctly written as `i &= i - 1`.
+
+An example of this lack of portability is any version of C++ prior to C++20. Before [C++20](https://en.cppreference.com/w/cpp/20), it was not guaranteed for signed integers to be implemented as [two's complement](https://en.wikipedia.org/wiki/Two%27s_complement), so the use of negative numbers ($-i$) alongside bitwise-and ($~\&~$) could ([theoretically](https://en.cppreference.com/w/cpp/language/types#Range_of_values)) lead to unexpected results.
+
The following implementation can be used like the other implementations, however it uses one-based indexing internally.
```{.cpp file=fenwick_sum_onebased}
@@ -335,6 +347,21 @@ struct FenwickTreeOneBasedIndexing {
}
};
```
+An equivalent implementation which utilizes the more portable formulas would have the following changes:
+
+```cpp
+ int sum(int idx) {
+ int ret = 0;
+ for (++idx; idx > 0; idx &= idx - 1)
+ ret += bit[idx];
+ return ret;
+ }
+
+ void add(int idx, int delta) {
+ for (++idx; idx < n; idx += idx - (idx & (idx - 1))
+ bit[idx] += delta;
+ }
+```
## Range operations
pFad - Phonifier reborn
Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies:
Alternative Proxy
pFad Proxy
pFad v3 Proxy
pFad v4 Proxy