I Introduction
For various InternetofThings (IoT) applications, a massive number of devices are expected to access wireless networks via random access or machinetomachine (M2M) communication protocols [1, 2]. In order to support such rapidly growing demands, an efficient random access protocol providing enhanced throughput should be essentially developed, which is operated in a distributed and decentralized manner to cope with excessive uplink accesses for future wireless networks [3, 4].
Traditionally, it has been well known that the tree (or splitting) algorithm [5] and its variants [6, 7, 8] provide improved throughputs. In particular, the firstcome firstserve (FCFS) splitting algorithm [6] has been known as the best one so far, which yields the maximum throughput of (packets/slot). Splitting algorithms in [5, 6, 7, 8] have been adopted for various wireless systems [9, 10, 12, 13, 14, 15], as the physical layer technologies have evolved. They have been particularly considered for a cellular system with smart antennas [9] and with successive interference cancellation (SIC) [10, 12]. In particular, for the systems with SIC, the splitting algorithms have been further investigated with several power control schemes incorporated in [16]. Note that in general, when more than one packets can be decoded in a slot with advanced signal processing techniques, it is called multipacket reception (MPR) capable channel and thus it encompasses the systems discussed in [9, 10, 12, 16].
In principle, SICenabled techniques can provide substantial throughput improvement of random access systems by utilizing collided packets in the decoding procedure after subtracting interference from decoded packets in a sequential manner. Because of such potential improvement, SIC techniques have been widely adopted and studied in the literature [10, 17, 18, 19, 20]. A treebased splitting algorithm has been developed in [10] in order to resolve collision efficiently by the help of a SICenabled receiver and its modified algorithms have been proposed robust to errorprone networks [12, 11], achieving the maximum throughput of (packets/slot). By allowing more refined feedback for collision resolve beyond treebased splitting algorithms or allowing a partially centralized control or cooperation between users [17, 18, 19, 20], it was shown that the maximum throughput can reach up to one, depending on the types of feedback information and network environment [19, 20]. Decentralized power control and allocation methods have been proposed in [21, 22, 23] for random access with capabilities of MRP and SIC in order to maximize the sum rate of random access systems by optimally allocating the transmit power of packets and their transmission rates.
Recently, more flexible and generalized random access or multiple access frameworks beyond the traditional slotbased user control paradigm have been considered to support massive M2M or IoT traffic demands in an efficient manner [27]. Coded random access incorporating with packet erasure codes or ratelesstype codes with SIC decoding has been designed in [24, 25, 26], which provide flexibility and robustness for network dynamics such as fading, bursty noise, and user interference, yielding improved throughput compared to the slotted random access systems. More recently, nonorthogonal random access (NOMA) has been studied as a promising solution to improve the spectral efficiency for future cellular systems, e.g., fifth generation (5G) systems [28, 29, 30]. For various different types of NOMA systems, SIC is the most wellknown detection techniques and actively studied to boost the spectral efficiency in NOMA systems [28, 29, 30].
In spite of recently developed random access or multiple access frameworks and advanced physical layer coding schemes, in this paper, we focus on slotted random access systems with SIC decoding and treebased splitting algorithm to revolve collision via simple feedback information. It is worthwhile to emphasize that the maximum throughput of (packets/slot) in [10] is still quite challenging to achieve in practical environment. There might be three major difficulties to be addressed: 1) SIC capabilities can be limited in practice due to imperfect physical layer processing or excessive complexity; 2) the backlog size is hard to obtain in realtime, which is needed for optimal transmission control; 3) the event of SIC failure can occur because of channel estimation error. Therefore, the primary aim of this paper to address those constraints and provide throughput close to 0.693 (packets/slot) in more practical random access environment. Our main contributions are summarized as follows:

A universal online learning method for estimating the current backlog size inspired by the concept of Bayes’ rule and an adaptive transmission control algorithm operated with the estimated backlog size are proposed.

A novel analytical methodology to derive the maximum throughput limit and the related optimal transmission probabilities are presented, considering practical constraints of bounded SIC capabilities and SIC failure events.

A practical random access system with SIC incorporated with a treebased splitting algorithm is proposed, achieving the maximum throughput limit by optimally controlling and adapting its transmission under bounded SIC capabilities and SIC failure events.
Even though we present the proposed online learning and adaptation method in a standard slotted random access system, the developed online learning framework proposed in this paper is universally applicable for other types of random access systems.
Ii System Model
We consider a wireless communication system in which an access point (AP) is located at the center of a cell and time is slotted of equal length for one packet transmission. Users arrive to the system randomly and each user is assumed to have one packet to transmit to the AP. Therefore, a packet and a user are indistinguishable and they are used interchangeably in this paper. Such a system scenario was frequently discussed in the literature and is also suitable for the machine type communications to support internet of things (IoT) where most of the machine devices transmit one small packet to the AP. It is worthy to note that we do not restrict the traffic arrival to be Poisson which is frequently assumed in the literature on random access systems. Due to the random arrival of the users, the number of active users in the system, i.e., the backlog size, also varies over time. Such a timevarying backlog size is the main difficulty of designing efficient protocols for random access systems.
The backlogged users randomly transmit their packets with a certain probability broadcasted by the AP and consequently the number of packets arrived at the AP also varies over time. When a single packet is received at the AP, it can decode the packet successfully. On the other hand, if the AP receives multiple packets simultaneously at a slot, it is possible to perform SIC for decoding. Considering the difficulty of implementing SIC practically at the physical layer, in this paper we abstract the SIC capability by . That is, if there are more than users’ simultaneous transmissions, the AP is not able to perform SIC due to some technical limitations of physical layer decoding and this slot is denoted as a collision slot. Alternatively, if the number of transmitting users is between , the AP initiates an SIC resolve procedure (SRP), denoted by SRP, from the next slot. Within SRP, only the transmitted users are allowed to retransmit their packets under the SICenabled collision resolve algorithm which will be introduced in Section V. Note that when initiating SRP, the AP is aware of the value which is the minimum information of performing SIC. With this information, SRP can be terminated when all the packets are successfully received. We denote by the duration for SRP and by its mean. For differentiation, we define the slots that are not in any SPR by the normal slots. We assume that the AP provides the feedback information of success, collision and, initiating and terminating the SRP at the end of each slot instantaneously.
Fig. 1 shows an example of packet transmissions over time when the SIC capability is . At Slots 2, 5, and 20, only one packet is transmitted which would be successfully received by the AP. At Slot 7, two packets are transmitted and the AP initiates an SRP which lasts for 3 slots, i.e., . In terms of decoding at the AP, as it receives only Packet 3 at Slot 10 and it is able to decode Packet 4 by performing SIC for the signal received previously at Slot 7. How those packets are retransmitted in the SRP will be detailed in Section V. Similarly, at Slot 13, an SRP is initiated and, in this case, the duration is 5. Finally, at Slot 19, as there are more than users’ transmissions, a collision occurs and the AP cannot initiate an SRP. We define the embedded points as the end of each idle slot, success slot, collision slot and the last slot of each SRP, while we do not count the idle and success slots occurred within any SRP for the embedded points. As a representative example, in Fig. 1, all the embedded points are tagged. The embedded points are always located in front of the normal slots.
For such a system, our primary purpose is to optimally control the transmissions of the backlogged users so that the system could provide maximum throughput as well as minimum access delay for each user. Note that as the communication of the system bases on the random access, the AP has no information about users, such as who have packets to transmit and what the current backlog size is. To achieve such a goal, we design a transmission control algorithm by clearing the following missions one by one.

How to design the SICenabled collision resolve algorithm so that the duration of SRP is minimized for each . Note that the users involved in SRP should transmit their packets in a distributed manner, i.e., the AP cannot schedule specific user’s transmission as it has no information about the users who have packets.

How to control the transmissions of the backlogged users in normal slots (other than the slots in an SRP) so that the system operates optimally while exploiting the potentials of the SICenabled collision resolve algorithm?
We should comment here that the random access systems with SIC was also analyzed in [10] and it showed that the maximum throughput could reach 0.693. However, [10] considered an extreme scenario of unbounded SIC capability of , i.e., no matter how many packets that are transmitted simultaneously, the AP is always able to perform SIC. This assumption also simplifies the optimal transmission control of the backlogged users: whenever an SRP is completed, all the backlogged users if exist in the system should transmit with probability one (which is called Gated Access in [10]). Unfortunately, such an operation no longer the optimal for the random access systems with bounded SIC capability . For example, if the backlog size is larger than at a slot, the transmission probability of one would always result in collisions in the consequent slots. Therefore, alternative transmission control algorithm should be designed when it comes to the bounded SIC capability.
We shall design the SICenabled collision resolve algorithm in Section V to solve the above mentioned first mission. In Section III, we start our discussion by assuming that the backlog size is known at each normal slot and investigate the optimal system performance. This discussion further motivates us to design a transmission control algorithm that can estimate the the backlog size online while also accommodating the advantages of the SICenabled collision resolve algorithm in Sections IV.
Iii System Optimization Assuming Known Backlog Size
We first investigate an ideal scenario where the backlog size is assumed to be known at each normal slot. In practice, such information cannot obtained in practice as the traffic arrivals are random. Investigation on such an ideal scenario enables us to observe a performance upper bound and also gives us useful insight for designing a transmission control algorithm.
Suppose that there are backlogged users in a normal slot and they transmit their packet with probability . Then, according to Renewal Theory, the expected service rate can be derived as
(1) 
where
denotes the binomial distribution with parameters
and . The numerator indicates the average number of packets successfully transmitted between two consecutive embedded points and the denominator indicates the corresponding duration.Once the value of is known, based on (III), we can find the optimal transmission probability that maximizes the service rate. If such an optimal probability could be controlled slot by slot, we can maximize the system throughput and minimize each user’s access delay. While we can optimize the system operation if the backlog size is known prior, it is hard to obtain such information in practice. Motivated by this, in the next Section IV, we propose a transmission control algorithm that has the capability of estimating the backlog size online. It is worthy to highlight that online estimation allows the algorithm to adapt to network dynamics.
Iv Proposed Transmission Control Algorithm
In order to estimate the backlog size in each slot, we adopt the concept of Bayes’ rule. As well known, based on a priori distribution, Bayes’ rule introduces a posteriori distribution given that an event is observed. Therefore, when estimating the backlog size, instead of giving an constant number for the estimated backlog size, Bayes’ rule estimates it through introducing a distribution. If we define the observed event as the number of users transmitted at each slot, we can update the estimated distribution of the backlog size slot by slot recursively. That is, the a posteriori
distribution of the current slot becomes the a priori distribution of the next slot. Unfortunately, to implement such a process, we need to calculate the distribution at each slot which accompanies high computational complexity. In order to reduce the computational burden that occurs at each slot, we assume the backlog size follows the Poisson distribution, which is uniquely characterized by its mean. Then, by simply updating the mean we can update its distribution. Such an assumption was also made in
[32] and was shown to be effective when applying to random access systems.Let the distribution of the backlog size to be Poisson with mean , i.e,
(2) 
Then, for the estimation of the backlog size, we need to update recursively slot by slot.
Iva Optimal Transmission Probability
For a given backlog size , the service rate is given in (III). If
is a random variable with distribution
, the average service rate is expressed as(3) 
where we have applied the fact that and denotes the average number of transmissions per slot. Once and are given, we can numerically search for the optimal that maximizes the above service rate . Table I summarizes the optimal and the corresponding maximum service rate , where the values of are obtained from Section V . Consequently, the optimal transmission probability for a given mean backlog size is derived as
(4) 
1  1  0.3678  1.3922 

2  1.378  0.5586  2.0458 
3  1.739  0.6352  2.7020 
4  2.060  0.6665  3.3833 
5  2.3762  0.6802  4.0681 
10  3.8734  0.6926  7.5626 
From Table I, we observe that if the SIC capability of 2 or 3 can be obtained, the service rate of 0.5586 or 0.6352 is achievable, which corresponds to 80.6% and 91.6% of the maximum throughput of 0.693 obtained with unbounded SIC capability [10]. Therefore, from the engineering viewpoint, one can expect that implementing the SIC capability of 3 at the physical layer is reasonably good. It is notable that higher SIC capability requires more complicated advanced signal processing techniques, which in general accompany higher implementation complexity and cost.
IvB Estimation on the Backlog Size
Suppose that the transmission probability is . Then, the probability that among users are transmitting is obtained as
(5) 
Moreover, given the a priori distribution of , the unconditional probability of transmitting users is be derived as
(6) 
If there are users’ simultaneous transmissions in a slot, by applying Bayes’ rule, we obtain
(7) 
which is the a posteriori distribution of the backlog size . Interestingly, it is still Poisson but with an updated mean by
(8) 
If the optimal transmission probability shown in (4) is applied, we have
(9) 
where is given in Table I. Now we are ready to discuss different channel observations such as idle, success, SRP, and collision.
IvB1 Idle event
An idle event occurs when . Then, the mean backlog size of the a posteriori distribution is obtained as
(10) 
IvB2 Success event
A success event occurs when . Consequently, the mean backlog size of the a posteriori distribution is
(11) 
Since one packet is successfully received by the at the AP after a success event, the estimation on the mean backlog size should be further subtracted by 1.
IvB3 SRP event
When , the AP initiates an SRP. For such an event, the mean backlog size of the a posteriori distribution is obtained as
(12) 
Since those packet would be eventually received by the AP at the end of the SRP, the estimation on the mean backlog size should be further subtracted by at the end of the SRP.
IvB4 Collision event
When , a collision occurs and the AP is not able to know the number of users involved in the collision. The probability that a collision occurs when there are backlogged users in the system is obtained as
(13) 
Consequently, the unconditional probability of a collision event is obtained as
(14) 
Then, by applying Bayes’ rule, we obtain the conditional probability of given a collision event by
(15) 
It is not difficult to check that the above distribution is no longer the Poisson while we still approximate the distribution by Poisson with its mean as
(16) 
If the optimal transmission probability of is applied, we have
(17) 
where and its value for each is listed in Table I.
IvC Algorithm Description
Algorithm 1 describes the overall procedure of our proposed algorithm. Note that the values of and are listed in Table I. Lines 3, 6, 11, 14 correspond to the estimation on the arrival rate where the previous arrival rate is given by multiplying the weight of . In Section VII, we set in simulation. Lines 4, 7, 12, 15 correspond to the update on the mean backlog size whose mathematical background was introduced in the previous subsection. Line 17 corresponds to the calculation on the optimal transmission probability.
V SIC Resolve Algorithm
The proposed transmission control algorithm, introduced in the previous Section IV, adapts the transmission probabilities for packet transmissions at each normal slot. For the design of the proposed algorithm, we applied the results of the average duration of SRP which is initiated by users’ simultaneous transmissions. Once an SRP starts, only the users involved in the previous transmissions retransmit their packets under the SIC resolve algorithm in the subsequent slots. If all the packets are successfully received by the AP, SRP is terminated. Section VA introduces the SIC resolve algorithm and Section VB analyzes the expected duration of each SRP, i.e., .
Va Algorithm Description
Let denote the set of the users that make the AP initiate an SRP where denotes the cardinality of the set . Then, the SRP is proceeded by invoking the function SIC_Resolve() recursively, which is described in Algorithm 2. is an output of the function SIC_Resolve() and it indicates the duration spent for the function.
At the first slot of SIC_Resolve(), the users in retransmit their packet with probability as in Line 6. This retransmission probability could be minimized for each as the AP knows the value of when invoking the function. Each user who has transmitted in the first slot considers itself as a member of a new set as in Line 7. If , the AP decodes the only transmitted user’s packet and additionally performs SIC based on the previously received signals. If the number of nontransmitted users is larger than 2, those users become a member of the set as in Line 14. If , i.e., the number of transmitted users is larger than 1, the SIC_Resolve is invoked for as in Line 17. After the completion of SIC_Resolve() if exists, the function SIC_Resolve runs again for if . The time duration for SIC_Resolve() is updated at Lines 8, 18, and 22.
VB Expected Duration of SIC Resolving Time
In this subsection, we analyze the expected SIC resolve time achievable by Algorithm 2.
VB1 Recursive relation for SIC resolving time
Let be the number of time slots for SIC resolving when users are initially involved. Obviously, . To derive the recursive relation between , denote as the number of time slots required for starting the inherent SIC resolve phase like Lines 17 or 21 in Algorithm 2, and as the number of users for the inherent SIC resolve phase. Notice that both and are random.
From the proposed SIC resolve algorithm, we have , , and for ,
(18) 
For better understanding of the relation between Algorithm 2 and , explanatory cases for resolving SIC started from three and four users are provided in Figs. 2 and 3, respectively. Specifically, in Fig. 2, three users (users 1 to 3) are initially involved and after random transmissions with probability , users 1 and 3 retransmit after three time slots, i.e., and in this case. Note that user 2 does not participate the remaining retransmission procedure because of the condition in Line 20 of Algorithm 2 is not satisfied. Then users 1 and 3 perform random transmissions with probability and user 3 retransmits alone after three time slots. Again, from the condition in Line 20 of Algorithm 2, user 1 does not participate the remaining retransmission procedure, which gives and in this case.
In Fig. 3, four users (users 1 to 4) are initially involved and after random transmissions with probability , user 1 retransmits after two time slots, i.e., and . Then users 2, 3, and 4 perform random transmissions with probability and users 2 and 4 retransmit after two time slots. In this case, and . Finally, users 2 and 4 perform random transmissions with probability and user 2 alone retransmits after four time slots, i.e., and . Note that users 2 and 4 retransmit after two time slots, but the corresponding received signal is removed at the AP and users 1 and 4 perform the same random transmissions, see the condition in Line 5 of Algorithm 2.
VB2 Expected Duration Analysis
Now consider how to derive the expected time required for resolving SIC, i.e., for .
From the definition of
, it follows the geometric distribution with success probability
and, as a result, we have(19) 
Also, for , we have
(20) 
which satisfies .
Therefore, , , and from (18),
(21) 
for . Then, can be recursively calculated from (VB2) by using (19) and (20).
Recall that is the retransmission probability of each involved user when users are involved for the SIC resolving, which can be optimized to minimize the expected SIC resolving time. Let be the set of optimal retransmission probabilities minimizing . Denote the minimum expected number of time slots for resolving SIC among users by . Then, it is given by
(22) 
Table II summarizes for . Through numerical calculations, we also observed that the setting of for all yields almost identical expected duration for resolving SIC as Table II. Therefore, from the engineering viewpoint, setting a constant probability of is a reasonably good choice.
2  2 

3  3.333 
4  4.761 
5  6.210 
10  13.426 
Vi Impact of SIC Failure
Although introducing SIC capability in the physical layer decoding can provide substantial throughput improvement in random access systems, SIC is quite challenging in practice, especially as the number of users involving in the SIC procedure increases due to imperfect channel estimation, error propagation during SIC process, etc. In this section, we will address such practical issues and how to modify the proposed online transmission control algorithm in Section IV, minimizing the performance degradation due to SIC failure.
Let be the probability of SIC failure in each of SIC steps. For simplicity, we assume that an event of SIC failure independently occurs with probability for each time slot involving more than one user are simultaneously transmitted.^{1}^{1}1The same analysis presented in Section VI is straightforwardly applicable for the case where the probability of SIC failure depends on the number of simultaneously transmitted users.
Obviously, the most efficient way to deal with such an SIC failure event is that the corresponding users involved in the event retransmit again except the users whose packets are already decoded, to construct the same received signal at the AP so that the AP can resume the rest of SIC procedure.
We can apply the same notion of the SIC resolving time to analyze the impact of SIC failure. Let be the number of time slots for SIC resolving when users are initially involved counting the retransmission slots due to SIC failure. Note that the analytical method introduced in Section VB is able to provide a simple universal tool for dealing with such probabilistic events. Define as the number of time slots required to have the received signal with no SIC failure at the AP excluding the originally received signal. Then, its mean is given by from the geometric distribution. Now define , which corresponds to the number of time slots for SIC resolving when users are initially involved assuming that there is no SIC failure for the received signal involving users.
Then, similar methods in Section VB are applicable to derive recursive relations between and . In particular, and . Also,
and . For ,
and .
Similarly, let be the minimum expected number of time slots for resolving SIC among users. Unlike the case of no SIC failure, is not optimal if . For instance, Table III summarizes and the corresponding for when .
2  3  0.5 

3  4.788  0.412 
4  6.633  0.343 
5  8.486  0.288 
10  17.802  0.163 
Now, from the derived , we can obtain the maximum service rate and the corresponding parameters required for applying the proposed online transmission control algorithm as the same manner in Section IV. For instance, Table IV summarizes and the corresponding parameters for online estimation when .
1  0.3678  1  1.3922 

2  0.4821  1.2580  2.1220 
3  0.5155  1.4700  2.8892 
4  0.5264  1.6380  3.6960 
5  0.5300  1.760  4.5461 
10  0.5316  1.8840  9.2964 
Vii Numerical Evaluation
In this section, we perform numerical evaluation to demonstrate the performance of the proposed online estimation and transmission control algorithm for random access with SIC. Recall the throughput improvement in Table I requires the backlog size in realtime, which is impossible for most random access systems. Furthermore, in practice, SIC is quite challenging especially as the number of users involving in the SIC procedure increases due to imperfect channel estimation, error propagation during SIC process, etc. In the following, we will address these practical issues and demonstrate that the proposed online estimation and adaptation framework can still provide such throughput improvement under more practical environment. Unless otherwise stated, we assume Poisson packet arrivals and the update weight in Algorithm 1 is set as for the online estimation.
Viia Online Estimation and Adaptation
In order to investigate the performance of the proposed algorithm in Section IV, Fig. 4 compares the average throughput achievable by the proposed algorithm with that of the ideal case where the backlog size is known. As seen in the figure, the maximum service rate supportable by the proposed algorithm is the same for that of the ideal case. In order to investigate the impact of the online estimation and transmission control closely, Fig. 5 plots the average access delay with respect to . The average delay rapidly increases as the arrival rate converges to the maximum service rate so that the system becomes unstable for this regime. Note that the average delay gap between the online estimation and the ideal case is less than four time slots for the stable regime of . The results demonstrate that the proposed online estimation and transmission control performs well by only learning from the historical information.
Not only for Poisson user arrivals, the proposed algoirthm is universally applicable for a board class of arrival distributions. As an explanatory case, we assume that for every 100 time slots the arrival rate is changed either the Poisson distribution of mean or zero with equal probability. Hence, the mean arrival rate is equal to . Figs. 6 and 7 plot the average throughput and delay of this nonPoisson arrival model, respectively. For comparison, we also plot the performances of the poisson arrival model. Because the considered nonPoisson model corresponds to on–off user arrivals, the delay performance is expected to worse that that of the Poisson model. As seen in the figures, the proposed algorithm still performs well even for nonPoisson user arrivals.
To investigate the realtime estimation performance of the proposed algorithm, Fig 8 compares the estimated backlog size (the parameter in Algorithm 1) with the actual backlog size by observing time slots when , where time episodes are averaged out to plot the figure. For each time episode, we set for and and for . The result shows that Algorithm 1 can acurately estimate the actual backlog size in realtime, providing an efficient transmission control for random access systems reflecting network dynamics.
ViiB Impact of SIC Failure
To evaluation the performance degradation due to SIC failure, Figs. 9 and 10 plot the average throughput and delay performance of the proposed algorithm when the probability of SIC failure is given by . It is worthwhile to mention that the maximum arrival rate supported by the proposed algorithm coincides with the analytical results of in Table IV. In spite of the performance degradation due to SIC failure, the results show that SIC is still beneficial for improving throughput compared to the known best achievable throughput of without applying SIC [6] even for the very restrictive environment of and . The overall delay performance shows similar tendency as the case of no SIC failure as in Fig. 5. Because of retransmissions of SICfailed packets, increasing the access delay is inevitable. Nonetheless, the proposed algorithm optimally adjusts the transmission probabilities during the SICenabled collision resolve procedure to minimize the retransmissions due to SIC failure in probabilistic sense.
Viii Conclusion
In this paper, we proposed the online estimation and adaptive transmission algorithm for slotted random access systems. The proposed algorithm has been designed to fully utilize the SIC capability of the AP by carefully adjusting its transmission via online estimated information. Thanks to such online estimation and adaptation features, it is universally applicable for general nonPoisson arrivals and also robust to limited SIC capability and SIC failure.
References
 [1] E. de Carvalho, E. Bjornson, J. H. Sorensen, P. Popovski, and E. G. Larsson, “ Random access protocols for massive MIMO,” IEEE Commun. Mag, vol. 55, no. 5, pp. 216222, May 2017.
 [2] R. Harwahyu, R. G. Cheng, C. H. Wei, and R. F. Sari, “Optimization of random access channel in NBIoT,” IEEE Internet Things J., vol. 5, no. 1, pp. 391402, Feb. 2018.
 [3] B. Singh, O. Tirkkonen, Z. Li, and M. A. Uusitalo, “Contentionbased access for ultrareliable low latency uplink transmissions,” IEEE Wireless Commun. Lett., vol. 7, no. 2, pp. 182185, Apr. 2018.
 [4] K. Kannan, and M. Devaraju, “QoS supported adaptive and multichannel MAC protocol in vehicular adhoc network,” Cluster Comput., pp. 113, Mar. 2018.
 [5] J. J. Capetanakis, “Tree algorithm for packet broadcast channel,” IEEE Trans. Inf. Theory, vol. IT25, no. 5, pp. 505515, Sep. 1979.
 [6] D. Bertsekas and R. Gallager, Data Networks, PrenticeHall, 1987.
 [7] M. Sidi and I. Cidon, “Splitting protocols in presence of capture,” IEEE Trans. Inf. Theory, vol. 31, no. 2, pp. 295301, Mar. 1985.
 [8] G. Fayolle, P. Flajolet, M. Hofri, and P. Jacquet, “Analysis of a stack algorithm for random multipleaccess communication,” IEEE Trans. Inf. Theory, vol. 31, no. 2, pp. 244254, Mar. 1985.
 [9] H. Jin and A. Acampora, “Performance of a treebased collision resolution algorithm in cellular systems with smart antennas,” IEEE Trans. Wireless Commun., vol. 6, no. 3, pp. 11431151, Mar. 2007.
 [10] Y. Yu and G. B. Giannakis, “Highthroughput random access using successive interference cancellation in a tree algorithm,” IEEE Trans. Inf. Theory, vol. 53, no. 12, pp. 46284639, Dec. 2007.
 [11] X. Wang, Y. Yu, and G. B. Giannakis, “Design and analysis of crosslayer tree algorithms for wireless random access,” IEEE Trans. Commun., vol. 55, no. 12, pp. 22532256, Dec. 2007.
 [12] X. Wang, Y. Yu, and G. B. Giannakis, “A robust highthroughput tree algorithm using successive interference cancellation,” IEEE Trans. Wireless Commun., vol. 7, no. 3, pp. 909919, Mar. 2008.
 [13] X. Qin and R. Berry, “Opportunistic splitting algorithms for wireless networks,” in Proc. IEEE INFOCOM, pp. 16621672, Mar. 2004.
 [14] D. Y. Kim, H. Nam, and H. Jin, “PseudoBayesian broadcasting algorithm for opportunistic splitting scheduling systems,” IEEE Trans. Veh. Technol., vol. 66, no. 6, pp. 54515455. Jun. 2017.
 [15] D. Xu and Y. Yao, “Splitting tree algorithm for decentralized detection in sensor networks,” IEEE Trans. Wireless Commun., vol. 12, no. 12, pp. 60246033, Dec. 2013.
 [16] A. D. Gore and A. Karandikar, “Powercontrolled FCFS splitting algorithm for wireless networks,” IEEE Trans. Veh. Technol., vol. 59, no. 2, pp. 842856, Feb. 2010.
 [17] Y. Hu, W. Yang, and Y. Cai, “Throughput analysis of slotted ALOHA with cooperative transmission using successive interference cancellation,” Science in China Series F: Information Sciences, vol. 52, no. 12, pp. 23542359, Dec. 2009.
 [18] M. Mollanoori and M. Ghaderi, “Uplink scheduling in wireless networks with successive interference cancellation,” IEEE Trans. Mobile Comput., vol. 13, no. 5, pp. 11321144, May 2014.
 [19] E. Paolini, G. Liva, and M. Chiani, “Graphbased random access for the collision channel without feedback: Capacity bound,” in Proc. IEEE ICC, Dec. 2011.
 [20] S.H. Lee, B. C. Jung, and S.W. Jeon, “Successive interference cancellation with feedback for random access networks,” IEEE Commun. Lett., vol. 21, no. 4, pp. 825828, Apr. 2017.
 [21] C. Xu, L. Ping, P. Wang, S. Chan, and X. Lin, “Decentralized power control for random access with successive interference cancellation,” IEEE J. Select. Areas Commun., vol. 31, no. 11, pp. 23872396, Nov. 2013.
 [22] H. Lin, K. Ishibashi, W.Y. Shin, and T. Fujii, “A simple random access scheme with multilevel power allocation,” IEEE Commun. Lett., vol. 19, no. 12, pp. 21182121, Dec. 2015.
 [23] Y. Li and Lin Dai, “Maximum sum rate of slotted aloha with successive interference cancellation,” IEEE Trans. Commun., vol. 66, no. 11, pp. 53855400, Nov. 2018.
 [24] E. Paolini, G. Liva, and M. Chiani, “Coded slotted ALOHA: A graphbased method for uncoordinated multiple access,” IEEE Trans. Inf. Theory, vol. 61, no. 12, pp. 68156832, Dec. 2015.
 [25] C. Stefanovic and P. Popovski, “ALOHA random access that operates as a rateless code,” IEEE Trans. Wireless Commun., vol. 61, no. 11, pp. 46534662, Nov. 2013.
 [26] L. Lu, L. You, and S. C. Liew “Networkcoded multiple access,” IEEE Trans. Mobile Comput., vol. 13, no. 12, pp. 28532869, Dec. 2014.
 [27] M. Shirvanimoghaddam, M. Dohler, and S. J. Johnson, “Massive nonorthogonal multiple access for cellular IoT: Potentials and limitations,” IEEE Commun. Mag, vol. 55, no. 9, pp. 5561, Sep. 2017.
 [28] H. Zhang, D.K. Zhang, W.X. Meng, and C. Li, “User pairing algorithm with SIC in nonorthogonal multiple access system,” in Proc. IEEE ICC, May 2016.
 [29] B. Ling, C. Dong, J. Dai, and J. Lin, “Multiple decision aided successive interference cancellation receiver for NOMA systems,” IEEE Wireless Commun. Lett. vol. 6, no. 4, pp. 498501, Aug. 2017.
 [30] K. Hicuchi and A. Benjebbour, “Nonorthogonal multiple access (NOMA) with successive interference cancellation for future radio access,” IEICE Trans. Commun., vol. E98?B, no. 3, pp. 403414, Mar. 2015.
 [31] B. S. Tsybakov and V. A. Mikhailov, “Random multiple packet access: Partandtry algorithm,” (in Russian) Probl. Pered. Inform., vol. 16, no. 4, pp. 65–79, Oct./Dec. 1980.
 [32] R. L. Rivest, “Network control by Bayesian broadcast,” IEEE Trans. Inf. Theory, vol. 33, no. 3, pp. 323328, Mar. 1987.