1 Introduction

Here, a tensor (or hypermatrix) refers to a multi-array of entries in some field, which can be viewed to be the coordinates of the classical tensor (as multilinear function) under an orthonormal basis. The eigenvalues of a tensor were introduced by Qi [17, 18] and Lim [14] independently. To find the eigenvalues of a tensor, Qi [17, 18] introduced the characteristic polynomial of a tensor, which is defined to be a resultant of a system of homogeneous polynomials. In general, there is not an explicit polynomial formula yet for resultants except in some very special cases; and many fundamental questions about resultants still remain open.

To our knowledge, there are currently three tools to compute a concrete resultant. The first one is Koszul complex, whose terms are given by the graded tensor product of a polynomial algebra and an exterior algebra, and the differential is built from objective polynomials in the resultant. The resultant is exactly equal to a certain characteristic of the related Koszul complex, see [1, 16]. The second one is the generalized trace, which is defined by Morozov and Shakirov [15]. Using the generalized traces and the Schur function, Cooper and Dutle [5] expressed the characteristic polynomial of a hypergraph and characterized several properties of the characteristic polynomial, including the characteristic polynomial for a single edge. Hu et al. [12] gave an expression of the characteristic polynomial of a tensor. Shao, Qi and Hu [19] gave a graph theoretic formula for the generalized trace. The third tool is the Poisson formula, which may provide an inductive computing method, see [11, Chapter 13, Theorem 1.2] or [13, Proposition 2.7]. For example, Cooper and Dutle computed the spectrum of the “all ones” tensors using the Poisson formula, see [6, Theorem 3]. We refer to [11, Chapter 13] and [7, Chapter 3] for an overview of calculation of resultants.

Recently, spectral hypergraph theory has been developed to explore connections between the structure of a uniform hypergraph and the eigenvalues of some related symmetric tensors. Cooper and Dutle [5] proposed the concept of the adjacency tensor for a uniform hypergraph. Shao et al. [20] proved that the adjacency tensor of a connected k-uniform hypergraph G has a symmetric H-spectrum if and only if k is even and G is odd-0bipartite. This result gives a method to check whether a connected even-uniform hypergraph is odd-bipartite or not. Fan et al. [9, 10] use the spectral symmetry or the eigenvariety associated with spectral radius of the adjacency tensor to characterize the structural properties of hypergraphs.

The characteristic polynomial of a hypergraph is defined to be the characteristic polynomial of its adjacency tensor. In this paper, we mainly aim to give a lower dimension formula to compute the characteristic polynomial of hypergraphs based on the Poisson formula and variants of the chip-firing game, and determine the characteristic polynomials of starlike hypergraphs including hyperstars and hyperpaths, recursively and explicitly.

For simplicity of notation, we denote \([n]=\{1, 2, \ldots , n\}\) and \([m, n]=\{m, m+1, \ldots , n\}\) for integers \(m<n\).

2 Preliminaries

In this section, we recall some basic notions and useful results on resultants and hypergraphs.

2.1 Resultants

Let \(F_1(x_1, \ldots , x_n), \ldots , F_n(x_1, \ldots , x_n)\) be n homogeneous polynomials over \(\mathbb {C}\) in variables \(x_1, \ldots , x_n\), where the degree of \(F_i\) is \(d_i>0\) for \(i \in [n]\). An important question is whether the system of equations

$$\begin{aligned} \left\{ \begin{array}{l} F_1(x_1, \ldots , x_n)=0,\\ \cdots \\ F_n(x_1, \ldots , x_n)=0 \end{array} \right. \end{aligned}$$
(2.1)

admits nontrivial solutions.

Generally, each \(F_i\) can be written as

$$\begin{aligned} F_i=\sum \limits _{|\alpha |=d_i} c_{i,\alpha }\mathbf{x}^\alpha , \end{aligned}$$

where \(\alpha =(i_1, \ldots , i_n)\), \(|\alpha |=i_1+\cdots +i_n\) and \(\mathbf{x}^\alpha =x_1^{i_1}\cdots x_n^{i_n}\). Note that the number of \(\alpha \)’s with \(|\alpha |=d\) is \(n+d-1 \atopwithdelims ()d\).

For each possible pair of indices \(i, \alpha \), we introduce a variable \(u_{i, \alpha }\). Then, given a polynomial \(P\in \mathbb {C}[u_{i, \alpha }: |\alpha |=d_i, i\in [n]]\), we let \(P(F_1, \ldots , F_n)\) denote the value obtained by replacing each variable \(u_{i, \alpha }\) in P with the corresponding coefficient \(c_{i, \alpha }\).

Theorem 2.1

[7, Chapter 3, Theorem 2.3] For fixed positive degrees \(d_1, \ldots , d_n\), there exists a unique polynomial \( \mathrm{Res} \in \mathbb {Z}[u_{i,\alpha }]\) satisfying the following properties:

  1. (i)

    If \(F_1, \ldots , F_n\in \mathbb {C}[x_1, \ldots , x_n]\) are homogeneous of degrees \(d_1, \ldots , d_n\), respectively, the system (2.1) has a nontrivial solution if and only if \(\mathrm{Res}(F_1, \ldots , F_n)=0\).

  2. (ii)

    \( \mathrm{Res}(x_1^{d_1}, \ldots , x_n^{d_n})=1\).

  3. (iii)

    \(\mathrm{Res}\) is irreducible, even regarded as a polynomial in \(\mathbb {C}[u_{i,\alpha }]\).

\(\mathrm{Res}(F_1, \ldots , F_n)\) is called the resultant of \(F_1, \ldots , F_n\). Resultants have important applications in algebraic geometry, algebraic combinatorics and spectral hypergraph theory. However, it is difficult to compute the resultant of general polynomials. Here, we list some useful properties and calculation methods of resultants which will be used in this paper.

Lemma 2.2

[5, Lemma 3.2] Let \(F_1, \ldots , F_n\in \mathbb {C}[x_1, \ldots , x_n]\) be homogeneous polynomials of degree \(d_1, \ldots , d_n\), respectively, and let \(G_1, \cdots , G_m\in \mathbb {C}[y_1, \ldots , y_m]\) be homogeneous polynomials of degree \(\delta _1, \ldots , \delta _m\), respectively. Then,

$$\begin{aligned} \mathrm {Res}(F_1, \ldots , F_n, G_1, \ldots , G_m)=\mathrm {Res}(F_1, \ldots , F_n)^{\prod \limits _{j=1}^m\delta _j} \mathrm {Res}(G_1, \ldots , G_m)^{\prod \limits _{i=1}^nd_i}. \end{aligned}$$

Lemma 2.3

[7, Chapter 3, Theorem 3.1] For a fixed \(j \in [n]\) and \(\lambda \in \mathbb {C}\),

$$\begin{aligned} \mathrm {Res}(F_1, \ldots , F_{j-1}, \lambda F_j, F_{j+1}, \ldots , F_n)=\lambda ^{d_1\cdots d_{j-1}d_{j+1} \cdots d_n} \mathrm {Res}(F_1, \ldots , F_n), \end{aligned}$$

where \(d_i\) is the degree of \(F_i\) for each \(i \in [n]\).

Next, we recall the Poisson formula. Given homogeneous polynomials \(F_1, \ldots , F_n\)\(\in \)\(\mathbb {C}[x_1, \ldots , x_n]\) of degree \(d_1, \ldots , d_n\), respectively, let

$$\begin{aligned} f_i(x_1, \ldots , x_{n-1}) = F_i(x_1, \ldots , x_{n-1}, 1), \ \ (1\le i\le n),\\ {\bar{F}}_{i}(x_1, \ldots , x_{n-1}) = F_i(x_1, \ldots , x_{n-1}, 0), \ \ (1\le i\le n-1). \end{aligned}$$

Observe that \(\bar{F}_1, \ldots , \bar{F}_{n-1}\) are still homogeneous in \(\mathbb {C}[x_1, \ldots , x_{n-1}]\) of degree \(d_1, \ldots , d_{n-1}\), respectively, but \(f_1,\ldots ,f_n\) are not homogeneous in general.

Lemma 2.4

[Poisson formula] Keep the above notation. If \(\mathrm {Res}(\bar{F}_1, \ldots , \bar{F}_{n-1})\ne 0\), then the quotient algebra \(A=\dfrac{\mathbb {C}[x_1, \ldots , x_{n-1}]}{\langle \,f_1, \ldots , f_{n-1}\,\rangle }\) has dimension \(d_1\cdots d_{n-1}\) as a vector space over \(\mathbb {C}\), where \(\langle \,f_1, \ldots , f_{n-1}\,\rangle \) is the ideal of the polynomial algebra \(C[x_1, \ldots , x_{n-1}]\) generated by \(f_1, \ldots , f_{n-1}\), and

$$\begin{aligned} \mathrm {Res}(F_1, \ldots , F_{n})=\mathrm {Res}({\bar{F}}_1, \ldots , {\bar{F}}_{n-1})^{d_n} \det (m_{f_{n}}: A \rightarrow A), \end{aligned}$$

where \(m_{f_n} : A \rightarrow A\) is the multiplication map given by \(f_n\).

Here, the above form of the Poisson formula follows from [7, Chapter 3, Theorem 3.4], which is different from the original one in [13].

2.2 Hypergraphs

A hypergraphH is a pair (VE), where V is the set of vertices and \(E\subset \mathcal {P}(V)\) is the set of edges. A hypergraph H is called k-uniform for an integer \(k\ge 2\) if for each \(\varepsilon \in E\), \(|\varepsilon |=k\). Clearly, a 2-uniform hypergraph is just a classical simple graph.

Definition 2.5

[5] Let \(H=(V, E)\) be a k-uniform hypergraph. The (normalized) adjacency tensor\(\mathcal {A}(H)=(a_{i_1\cdots i_k})_{i_1, \cdots , i_k\in V}\) is defined by

$$\begin{aligned} a_{i_1\cdots i_k}= \left\{ \begin{array}{ll} \dfrac{1}{(k-1)!}, &{} \hbox {if~} \{i_1,\ldots , i_k\}\in E, \\ 0, &{} \hbox {otherwise}. \end{array} \right. \end{aligned}$$

For convenience, we use the following notation. Let V be a finite set and m a positive integer. For each \(\mathbf e =(i_1, \ldots , i_m)\in V^{m}\) and \(\mathbf c =(c_1, \ldots , c_m)\in \mathbb {N}^m\), we denote \(\mathbf{x}_\mathbf{e }^\mathbf{c }=x_{i_1}^{c_1}\cdots x_{i_m}^{c_m}\). We also write \(\mathbf{x}_\mathbf{e }^{\mathbb {1}}\) as \(\mathbf{x}_\mathbf e \), where \(\mathbb {1}=(1, \ldots , 1)\in \mathbb {N}^m\). If \(V=[n]\), \(\mathbf c =(c_1, \ldots , c_n) \in \mathbb {N}^n\), we write \(\mathbf{x}_{[n]}^\mathbf{c }\) as \(\mathbf{x}^\mathbf{c }\).

The eigenvalues of a tensor were introduced by Qi [17, 18] and Lim [14] independently. The adjacency tensor of a uniform hypergraph was introduced by Cooper and Dutle [5]. Here, we briefly give the definition of eigenvalues of uniform hypergraphs based on the above.

Definition 2.6

[5, 17] Let \(H=(V, E)\) be a k-uniform hypergraph and \(\mathcal {A}=(a_{i_1\cdots i_k})\) be the adjacency tensor of H. For some \(\lambda \in \mathbb {C}\), if there exists a nonzero vector \(\mathbf{x}\in \mathbb {C}^{|V|}\) such that for each \(j \in V\),

$$\begin{aligned} \sum \limits _{i_2,i_3,\ldots ,i_k \in V} a_{ji_2 i_3 \ldots i_k } x_{i_2}x_{i_3} \cdots x_{i_k}=\lambda x_j^{k-1}, \end{aligned}$$

or equivalently, for each \(v\in V\),

$$\begin{aligned} \sum \limits _{v\in \varepsilon \in E} \mathbf{x}_{\varepsilon \backslash \{v\}}=\lambda x_v^{k-1}, \end{aligned}$$

then \(\lambda \) is called an eigenvalue of H.

For each \(v\in V\), define

$$\begin{aligned} F_v=\lambda x_v^{k-1}-\sum \limits _{v\in \varepsilon \in E} \mathbf{x}_{\varepsilon \backslash \{v\}}. \end{aligned}$$

The polynomial

$$\begin{aligned} \phi _H(\lambda )=\mathrm {Res}(F_v: v\in V) \end{aligned}$$

in the indeterminant \(\lambda \) is called the characteristic polynomial of H. Consequently, \(\lambda \) is an eigenvalue of H if and only if \(\phi _H(\lambda )=0\).

2.3 The dollar game on graphs

Let \(G=(V,E)\) be a simple graph. Recall that a configuration\({\mathbf {c}}\) on G means a function \({\mathbf {c}}: V \rightarrow \mathbb {N}\), which can be understood as a pile of \({\mathbf {c}}(v)\) tokens (chips, or dollars) at each vertex v. A dollar game on G starts from a configuration \({\mathbf {c}}\). At each step of the game, a vertex v is fired, that is, dollars move from v to its adjacent vertices, one dollar going along each edge incident to v. Fix a vertex w of G, called the bank vertex. A vertex v other than w can be fired if and only if \({\mathbf {c}}(v)\ge \deg (v)\), where \(\deg (v)\) is the degree of the vertex v. The bank vertex w is allowed to go into debt such that w can be fired if and only if no other firing is possible.

Suppose that \(\mathcal {X}\) is a non-empty finite sequence of (not necessarily distinct) vertices of G, such that starting from a configuration \({\mathbf {c}}\), the vertices can be fired in the order of \(\mathcal {X}\). If v occurs x(v) times, we shall refer to x as the representative vector for \(\mathcal {X}\). The configuration \({\mathbf {c}}'\) after the sequence of firing \(\mathcal {X}\) is given by

$$\begin{aligned} {\mathbf {c}}'={\mathbf {c}}-Lx, \end{aligned}$$

where L is the Laplacian matrix of G.

The dollar game on graphs was introduced by Biggs [2] and is a variant of the chip-firing game. It is often described in terms of “snowfalls” and “avalanches” in the literature. A configuration \({\mathbf {c}}\) is said to be stable if \(0\le {\mathbf {c}}(v)<\deg (v)\) for every \(v\ne w\). A sequence of firing is w-legal if and only if each occurrence of a vertex \(v\ne w\) follows a configuration t with \(t(v)\ge \deg (v)\) and each occurrence of w follows a stable configuration. A configuration \({\mathbf {c}}\) on G is said recurrent if there is a w-legal sequence for \({\mathbf {c}}\) which leads to the same configuration. A critical configuration \({\mathbf {c}}\) means that \({\mathbf {c}}\) is both stable and recurrent. We refer to [2] for more details.

Lemma 2.7

[2, Theorem 6.2] If G is a connected graph, then the number of critical configurations is equal to the number of spanning trees of G.

Example 2.8

Let \(K_k\) be a complete graph on k vertices. Then, the number of critical configurations is \(k^{k-2}\).

3 Poisson formula for characteristic polynomials of hypergraphs

3.1 Poisson formula for hypergraphs

Let \(H=(V, E)\) be a k-uniform hypergraph. Recall that the characteristic polynomial of H is defined as

$$\begin{aligned} \phi _H(\lambda )=\mathrm {Res}(F_v: v\in V),\end{aligned}$$

where

$$\begin{aligned} F_v=\lambda x_v^{k-1}-\sum \limits _{v\in \varepsilon \in E} \mathbf{x}_{\varepsilon \backslash \{v\}}\in \mathbb {C}[x_v: v\in V]. \end{aligned}$$

In order to use Poisson formula for the resultant \(\mathrm {Res}(F_v: v\in V)\), we need fix a vertex w in V. Denote by \(E_w\) the set of all edges containing the vertex w and \(\varepsilon _{\widehat{w}}=\varepsilon \backslash \{w\}\) for each \(\varepsilon \in E_w\). Then, we have

$$ f_w=\lambda -\sum \limits _{\varepsilon \in E_w} \mathbf{x}_{\varepsilon _{\widehat{w}}},$$

and for all \(v\in V\backslash \{w\}\),

$$\begin{aligned}&f_v= \lambda x_v^{k-1}-\sum \limits _{v\in \varepsilon \in E\backslash E_w}\mathbf{x}_{\varepsilon \backslash \{v\}}-\sum \limits _{v\in \varepsilon \in E_w} \mathbf{x}_{\varepsilon _{\widehat{w}}\backslash \{v\}}, \nonumber \nonumber \\&\bar{F}_v= \lambda x_v^{k-1}-\sum \limits _{v\in \varepsilon \in E\backslash E_w}\mathbf{x}_{\varepsilon \backslash \{v\}}, \end{aligned}$$
(3.1)

where \(f_w\) and \(f_v\) are obtained from \(F_w\) and \(F_v\) by taking \(x_w=1\), respectively, and \(\bar{F}_v\) is obtained from \(F_v\) by taking \(x_w=0\).

Deleting the vertex w in V and the edges in \(E_w\), one can obtain a sub-hypergraph \(\widehat{H}=(\widehat{V}, \widehat{E})\). To be precise, \(\widehat{V}=V\backslash \{w\}\) and \(\widehat{E}=E\backslash E_w\).

Lemma 3.1

Retain the above notation. Then,

$$\begin{aligned} \phi _H(\lambda )=\phi _{\widehat{H}}(\lambda )^{k-1}\det (m_{f_w}: A \rightarrow A), \end{aligned}$$
(3.2)

where A is the quotient algebra \(\dfrac{\mathbb {C}[x_v: v\in \widehat{V}]}{\langle \,f_v: v\in \widehat{V}\,\rangle }\) and \(m_{f_w}\) is the multiplication map of A given by \(f_w\).

Proof

By Lemma 2.4, the characteristic polynomial of H is

$$\begin{aligned} \phi _H(\lambda )=\mathrm {Res}(\bar{F}_v: v\in \widehat{V})^{k-1}\det (m_{f_w}: A \rightarrow A).\end{aligned}$$

Considering the sub-hypergraph \(\widehat{H}\) of H, by Eq. (3.1), we have

$$\begin{aligned} \phi _{\widehat{H}}(\lambda )=\mathrm {Res}(\bar{F}_v: v\in \widehat{V}).\end{aligned}$$

The result follows. \(\square \)

By definition, the algebra A is \((k-1)^{r-1}\)-dimensional as a vector space over \(\mathbb {C}\) where r is the number of vertices of H. In general, it is difficult to compute the determinant \(\det (m_{f_w}: A \rightarrow A)\). However, we can give a description for some special cases.

3.2 Hypergraphs with a cut vertex

Let \(H=(V, E)\) be a k-uniform connected hypergraph and \(w\in V\). Denote \(E_{\widehat{w}}=\{\varepsilon _{\widehat{w}}\mid \varepsilon \in E_w\}\). Deleting the vertex w, we can get a (non-uniform) hypergraph \(\widetilde{H}=(\widetilde{V}, \widetilde{E})\), with \(\widetilde{V}=\widehat{V}=V\backslash \{w\}\) and \(\widetilde{E}=(E\backslash E_w)\cup E_{\widehat{w}}\). Recall the vertex w is called a cut vertex if \(\widetilde{H}\) is not connected; see Fig. 1. Suppose that w is a cut vertex and \(\widetilde{H}_1=(\widetilde{V}_1, \widetilde{E}_1), \cdots , \widetilde{H}_n=(\widetilde{V}_n, \widetilde{E}_n)\) (\(n\ge 2\)) are the connected components of \(\widetilde{H}\). For each \(i\in [n]\), we set \(V_i=\widetilde{V}_i\), \(E_i=\widetilde{E}_i\backslash E_{\widehat{w}}\), and then obtain a sub-hypergraph \(H_i=(V_i, E_i)\) of H. Note that each \(H_i\) is a k-uniform hypergraph and may not be connected.

Fig. 1
figure 1

A k-uniform hypergraph H with a cut vertex w

For each \(i\in [n]\), we denote

$$\begin{aligned} E_w^i=\{\varepsilon \in E_w \mid \varepsilon \cap V_i\ne \varnothing \}. \end{aligned}$$

By definition, we have

$$\begin{aligned} \phi _H(\lambda )=\mathrm {Res}(F_v: v\in V), \end{aligned}$$

where

$$\begin{aligned} F_w= \lambda x_w^{k-1}-\sum \limits _{i=1}^n\sum \limits _{\varepsilon \in E_w^i} \mathbf{x}_{\varepsilon _{\widehat{w}}}, \end{aligned}$$

and for all \(v_i\in V_i\), \(i \in [n]\),

$$\begin{aligned} F_{v_i}=\lambda x_{v_i}^{k-1}-\sum \limits _{v_i\in \varepsilon \in E_i} \mathbf{x}_{\varepsilon \backslash \{v_i\}}-\sum \limits _{v_i\in \varepsilon \in E_w^i} \mathbf{x}_{ \varepsilon \backslash \{v_i\}}, \end{aligned}$$

Therefore,

$$\begin{aligned} f_w= \lambda -\sum \limits _{i=1}^n\sum \limits _{\varepsilon \in E_w^i} \mathbf{x}_{\varepsilon _{\widehat{w}}}, \end{aligned}$$

and for all \(v_i\in V_i\), \(i \in [n]\),

$$\begin{aligned}&f_{v_i} = \lambda x_{v_i}^{k-1}-\sum \limits _{v_i\in \varepsilon \in E_i} \mathbf{x}_{\varepsilon \backslash \{v_i\}}-\sum \limits _{v_i\in \varepsilon \in E_w^i} \mathbf{x}_{\varepsilon _{\widehat{w}} \backslash \{v_i\}}, \\&\bar{F}_{v_i} = \lambda x_{v_i}^{k-1}-\sum \limits _{v_i\in \varepsilon \in E_i} \mathbf{x}_{\varepsilon \backslash \{v_i\}}. \end{aligned}$$

Let A be the quotient algebra \(\dfrac{\mathbb {C}[x_v: v\in \widehat{V}]}{\langle \,f_v: v\in \widehat{V}\,\rangle }\), and let \(m_{f_w}: A \rightarrow A\) be the multiplication map given by \(f_w\). Since \(V_1, \ldots , V_n\) form a partition of \(\widehat{V}\) and for each \(v_i\in V_i, f_{v_i}\in \mathbb {C}[x_v: v\in V_i]\), each monomial of A can be expressed as \(\mathbf{x}_{\widehat{V}}^{{\mathbf {c}}}=\mathbf{x}_{V_1}^{{\mathbf {c}}_1}\mathbf{x}_{V_2}^{{\mathbf {c}}_2} \cdots \mathbf{x}_{V_n}^{{\mathbf {c}}_n}\) by considering \(f_v=0\) inside A for \(v \in V_i\) and \(i \in [n]\). So we have \(A=A_1 \otimes \cdots \otimes A_n\), where

$$\begin{aligned} A_i= \dfrac{\mathbb {C}[x_v: v\in V_i]}{\langle \,f_v: v\in V_i\,\rangle }, \; i \in [n]. \end{aligned}$$

We define and denote \(m_{i, w}: A_i \rightarrow A_i\) the linear map given by the multiplication by \(\sum \limits \nolimits _{\varepsilon \in E_w^i} \mathbf{x}_{\varepsilon _{\widehat{w}}}\) for each \(i\in [n]\). Then,

$$\begin{aligned} m_{f_w}=\lambda \mathrm{Id}_{A} -\sum \limits _{i=1}^n \mathrm{Id}_{A_1}\otimes \cdots \otimes \mathrm{Id}_{A_{i-1}}\otimes m_{i,w} \otimes \mathrm{Id}_{A_{i+1}}\otimes \cdots \otimes \mathrm{Id}_{A_n}, \end{aligned}$$

where \(\mathrm{Id}\) denotes the identity map on certain vector space.

Assumption 1

For each \(i\in [n]\), there exists an ordered \(\mathbb {C}\)-basis \(\mathbf{x}^{\alpha _{i, 1}}, \ldots , \mathbf{x}^{\alpha _{i, d_i}}\) for \(A_i\) such that the matrix of \(m_{i, w}\) with respect to this basis is a lower triangular matrix with the diagonal entry \(\alpha _{i, j_i}\), \( j_i \in [d_i]\), where \(d_i=(k-1)^{r_i}\) and \(r_i=|V_i|\).

Under Assumption 1, \(\{\mathbf{x}^{\alpha _{1, j_1}}\cdots \mathbf{x}^{\alpha _{n, j_n}}\mid j_i\in [d_i], i\in [n]\}\) with the left lexicographic order is a basis for A such that the matrix of \(m_{f_w}\) with respect to this basis is still a lower triangular matrix with diagonal entries \(\lambda -\sum \limits _{1\le i\le n}\alpha _{i, j_i}\) for \(j_i \in [d_i]\) and \(i \in [n]\). In this situation, we have

$$\begin{aligned} \det (m_{f_w}: A \rightarrow A)={\mathop {\mathop {\prod }\limits _{1\le j_i\le d_i}}\limits _{1\le i\le n}} \left( \lambda -\sum \limits _{i=1}^n \alpha _{i, j_i} \right) . \end{aligned}$$
(3.3)

Corollary 3.2

Let H be a k-uniform hypergraph with a cut vertex w. Then, under Assumption 1,

$$\begin{aligned} \phi _H(\lambda )=\prod \limits _{i=1}^n\phi _{H_i}(\lambda )^{(k-1)^{1+\sum \limits _{j\ne i}r_j}} {\mathop {\mathop {\prod }\limits _{1\le j_i\le d_i}}\limits _{1\le i\le n}} \left( \lambda -\sum \limits _{i=1}^n \alpha _{i, j_i} \right) . \end{aligned}$$

Proof

By definition, the characteristic polynomial of \(H_i=(V_i, E_i)\) is

$$\begin{aligned} \phi _{H_i}(\lambda )=\mathrm {Res}(\bar{F}_{v_i}: v_i\in V_i). \end{aligned}$$

By Lemma 2.2, we have

$$\begin{aligned} \mathrm {Res}(\bar{F}_v: v\in \widehat{V})= \prod \limits _{i=1}^n \mathrm {Res}(\bar{F}_{v_i}: v_i\in V_i)^{(k-1)^{\sum \limits _{j\ne i} |V_j|}} = \prod \limits _{i=1}^n \phi _{H_i}(\lambda )^{(k-1)^{\sum \limits _{j\ne i} r_j}}. \end{aligned}$$

From Lemma 3.1, it follows that

$$\begin{aligned} \phi _H(\lambda )= & {} \mathrm {Res}(\bar{F}_v: v\in \hat{V})^{k-1} \det (m_{f_w}: A \rightarrow A)\\= & {} \prod \limits _{i=1}^n\phi _{H_i}(\lambda )^{(k-1)^{1+\sum \limits _{j\ne i}r_j}}{\mathop {\mathop {\prod }\limits _{j_i \in [d_i]}}\limits _{i \in [n]}} \left( \lambda -\sum \limits _{i=1}^n \alpha _{i, j_i} \right) . \end{aligned}$$

\(\square \)

3.3 Hypergraphs with a cored vertex

Let \(H=(V, E)\) be a k-uniform hypergraph. Recall that a vertex \(w \in V\) is called a cored vertex if it is contained in only one edge; see Fig. 2. Deleting the cored vertex w and the edge \(\varepsilon _w\) containing w, one can obtain a sub-hypergraph \(\widehat{H}=(\widehat{V}, \widehat{E})\) with \(\widehat{V}=V\backslash \{w\}\) and \(\widehat{E}=E\backslash \{\varepsilon _w\}\). Then,

$$\begin{aligned}&F_w = \lambda x_w^{k-1}-\mathbf{x}_{\varepsilon _{\widehat{w}}}, \\&F_v = \lambda x_v^{k-1}-\sum \limits _{v\in \varepsilon \in \widehat{E}} \mathbf{x}_{\varepsilon \backslash \{v\}}-\sum \limits _{v\in \varepsilon _w} \mathbf{x}_{\varepsilon _w\backslash \{v\}}, \; v\ne w. \end{aligned}$$

Moreover, \(f_w = \lambda -\mathbf{x}_{\varepsilon _{\widehat{w}}}\), and for all \(v\ne w\),

$$\begin{aligned}&f_v = \lambda x_v^{k-1}-\sum \limits _{v\in \varepsilon \in \widehat{E}} \mathbf{x}_{\varepsilon \backslash \{v\}}-\sum \limits _{v\in \varepsilon _w} \mathbf{x}_{\varepsilon _{\widehat{w}}\backslash \{v\}},\\&\bar{F}_v = \lambda x_v^{k-1}-\sum \limits _{v\in \varepsilon \in \widehat{E}} \mathbf{x}_{\varepsilon \backslash \{v\}}. \end{aligned}$$
Fig. 2
figure 2

A hypergraph H with a cored vertex w

Corollary 3.3

Let H be a k-hypergraph with a cored vertex w. Retain the above notation. Then,

$$\begin{aligned} \phi _H(\lambda )=\phi _{\widehat{H}}(\lambda )^{k-1} \det (m_{f_w}: A\rightarrow A), \end{aligned}$$

where A is the quotient algebra \(\dfrac{\mathbb {C}[x_v: v\in \widehat{V}]}{\langle \,f_v: v\in \widehat{V}\,\rangle }\) and \(m_{f_w}\) is the multiplication map of A given by \(f_w\).

4 Hyperpaths

Let \(P_n^k\) be a k-uniform hyperpath with n edges or of length n, which has vertices labeled as \(0, 1, \ldots , r=n(k-1)\) from left to right as in Fig. 3, and edges \(\varepsilon _t=\{(k-1)(t-1), (k-1)(t-1)+1, \ldots , (k-1)t\}\) for \( t \in [n]\).

Fig. 3
figure 3

A k-uniform hyperpath \(P_n^k\) with n edges

In this section, we will give a recursive formula for the characteristic polynomials of hyperpaths. By the Poisson formula introduced in Lemma 3.1, it suffices to compute the related determinant. For this, we introduce the following model of the dollar game on hypergraphs.

4.1 The dollar game on hypergraphs and the firing graph

We now define a dollar game on a hypergraph, considered as a variant of the dollar game on a graph. Let \(H=(V, E)\) be a k-uniform hypergraph with a specified bank vertex w. A function \({\mathbf {c}}: V \rightarrow \mathbb {N}\) is called a configuration onH. A dollar game starts from a configuration \({\mathbf {c}}\). At each of step of the game, a vertex v is fired on an edge \(\varepsilon \) containing v, that is, the vertex v decreases \(k-1\) dollars, and each of the vertices of \(\varepsilon \) other than v increases 1 dollar, which yields a new configuration on H. Note that such v will be fired on all edges containing v, one edge each step, which yields \(\deg (v)\) new configurations. A vertex v other than w can be fired if and only if \({\mathbf {c}}(v)\ge k-1\). The bank vertex w can be fired regardless of the value \({\mathbf {c}}(w)\), but only when no other vertex can be fired. We say that a configuration \({\mathbf {c}}\) is stable if \(0 \le {\mathbf {c}}(v)<k-1\) for all \(v \ne w\). The above setting of the dollar game on hypergraphs is different from that of the dollar game on simple graphs in Sect. 2.3, but will be useful for our discussion.

Let “\(\le \)” be a total ordering on the set \(\widehat{V}=V\backslash \{w\}\). Let \({\mathbf {c}}\) be a configuration on H. The weight of \({\mathbf {c}}\) is defined and denoted to be \(\omega ({\mathbf {c}})=\sum _{v\in \widehat{V}} {\mathbf {c}}(v)\). We define the left anti-lexicographical order\(\prec \) for all configurations on H. To be precise, for any configurations \({\mathbf {c}}\) and \({\mathbf {c}}'\), \({\mathbf {c}}\prec {\mathbf {c}}' \) if and only if either \(\omega ({\mathbf {c}})<\omega ({\mathbf {c}}')\), or \(\omega ({\mathbf {c}})=\omega ({\mathbf {c}}')\), \({\mathbf {c}}(i)={\mathbf {c}}'(i)\) for any \(1\le i\le t-1\) and \({\mathbf {c}}(t)>{\mathbf {c}}'(t)\) for some t.

Based on the above discussion, we now define a directed graph, called a firing graph\(\mathcal {G}({\mathbf {c}}_0)\) associated with a stable configuration \({\mathbf {c}}_0\) on the hypergraph H, which is closely related to the dollar game on H starting from \({\mathbf {c}}_0\). Here, a vertex of \(\mathcal {G}({\mathbf {c}}_0)\) is a configuration and a directed edge will be called an arrow.

Step 1 Initially, set \(\mathcal {V}_0=\{{{\mathbf {c}}}_0\}\) and \(\mathcal {E}_0=\varnothing \).

Step 2 For each edge \(\varepsilon \) containing w, let \(\bar{{\mathbf {c}}}_{0,\varepsilon }\) be a configuration given by

$$\begin{aligned} \bar{{\mathbf {c}}}_{0,\varepsilon }(v)= \left\{ \begin{array}{ll} {{\mathbf {c}}}_0(v)+1, &{} \mathrm{if} v\in \varepsilon \backslash \{w\}, \\ {{\mathbf {c}}}_0(v), &{} \hbox {otherwise}. \end{array} \right. \end{aligned}$$

Put \(\mathcal {V}_1=\mathcal {V}_0\cup \{\bar{\mathbf {c}}_{0,\varepsilon }: w \in \varepsilon \}\) and \(\mathcal {E}_1=\{{\mathbf {c}}_0\xrightarrow {\alpha _{w, \varepsilon }}\bar{\mathbf {c}}_{0,\varepsilon }: w \in \varepsilon \}\), where the arrow \(\alpha _{w, \varepsilon }\) means that the configuration \(\bar{\mathbf {c}}_{0,\varepsilon }\) is obtained from \({\mathbf {c}}_0\) by firing the bank vertex w on the edge \(\varepsilon \).

Step 3 If all of configurations in \(\mathcal {V}_i\backslash \mathcal {V}_{i-1}\) are stable, then we define \(\mathcal {V}(\mathcal {G}({\mathbf {c}}_0))=\mathcal {V}_i\) and \(\mathcal {E}(\mathcal {G}({\mathbf {c}}_0))=\mathcal {E}_i\). Otherwise, for each non-stable configuration \({\mathbf {c}}\in \mathcal {V}_i\backslash \mathcal {V}_{i-1}\), we choose the vertex

$$u_{{\mathbf {c}}}:=\max \{v\in \widehat{V} \mid {\mathbf {c}}(v)\ge k-1\}.$$

For each edge \(\varepsilon \in E\) containing \(u_{{\mathbf {c}}}\), we define an arrow \({\mathbf {c}}\xrightarrow {\alpha _{u_{{\mathbf {c}}}, \varepsilon }} \bar{\mathbf {c}}_{\varepsilon }\), where the configuration \(\bar{\mathbf {c}}_\varepsilon \) is given by

$$\begin{aligned} \bar{\mathbf {c}}_\varepsilon (v)= \left\{ \begin{array}{ll} {\mathbf {c}}(v)-(k-1), &{} \hbox {if~}v=u_{{\mathbf {c}}}, \\ {\mathbf {c}}(v)+1, &{} \hbox {if~} v\in \varepsilon \backslash \{u_{{\mathbf {c}}}\}, \\ {\mathbf {c}}(v), &{} \hbox {otherwise}, \end{array} \right. \end{aligned}$$

and

$$\begin{aligned}&\mathcal {V}_{i+1} = \mathcal {V}_i\cup \{\bar{\mathbf {c}}_\varepsilon \mid u_{{\mathbf {c}}} \in \varepsilon \in E, {\mathbf {c}}\ \hbox {is not stable in}\ \mathcal {V}_i\backslash \mathcal {V}_{i-1}\},\\&\mathcal {E}_{i+1} = \mathcal {E}_i\cup \{{\mathbf {c}}\xrightarrow {\alpha _{u_{{\mathbf {c}}}, \varepsilon }} \bar{\mathbf {c}}_{\varepsilon } \mid u_{{\mathbf {c}}}\in \varepsilon \in E, {\mathbf {c}}\ \hbox {is not stable in }\ \mathcal {V}_i\backslash \mathcal {V}_{i-1}\}. \end{aligned}$$

Note that \(\bar{\mathbf {c}}_\varepsilon \) may have been in \(\mathcal {V}_i\). If \(\bar{\mathbf {c}}_\varepsilon \notin \mathcal {V}_i\), we say that \(\bar{\mathbf {c}}_\varepsilon \) is obtained from \({\mathbf {c}}\) by firing \(u_{{\mathbf {c}}}\) on \(\varepsilon \).

Step 3 tells us the firing rule, that is, which vertex will be fired at the next step among all non-stable vertices other than the bank vertex. We use the maximal principle, the non-stable vertex with the maximal order will be firstly fired. In addition, from the construction of \(\mathcal {G}({\mathbf {c}}_0)\), we have \(0\le \omega ({\mathbf {c}}) \le \omega ({\mathbf {c}}_0)+k-1\) for any \({\mathbf {c}}\in \mathcal {V}(\mathcal {G}({\mathbf {c}}_0))\), which implies that \(\mathcal {G}({\mathbf {c}}_0)\) is a finite directed graph.

Example 4.1

Let \(P^3_3\) be the 3-uniform hyperpath with 3 edges as in Fig. 3 by taking \(n=3\) and \(k=3\). Let 0 be the bank vertex of \(P^3_3\), and let \({\mathbf {c}}_0=(*,1,1,1,1,0,0)\). Here, the value of the bank vertex 0 is omitted in each configuration of the dollar game starting from \({\mathbf {c}}_0\). The firing graph \(\mathcal {G}({\mathbf {c}}_0)\) is drawn in Fig. 4, where each vertex within a circle means that it will be fired at the next step.

Fig. 4
figure 4

Firing graph of \(P_3^3\) associated with \({\mathbf {c}}_0\)

Generally speaking, it is difficult to obtain the firing graphs \(\mathcal {G}({\mathbf {c}}_0)\) for all stable configurations \({\mathbf {c}}_0\) on a general hypergraph. However, for some special classes of hypergraphs, e.g., hyperpaths, we can characterize the structure of \(\mathcal {G}({\mathbf {c}}_0)\).

4.2 Firing graphs of hyperpaths

In this part, we characterize the structure of the firing graph of \(P_n^k\) in Fig. 3 associated with the arbitrary fixed stable configuration, where the vertex 0 is the bank vertex. For a configuration \({\mathbf {c}}\) on \([0,n(k-1)]\), we denote \({\mathbf {c}}=({\mathbf {c}}(0),{\mathbf {c}}^1, \ldots , {\mathbf {c}}^n)\), where \({\mathbf {c}}^i\) is the restriction of \({\mathbf {c}}\) on \(\widehat{\varepsilon }_i:=\varepsilon _i \backslash \{(i-1)(k-1)\}\), i.e.,

$$\begin{aligned} {\mathbf {c}}^i=({\mathbf {c}}((i-1)(k-1)+1), \ldots , {\mathbf {c}}(i(k-1))), i \in [n]. \end{aligned}$$

Let \(\tilde{\mathbf {c}}^i\) be the restriction of \({\mathbf {c}}\) on \(\varepsilon _i\), i.e.,

$$\begin{aligned} \tilde{\mathbf {c}}^i=({\mathbf {c}}((i-1)(k-1)),{\mathbf {c}}((i-1)(k-1)+1), \ldots , {\mathbf {c}}(i(k-1))), i \in [n]. \end{aligned}$$

Then, for each \(i \in [n]\), \(\tilde{\mathbf {c}}^i\) can be considered as a stable configuration on the complete graph \(K_k\) with the vertex set \(\varepsilon _i\) and the bank vertex \((i-1)(k-1)\).

Let \(\mathcal {S}_k\) be the set of all stable configurations on the complete graph \(K_k\), and \(\mathcal {C}_k\) be the set of all critical configurations on \(K_k\). Denote by \(\mathcal {B}\) the set of all stable configurations on \(P_n^k\), and for \(s\in [0, n-1]\)

$$\begin{aligned} \mathcal {B}_s=\{({\mathbf {c}}^1, \ldots , {\mathbf {c}}^n)\mid \tilde{\mathbf {c}}^i\in \mathcal {C}_k \hbox {~for~} i \in [s], \tilde{\mathbf {c}}^{s+1}\notin \mathcal {C}_k, \hbox {~and~} \tilde{\mathbf {c}}^i\in \mathcal {S}_k \hbox {~for~} i \in [s+1, n]\}, \end{aligned}$$

and \(\mathcal {B}_n=\{({\mathbf {c}}^1, \ldots , {\mathbf {c}}^n)\mid \tilde{\mathbf {c}}^i\in \mathcal {C}_k, i\in [n]\}\). Note that the bank vertex in the graph dollar game for every edge other than the first edge inherits its stability with respect to hypergraph dollar game from the stability in the graph dollar game for the previous edge. So, \(\mathcal {B}\) is the disjoint union of \(\mathcal {B}_0, \ldots , \mathcal {B}_n\).

Lemma 4.2

Suppose that \({\mathbf {c}}_0\) is a stable configuration on \(P_n^k\), the labeled k-uniform hyperpath described in the opening of Sect. 4, with 0 as the bank vertex, and \(\mathcal {G}({\mathbf {c}}_0)\) is the firing graph of \(P_n^k\) associated with \({\mathbf {c}}_0\). Then,

  1. (i)

    for any \(t\in [n]\), if \({\mathbf {c}}(t(k-1))=k-1\), then \({\mathbf {c}}(v)<k-1\) for any \(v>t(k-1)\).

  2. (ii)

    for any \(t\in [n]\), the vertex \(t(k-1)\) is fired at \({\mathbf {c}}\in \mathcal {V}(\mathcal {G}({\mathbf {c}}_0))\) if and only if \({\mathbf {c}}(t(k-1))=k-1\).

  3. (iii)

    for any configuration \({\mathbf {c}}\in \mathcal {V}(\mathcal {G}({\mathbf {c}}_0))\), \(\omega ({\mathbf {c}})\ge \omega ({\mathbf {c}}_0)\), where the equality holds if and only if the vertices \(1, \ldots , k-1\) have been fired on \(\varepsilon _1\) exactly once along any directed path from \({\mathbf {c}}_0\) to \({\mathbf {c}}\).

Proof

(i) Assume to the contrary, there exists a configuration \({\mathbf {c}}\) such that \({\mathbf {c}}(t(k-1))=k-1\), and \({\mathbf {c}}(v)\ge k-1\) for some \(v>t(k-1)\). According to the construction of firing graphs, there exists a directed path \(P({\mathbf {c}}_0, {\mathbf {c}})\) from \({\mathbf {c}}_0\) to \({\mathbf {c}}\) (without repeated vertices on the path). We also assume that \({\mathbf {c}}\) is the first occurrence of the configuration with the above property, so that along the path \(P({\mathbf {c}}_0, {\mathbf {c}})\), except \({\mathbf {c}}\) no configurations of \(P({\mathbf {c}}_0, {\mathbf {c}})\) hold the above property.

Case 1: \(t(k-1)\) is not fired before \({\mathbf {c}}\) along the path \(P({\mathbf {c}}_0, {\mathbf {c}})\). Then, \({\mathbf {c}}(v)={\mathbf {c}}_0(v)<k-1\) for any \(v>t(k-1)\) since such v cannot gain more dollars before \(t(k-1)\) is fired, a contradiction.

Case 2: \(t(k-1)\) is fired at some configuration say \({\mathbf {c}}'_{-1}\) before \({\mathbf {c}}\) along the path \(P({\mathbf {c}}_0, {\mathbf {c}})\). We may assume \(t(k-1)\) is last fired at \({\mathbf {c}}'_{-1}\) before \({\mathbf {c}}\). Consider the directed path \(P({\mathbf {c}}_0, {\mathbf {c}})\):

Note that \({\mathbf {c}}'_{-1}(t(k-1))\ge k-1\) and \({\mathbf {c}}'_{-1}(v)<k-1\) for any \(v>t(k-1)\) by the firing rule. If \({\mathbf {c}}(t(k-1)) \ge k\), there exists a configuration \({\mathbf {c}}'\) before \({\mathbf {c}}'_{-1}\) such that \({\mathbf {c}}'(t(k-1)) = k-1\) and \({\mathbf {c}}'(v) \ge k-1\) for some \(v > t(k-1)\), a contradiction to the assumption about \({\mathbf {c}}\). So \({\mathbf {c}}'_{-1}(t(k-1))= k-1\), and \({\mathbf {c}}'_0(t(k-1))=0\).

From \({\mathbf {c}}'_0\) to \({\mathbf {c}}\), the vertex \(t(k-1)\) obtains \(k-1\) dollars by firing some vertices on \(\varepsilon _t\) or \(\varepsilon _{t+1}\). If \(t(k-1)\) obtains one dollar by firing some vertex on \(\varepsilon _t\) at configuration say \({\mathbf {c}}'_s\), then \({\mathbf {c}}'_s(v)<k-1\) for any \(v>t(k-1)\) by the firing rule. Furthermore, such v cannot gain more dollars from \({\mathbf {c}}'_s\) to \({\mathbf {c}}\) as \(t(k-1)\) would not be fired before \({\mathbf {c}}\). So, \({\mathbf {c}}(v)={\mathbf {c}}'_s(v)<k-1\) for any \(v>t(k-1)\), a contradiction.

Otherwise, \(t(k-1)\) obtains \(k-1\) dollars only by firing vertices on \(\varepsilon _{t+1}\). Then,

$$\begin{aligned} {\mathbf {c}}'_0(v)= \left\{ \begin{array}{ll} 0, &{} \hbox {~if~} v=t(k-1)\\ {\mathbf {c}}'_{-1}(v)+1 , &{} \hbox {~if~} v \in [t(k-1)+1,(t+1)(k-1)], \\ {\mathbf {c}}'_{-1}(v), &{} \hbox {~if~} v\ge (t+1)(k-1)+1. \end{array} \right. \end{aligned}$$

It follows that the number of configurations in \(\{{\mathbf {c}}'_i\mid \alpha _{u_i, *}=\alpha _{u_i, \varepsilon _{t+1}}, i=0, 1, \ldots , m\}\) is \(k-1\) as \({\mathbf {c}}(t(k-1))=k-1\). Observe that for any \(v \in [t(k-1)+1,(t+1)(k-1)-1]\),

$$\begin{aligned} {\mathbf {c}}'_i(v)\le {\mathbf {c}}'_0(v)+(k-2)=c'_{-1}(v)+(k-1)<2k-2. \end{aligned}$$

It forces that each vertex in \([t(k-1)+1,(t+1)(k-1)-1]\) is fired exactly once on the edge \(\varepsilon _{t+1}\) from \({\mathbf {c}}'_0\) and \({\mathbf {c}}\), implying that the vertex \((t+1)(k-1)\) is also fired exactly once on \(\varepsilon _{t+1}\).

Let \(u_j\) be the vertex in \([t(k-1)+1,(t+1)(k-1)]\) last fired on \(\varepsilon _{t+1}\). Then, \({\mathbf {c}}'_j(v) < k-1\) for any \(v >(t+1)(k-1)\) by the firing rule. Also, we have \({\mathbf {c}}'_{{j+1}}(t(k-1))=k-1\), \({\mathbf {c}}'_{{j+1}}(v) < k-1\) for any \(v \in [t(k-1)+1,(t+1)(k-1)]\) by the assumption on \(u_j\), and \({\mathbf {c}}'_{j+1}(v) < k-1\) for any \(v >(t+1)(k-1)\) by the firing rule, which implies that \({\mathbf {c}}'_{j+1}(v) < k-1\) for any \(v >t(k-1)\). As \(t(k-1)\) is last fired at \({\mathbf {c}}'_{-1}\), we have \({\mathbf {c}}'_{{j+1}}={\mathbf {c}}\), also a contradiction.

(ii) The sufficiency of (ii) follows by the result (i) and the firing rule. For the necessity of (i), if \(t(k-1)\) is fired, then \({\mathbf {c}}(t(k-1)) \ge k-1\). If \({\mathbf {c}}(t(k-1)) \ge k\), there exists a configuration \({\mathbf {c}}'\) before \({\mathbf {c}}\) such that \({\mathbf {c}}'(t(k-1)) = k-1\) and \({\mathbf {c}}'(v) \ge k-1\) for some \(v > t(k-1)\), a contradiction to (i).

(iii) By definition, we have \(\omega (\bar{\mathbf {c}}_0)=\omega ({\mathbf {c}}_0)+k-1\). Observe that for any arrow \({\mathbf {c}}\xrightarrow {\alpha _{u, \varepsilon _t}} {\mathbf {c}}'\), where \(u \ne 0\), \(\omega ({\mathbf {c}}')=\omega ({\mathbf {c}})-1\) if \(\varepsilon _t=\varepsilon _1\), and \(\omega ({\mathbf {c}}')=\omega ({\mathbf {c}})\) if \(\varepsilon _t\ne \varepsilon _1\). Therefore, we need only show that no vertex is fired on \(\varepsilon _1\) more than one time in any directed path starting from \({\mathbf {c}}_0\). Let

be such a walk in \(\mathcal {G}({\mathbf {c}}_0)\). If there exists a vertex in \([k-1]\) be fired more than one time on \(\varepsilon _1\), we may assume that \(u_t\) is the first occurrence of the vertex satisfying \(\alpha _{u_t, \varepsilon _1}=\alpha _{u_\ell , \varepsilon _1}\) for some \(\ell < t\). If \(1\le u_\ell = u_t\le k-2\), by our assumption, we have

$$\begin{aligned} {\mathbf {c}}_t(u_t)\le \bar{\mathbf {c}}_0(u_t)-(k-1)+(k-2)={\mathbf {c}}_0(u_t)<k-1, \end{aligned}$$

since there are at most \(k-2\) vertices different from \(u_t\) in \(\varepsilon _1\backslash \{0\}\), each fired at most once, a contradiction.

If \(u_\ell =u_t=k-1\), we have \({\mathbf {c}}_\ell (v)<k-1\) for any \(v>k-1\) by the firing rule, \(u_i \in [k-2]\) for each \(i \in [\ell +1,t-1]\) and each of them is fired only one time by our assumption on \(u_t\). By the part (i), we have \({\mathbf {c}}_\ell (k-1)=k-1\). Therefore,

$$\begin{aligned} {\mathbf {c}}_t(k-1)\le {\mathbf {c}}_\ell (k-1)-(k-1)+(k-2)=k-2. \end{aligned}$$

This is also a contradiction. \(\square \)

Lemma 4.3

Let \(P_1^k\) be the k-uniform hyperpath on vertices \(0, 1, \ldots , k-1\) with a single edge \(\varepsilon \). Let \({\mathbf {c}}_0\) be a stable configuration on \(P_1^k\) with 0 as the bank vertex. Then, \({\mathbf {c}}_0\) is a critical configuration if and only if the firing graph \(\mathcal {G}({\mathbf {c}}_0)\) is a directed cycle of length k, and \({\mathbf {c}}_0\) is a non-critical configuration if and only if \(\mathcal {G}({\mathbf {c}}_0)\) is a directed path of length less than k.

Proof

By the firing rule, as \(P_1^k\) contains only one edge \(\varepsilon \), each firing is on the edge \(\varepsilon \) so that every vertex (configuration) of \(\mathcal {G}({\mathbf {c}}_0)\) has out-degree 1. So we have a sequence of finite configurations: \({\mathbf {c}}_0 \bar{\mathbf {c}}_0 {\mathbf {c}}_1 \cdots {\mathbf {c}}_m\), starting from \({\mathbf {c}}_0\) and ending at \({\mathbf {c}}_m\). Note \(\omega (\bar{\mathbf {c}}_0)=\omega ({\mathbf {c}}_0)+k-1\), and for any \({\mathbf {c}}\xrightarrow {\alpha _{u, \varepsilon _t}} {\mathbf {c}}'\), where \(u \ne 0\), \(\omega ({\mathbf {c}}')=\omega ({\mathbf {c}})-1\). So, starting from \(\bar{\mathbf {c}}_0\), the weight of configurations in the sequence is strictly decreasing.

If \({\mathbf {c}}_m={\mathbf {c}}_0\), or equivalently \({\mathbf {c}}_0\) is critical, then \(\mathcal {G}({\mathbf {c}}_0)\) is a directly cycle. As \(\omega ({\mathbf {c}}_m)=\omega ({\mathbf {c}}_0)\), by Lemma 4.2(iii), each vertex of \(\varepsilon \) is fired exactly once, implying that \(\mathcal {G}({\mathbf {c}}_0)\) is a directed cycle of length k. Otherwise, \({\mathbf {c}}_0\) is non-critical, and \(\mathcal {G}({\mathbf {c}}_0)\) is a directed path. We assert that not all of vertices in \(\varepsilon \) is fired; otherwise, by the proof of Lemma 4.2(iii), no vertex in \(\varepsilon \) is fired more than once so that each of them is fired exactly once. In this case, we would get \({\mathbf {c}}_m={\mathbf {c}}_0\), a contradiction. So \(\mathcal {G}({\mathbf {c}}_0)\) is a directed path of length less than k. \(\square \)

Remark 4.4

In Lemma 4.3, the dollar game on the hypergraph \(P_1^k\) (a single edge \(\varepsilon \)) is equivalent to considering the standard dollar game on a complete (simple) graph \(K_k\) on vertices \(0, 1, \ldots , k-1\) with 0 as the bank vertex. For the dollar game on simple graphs, the order of the firing can be arbitrary. However, here for dollar game on \(K_k\), the firing order obeys the maximal principle; see Sect. 4.1.

Next, we describe the structure of the firing graph of \(P_n^k\) associated with any stable configuration.

Proposition 4.5

Let \(P_n^k\) be the k-uniform hyperpath as in Fig. 3. Then, for a given stable configuration \({\mathbf {c}}_0\in \mathcal {B}_s\) on \(P_n^k\) with 0 as the bank vertex, where \(s\in [n]\), the firing graph \(\mathcal {G}({\mathbf {c}}_0)\) of \(P_n^k\) associated with \({\mathbf {c}}_0\) has the structure as in Fig. 5, where \(\mathcal {G}'\) is the subgraph induced by the configurations obtained by first firing the vertex \(s(k-1)\) on the edge \(\varepsilon _{s+1}\), each directed cycle in \(\mathcal {G}({\mathbf {c}}_0)\backslash \mathcal {G}'\) has length k, such that

  1. (i)

    \(\mathcal {G}'\) is not empty if and only if \(s<n\);

  2. (ii)

    if \(s<n\), there are no arrows from \(\mathcal {G}'\) to \(\mathcal {G}({\mathbf {c}}_0)\backslash \mathcal {G}'\);

  3. (iii)

    if \(s<n\), \({\mathbf {c}}_0\prec {\mathbf {c}}'\) for any \({\mathbf {c}}'\in \mathcal {G}'\).

Fig. 5
figure 5

Firing graph of \(P_n^k\) associated with \({\mathbf {c}}_0\in \mathcal {B}_s\)

Proof

(i) Let \({\mathbf {c}}_0\in \mathcal {B}_s\) be a fixed stable configuration on \(P_n^k\), where \(s\in [n]\). Recall that \(\tilde{\mathbf {c}}_0^i\) is a stable configuration on \(K_k\) for each \(i\in [n]\), where \(\tilde{\mathbf {c}}_0^i\) is defined at the beginning of this section. By definition, \({\mathbf {c}}_0\in \mathcal {B}_s\) means that \(\tilde{\mathbf {c}}_0^i\in \mathcal {C}_k\) for \(1\le i\le s\) and \(\tilde{\mathbf {c}}^{s+1}_0\in \mathcal {S}_k\backslash \mathcal {C}_k\). By Remark 4.4, to fire the vertices of \(\varepsilon _i\) on the edge \(\varepsilon _i\), it is equivalent to consider the dollar game on \(K_k\) on vertices of \(\varepsilon _i\) with the bank vertex \((i-1)(k-1)\), as only when \((i-1)(k-1)\) is firstly fired, the other vertices of \(\varepsilon _i\) can gain more dollars. By Lemma 4.3, as \(\tilde{\mathbf {c}}_0^1\) is a critical configuration, we get a directed cycle of length k by firing the vertices \(0,1, \ldots , k-1\) on \(\varepsilon _1\) and \(k-1\) also \(\varepsilon _2\), where \({\mathbf {c}}_1\) is the configuration such that \(k-1\) is a non-stable vertex with the maximal order in \({\mathbf {c}}_1\) ; see Fig. 6.

Fig. 6
figure 6

Directed cycle obtained by firing the vertices in \(\varepsilon _1\)

Similarly, since \(\tilde{\mathbf {c}}_0^2, \ldots , \tilde{\mathbf {c}}_0^s\) are critical stable configurations, each \(\varepsilon _i\) yields a directed cycle of length k by firing the vertices of \(\varepsilon _i\) on \(\varepsilon _i\), \(i \in [s]\). Hence, we get s directed cycles, which implies that \(\mathcal {G}'\) is not empty if and only if \(s <n\).

(ii) Denote

$$\mathcal {V}'=\{{\mathbf {c}}'\in \mathcal {G}({\mathbf {c}}_0)\mid {\mathbf {c}}' \hbox {~is obtained from~} {\mathbf {c}}_s \hbox {~by first firing~} s(k-1) \hbox {~on~} \varepsilon _{s+1} \},$$

and \(\omega '({\mathbf {c}})=\sum \limits _{v=s(k-1)+1}^{n(k-1)} {\mathbf {c}}(v)\) for any \({\mathbf {c}}\in \mathcal {G}({\mathbf {c}}_0)\). We claim that \(\omega '({\mathbf {c}}')>\omega '({\mathbf {c}})\) for any \({\mathbf {c}}'\in \mathcal {V}'\) and \({\mathbf {c}}\in \mathcal {V}(\mathcal {G}({\mathbf {c}}_0))\backslash \mathcal {V}'\).

In fact, we choose a configuration \({\mathbf {c}}'_m \in \mathcal {V}'\) such that \(\omega '({\mathbf {c}}'_m)=\min \{\omega '({\mathbf {c}}')\mid {\mathbf {c}}'\in \mathcal {V}'\}\), and a directed path \(P({\mathbf {c}}_s, {\mathbf {c}}'_m)\) from \({\mathbf {c}}_s\) to \({\mathbf {c}}'_m\) as follows:

We may also assume that \(\omega '({\mathbf {c}}'_i)>\omega '({\mathbf {c}}'_m)\) for any \(i \in [m-1]\). By the firing rule, we have \(\alpha _{u'_{m-1}, *}\ne \alpha _{u'_{m-1}, \varepsilon _j}\) for any \(j \in [s]\); otherwise \(\omega '({\mathbf {c}}'_{m-1})=\omega '({\mathbf {c}}'_m)\). So there exists some configurations among \({\mathbf {c}}'_1,\ldots ,{\mathbf {c}}'_{m-1}\) fired on some edge \(\varepsilon _t\) for \(t \ge s+1\).

If there exists \(\alpha _{u'_i, *}=\alpha _{s(k-1), \varepsilon _{s+1}}\) for some \(i \in [m-1]\), letting \(u'_i\) be the last one for the vertex \(s(k-1)\) fired on \(\varepsilon _{s+1}\), then it suffices to consider the restriction of the configurations on \([s(k-1), n(k-1)]\), which is equivalent to considering the firing graph of \(P^k_{n-s}\) associated with the configuration \(({\mathbf {c}}'_i(s(k-1)), \ldots , {\mathbf {c}}'_i(n(k-1)))\) with \(s(k-1)\) as the bank vertex, where \(P^k_{n-s}\) is the sub-hypergraph of \(P_n^k\) consisting of the edges \(\varepsilon _{s+1}, \cdots , \varepsilon _n\). By Lemma 4.2 (iii), we have \(\omega '({\mathbf {c}}'_m)\ge \omega '({\mathbf {c}}'_i)\), a contradiction. So, from \({\mathbf {c}}'_1\) to \({\mathbf {c}}'_m\), \(s(k-1)\) is no longer fired on \(\varepsilon _{s+1}\).

By the above discussion, in order to compute the \(\omega '({\mathbf {c}}'_m)\), it suffices to consider the firing graph of \(P^k_{n-s}\) associated with the configuration \(({\mathbf {c}}_s(s(k-1)), \ldots , {\mathbf {c}}_s(n(k-1)))\) with \(s(k-1)\) as the bank vertex. Note that \({\mathbf {c}}_s(v)={\mathbf {c}}_0(v)=\hat{\mathbf {c}}(v)\) for any \(v >s(k-1)\) and any \(\hat{\mathbf {c}}\) in the front s directed cycles of \(\mathcal {G}({\mathbf {c}}_0)\). By Lemma 4.2 (iii) again, we have

$$\omega '({\mathbf {c}}'_m)>\omega '({\mathbf {c}}_s)=\omega '({\mathbf {c}}_0)=\omega '(\hat{\mathbf {c}}),$$

since \(\tilde{\mathbf {c}}_0^{s+1}=({\mathbf {c}}_0(s(k-1)), \ldots , {\mathbf {c}}_0((s+1)(k-1)))\) is a non-critical. So the claim follows, and the subgraph \(\mathcal {G}'\) induced by \(\mathcal {V}'\) contains no configurations in the front s directed cycles of \(\mathcal {G}({\mathbf {c}}_0)\).

(iii) We will prove the assertion by induction on s. If \(s=0\), then Lemma 4.2(iii), \(\omega ({\mathbf {c}}')>\omega ({\mathbf {c}}_0)\) for any \({\mathbf {c}}'\) in \(\mathcal {G}'\) as \(\tilde{\mathbf {c}}^1_0\) is not critical, implying \({\mathbf {c}}_0 \prec {\mathbf {c}}'\). If \(s\ge 1\), also by Lemma 4.2 (iii), we have \(\omega ({\mathbf {c}}')\ge \omega ({\mathbf {c}}_0)\), and \(\omega ({\mathbf {c}}')=\omega ({\mathbf {c}}_0)\) if and only if the vertices \(1, \ldots , k-1\) have been fired on \(\varepsilon _1\) exactly once along any directed path from \({\mathbf {c}}_0\) to \({\mathbf {c}}'\). So it suffices to consider the case of \(\omega ({\mathbf {c}}')=\omega ({\mathbf {c}}_0)\).

Suppose that \(P({\mathbf {c}}_0, {\mathbf {c}}')\) is a directed path from \({\mathbf {c}}_0\) to \({\mathbf {c}}'\) as follows:

As \(\omega ({\mathbf {c}}')=\omega ({\mathbf {c}}_0)\), we have a unique \(\sigma (v) \in [t]\) such that \(\alpha _{u_{\sigma (v)}, *}=\alpha _{v, \varepsilon _1}\) for each \(v \in [k-1]\) by Lemma 4.2 (iii). It follows that

$$\begin{aligned} {\mathbf {c}}'(v)={\mathbf {c}}_0(v) \hbox {~for~} v \in [k-2]. \end{aligned}$$
(4.1)

By the construction of \(\mathcal {G}({\mathbf {c}}_0)\), there exist \(1\le i_1<j_1\le t\) such that \(\alpha _{u_{i_1}, *}=\alpha _{k-1, \varepsilon _2}\) and \(\alpha _{u_{j_1}, *}=\alpha _{k-1, \varepsilon _1}\). We also have

$$\begin{aligned} {\mathbf {c}}_0(v)={\mathbf {c}}_{i_1}(v), \; {\mathbf {c}}_{j_1}(v)={\mathbf {c}}'(v), \hbox {~for any~} v > k-1. \end{aligned}$$
(4.2)

If there exists an \(\ell _1\) such that \(i_1<\ell _1<j_1\) and \(1\le u_{\ell _1}\le k-2\), then by Lemma 4.2 (ii),

$$\begin{aligned}&{\mathbf {c}}_{{i_1}}(k-1) = {\mathbf {c}}_0(k-1)+1+m_1=k-1,\\&{\mathbf {c}}_{{j_1}}(k-1) = {\mathbf {c}}_{{\ell _1}}(k-1)+m_2=k-1, \end{aligned}$$

where \(m_1\) (respectively, \(m_2\)) is the number of vertices of \([k-2]\) fired on \(\varepsilon _1\) from \({\mathbf {c}}_1\) to \({\mathbf {c}}_{{i_1}-1}\) (respectively, from \({\mathbf {c}}_{{\ell _1}}\) to \({\mathbf {c}}_{{j_1}-1}\)). By the firing rule, \({\mathbf {c}}_{\ell _1}(k-1)<k-1\). So, after \({\mathbf {c}}_{{j_1}}\), only the vertices of \([k-2]\), \(k-1-(m_1+m_2+1)\) vertices in total, can be fired (certainly on \(\varepsilon _1\)). We have

$$\begin{aligned} \mathbf c '(k-1)&=k-1- (m_1+m_2+1)\\ {}&=\mathbf c _0(k-1)-((k-1)-\mathbf c _{u_{\ell _1}}(k- 1))<\mathbf c _0(k-1), \end{aligned}$$

which implies that \({\mathbf {c}}_0\prec {\mathbf {c}}'\).

Otherwise, we consider the dollar game on \(P^k_{n-1}\), where \(P^k_{n-1}\) is the sub-hypergraph of \(P_n^k\) consisting of the edge \(\varepsilon _2, \ldots , \varepsilon _n\) with \(k-1\) as the bank vertex. The game starts from \(\hat{\mathbf {c}}_{i_1}\) and produces a sequence of configurations: \(\hat{\mathbf {c}}_{i_1}, \hat{\mathbf {c}}_{i_1+1},\ldots , \hat{\mathbf {c}}_{j_1}\), where \(\hat{\mathbf {c}}_{p}\) is a restriction of \({\mathbf {c}}_p\) on the vertices of \(P^k_{n-1}\) for \(p \in [i_1,j_1]\). By the firing rule, as no vertices are fired on \(\varepsilon _1\), for \(p \in [i_1,j_1], \omega ({\mathbf {c}}_p)\)’s are constant, and \({\mathbf {c}}_p(v)\)’s are constant for each \(v\in [k-2]\). Let \(\hat{\omega }(\hat{\mathbf {c}}_p)=\sum _{v=k}^{n(k-1)}{\mathbf {c}}_p(v)\). As \({\mathbf {c}}_{{i_1}}(k-1)={\mathbf {c}}_{j_1}(k-1)=k-1\), we have \(\hat{\omega }(\hat{\mathbf {c}}_{i_1})=\hat{\omega }(\hat{\mathbf {c}}_{j_1})\), which implies that each vertex on \(\varepsilon _2\) rather than \(k-1\) will be fired exactly once by Lemma 4.2 (iii). So

$$\begin{aligned} {\mathbf {c}}'(k-1)={\mathbf {c}}_0(k-1). \end{aligned}$$
(4.3)

By induction, we have \(\hat{\mathbf {c}}_{i_1} \prec \hat{\mathbf {c}}_{j_1}\), and there exists a vertex \(u>k-1\) such that \(\hat{\mathbf {c}}_{i_1}(v)= \hat{\mathbf {c}}_{j_1}(v)\) for \(k-1 \le v \le u\) and \(\hat{\mathbf {c}}_{i_1}(u)> \hat{\mathbf {c}}_{j_1}(u)\). So, combining Eqs. (4.1), (4.3) and (4.2), we have \({\mathbf {c}}_0(v)= {\mathbf {c}}'(v)\) for \(1 \le v \le u\) and \({\mathbf {c}}_0(u)> {\mathbf {c}}'(u)\), implying \({\mathbf {c}}_0 \prec {\mathbf {c}}'\). \(\square \)

4.3 Formulas from firing graph

Our goal is to compute the determinant of the multiplication map by \(f_w\) of \(A=\dfrac{\mathbb {C}[x_v: v\in \widehat{V}]}{\langle \,f_v: v\in \widehat{V}\,\rangle }\) associated with the hyperpath \(P^k_n\) as in Fig. 3, where w is taken to be the vertex 0. We focus on the firing graph \(\mathcal {G}({\mathbf {c}}_0)\) of the structure shown in Fig. 5. If a configuration \({\mathbf {c}}\) in \(\mathcal {G}({\mathbf {c}}_0)\) refers to a homogeneous polynomial \(\mathbf{x}^{{\mathbf {c}}}=\prod \limits _{{{\mathbf {c}}}\in \widehat{V}} x_v^{{\mathbf {c}}(v)}\) by ignoring the bank vertex, then we have

$$\begin{aligned} m_{f_0}(\mathbf{x}^{{\mathbf {c}}_0})=\lambda \mathbf{x}^{{\mathbf {c}}_0}-\mathbf{x}^{\widehat{\varepsilon }_1}\mathbf{x}^{{\mathbf {c}}_0}=\lambda \mathbf{x}^{{\mathbf {c}}_0}-\mathbf{x}^{\bar{\mathbf {c}}_0}. \end{aligned}$$

If there is a non-stable vertex v in a configuration \({\mathbf {c}}\) with maximal order, then v will be fired on each edge containing v, which yields configurations \({\mathbf {c}}'\) and arcs \(({\mathbf {c}},{\mathbf {c}}')\). As \(f_v=0\) inside A, \(x_v^{k-1}=\frac{1}{\lambda } (\sum _{v \in \varepsilon \in E\backslash E_w} \mathbf{x}_{\varepsilon \backslash \{v\}}+ \sum _{v \in \varepsilon \in E_w}\mathbf{x}_{e_{\hat{w}}\backslash \{v\}})\). So we have

$$\begin{aligned} {\mathbf{x}}^{{\mathbf {c}}}=\sum _{({{\mathbf {c}}}, {{\mathbf {c}}}')\in \mathcal {E}(\mathcal {G}({\mathbf {c}}_0))} \dfrac{1}{\lambda }\mathbf{x}^{{\mathbf {c}}'} \end{aligned}$$

for any \({\mathbf {c}}\ne {\mathbf {c}}_0\). Intuitively, a firing of the bank vertex means a multiplication of \(f_0\), and a firing of other vertex v means a substitution of \(x_v^{k-1}\) by using \(f_v=0\).

In view of this, we consider the weighted directed graph, still denoted by \(\mathcal {G}({\mathbf {c}}_0)\), by assigning the weight 1 on the arrow \(\alpha _{0,\varepsilon _1}\) and the weight \(\frac{1}{\lambda }\) on the others. In order to obtain more explicit formulas, we need simplify the weighted firing graph \(\mathcal {G}({\mathbf {c}}_0)\). Observe that the directed graph \(\mathcal {G}({\mathbf {c}}_0)\) may contain non-stable configurations. Next, we will erase all non-stable configurations by modifying the weight until all of vertices are stable configurations.

Define a function \(g^i(x)\) in indeterminant \(\lambda \) recursively:

$$\begin{aligned} g^{-1}(x)=0, g^{0}(x)=1, g^{1}(x)=g(x)=\dfrac{1}{1-\frac{x}{\lambda ^k}}, g^{i}(x)=g^{i-1}(g(x)) \hbox {~for~} i \ge 2. \end{aligned}$$
(4.4)

Lemma 4.6

Let \({\mathbf {c}}_0\) be a configuration in \(\mathcal {B}_s\), \(0\le s \le n\). Then, we have

$$\begin{aligned} {\mathbf{x}}^{\bar{{\mathbf {c}}}_0}=\dfrac{g^{s-1}(1)}{\lambda ^{k-1}}\mathbf{x}^{{{\mathbf {c}}}_0}+\sum _{{{\mathbf {c}}}'\in \mathcal {S}(\mathcal {G}')} h_{{\mathbf {c}}'}(\lambda )\mathbf{x}^{{\mathbf {c}}'}. \end{aligned}$$

where \(\bar{\mathbf {c}}_0(v)={\left\{ \begin{array}{ll} {\mathbf {c}}_0(v)+1, &{} 1\le v\le k-1,\\ {\mathbf {c}}_0(v), &{}k\le v\le n(k-1)\end{array}\right. }\), \(g^i(x)\) is defined as in (4.4), \(h_{{\mathbf {c}}'}(\lambda )\) is a function in \(\lambda \) for each \({\mathbf {c}}'\), and \(\mathcal {S}(\mathcal {G}')\) is the set of all stable configurations in \(\mathcal {G}'\) defined in Proposition 4.5.

Proof

Clearly, for \(s=0\),

$$\begin{aligned} \mathbf{x}^{\bar{\mathbf {c}}_0}=\sum _{{\mathbf {c}}'\in \mathcal {S}(\mathcal {G}')} h_{{\mathbf {c}}'}(\lambda )\mathbf{x}^{{\mathbf {c}}'}=\frac{g^{s-1}(1)}{\lambda ^{k-1}}\mathbf{x}^{{\mathbf {c}}_0}+\sum _{{\mathbf {c}}'\in \mathcal {S}(\mathcal {G}')} h_{{\mathbf {c}}'}(\lambda )\mathbf{x}^{{\mathbf {c}}'}, \end{aligned}$$

and for \(s=1\),

$$\begin{aligned} \mathbf{x}^{\bar{\mathbf {c}}_0}=\dfrac{1}{\lambda ^{k-1}}\mathbf{x}^{{\mathbf {c}}_0}+\sum _{{\mathbf {c}}'\in \mathcal {S}(\mathcal {G}')} h_{{\mathbf {c}}'}(\lambda )\mathbf{x}^{{\mathbf {c}}'}. \end{aligned}$$

For \(s\ge 2\), from the subgraph of \(\mathcal {G}({\mathbf {c}}_0)\) in Fig. 7, we get

$$\begin{aligned} \mathbf{x}^{{\mathbf {c}}_{s-1}}=\dfrac{1}{\lambda }\mathbf{x}^{\bar{\mathbf {c}}_{s-1}}+\dfrac{1}{\lambda ^k}\mathbf{x}^{{\mathbf {c}}_{s-1}}+\dfrac{1}{\lambda ^{l_s}}\mathbf{x}^{{\mathbf {c}}'_1}, \end{aligned}$$

where \(l_s\) is the length of the directed path from \({\mathbf {c}}_{s-1}\) to \({\mathbf {c}}'_1\). It follows that

$$\begin{aligned} \mathbf{x}^{{\mathbf {c}}_{s-1}}=g(1)\dfrac{1}{\lambda } \mathbf{x}^{\bar{\mathbf {c}}_{s-1}}+g(1)\dfrac{1}{\lambda ^{l_s}}\mathbf{x}^{{\mathbf {c}}'_1}. \end{aligned}$$

Therefore, we can erase the sth directed cycle by adding an arrow \(({\mathbf {c}}_{s-1}, {\mathbf {c}}'_1)\) with weight \(\dfrac{g(1)}{\lambda ^{l_s}}\) and replacing the weight of \(({\mathbf {c}}_{s-1}, \bar{\mathbf {c}}_{s-1})\) by \(\dfrac{g(1)}{\lambda }\); see Fig. 8.

Fig. 7
figure 7

Weighted sth direct cycle by firing the vertices in \(\varepsilon _s\)

Fig. 8
figure 8

Modified weighted \((s-1)\)th direct cycle by erasing the sth directed cycle

Let \(l_j\) be the length of the directed path from \({\mathbf {c}}_{j-1}\) to \({\mathbf {c}}_j\) for \(2\le j\le s-1\). From this modified weighted subgraph, we immediately get

$$\begin{aligned} \mathbf{x}^{{\mathbf {c}}_{s-2}}=\dfrac{1}{\lambda }\mathbf{x}^{\bar{\mathbf {c}}_{s-2}}+g(1)\dfrac{1}{\lambda ^k}\mathbf{x}^{{\mathbf {c}}_{s-2}}+g(1)\dfrac{1}{\lambda ^{l_{s-1}+l_s}}\mathbf{x}^{{\mathbf {c}}'_1}, \end{aligned}$$

and

$$\begin{aligned} \mathbf{x}^{{\mathbf {c}}_{s-2}}=g^{2}(1)\dfrac{1}{\lambda } \mathbf{x}^{\bar{\mathbf {c}}_{s-2}}+g^{2}(1)g(1)\dfrac{1}{\lambda ^{l_{s-1}+l_s}}\mathbf{x}^{{\mathbf {c}}'_1}. \end{aligned}$$

Assume that \(\mathbf{x}^{{\mathbf {c}}_{s-i}}=g^{i}(1)\dfrac{1}{\lambda } \mathbf{x}^{\bar{\mathbf {c}}_{s-i}}+\dfrac{g^{i}(1)\cdots g^{2}(1)g(1)}{\lambda ^{l_{s-i+1}+l_{s-i+2}\cdots +l_s}}\mathbf{x}^{{\mathbf {c}}'_1}\) for \(i\ge 1\). Then, from the modified weighted subgraph (Fig. 9), we have

$$\begin{aligned} \mathbf{x}^{{\mathbf {c}}_{s-(i+1)}}= \dfrac{1}{\lambda } \mathbf{x}^{\bar{\mathbf {c}}_{s-(i+1)}} +g^{i}(1)\dfrac{1}{\lambda ^k} \mathbf{x}^{{\mathbf {c}}_{s-(i+1)}} +\dfrac{g^{i}(1)\cdots g^{2}(1)g(1)}{\lambda ^{l_{s-i}+l_{s-i+1}+\cdots +l_s}}\mathbf{x}^{{\mathbf {c}}'_1}, \end{aligned}$$

and hence

$$\begin{aligned} \mathbf{x}^{{\mathbf {c}}_{s-(i+1)}}= g^{i+1}(1)\dfrac{1}{\lambda } \mathbf{x}^{\bar{\mathbf {c}}_{s-(i+1)}} +\dfrac{g^{i+1}(1)\cdots g^{2}(1)g(1)}{\lambda ^{l_{s-i}+l_{s-i+1}+\cdots +l_s}}\mathbf{x}^{{\mathbf {c}}'_1}. \end{aligned}$$
Fig. 9
figure 9

Modified weighted \((s-i)\)th direct cycle by erasing the \((s-i+1)\)th, \(\ldots \), sth directed cycles

Taking \(i=s-2\), we get

$$\begin{aligned} \mathbf{x}^{{\mathbf {c}}_{1}}= g^{s-1}(1)\dfrac{1}{\lambda } \mathbf{x}^{\bar{\mathbf {c}}_{1}} +\dfrac{g^{s-1}(1)\cdots g^{2}(1)g(1)}{\lambda ^{l_{2}+l_{3}+\cdots +l_s}}\mathbf{x}^{{\mathbf {c}}'_1}. \end{aligned}$$

It follows that

$$\begin{aligned} \mathbf{x}^{\bar{\mathbf {c}}_0}=\dfrac{g^{s-1}(1)}{\lambda ^{k-1}}{\mathbf{x}}^{{\mathbf {c}}_0}+\prod _{i=1}^{s-1} g^{i}(1) \dfrac{1}{\lambda ^{\sum _{j=1}^{s}l_j}}{\mathbf{x}}^{{\mathbf {c}}'_1}, \end{aligned}$$

where \(l_1\) is the length of the directed path from \(\bar{\mathbf {c}}_0\) to \({\mathbf {c}}_1\). By Proposition 4.5 (i), we know that there is not an arrow from \(\mathcal {G}'\) to \(\mathcal {G}\) and therefore \(\mathbf{x}^{{\mathbf {c}}'_1}=\sum _{{\mathbf {c}}'\in \mathcal {S}(\mathcal {G}')}\bar{h}_{{\mathbf {c}}'}(\lambda ) \mathbf{x}^{{\mathbf {c}}'}\) for some function \(\bar{h}_{{\mathbf {c}}'}(\lambda )\) in \(\lambda \). It follows that

$$\begin{aligned} \mathbf{x}^{\bar{\mathbf {c}}_0}=\dfrac{g^{s-1}(1)}{\lambda ^{k-1}}\mathbf{x}^{{\mathbf {c}}_0}+\sum _{{\mathbf {c}}'\in \mathcal {S}(\mathcal {G}')} h_{{\mathbf {c}}'}(\lambda )\mathbf{x}^{{\mathbf {c}}'}, \end{aligned}$$

where \(h_{{\mathbf {c}}'}(\lambda )=\prod _{i=1}^{s-1} g^{i}(1) \dfrac{1}{\lambda ^{\sum _{j=1}^{s}l_j}}\bar{h}_{{\mathbf {c}}'}(\lambda )\). \(\square \)

4.4 The characteristic polynomial of a hyperpath

We will give a recursive formula of the characteristic polynomial of hyperpaths. Observe the number of stable configurations on \(K_k\) is \((k-1)^{k-1}\), and the number of critical configurations on \(K_k\) is \(k^{k-2}\) by Example 2.8. So \(\mathcal {B}_s\) contains exactly \(\mu _{n,k}(s)\) elements, which is defined as follows:

$$\begin{aligned} \mu _{n,k}(s)= \left\{ \begin{array}{ll} k^{s(k-2)}((k-1)^{k-1}-k^{k-2})(k-1)^{(n-s-1)(k-1)}, &{} s \in [0, n-1],\\ k^{n(k-2)}, &{} s=n. \end{array} \right. \end{aligned}$$
(4.5)

Theorem 4.7

Let \(\phi ^P_{n, k}(\lambda )\) be the characteristic polynomial of the hyperpath \(P_n^k\), where \(n\ge 2\). Then,

$$\begin{aligned} \phi ^P_{n, k}(\lambda )=\lambda ^{(k-2)(k-1)^{n(k-1)}}\prod _{s=0}^{n}\left( \lambda -\dfrac{g^{s-1}(1)}{\lambda ^{k-1}}\right) ^{\mu _{n,k}(s)}\phi ^P_{n-1, k}(\lambda )^{(k-1)^{k-1}} \end{aligned}$$

where \(g^i(x)\) is defined in (4.4) and \(\mu _{n,k}(s)\) is defined in (4.5).

Proof

By definition, the characteristic polynomial of \(P_n^k\) is

$$\begin{aligned} \phi ^P_{n, k}(\lambda )=\mathrm {Res}(F_0, F_1, \ldots , F_{r}), \end{aligned}$$

where

$$\begin{aligned} F_i(x_0, x_1, \ldots , x_r)=\lambda x_i^{k-1}-\sum \limits _{i\in \varepsilon \in E} \mathbf{x}^{\varepsilon \backslash \{i\}} \end{aligned}$$

is a homogeneous polynomial in variables \(x_0, x_1, \ldots , x_r\) of degree \(k-1\) for \(i=0, 1, \ldots , r=n(k-1)\). Clearly,

$$\begin{aligned} f_i(x_1, \ldots , x_r)= \left\{ \begin{array}{ll} \lambda -x_1\cdots x_{k-1}, &{} \hbox {~if~}i=0,\\ \lambda x_i^{k-1}-x_1\cdots x_{i-1}x_{i+1}\cdots x_{k-1}, &{} \hbox {~if~} i \in [k-2],\\ \lambda x_{k-1}^{k-1} -x_1\cdots x_{k-2}-x_k\cdots x_{2(k-1)}, &{} \hbox {~if~} i=k-1,\\ F_i(x_0, x_1, \ldots , x_r), &{} \hbox {~otherwise}, \end{array} \right. \end{aligned}$$

and

$$\begin{aligned} \bar{F}_i(x_1, \ldots , x_r)= \left\{ \begin{array}{ll} \lambda x_i^{k-1}, &{} \hbox {~if~} i \in [k-2],\\ \lambda x_{k-1}^{k-1} -x_k\cdots x_{2(k-1)}, &{} \hbox {~if~}i=k-1,\\ F_i(x_0, x_1, \ldots , x_r), &{} \hbox {~otherwise}. \end{array} \right. \end{aligned}$$

By Lemmas 2.2 and 2.3, we have

$$\begin{aligned} \phi ^P_{n, k}(\lambda )= & {} \mathrm {Res}(\bar{F}_1, \ldots , \bar{F}_r)^{k-1} \det (m_{f_0})\\= & {} \mathrm {Res}(\lambda x_1^{k-1}, \ldots , \lambda x_{k-2}^{k-1}, \bar{F}_{k-1}, \ldots , \bar{F}_r)^{k-1} \det (m_{f_0})\\= & {} \left( \mathrm {Res}(\lambda x_1^{k-1}, \ldots , \lambda x_{k-2}^{k-1})^{(k-1)^{(n-1)(k-1)+1}}\mathrm {Res}(\bar{F}_{k-1}, \ldots , \bar{F}_r)^{(k-1)^{k-2}}\right) ^{k-1}\\&\times \det (m_{f_0})\\= & {} \lambda ^{(k-2)(k-1)^{n(k-1)}}\phi ^P_{n-1, k}(\lambda )^{(k-1)^{k-1}}\det (m_{f_0}), \end{aligned}$$

where \(m_{f_0}\) is the multiplication map of the quotient algebra

$$\begin{aligned} A=\mathbb {C}[x_1, \ldots , x_r]/\langle f_1, \ldots , f_r\rangle \end{aligned}$$

given by \(m_{f_0}(x_1^{i_1}\cdots x_r^{i_r})=\lambda x_1^{i_1} \cdots x_r^{i_r}-x_1^{i_1+1}\cdots x_{k-1}^{i_{k-1}+1}x_{k}^{i_k} \cdots x_{r}^{i_r}\).

We choose a \(\mathbb {C}\)-basis \(\mathbf {B}=\{\mathbf{x}^{{\mathbf {c}}}{\mid } {{\mathbf {c}}}: [r]\rightarrow [k-2]\}\) for A, where \(\mathbf{x}^{{\mathbf {c}}}=x_1^{{{\mathbf {c}}}(1)}\cdots x_r^{{{\mathbf {c}}}(r)}\). In fact, for any \(\mathbf{x}^{{\mathbf {c}}} \in \mathbf {B}\), \(\tilde{\mathbf {c}}=(\tilde{\mathbf {c}}(0),{\mathbf {c}}(1),\ldots ,{\mathbf {c}}(r))\) can be viewed as a stable configuration on \(P_n^k\), where 0 is the bank vertex whose value can be omitted. We denote by \(\mathcal {B}\) the set of all stable configurations on \(P_n^k\). Clearly, there is a one-to-one correspondence between \(\mathbf {B}\) and \(\mathcal {B}\) by ignoring the bank vertex, and the left anti-lexicographical ordering on \(\mathcal {B}\) gives a total ordering on \(\mathbf {B}\). To be precise, \(x_1^{i_1}\cdots x_r^{i_r}\prec x_1^{j_1}\cdots x_r^{j_r}\) if and only if \(\sum \limits _{t=1}^r i_t<\sum \limits _{t=1}^r j_t\), or \(\sum \limits _{t=1}^r i_t=\sum \limits _{t=1}^r j_t\), \(i_s=j_s\) for \(1\le s\le t-1\) and \(i_t>j_t\) for some \(t\in [r]\).

Retaining the notation in Sect. 4.2, we observe that \(\mathcal {B}\) is exactly the disjoint union of \(\mathcal {B}_0, \mathcal {B}_1, \ldots , \mathcal {B}_n\). For each \(s \in [n]\), the number of configurations in \(\mathcal {B}_s\) is exactly \(\mu _{n,k}(s)\). Let \({\mathbf {c}}_0\) be a configuration in \(\mathcal {B}_s\). By Lemma 4.6, we know

$$\begin{aligned} m_{f_0}(\mathbf{x}^{{{\mathbf {c}}}_0})=\mathbf{x}^{\bar{{\mathbf {c}}}_0}=\left( \lambda -\dfrac{g^{s-1}(1)}{\lambda ^{k-1}}\right) \mathbf{x}^{{\mathbf {c}}_0}-\sum _{{\mathbf {c}}'\in \mathcal {S}(\mathcal {G}')} h_{{\mathbf {c}}'}(\lambda )\mathbf{x}^{{\mathbf {c}}'}, \end{aligned}$$

where \(\mathcal {G}'\), \(\mathcal {S}(\mathcal {G}')\) and \(h_{{\mathbf {c}}'}(\lambda )\) are defined in Proposition 4.5 or Lemma 4.6. On the other hand, by Proposition 4.5 (ii), we have \({\mathbf {c}}_0\prec {\mathbf {c}}'\) for any \({\mathbf {c}}'\) in \(\mathcal {G}'\). It follows that the matrix of \(m_{f_0}\) associated with the ordered basis \(\mathbf {B}\) is a lower triangle matrix with \(\frac{g^{s-1}(1)}{\lambda ^{k-1}}\) appearing on the diagonal exactly \(\mu _{n,k}(s)\) times for \(s \in [0,n]\). So

$$\begin{aligned} \det (m_{f_0})=\prod _{s=0}^{n}\left( \lambda -\dfrac{g^{s-1}(1)}{\lambda ^{k-1}}\right) ^{\mu _{n,k}(s)}, \end{aligned}$$

and

$$\begin{aligned} \phi ^P_{n, k}(\lambda )= \lambda ^{(k-2)(k-1)^{n(k-1)}}\prod _{s=0}^{n}\left( \lambda -\dfrac{g^{s-1}(1)}{\lambda ^{k-1}}\right) ^{\mu _{n,k}(s)} \phi ^P_{n-1, k}(\lambda )^{(k-1)^{k-1}}. \end{aligned}$$

\(\square \)

Example 4.8

By Theorem 4.7, we can deduce the characteristic polynomial of \(P_n^k\) inductively from an edge. The following are the characteristic polynomials \(\phi ^P_{n, k}(\lambda )\) of \(P_n^k\) for some specified n and k.

\(\phi ^P_{1,3}(\lambda )=\lambda ^3(\lambda ^3-1)^3\) of degree 12;

\(\phi ^P_{2, 3}(\lambda )=\lambda ^{35}(\lambda ^3-1)^6(\lambda ^3-2)^9\) of degree 80;

\(\phi ^P_{3, 3}(\lambda )=\lambda ^{151}(\lambda ^3-1)^{27}(\lambda ^3-2)^{18}(\lambda ^6-3\lambda ^3+1)^{27}\) of degree 448;

\(\phi ^P_{4, 3}(\lambda )=\lambda ^{891}(\lambda ^3-1)^{201}(\lambda ^3-2)^{81}(\lambda ^3-3)^{81}(\lambda ^6-3\lambda ^3+1)^{54}\) of degree 2304;

\(\phi ^P_{1, 4}(\lambda )=\lambda ^{44}(t^4-1)^{16}\) of degree 108;

\(\phi ^P_{2, 4}(\lambda )=\lambda ^{2671}(\lambda ^4-1)^{352}(\lambda ^4-2)^{256}\) of degree 5103;

\(\phi ^P_{3, 4}(\lambda )=\lambda ^{95774}(\lambda ^4-1)^{11440}(\lambda ^4-2)^{5632}(\lambda ^8-3\lambda ^4+1)^{4096}\) of degree 196830.

5 Starlike hypergraphs

In this section, we will deal with a class of k-uniform hypergraphs, called starlike hypergraph and denoted by \(S^k_{n_1,\ldots ,n_m}\), which is obtained from m hyperpaths of length \(n_1,\ldots ,n_m\) by sharing a common vertex w, where \(m \ge 1\) and \(n_i \ge 1\) for \(i \in [m]\); see Fig. 10. When \(m=1\), it is a hyperpath \(P_{n_1}^k\). When \(n_i=1\) for each \(i \in [m]\), it is called a hyperstar with m edges and denoted by \(S^k_m\).

Fig. 10
figure 10

Starlike hypergraph \(S^k_{n_1,n_2,\ldots ,n_m}\)

Theorem 5.1

Let \(S^k_{n_1,\ldots ,n_m}\) be the k-uniform starlike hypergraph with exactly t hyperpaths having length 1, where \(0 \le t \le m\). Then, the characteristic polynomial \(\phi _{n_1,\ldots ,n_m;k}^{SH}(\lambda )\) of \(S^k_{n_1,\ldots ,n_m}\) is

$$\begin{aligned}&\lambda ^{(m(k-2)+t)(k-1)^{\sum \limits _{j=1}^m r_j}} {\mathop {\mathop {\prod \limits }\limits _{i \in [m]}}\limits _{n_i>1}} \phi ^P_{n_i-1,k}(\lambda )^{(k-1)^{k-1+\sum \limits _{j \ne i}r_j}}\\&\qquad \times \,{\mathop {\mathop {\prod \limits }\limits _{s_i \in [0,n_i]}}\limits _{1 \le i \le m}}\left( \lambda -\sum \limits _{i=1}^{m}\dfrac{g^{s_i-1}(1)}{\lambda ^{k-1}}\right) ^{\prod \limits _{i=1}^m \mu _{n_i,k}(s_i)}, \end{aligned}$$

where \(r_i=n_i(k-1)\) for \(i \in [m]\), \(\phi ^P_{n,k}(\lambda )\) denotes the characteristic polynomial of \(P_n^k\), \(g^s(x)\) and \(\mu _{n,k}(s)\) are defined in (4.4) and (4.5), respectively.

Proof

For each \(i \in [m]\), let \(\varepsilon _i\) be the edge of \(P_{n_i}^k\) containing w, and \(\widehat{\varepsilon }_i=\varepsilon _i \backslash \{w\}\). Let \(H_i\) be sub-hyperpath of \(P_{n_i}^k\) by deleting the vertex w and the edge \(\varepsilon _i\) for \(i \in [m]\), which is a hyperpath \(P_{n_i-1}^k\) together with \(k-2\) isolated vertices if \(n_i>1\) or \(r_i(=k-1)\) isolated vertices otherwise. If \(n_i>1\), by Lemmas 2.2 and 2.3, we get

$$\begin{aligned} \phi _{H_i}(\lambda )=\lambda ^{(k-2)(k-1)^{r_i-1}}\phi _{n_i-1,k}^P(\lambda )^{(k-1)^{k-2}}; \end{aligned}$$
(5.1)

otherwise,

$$\begin{aligned} \phi _{H_i}(\lambda )=\lambda ^{(k-1)(k-1)^{r_i-1}}. \end{aligned}$$
(5.2)

So, combining (5.1) and (5.2), we have

$$\begin{aligned} \prod \limits _{i=1}^m \phi _{H_i}^{(k-1)^{1+\sum \limits _{j \ne i} r_j}}=\lambda ^{(m(k-2)+t)(k-1)^{\sum \limits _{j=1}^m r_j}} {\mathop {\mathop {\prod \limits }\limits _{i \in [m]}}\limits _{n_i>1} \phi ^P_{n_i-1,k}(\lambda )^{(k-1)^{k-1+\sum \limits _{j \ne i}r_j}}}. \end{aligned}$$
(5.3)

Let V be the vertex set of \(S^k_{n_1,\ldots ,n_m}\), and let \(V_i\) be the vertex set of \(P_{n_i}^k\) for \(i \in [m]\). Note that w is a cut vertex of \(S^k_{n_1,\ldots ,n_m}\) and recall the discussion in Sect. 3.2. The quotient algebra \(A=\dfrac{\mathbb {C}[x_v: v \in V\backslash \{w\}]}{\langle \,f_v:v \in V\backslash \{w\}\,\rangle }\) holds

$$\begin{aligned} A =A_1 \otimes \cdots \otimes A_m, \end{aligned}$$

where \(A_i=\dfrac{\mathbb {C}[x_v: v \in V_i\backslash \{w\}]}{\langle \,f_v:v \in V_i\backslash \{w\}\,\rangle }\) for \(i \in [m]\), and

$$\begin{aligned} f_v=\lambda x_v^{k-1}-\sum \limits _{v \in \varepsilon \in E(P_{n_i}^k)}x_{\varepsilon \backslash \{v\}}-\sum \limits _{v \in \widehat{\varepsilon }_i}x_{\widehat{\varepsilon }_i \backslash \{v\}}, \; v \in V_i\backslash \{w\}, i \in [m]. \end{aligned}$$

Let

$$\begin{aligned} f_w=\lambda -\sum \limits _{i=1}^m x_{\widehat{\varepsilon }_i} \end{aligned}$$

and \(m_{f_w}: A \rightarrow A\) be the linear map given by multiplication by \(f_w\). Then,

$$\begin{aligned} m_{f_w}=\lambda \mathrm{Id}_A-\sum \limits _{i=1}^m \mathrm{Id}_{A_1} \otimes \cdots \otimes \mathrm{Id}_{A_{i-1}} \otimes m_{i,w} \otimes \mathrm{Id}_{A_{i+1}} \otimes \cdots \otimes \mathrm{Id}_{A_m}, \end{aligned}$$

where \(m_{i,w}: A_i \rightarrow A_i\) is the linear map given by multiplication by \(\mathbf{x}_{\widehat{\varepsilon }_i}\) for \(i \in [m]\).

Fix an \(i \in [m]\). Let \(\mathbf {B}=\{\mathbf{x}^{{\mathbf {c}}}|{\mathbf {c}}: V_i\backslash \{w\} \rightarrow [k-2]\}\) be a basis of \(A_i\), where \(\mathbf{x}^{{\mathbf {c}}}=\prod _{v \in V_i\backslash \{w\}}x_v^{{\mathbf {c}}(v)}\), and let \(\mathcal {B}\) be the set of stable configuration of \(P^k_{n_i}\) with w as the bank vertex. As discussed in Theorem 4.7, there is a one-to-one correspondence between \(\mathbf {B}\) and \(\mathcal {B}\) by ignoring the bank vertex. We also have a left anti-lexicographical ordering \(\prec \) on \(\mathbf {B}\) arising from the order of \(\mathcal {B}\). Retaining the notation in Sect. 4.2, \(\mathcal {B}\) is the disjoint union of \(\mathcal {B}_0,\mathcal {B}_1,\ldots ,\mathcal {B}_{n_i}\), and the number of configurations in \(\mathcal {B}_s\) is \(\mu _{n_i,k}(s)\) for each \(s \in [0,n_i]\).

For each \({\mathbf {c}}_0 \in \mathcal {B}_s\), by Lemma 4.6,

$$m_{i,w}(\mathbf{x}^{{\mathbf {c}}_0})=\mathbf{x}^{\bar{\mathbf {c}}_0}=\dfrac{g^{s-1}(1)}{\lambda ^{k-1}}\mathbf{x}^{{\mathbf {c}}_0}+\sum _{{\mathbf {c}}' \in \mathcal {S}(\mathcal {G}')}h_{{\mathbf {c}}'}(\lambda )\mathbf{x}^{{\mathbf {c}}'},$$

where \(\mathcal {G}'\), \(\mathcal {S}(\mathcal {G}')\) and \(h_{{\mathbf {c}}'}(\lambda )\) are defined in Proposition 4.5 or Lemma 4.6. Note that \({\mathbf {c}}_0 \prec {\mathbf {c}}'\) for any \(c' \in \mathcal {G}'\). So the matrix of \(m_{i,w}\) associated with the basis \(\mathbf {B}\) under the above order is a lower triangle matrix with \(\frac{g^{s-1}(1)}{\lambda ^{k-1}}\) appearing on the diagonal exactly \(\mu _{n_i,k}(s)\) times for \(s \in [0,n_i]\).

By the above discussion,

$$\begin{aligned} \det (m_{f_w})={\mathop {\mathop {\prod \limits }\limits _{s_i \in [0,n_i]}}\limits _{1 \le i \le m}}\left( \lambda -\sum \limits _{i=1}^{m}\dfrac{g^{s_i-1}(1)}{\lambda ^{k-1}}\right) ^{\prod \limits _{i=1}^m \mu _{n_i,k}(s_i)}. \end{aligned}$$
(5.4)

The result follows by Corollary 3.2 and the equalities (5.3) and (5.4). \(\square \)

Taking \(n_i=1\) for \(i \in [m]\) in Theorem 5.1, we get the characteristic polynomial of hyperstar \(S_m^k\).

Corollary 5.2

Let \(\phi _{n,k}^S\) be the characteristic polynomial of the k-uniform hyperstar \(S_m^k\) with m edges. Then,

$$\begin{aligned} \phi _{m,k}^S(\lambda )=\lambda ^{r(k-1)^r} \prod \limits _{p=0}^m \left( \lambda -\dfrac{p}{\lambda ^{k-1}}\right) ^{{m \atopwithdelims ()p} k^{(k-2)p}\left( (k-1)^{{k-1}}-k^{k-2}\right) ^{n-p} } \end{aligned}$$

where \(r=m(k-1)\).

Corollary 5.3

[5, Theorem 4.3] Let E be the k-uniform hypergraph with k vertices and a single edge. Then,

$$\begin{aligned} \phi _E(\lambda )=\lambda ^{k(k-1)^{k-1}-k^{k-1}}(\lambda ^k-1)^{k^{k-2}}. \end{aligned}$$

Corollary 5.4

Let \(S_2^k\) be a k-uniform hyperstar with two edges or hyperpath with two edges. Then,

$$\begin{aligned} \phi _{S_2^k}(\lambda )=\lambda ^{\mu _k}(\lambda ^k-1)^{2k^{k-2}\left( (k-1)^{k-1}-k^{k-2}\right) }(\lambda ^k-2)^{k^{2(k-2)}}, \end{aligned}$$

where \(\mu _k=(2k-1)(k-1)^{2(k-1)}-2k^{k-1}(k-1)^{k-1}+k^{2k-3}\).

Corollary 5.5

The characteristic polynomial of the starlike hypergraph \(S^k_{1,1,2}\) is

$$\begin{aligned} \phi _{1,1,2;k}^{SH}(\lambda )&= \lambda ^{(4k-3)(k-1)^{4(k-1)}-(k^{k-1}+3k^{k-2})(k-1)^{3(k-1)}+3k^{2(k-2)}(k-1)^{2(k-1)}-k^{3(k-2)}(k-1)^{k-1}}\\&\cdot (\lambda ^k-1)^{k^{k-2}(k-1)^{3(k-1)}}\\&\cdot \left( \lambda -\frac{1}{\lambda ^{k-1}}\right) ^{k^{k-2}(3(k-1)^{k-1}-k^{k-2})((k-1)^{k-1}-k^{k-2})^2}\\&\cdot \left( \lambda -\frac{2}{\lambda ^{k-1}}\right) ^{k^{2(k-2)}(3(k-1)^{k-1}-2k^{k-2})((k-1)^{k-1}-k^{k-2})}\\&\cdot \left( \lambda -\frac{3}{\lambda ^{k-1}}\right) ^{k^{3(k-2)}((k-1)^{k-1}-k^{k-2})}\\&\cdot \left( \lambda -\frac{\lambda }{\lambda ^k-1}\right) ^{k^{2(k-2)}((k-1)^{k-1}-k^{k-2})^2}\\&\cdot \left( \lambda -\frac{1}{\lambda ^{k-1}}-\frac{\lambda }{\lambda ^k-1}\right) ^{2k^{3(k-2)}((k-1)^{k-1}-k^{k-2})}\\&\cdot \left( \lambda -\frac{2}{\lambda ^{k-1}}-\frac{\lambda }{\lambda ^k-1}\right) ^{k^{4(k-2)}}. \end{aligned}$$

In particular, the characteristic polynomial \(\phi _{1,1,2;3}^{SH}(\lambda )\) of \(S^3_{1,1,2}\) is

$$\begin{aligned} \lambda ^{980}(\lambda ^3-1)^{75}(\lambda ^3-2)^{54}(\lambda ^3-3)^{27}(\lambda ^4-2\lambda )^9(\lambda ^6-3\lambda ^3+1)^{54}(\lambda ^6-4\lambda ^3+2)^{81} \end{aligned}$$

of degree 2294.

6 Conclusion

We give an explicit and recursive formula for the characteristic polynomial of the adjacency tensor of a starlike hypergraph, which is a resultant of a system of polynomials related to the structure of the hypergraph. Variants of the chip-firing game on simple graphs such as dollar game on simple graphs or hypergraphs are applied to analyze these resultants. As such, we provide a combinatorial method for computing resultants, which will have potential use for commutative algebra, algebraic geometry and physical fields.

We note that there are many numerical methods and algorithms for computing partial (real or extreme) eigenvalues of a general (symmetric) tensor; see e.g., Chang et al. [3], Chen et al. [4], Cui et al. [8] and the references therein. We also note the starlike hypergraph is a power hypergraph \(G^k\), which is obtained from a starlike simple graph G by adding \(k-2\) vertices to each of its edges. Zhou et al. [21] proved that if \(\lambda \) is a nonzero eigenvalue of G or any subgraph of G, then \(\lambda ^{\frac{2}{k}}\) is an eigenvalue of \(G^k\). In fact, the nonzero eigenvalues of \(G^k\) are exactly those eigenvalues arising from G in the above way. Our method allows for calculating these same eigenvalues, with the added benefit of determining the algebraic multiplicity of each eigenvalue.