gh-35757: Implementing the type B Ish arrangement.
    
<!-- Please provide a concise, informative and self-explanatory title.
-->
<!-- Don't put issue numbers in the title. Put it in the Description
below. -->
<!-- For example, instead of "Fixes #12345", use "Add a new method to
multiply two integers" -->

We implement the type B Ish arrangement following a recent preprint by
Tran and Tsujie arXiv: 2304.12022.

### 📚 Description

<!-- Describe your changes here in detail. -->
<!-- Why is this change required? What problem does it solve? -->
<!-- If this PR resolves an open issue, please link to it here. For
example "Fixes #12345". -->
<!-- If your change requires a documentation PR, please link it
appropriately. -->

### 📝 Checklist

<!-- Put an `x` in all the boxes that apply. It should be `[x]` not `[x
]`. -->

- [x] The title is concise, informative, and self-explanatory.
- [x] The description explains in detail what this PR is about.
- [x] I have linked a relevant issue or discussion.
- [x] I have created tests covering the changes.
- [x] I have updated the documentation accordingly.

### ⌛ Dependencies

<!-- List all open PRs that this PR logically depends on
- #12345: short description why this is a dependency
- #34567: ...
-->

<!-- If you're unsure about any of these, don't hesitate to ask. We're
here to help! -->
    
URL: #35757
Reported by: Travis Scrimshaw
Reviewer(s): Frédéric Chapoton

Author: Release Manager

src/doc/en/reference/references/index.rst

@@ -6120,6 +6120,10 @@ REFERENCES:
 .. [TOPCOM] \J. Rambau, TOPCOM
             <http://www.rambau.wm.uni-bayreuth.de/TOPCOM/>.
 
+.. [TT2023] Tan Nhat Tran and Shuhei Tsujie.
+            *A type B analog of Ish arrangement*.
+            Preprint, :arxiv:`2304.12022`. (2023)
+
 .. [TTWL2009] Trebst, Troyer, Wang and Ludwig, A short introduction to
               Fibonacci anyon models, :arxiv:`0902.3275`.
 

src/sage/geometry/hyperplane_arrangement/library.py

@@ -420,6 +420,85 @@ def Ish(self, n, K=QQ, names=None):
         A.characteristic_polynomial.set_cache(charpoly)
         return A
 
+    def IshB(self, n, K=QQ, names=None):
+        r"""
+        Return the type B Ish arrangement.
+
+        INPUT:
+
+        - ``n`` -- integer
+        - ``K`` -- field (default:``QQ``)
+        - ``names`` -- tuple of strings or ``None`` (default); the
+          variable names for the ambient space
+
+        OUTPUT:
+
+        The type `B` Ish arrangement, which is the set of `2n^2` hyperplanes
+
+        .. MATH::
+
+            \{ x_i \pm x_j = 0 : 1 \leq i < j \leq n \}
+            \cup
+            \{ x_i = a : 1 \leq i\leq n, \quad i - n \leq a \leq n - i + 1 \}.
+
+        EXAMPLES::
+
+            sage: a = hyperplane_arrangements.IshB(2)
+            sage: a
+            Arrangement of 8 hyperplanes of dimension 2 and rank 2
+            sage: a.hyperplanes()
+            (Hyperplane 0*t0 + t1 - 1,
+             Hyperplane 0*t0 + t1 + 0,
+             Hyperplane t0 - t1 + 0,
+             Hyperplane t0 + 0*t1 - 2,
+             Hyperplane t0 + 0*t1 - 1,
+             Hyperplane t0 + 0*t1 + 0,
+             Hyperplane t0 + 0*t1 + 1,
+             Hyperplane t0 + t1 + 0)
+            sage: a.cone().is_free()
+            True
+
+        .. PLOT::
+            :width: 500 px
+
+            a = hyperplane_arrangements.IshB(2)
+            sphinx_plot(a.plot())
+
+        ::
+
+            sage: a = hyperplane_arrangements.IshB(3); a
+            Arrangement of 18 hyperplanes of dimension 3 and rank 3
+            sage: a.characteristic_polynomial()
+            x^3 - 18*x^2 + 108*x - 216
+            sage: b = hyperplane_arrangements.Shi(['B', 3])
+            sage: b.characteristic_polynomial()
+            x^3 - 18*x^2 + 108*x - 216
+
+        TESTS::
+
+            sage: a.characteristic_polynomial.clear_cache()
+            sage: a.characteristic_polynomial()
+            x^3 - 18*x^2 + 108*x - 216
+
+        REFERENCES:
+
+        - [TT2023]_
+        """
+        H = make_parent(K, n, names)
+        x = H.gens()
+        hyperplanes = []
+        for i in range(n):
+            for j in range(i+1, n):
+                hyperplanes.append(x[i] - x[j])
+                hyperplanes.append(x[i] + x[j])
+            for a in range(i+1-n, n-i+1):
+                hyperplanes.append(x[i] - a)
+        A = H(*hyperplanes)
+        x = polygen(QQ, 'x')
+        charpoly = (x - 2*n) ** n
+        A.characteristic_polynomial.set_cache(charpoly)
+        return A
+
     def linial(self, n, K=QQ, names=None):
         r"""
         Return the linial hyperplane arrangement.