Range Sum Query 2D - Mutable

Design a data structure that supports point updates and returns the sum of any submatrix in a 2D matrix, optimizing for many interleaved update and sumRegion operations. Expect to use a dynamic 2D structure (e.g., 2D BIT or 2D segment tree) to achieve much better than O(mn) per operation — typically logarithmic time per update/query.