
From the equations in (\[E:opt-p\]), we can solve the optimal $\mathbf{p}^\star\triangleq [p_1^\star, \dots, p_{K_P}^\star]^T$ as $$\label{E:PA}
  \mathbf{p}^\star = \mathbf{M}^{-1} \boldsymbol{\zeta},$$ where $\boldsymbol{\zeta} = [\zeta_1(\mathbf{W}_f^\star), \dots, \zeta_{K_P}(\mathbf{W}_f^\star)]^T$, and $\mathbf{M}\in\mathbb{C}^{{K_P}\times {K_P}}$ is defined as $$[\mathbf{M}]_{kj} = \left\{
      \begin{array}{ll}
        \frac{1}{\gamma_k}|\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{d,k}^\star|^2, & k = j,\\
         - |\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{t,j}^\star|^2, & k \neq j.
      \end{array}
  \right.$$

Finally, we summarize the proposed low-complexity algorithm to solve the original optimization problem (\[E:problem-MU\]) in general multiuser case in Table \[tab:Distributed-P-B-allocation\].

------------------------------------------------------------------------

[aaa]{} **Initialization:** Set $i=0$, $\underline{R_{sum}}=0$, and $\overline{R_{sum}} = \sum_{k=1}^{K_P} \log(1 +\mathtt{SINR}_k^{ub})$, where $\mathtt{SINR}_k^{ub} = \frac{P_0\|\mathbf{h}_{Pk}\|^2}{\sigma_n^2}$ is an upper bound of the SINR of ${\text{PUE}}_k$, which is achieved when the ICI disappears and only ${\text{PUE}}_k$ is served by the PBS.\
**Bisection Iteration:** At the $i$-th iteration, set $i \leftarrow i + 1$.

-   Compute $R_{sum}^0 = \frac{\underline{R_{sum}} + \overline{R_{sum}}}{2}$.

-   Given $R_{sum}^0$, compute $\{\bar{\lambda}_k
From the equations in (\[E: opt - p\ ]), we can solve the optimal $ \mathbf{p}^\star\triangleq [ p_1^\star, \dots, p_{K_P}^\star]^T$ as $ $ \label{E: PA } 
   \mathbf{p}^\star = \mathbf{M}^{-1 } \boldsymbol{\zeta},$$ where $ \boldsymbol{\zeta } = [ \zeta_1(\mathbf{W}_f^\star), \dots, \zeta_{K_P}(\mathbf{W}_f^\star)]^T$, and $ \mathbf{M}\in\mathbb{C}^{{K_P}\times { K_P}}$ is defined as $ $ [ \mathbf{M}]_{kj } = \left\ { 
       \begin{array}{ll } 
         \frac{1}{\gamma_k}|\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{d, k}^\star|^2, & k = j,\\ 
          - |\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{t, j}^\star|^2, & k \neq j. 
       \end{array } 
   \right.$$ 

 ultimately, we sum up the proposed low - complexity algorithm to resolve the original optimization trouble (\[E: problem - MU\ ]) in general multiuser case in Table   \[tab: circulate - P - B - allocation\ ]. 

 ------------------------------------------------------------------------ 

 [ aaa ] { } * * low-level formatting :* * rig $ i=0 $, $ \underline{R_{sum}}=0 $, and $ \overline{R_{sum } } = \sum_{k=1}^{K_P } \log(1 + \mathtt{SINR}_k^{ub})$, where $ \mathtt{SINR}_k^{ub } = \frac{P_0\|\mathbf{h}_{Pk}\|^2}{\sigma_n^2}$ is an upper bound of the SINR of $ { \text{PUE}}_k$, which is achieve when the ICI disappears and only $ { \text{PUE}}_k$ is served by the PBS.\ 
 * * Bisection Iteration :* * At the $ i$-th iteration, set $ i \leftarrow i + 1$. 

 -    Compute $ R_{sum}^0 = \frac{\underline{R_{sum } } + \overline{R_{sum}}}{2}$. 

 -    give $ R_{sum}^0 $, compute $ \{\bar{\lambda}_k

Frol the equations in (\[E:opt-p\]), we can solve tkw optikal $\mafhbf{p}^\staf\triangleq [p_1^\star, \dots, p_{K_P}^\stac]^T$ aw $$\labtj{E:PA}
  \mathbf{p}^\star = \mathbf{M}^{-1} \holdsymbil{\zeua},$$ where $\boldsymbol{\zeta} = [\dzta_1(\mafmbf{W}_f^\vvar), \dots, \zeta_{K_P}(\kathbf{W}_f^\stdr)]^T$, and $\mathbf{K}\iv\mcthbb{C}^{{K_P}\times {K_P}}$ is defined as $$[\mathbs{M}]_{kj} = \lrfh\{
      \begin{artay}{ll}
        \frac{1}{\gamma_k}|\mathbf{h}_{Pk}^H\tilde{\mathbf{s}}_{d,k}^\star|^2, & k = j,\\
         - |\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{t,u}^\stag|^2, & k \neq j.
      \enf{array}
  \rigyt.$$

Figqlly, we summxrize the krmposed low-domplexity algorithm to solve tfe ornginal optinizatlmn problem (\[X:problvm-MU\]) in general multigser caxe in Table \[tab:Cisvribyted-P-B-allocation\].

------------------------------------------------------------------------

[aaa]{} **Iiitialization:** Set $i=0$, $\onderline{R_{vuj}}=0$, and $\overline{R_{sun}} = \sum_{k=1}^{N_P} \lmg(1 +\mxrhtg{SIHR}_l^{ug})$, wherf $\methtt{SINR}_k^{ug} = \frac{P_0\|\matybf{h}_{Pk}\|^2}{\sigma_n^2}$ is an ip[vt bound of tge SINW jf ${\text{PUE}}_k$, which is achieved when the PCI sisappears and only ${\texr{PUE}}_k$ is served by thg PBS.\
**Biseceion Iteration:** At the $i$-th iteration, set $i \leftarrmw i + 1$.

-   Eimputd $R_{dum}^0 = \frac{\underline{R_{sum}} + \overline{R_{sum}}}{2}$.

-   Given $W_{auk}^0$, bompute $\{\bar{\lambda}_h
From the equations in (\[E:opt-p\]), we can optimal [p_1^\star, \dots, as $$\label{E:PA} \mathbf{p}^\star = \dots, \zeta_{K_P}(\mathbf{W}_f^\star)]^T$, and {K_P}}$ is defined $$[\mathbf{M}]_{kj} = \left\{ \begin{array}{ll} \frac{1}{\gamma_k}|\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{d,k}^\star|^2, & = j,\\ - |\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{t,j}^\star|^2, & k \neq j. \end{array} \right.$$ Finally, we summarize proposed low-complexity algorithm to solve the original optimization problem (\[E:problem-MU\]) in general multiuser in \[tab:Distributed-P-B-allocation\]. [aaa]{} Set $i=0$, $\underline{R_{sum}}=0$, and $\overline{R_{sum}} = \sum_{k=1}^{K_P} \log(1 +\mathtt{SINR}_k^{ub})$, where $\mathtt{SINR}_k^{ub} = \frac{P_0\|\mathbf{h}_{Pk}\|^2}{\sigma_n^2}$ is an upper of the SINR of ${\text{PUE}}_k$, which is achieved the ICI disappears and ${\text{PUE}}_k$ is served by the **Bisection At the iteration, $i i + 1$. Compute $R_{sum}^0 = \frac{\underline{R_{sum}} + \overline{R_{sum}}}{2}$. - Given $R_{sum}^0$, compute $\{\bar{\lambda}_k

From the equations in (\[E:opt-p\]), we Can solve thE optiMal $\MatHbF{p}^\stAr\trIangleq [p_1^\star, \doTS, p_{K_P}^\Star]^T$ as $$\label{E:PA}
  \mathbf{p}^\Star = \mAtHBf{M}^{-1} \bOLdSymboL{\zeta},$$ whERe $\BOLdsYmBoL{\zeTa} = [\ZEtA_1(\mathBf{W}_F^\star), \doTs, \zeta_{K_P}(\maThbF{W}_F^\star)]^T$, and $\matHBf{m}\in\mathbb{C}^{{k_P}\tImes {K_P}}$ is defiNed As $$[\mathBf{m}]_{kj} = \LEft\{
      \beGin{Array}{Ll}
        \frac{1}{\GAmma_k}|\mAthbf{h}_{Pk}^H\TiLDe{\mathBF{w}}_{d,k}^\staR|^2, & K = J,\\
         - |\mAthbF{h}_{Pk}^H\tilde{\mathbf{w}}_{T,J}^\sTAr|^2, & k \neq j.
      \end{arraY}
  \right.$$

fiNAlLY, We sUmmArize the prOpOsed lOW-compleXItY ALGorIThm to solve the Original optIMizAtion pRoBleM (\[e:problEm-MU\]) iN gENerAl multiuser Case In Table \[taB:DistrIButed-P-B-ALlocatiOn\].

------------------------------------------------------------------------

[aaa]{} **INitIalIzatIOn:** seT $i=0$, $\uNdERliNE{R_{Sum}}=0$, ANd $\oVerline{R_{SuM}} = \sUm_{k=1}^{K_P} \Log(1 +\mATHTT{SINr}_k^{uB})$, wheRe $\matHtt{SINR}_k^{ub} = \fraC{P_0\|\mAthbF{H}_{Pk}\|^2}{\Sigma_N^2}$ is an UppeR bOund oF the SInR of ${\tExT{PUE}}_k$, which is achIeveD when the IcI dIsAppEaRs and ONly ${\texT{PUe}}_k$ iS served By the PBs.\
**bisEcTION ITeration:** At the $i$-th itErATIoN, set $i \lefTarrow I + 1$.

-   coMpUTe $R_{sum}^0 = \frAc{\UndErliNE{r_{sum}} + \oVerlINe{r_{sum}}}{2}$.

-   GiveN $R_{sum}^0$, cOMpUtE $\{\bar{\lamBdA}_k

From the equations in (\[ E:opt-p\]) , wecan  so lv e th e op timal $\mathbf { p}^\ star\triangleq [p_1^\s tar,\d o ts,p _{ K_P}^ \star]^ T $a s  $$ \l ab el{ E: P A} 
  \m ath bf{p}^\ star = \ma thb f{ M}^{-1} \bol d sy mbol{\zeta },$ $ where $\bo lds ymbol{ \z eta }  = [\ zet a_1(\ mathbf { W}_f^\ star), \d ot s , \zet a _{K_P}( \ m at hbf{ W}_f^\star)]^T$,a nd $\mathbf{M}\in \mathb b{ C }^ { { K_P }\t imes {K_P} }$  is d e fined a s  $ $ [ \ mat h bf{M}]_{kj} =  \left\{
   \begin {a rra y }{ll}     \fr ac{1}{\gamm a_k} |\mathbf{ h}_{Pk } ^H\tild e {\mathb f{w}}_ {d, k}^ \sta r |^ 2,  &k=  j, \ \      - |\ma th bf {h}_{ Pk}^ H \ t i lde{ \ma thbf {w}}_ {t,j}^\star|^ 2,& k\ neq  j.
    \end {a rray} 
  \ri ght.$ $
Finally, we su mmar ize the p rop os edlo w-com p lexity  al gor ithm to  solvet heor i g i na l optimization pro bl e m  ( \[E:prob lem-MU \ ])  i n  general  m ult iuse r casein T a bl e \[tab: Distri b ut ed -P-B-al lo cation \] .
--- ----- - ---- ------ -------- ----- - -------------- - ------------- - -- - - -- - ---
[a aa]{} **Ini tial i zati on:* *  S et$ i=0$,  $\un de r li n e{R_{sum}}=0$, and$\ overli ne{R_ {sum}} = \sum _{k=1}^{K_ P } \log(1 + \mat h tt { SINR}_k^{ub})$ , whe re $\matht t {SINR}_k ^{ub}  = \frac {P_0\|\ma t h bf{h}_{P k}\ |^2 }{\ sig m a _n ^2}$ is an up p e r bo un d of th e S INR of${\ tex t{P UE} }_ k$, which  is achi ev ed  w he n t he IC I  disappe ar s a nd  on ly ${ \ text{P UE}}_ k$ i sse r ved  by the PB S . \
** Bi se ctio n I te ratio n:** Atthe $i$ -th itera tio n , se t$i  \lefta rrow i + 1$.
-    Compute  $ R_{ sum}^0 =  \frac{\ underline{R_{sum}} + \o v erline{ R_{ sum}} }{2} $.

-   G ive n $R_{ sum } ^0$, c ompute  $\{\ ba r{\ l a mbda} _ k 

From the_equations in_(\[E:opt-p\]), we can solve_the optimal_$\mathbf{p}^\star\triangleq_[p_1^\star, \dots,_p_{K_P}^\star]^T$_as $$\label{E:PA}
 _\mathbf{p}^\star = \mathbf{M}^{-1}_\boldsymbol{\zeta},$$ where $\boldsymbol{\zeta} =_[\zeta_1(\mathbf{W}_f^\star), \dots, \zeta_{K_P}(\mathbf{W}_f^\star)]^T$,_and_$\mathbf{M}\in\mathbb{C}^{{K_P}\times {K_P}}$ is defined as $$[\mathbf{M}]_{kj} = \left\{
      \begin{array}{ll}
__ ___  \frac{1}{\gamma_k}|\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{d,k}^\star|^2, & k_= j,\\
   _ _   - |\mathbf{h}_{Pk}^H\tilde{\mathbf{w}}_{t,j}^\star|^2, & k \neq_j.
_  _  \end{array}
  \right.$$

Finally, we summarize the proposed_low-complexity algorithm to solve the original_optimization problem (\[E:problem-MU\])_in_general_multiuser case in Table \[tab:Distributed-P-B-allocation\].

------------------------------------------------------------------------

[aaa]{}_**Initialization:** Set $i=0$, $\underline{R_{sum}}=0$, and $\overline{R_{sum}}_= \sum_{k=1}^{K_P} \log(1 +\mathtt{SINR}_k^{ub})$, where $\mathtt{SINR}_k^{ub}_= \frac{P_0\|\mathbf{h}_{Pk}\|^2}{\sigma_n^2}$ is an upper bound of_the SINR of ${\text{PUE}}_k$, which is_achieved when the ICI disappears_and only_${\text{PUE}}_k$ is served by the_PBS.\
**Bisection Iteration:** At_the $i$-th_iteration, set $i_\leftarrow i + 1$.

-  _Compute $R_{sum}^0 =_\frac{\underline{R_{sum}} + \overline{R_{sum}}}{2}$.

-   Given_$R_{sum}^0$,_compute $\{\bar{\lambda}_k