Entangled N-photon states for fair and optimal social decision making

Situations involving competition for resources among entities can be modeled by the competitive multi-armed bandit (CMAB) problem, which relates to social issues such as maximizing the total outcome and achieving the fairest resource repartition among individuals. In these respects, the intrinsic randomness and global properties of quantum states provide ideal tools for obtaining optimal solutions to this problem. Based on the previous study of the CMAB problem in the two-arm, two-player case, this paper presents the theoretical principles necessary to find polarization-entangled N-photon states that can optimize the total resource output while ensuring equality among players. These principles were applied to two-, three-, four-, and five-player cases by using numerical simulations to reproduce realistic configurations and find the best strategies to overcome potential misalignment between the polarization measurement systems of the players. Although a general formula for the N-player case is not presented here, general derivation rules and a verification algorithm are proposed. This report demonstrates the potential usability of quantum states in collective decision making with limited, probabilistic resources, which could serve as a first step toward quantum-based resource allocation systems.


Scientific Reports
| (2020) 10:20420 | https://doi.org/10.1038/s41598-020-77340-3 www.nature.com/scientificreports/ Entanglement is a fundamental property of composite systems in quantum physics that has attracted interest in game theory for solving deterministic problems with payoff matrices by finding Nash equilibrium 19 in competitive situations [20][21][22] . Some reports have also highlighted the benefits of using entanglement and a quantum approach in machine learning [23][24][25] . However, as the CMAB problem is probabilistic and the same kinds of solution algorithms cannot be followed, new approaches are needed.
We previously numerically and experimentally demonstrated that the polarization degree of freedom of photons can be used to allocate the rewards from two machines to two users efficiently by using polarizationentangled photon pair states 26 . In the present report, we extended this work to the use of N polarized photons in an N-photon quantum superposition state to solve the CMAB in more general situations with N ≥ 2 and two choices. The goal is to maximize both the reward output from the machines and the equality, or fairness, of the repartition of rewards among players, given that their possible actions are limited to a rotation of their own polarization measurement basis at every turn. More precisely, the intent is to obtain N-photon quantum superposition states with such properties that the optimal strategy for each player leads to an optimal situation for all players in terms of total outcome and fair repartition.
In this report, after defining the general situation under study, we identify the set of constraints on N-photon quantum superposition states to optimize the reward allocation for any N. From there, we present the derivation of the corresponding states for three, four, and five users and a comparison of the corresponding properties for the MAB solution. Finally, we discuss the convergence estimation of error-correction protocols with these states.

Results
Definitions and formulation of the problem. Problem description and hypotheses. The MAB problem is used to reproduce situations in which the available choices have outcomes with different probabilities and/or reward amounts, which are unknown to the user(s). In the CMAB problem, several users simultaneously face the same choices, also referred to as "machines" in the following, with or without communication with each other. If several users select the same machine, the potential reward is split between them. Figure 1 represents a typical example case of internet connection through relays, where the concepts of competition and limited communication between users are clearly involved.
According to the situation in Fig. 1, we defined four rules for our study: (i) Each machine has a fixed reward in case of success, equal to 1 for all, and 0 otherwise; (ii) If a machine is selected by k users and gives a reward, each user receives an individual reward of 1/k; (iii) The users know the reward amount given by each machine in case of success, as well as the total number of users, but not the probability of reward of the machines; (iv) The users do not communicate with each other.
Rule (i) is only for calculational convenience and does not cause loss of generality, as arbitrary and different rewards between machines produce the same results. The main requirements for Rule (ii) are that the total rewards of the machines remain the same regardless of how many users select them, and that some users are not favored over others. Without Rule (iii), users would have no means of deducing anything about the decisions of the other users and their consequences; thus, it is necessary to elaborate the optimization algorithms. Rule (iv) restricts the kinds of algorithms usable here and emulates actual network architecture constraints.
Definitions and formalism. We consider now the CMAB with N ≥ 2 and two machines, A and B. Let x i,A and x i,B be the reward amounts for turn i of machines A and B, respectively, and let r i,j be the reward amount received In every run, each player selects one machine to use to try to connect to the internet, which does so with a (unknown) probability P 1,2 . In case of connection, each player receives the same, shared bandwidth as the other players who selected the same machine. www.nature.com/scientificreports/ by user j ∈ 1, N in turn i . For a given turn i , user j makes a choice between A and B, with its reward obeying Rule (ii) such that Here ε i,j,k corresponds to the actual selection of machine k in turn i by user j , with k=A,B ε i,j,k = 1 . After a given number of turns n t , we define R j and R as R j is the accumulated reward for user j, and R is the total accumulated reward. Since N ≥ 2, it is evident that the optimal strategy involves having each machine selected by at least one user in every turn, so that the total accumulated reward R is maximized.
Fairness. Optimizing the total accumulated reward R does not imply anything about the repartitioning of the total reward among users. In the simplest case of two users and two machines with constant, different probabilities of giving the same reward, having the users always select the same machines would result in a heavily unbalanced situation. In contrast, if both users recognize the best machine and always select it simultaneously, they receive the same reward amount, although the total accumulated reward R decreases.
To quantify this effect, a metric is needed that indicates in a straightforward manner whether or not the current distribution is fair between users. This notion of fairness has been particularly studied in socioeconomics, in the use of income inequality indices such as the Gini index 27 , Hoover index 28 , and Theil index 29-31 , like the information entropy used in telecommunication. Telecommunication technology also relies on these metrics as figures of merit for resource sharing, with several definitions having been proposed and studied [32][33][34] .
In this study, we decided to use the Jain index I J to estimate the fairness between users for a given repartition R j j∈1,N after a given number of trials, where I J is given by The choice of this metric is justified by several important properties, following the discussion by 34 : • The metric should be continuous with respect to the reward variables; • The definition should not vary with N, explaining the normalization by N in the denominator; • Its value interval is between 0 and 1; • It is symmetric and sensitive to the variation in the reward of any user.
The first two properties are essential to establish a metric that is valid for any N and any reward distribution. Concerning the Jain index specifically, its main advantage is its intuitiveness: if only one user receives a reward while N − 1 others do not, then I J = 1 N , and for k equally rewarded users and N − k left-out ones, I J = k N . This aspect explains why I J is already used widely 34 in telecommunication to estimate the fairness of bandwidth allocation between channels.
Finally, fairness can be understood in two inequivalent ways: • Fairness of the reward accumulated so far, or • Fairness of the expected value of rewards to be obtained in the next trial (regardless of the accumulated reward differences).
The first case corresponds to active fairness, in which it is attempted to correct any uneven accumulated repartition in every step. The second case is characterized by a passive pattern, in which it is attempted to make the subsequent trials fair, without trying to correct unfairness from past events. In the example of telecommunication, an algorithm following the first case will attempt to compensate for differences in the accumulated data amount in every trial, while the second one will only attempt to achieve an equal data transmission rate on average between users for the subsequent trials. Since the second situation corresponds to a sufficient condition for the first situation after a transition period, we will focus here only on passive or instantaneous fairness. Note that the individual rewards in every step do not need to be equal for the situation to be fair: it only matters that the average expected value is the same for all users.
Targets and performance. With these tools, we can formulate the problem to be solved in this CMAB situation.
How can the maximum total reward available for all users be obtained, while simultaneously guaranteeing instantaneous fairness between users? Considering this objective, we used a modified version of the Jain index that involves both fairness and total reward performance. Let X k (n t ) = n t i=1 x i,k be the accumulated reward given (3) www.nature.com/scientificreports/ by machine k after n t trials and X(n t ) = X A (n t ) + X B (n t ) be the total reward available to users. The modified index is defined as This quantity corresponds to the fairness index considered by the global performance, which we call the pondered index. It possesses the same advantages as the Jain index while also evaluating the efficiency of the global strategy: if all users make the same choice, the repartition will be fair with poor performance, while a single user who always selects the best machine alone will generate a strongly unbalanced repartition. Thus, I p ∈ [0, 1] and is equal to 1 if and only if the repartition is fair and the total available reward is obtained by the users. The strategies discussed herein are all intended to make I p as close to 1 as possible. N = 2 case and constraints on the required quantum state. This section discusses the N = 2 case from the perspective of the fairness discussion introduced before and with the use of quantum superposition states employing the polarization of photons. The purpose is to apply the formalism to a well-known case so that situations with more than two users can be easily treated.
This situation is the simplest case of users having to share limited and uncertain resources among them, without direct communication between them. One possibility is to ensure, by construction, equality among users in the case of a sub-optimal reward produced by a selfish strategy. To maximize the total reward, each machine should be selected by one user in every step, and to maximize fairness, each user should select one machine as often as the other. In other words, users should avoid conflicts of decision between them. This principle leads to the difficult task of guaranteeing that users never select the same machine simultaneously when they cannot communicate with each other.
In our previous work 26 , we studied the situation in which users receive photon pairs in a quantum superposition of polarization states, while only being able to rotate a half-waveplate to optimize their rewards. Comparisons with other strategies such as independent random decision making and situations with selfish users are also provided. While the latter part will not be discussed here, we will describe the quantum state used and its properties in relation to the two-user CMAB.
Derivation of the state. We demonstrated previously 26 that polarization-entangled photon pairs can solve this issue when they are in state |ψ(φ)� 2 given by where φ is a real number and |H� and |V � are the horizontal and vertical linear polarization states of single photons, respectively. In that setup, each user has a half waveplate, which transforms any input state according to the angle θ/2 of its fast axis with respect to the |H� direction of the photon source, following Generally speaking, state |ψ(φ)� 2 only works if users have the same polarization measurement base {|H�, |V �} , or equivalently if θ = 0 [90 • ] with respect to the photon source (there is then no mixing of states while projecting onto the orthogonal states of a polarizing beam splitter cube). In contrast, for φ = π, which is invariant under any simultaneous rotation of the polarization measurement bases of the users. In this case, there is no need for direct communication between users to achieve an optimal situation, since only the relative angle between the measurement bases of the users matters.
Properties of the state. Next, assuming the state |ψ� 2 is sent, we discuss the evolution of fairness and total performance of the users when they modify their polarization bases independently by rotating their half-waveplates. It is assumed that this state is well controlled and reproducible and that the detection efficiency is 100%; in practice, limited detection efficiency can be overcome by post-selection or coincidence schemes, at the expense of lower fidelity of the state. Figure 2 shows the evolution of fairness and total performance with the rotation angles of the waveplates of both users after 1000 trials, averaged over 20 repetitions. Figure 2b shows that the fairness remains constant regardless of the combination of angles, which means that both players always receive the same accumulated reward regardless of their own angles. In addition, Fig. 2a demonstrates that the total reward depends only on the relative angle between the waveplates of the users, in accordance with the global rotational invariance of the state. From the perspective of the users, improving the reward of one user is thus equivalent to improving the rewards of both users equally. Finally, Fig. 2c presents the modified Jain index I p , which is equivalent here to the evolution of the total reward due to the perfect fairness in every case between users.
(4) www.nature.com/scientificreports/ Realignment algorithm convergence. This section discusses the determination of how users can actually proceed to have their polarization measurement bases aligned with each other (modulo 180°). According to Fig. 2, for any given tilt angle of one user, the other one can find an optimal position by tuning only its own waveplate. In other words, the action of one user is sufficient to find an optimal configuration for both.
In the previous study 26 , only one user attempted to find a correct angle configuration in the realignment algorithm, without any information about the angle of the other user. However, this strategy implies agreement regarding which user should act, because if both users act simultaneously they may not find an equilibrium situation.
These observations motivated us to implement a new, simpler, and scalable strategy, which involves selecting a random waveplate position when too many conflict situations are recorded over time. Given that Rules (i), (ii), and (iii) are verified, all users can recognize simultaneously when there is total conflict of decision and thus when realignment is needed, based on the reward amounts they receive (1/N in the case of success and full conflict). The exact algorithm is presented in section 1 of the Supplementary information.
To compare the performances of the two algorithms, we numerically studied a set of 100 random initial angles for the measurement basis of both users, with any value between 0° and 360° in increments of 5°, from which one player performed either the algorithm corresponding to the autonomous polarization-basis alignment of our previous study 26 (under assumption (II) in that report) or the proposed memory-based random algorithm. Briefly, the first algorithm supposes that no information is available to either user about the measurement basis position of the other user, and only one user adjusts its basis tilt in small increments until conflict is avoided. Each set of initial angles is repeated 20 times, and the averaged I p is estimated from the angle configuration in several time steps. We used identical parameters for both algorithms, with a memory capacity of eight rewardgiving events, a conflict event threshold of 2 before the angle was changed, and an angle increment of 5° for the previous algorithm. These parameters correspond to a trade-off between sensitivity to conflicts and robustness against possible errors. Figure 3 shows that both algorithms converge toward a stable, optimal configuration after several hundred plays, with the proposed algorithm converging faster. Using a coarser angle increment in the previous algorithm would probably make it converge faster, with the risk of skipping an optimal situation. As such, the proposed algorithm performs well for any initial configuration, and we chose to use it again for realignment estimation for higher numbers of players.  Constraints on the quantum state. Several criteria on the quantum state to be sent to the users can be introduced following the observations made in the N = 2 case: • Invariance of the measurement probabilities with respect to simultaneous rotation of the measurement bases of all users; • Symmetry of the state with respect to all other users; • No term with all users selecting the same choice simultaneously.
The first condition comes from the requirement that the state be device-independent; otherwise, the problem depends on the choice of the measurement basis of the photon source. The second condition means that with an odd number of users, such as three, the terms should have equal probabilities between an unbalanced situation (such as two users on A and one user on B) and its mirror (two users on B and one on A), as well as equal probabilities between all permutations of these terms. Finally, the last condition is associated with the search for the maximum total reward, which requires that all choices be attempted in every trial.
Derivation of the quantum state. Let us start from the three conditions given in the previous section to discuss whether or not there exists a family of compatible states for N = 3. This derivation is intended to demonstrate how these conditions constrain the form of the final quantum state obtained.
Global rotational invariance is best described using the density matrix formulation for the quantum states. Let |ψ� 3 be the entangled state under study and ρ 3 the associated density matrix. Following the last criteria given in the previous section and using the polarization basis {|H�, |V �} of the photon source, the state in the threephoton Hilbert space can be written as in its canonical basis, where the coefficients are all complex and include the normalization condition. From this equation, the necessary conditions on the coefficients can be derived such that the invariance of the state probability measurements does not vary under a simultaneous rotation θ ∈ R of all measurement bases, denoted hereafter as |θ�, θ ′ in the photon source reference frame. This condition involves the global rotation operator R 3 (θ) , defined in the three-photon Hilbert space H from the two-dimensional operator r(θ) as Then, the condition of invariance of the measurement probabilities under global rotation can be expressed as After applying (10) to the canonical basis of the three-photon Hilbert space, the following non-redundant conditions can be obtained: From (11) and (13), one can deduce that the complex cubic roots of 1 are involved. In addition, as the global phase of the state is irrelevant, there are 11 unknown variables to find (the amplitudes and phases of the parameters), with only 10 equations ((13) summarizes six equations). Four states can solve the equations, corresponding to the permutation of the phase differences and the choice of ±i in (12): All states have the same properties that are relevant to the present analysis, with symmetry in the rotation parameter space of each player. In the following, we discuss only the state 9) R 3 (θ) =r(θ) ⊗r(θ) ⊗r(θ),r(θ) = cos θ sin θ − sin θ cos θ .
(10) ∀θ ∈ R, ∀|φ� ∈ H, �φ|R(θ) † ρ 3R3 (θ)|φ� = �φ|ρ 3 |φ�. www.nature.com/scientificreports/ State properties. As discussed about the constraints, the total reward for all users is maximized when at least one user selects each machine in every turn. Thus, the total reward is reduced in conflict situations in which all users select one machine. Figure 4 shows three-dimensional scatter plots of the total reward, fairness, and pondered index as functions of the angle of the polarization measurement basis of each user relative to the photon source basis, with 5° discretization. The color and size of each sphere are functions of the variable studied. The fairness exhibits the greatest difference from the N = 2 case here, as there exist situations in which both total reward and fairness are sub-optimal. The total rewards still vary considerably over the parameter space, and the combination of the variations explains the more significant variation of the pondered index with respect to the N = 2 case.
Note that in addition to π-periodic optimal situations and invariance with respect to simultaneous rotation, additional optimal situations can be found for specific combinations of angles, indicated by the red circles in Fig. 4c and corresponding to rotations of π/3 or 2π/3 . This characteristic corresponds to the fact that permutations of the phases of several terms in |ψ� 3 yield the same performance. This feature leads to a parameter space in which many combinations of angles between the three users are favorable.

Realignment algorithm convergence.
In accordance with the objective of achieving device-independent behavior, the users should be able to find the optimal situation without needing to communicate with each other. Supposing that their measurement bases are misaligned with each other, there should be a mean of retrieving a favorable situation as quickly as possible. However, unlike in the N = 2 case, the actions of one user alone are not sufficient to reach the optimal angle combination in the general case, as can be seen in Fig. 4 along the tuning parameter axis of each user taken separately. Consequently, the random exploration strategy introduced in Realignment algorithm convergence of the N = 2 case is needed.
To estimate the efficiency and speed of convergence of this strategy, we selected 100 initial random angle combinations and repeated the calculations 20 times for each combination to obtain a statistical average. At different time steps, the angle combination was taken and attempted 1000 times to estimate the instantaneous I p . Figure 5a shows the evolution of I p over time when one, two, or three users follow this strategy while the other(s) remain fixed, still with 5° discretization of the available angles for computational purposes. It is evident that when at least two users react to conflicts, the average I p increases steadily until it reaches more than 0.98 after 10,000 events, while it remains bounded at 0.95 on average when only one user moves. At least two users are thus necessary to reach the optimal situation.
We also tested the evolutionary stability 35 of this strategy by checking whether users are at a disadvantage or an advantage if they use this strategy or if they remain passive, relying on other users to improve the situation. Figure 5b depicts the estimated reward given the current angle configuration for each user when only user  www.nature.com/scientificreports/ 1 is moving, while Fig. 5c shows the situation when users 1 and 2 react to conflicts, but not user 3. As can be seen, a passive user always receives a lesser or equal reward compared to the active users, while all users always benefit from attempting to improve the situation. It is thus beneficial for users to follow this algorithm rather than remaining passive, which confirms the evolutionary stability of the algorithm. Besides, if any user tries to increase its own outcome, other users can easily recognize it and get even higher reward than the selfish user, thus pushing users to play for the common best situation. N = 4 case. We have already identified the differences between odd and even numbers of users, mainly the necessary symmetry between terms that gives unbalanced rewards, such as when there are two users for machine A and one user for machine B, as well as similarity, such as invariance by simultaneous rotation of the measurement bases of the users. This part will focus on the N = 4 case to identify the properties that are preserved when the situation is scaled up according to the number of users.
Differences from the N = 2 and N = 3 cases. The same constraints apply to the state to be obtained with respect to the N = 3 case: invariance with respect to simultaneous rotation of the polarization measurement bases, symmetry between users, and avoidance of terms in which all users select the same machine simultaneously. However, this last constraint is fulfilled by two types of terms: • Asymmetric terms in which three users select one machine and one user selects the other, such as |HHHV � and |VVVH�, • Symmetric terms in which users are split equally between machines, such as |HHVV � and its permutations.
A priori, both kinds of terms are usable here as long as their relative amplitudes are equal between permutations of the photons of the users; other considerations are required to determine the expression of the state.
Derivation of the quantum state. The state can be derived using the same technique as for N = 3, starting from global rotation invariance of the density matrix of the target state. We chose to study separately the asymmetric and symmetric terms for clarity.
Let |S� 4 and |A� 4 respectively be the symmetric and asymmetric target states in the corresponding Hilbert space, with all coefficients being complex: Then, the invariance by simultaneous rotation of any angle θ using the operator R 4 (θ) is applied. For the symmetric states |S� 4 , the conditions on the complex coefficients for invariance under operator R 4 (θ) follow: As before, the global phase can be fixed at 0 without loss of generality, which gives 10 equations for 11 unknowns. The constraints are the same as for the N = 3 case, with the difference that there is no ±i factor anymore between mirror terms. There are two possible states of this kind: In the following, we will discuss only the case z = e 2iπ 3 , as both states produce the same results in this analysis. For |A� 4 , the global rotation invariance gives the conditions Excluding the global phase, there are 15 unknowns for 13 equations, leaving the choice of one relative phase between a terms and the permutation order (for a given phase difference between, say, a 1 and a 2 , we can have either a 3 = −a 1 or a 4 = −a 1 ). This characteristic gives the family of states (22) a 1 + a 2 + a 3 + a 4 = 0 State properties. As before, the objective is to study the variations of the performances of the states when the users vary their waveplate angles. As a four-dimensional representation of the parameter space is inconvenient to display and comprehend, we used the global rotation invariance and fixed the angle of one user at 0°, while the three other users varied their angles. The other configurations can be deduced by translation of the threedimensional figure along the three-dimensional diagonal. Figure 6 shows the total reward, fairness, and I p for the symmetric state |S� 4 when users 1, 2, and 3 tilt their waveplates, with 5° discretization, and user 4 remains fixed at 0. While the fairness remains almost optimal in any configuration, the conflict rate (and, thus, the total reward) changes significantly while only particular states give optimal situations, typically with a period of 180°. Between them, conflicts and/or fairness are suboptimal, and the density of optimal combinations in the four-dimensional space is relatively low. Regarding realignment, at least three users generally need to tune their half-waveplates to achieve an optimal situation. Figure 6c shows that for any fixed angle, for example, for user 3, there is no optimal configuration unless users 3 and 4 have the same angle (here, 0°). This finding is confirmed by the simulation results presented in Sect. 2 of the Supplementary information.
All of the asymmetric states |A� 4 (φ) show fairness greater than 0.995 for any angle configuration. Thus, the variations of I p only depend on the total reward. To understand the influence of φ on the state performance, we studied the following states: Figure 7a-c show only I p for the asymmetric states |A� 4 (0),|A� 4 + π 2 , and |A� 4 (π) , respectively, with the same user configuration as for |S� 4 . For φ = 0 and φ = π , there exist several planes of optimal situations, passing   www.nature.com/scientificreports/ through points (0 [π], 0 [π], 0 [π]) that are not on the same edge of the cube of length π . However, for φ = π 2 , the optimal situations are only on lines corresponding to the intersections of the planes visible in (a) and (c).
Optimal planes such as those in Fig. 7 can be obtained for any permutation of the terms in (26) as long as φ = 0 [π] . The locations of these planes can be explained by the choice of phase terms in (26) and the sign of the phase difference in the second group of terms. Different signs give different planes in the parameter, although they always pass through points (0 [π], 0 [π], 0 [π]) that are not on the same edge of the cube of length π for a given plane. Furthermore, if user 4 has any fixed angle θ 4 , the pattern in the parameter space of the three remaining users will simply be shifted, and the new optimal planes will pass through points (θ 4 [π], θ 4 [π], θ 4 [π]).
Regarding the realignment procedure, in all cases and for any combination of fixed angles for two users, such as users 3 and 4, there always exists an optimal situation accessible to the two remaining users, as can be seen in Fig. 7. Thus, two users are sufficient to obtain optimal conditions, even though the other users remain fixed. In addition, fairness is always achieved; thus, it is of the best interest of each user to improve the situation for everyone, which is directly linked to the situation of the given user. This finding is also confirmed by the simulation results corresponding to this case, as presented in section 2 of the Supplementary information. Besides, when φ = 0 [π] , one user can reach an optimal situation from any initial 4-angle configuration by itself, as all lines parallel to x, y or z axes intersect an optimal plane: in those particular cases, the action of one user is enough to realign the whole system.

Discussion
In the previous sections, we demonstrated the application of quantum superposition of the polarization states of two, three, and four users who must decide between two machines. Explicit expressions of the optimal states in each case were found, as well as expressions of the optimal states in the N = 5 case, whose derivation is presented in Section 3 of the Supplementary information. In addition, we developed a script using Wolfram Mathematica, presented in Section 4 of the Supplementary information, that can identify whether or not a given N-photon state is symmetric and invariant under simultaneous rotation of all users. A generalization of the theoretical work presented here for any N would be an interesting development, as entangled photon states with more than 10 photons have already been generated experimentally. Based on the empirical observations, we hypothesize that complex N-roots of 1 are involved in the phase coefficients for the optimal state formulation, regardless of N.
In addition, we observed differences between odd and even N: while maximum fairness is achieved for even N regardless of the combination of user waveplate angles, odd N does not exhibit this property. This fact remains to be explained.
As for any resource sharing system, security is a major concern similarly with all telecommunication systems. The system described in this article relies on two fundamental features: • The properties are device independent, so any new user can be added with initial alignment verification.
• Equality among users and resource allocation performance can be achieved without having to rely blindly on a remote central entity.
Regarding the first point, device independence is guaranteed by the rotation-independent property of the quantum states used in this study. Regarding the second point, the central entity corresponds to the N-photon source providing for all users. For this objective to be achieved, the users need a way to check whether the state received indeed corresponds to the optimal state expected (given the total number of simultaneous users and the expected total throughput). Such protocols have already been developed and verified for entangled photon pairs in the case of quantum key distribution [36][37][38] , including in a star-type network between an arbitrary number of users 39 , and several studies have proven the feasibility for more than two parties [40][41][42] , although limitations have been identified for large N 43 . All this may be of particular interest to secure applications of quantum information and resource sharing such as voting strategies using quantum states 44 .
Last but not least, this analysis relies on means of producing fully tunable quantum superposition of N photons with their polarizations as the entangled degree of freedom. Recent works have shown ways to obtain a Greenberger-Horne-Zeilinger state for up to 10 photons 45 , while experiments have already succeeded in generating any 6-photon state 46 . As such, the production of N-photon states is less a theoretical issue than a technological one.

Conclusion
Solving the CMAB problem implies maximizing the total outcome from given choices, as well as ensuring fair repartition of this outcome among all users. In this article, we theoretically and numerically demonstrated that carefully chosen polarization-entangled N-photon states can solve this problem for at least five players who must choose simultaneously between two choices or machines. Different behaviors, such as guaranteed fairness of outcome repartition only for even N, were identified between even and odd numbers of players, due to the fundamental properties of the polarization degree of freedom of photons. Nonetheless, in every case, the properties of these states imply that the best strategy for each player is to attempt actively to reach a measurement configuration that corresponds to a common optimum for all players simultaneously, thus achieving evolutionary stability. These quantum states are derived from a few basic rules, such as global invariance by rotation of the measurement bases of the players and symmetry of the state by permutation of the players and/or choices. Although a detailed derivation of the favorable states for any N has not been performed, we developed an algorithm that can verify whether any entangled state follows the set of rules described in this report. In addition, we demonstrated that, under basic assumptions, a fairly simple one-sided algorithm enables the users to correct any misalignment from an optimal configuration autonomously, without the need to trust a central entity. By