Item Listing Optimization for E-commerce Websites based on Diversity
Naoki Nishimura Kotaro Tanahashi Koji Suganuma Masamichi J. Miyama Masayuki Ohzeki

Application of Clumpy Torus Model to Broadband X-ray Spectra of
Two Seyfert 1 Galaxies: IC 4329A and NGC 7469
    Toshihiro Kawaguchi
    December 30, 2023
=============================================================================================================





For e-commerce websites, deciding the manner in which items are listed on webpages is an important issue because it can dramatically affect item sales.
One of the simplest strategies of listing items to improve the overall sales is to do so in a descending order of sales or sales numbers.
However, in lists generated using this strategy, items with high similarity are often placed consecutively.
In other words, the generated item list might be biased toward a specific preference.
Therefore, this study employs penalties for items with high similarity being placed next to each other in the list and transforms the item listing problem to a quadratic assignment problem (QAP).
The QAP is well-known as an NP-hard problem that cannot be solved in polynomial time.
To solve the QAP, we employ quantum annealing (QA), which exploits the quantum tunneling effect to efficiently solve an optimization problem.
In addition, we propose a problem decomposition method based on the structure of the item listing problem because the quantum annealer we use (i.e., D-Wave 2000Q) has a limited number of quantum bits.
Our experimental results indicate that we can create an item list that considers both sales and diversity.
In addition, we observe that using the problem decomposition method based on a problem structure can lead to a better solution with the quantum annealer in comparison with the existing problem decomposition method.
Keywords: item listing, e-commerce, quadratic assignment problem, quantum annealing, D-Wave, problem decomposition




§ INTRODUCTION

Several companies have recently started operating e-commerce websites to sell their items and services to the public considering the widespread use of the internet.
For these companies, deciding on the order in which items are listed on their website's pages is important because this ordering has the potential to dramatically affect the sales of their items or services.
Figure <ref> shows a snapshot of a hotel reservation website. 
This is an example of the items being listed on an e-commerce page.
On this website, hotels at different locations are listed in the order of popularity from top to bottom.
These sorted items for display on webpages are collectively referred to as an item list.



To improve sales on e-commerce websites, placing items in the descending order of sales or sales numbers is a simple strategy for determining the list order of items <cit.>.
In addition, the sales of an item can be estimated when it is placed at different positions in the item list and the position of each product is determined to maximize the total sales estimate.
In particular, if s_ij is the estimated sales of an item i ∈ I when it is placed in a position j ∈ J, then the total sales for all items can be maximized by solving the following integer programming problem <cit.>:


    [                                 ∑_i ∈ I∑_j ∈ J s_ijx_ij    ∑_i ∈ Ix_ij = 1,  j ∈ J,    ∑_j ∈ Jx_ij = 1,  i ∈ I, x_ij∈{0,1},  i ∈ I,  j ∈ J. ]

where x_ij is a binary variable that indicates whether or not to assign item i to position j. 
The abovementioned constraints ensure that only one item is allocated to each position, and only one position is allocated to each item.
In this study,


    S( x)=∑_i ∈ I∑_j ∈ Js_ijx_ij

is referred to as the sales term for x=(x_11,x_12,⋯).
This problem can be interpreted as a network flow problem, and there exists an efficient technique to solve such a problem in polynomial time.
Furthermore, the solution obtained by solving this network flow problem with x_ij∈ [0,1] coincides with the solution of the abovementioned integer programming problem <cit.>.

However, in the case of the list of items generated using such a strategy, the relationship between the different objects is ignored because the sales of each item s_ij are considered independently.
For example, let us assume that customers visit an e-commerce website and browse the page of a particular item group.
If the relationships among different items are not considered while placing items in an item list, several items with high similarities can possibly be placed close to each other, thereby reducing the value of the item list for customers in terms of item diversity.
Considering this, several attempts have been made to include item diversity in item recommendation lists for users to ensure that they find the recommendation lists useful <cit.>.
In these previous studies, the measures of diversity in the item recommendation lists for customers were improved by solving the maximum matching problem of the bipartite graph obtained after Top-N recommendation.

The present study includes diversity in the item list targeting the entire user base using a proposed method of generating an item list by imposing penalties when items with a high similarity are placed at adjacent locations.
Thus, considering both sales and diversity, the item list generation problem can be formulated as a quadratic assignment problem (QAP) as detailed below.

We employ quantum annealing (QA) herein to solve the QAP <cit.>.
An optimization problem formulated with discrete variables can be efficiently solved using the Ising model or a quadratic unconstrained binary optimization problem (QUBO) because of the introduction of the quantum tunneling effect by QA.
Currently, the protocol of QA is artificially realized in an actual quantum device known as a quantum annealer <cit.>.
The quantum annealer has been tested for numerous applications, including portfolio optimization <cit.>, protein folding simulation  <cit.>, molecular similarity problem <cit.>, computational biology <cit.>, job-shop scheduling <cit.>, traffic optimization <cit.>, election forecasting <cit.>, machine learning <cit.>, and for automated guided vehicles in plants <cit.>.
In addition, several other studies have been conducted to efficiently solve various problems using the quantum annealer <cit.>.

In particular, our problem can be solved using the quantum annealer by formulating our QAP as a QUBO.
However, a QUBO for such a large number of items cannot be directly solved in one instance with the current state-of-the-art quantum annealer, namely D-Wave 2000Q, because the D-Wave 2000Q employs the chimera graph. 
The physical qubits available on D-Wave 2000Q are less than 2048 because the qubits might have defects. 
In addition, the connection between the physical qubits is sparse and limited on the chimera graph.
Thus, several embedding techniques have been proposed; however, the number of logical qubits available to represent the optimization problems to be solved is drastically reduced <cit.>.
To address this issue, a heuristic method has been proposed to solve a large-sized problem using a limited number of hardware bits.
D-Wave Systems, which is the manufacturer of D-Wave 2000Q, has developed an open-source software  <cit.> that solves a large-sized problem by dividing it into small subproblems.
However, the decomposed QAP might not necessarily lead to feasible solutions because  selects the subset of variables in the order of the energy impact of each variable for division of a problem.
Furthermore, in the literature <cit.>, the division of an original problem into subproblems based on its structure is a promising method to efficiently solve large-sized problems using D-Wave 2000Q.
Considering this, we propose herein a method of obtaining better objective values for the QAP problem compared to those of the existing method <cit.> in the same calculation time by decomposing the problem based on the set of items and positions that they can be assigned to in an item list.
In addition, we assess the performance of the proposed method using the actual access log of a hotel reservation website.

The primary contributions of our study are summarized as follows:


  * We propose a method of creating item lists on an e-commerce website as a QAP considering the sales and diversity of the items.

  * We convert the QAP to the QUBO to solve the abovementioned problem with D-wave 2000Q.

  * We propose a decomposition technique exploiting the structure of the item list.





§ MODELS



 §.§ Formulating the Item Listing Optimization Problem as a QAP

We introduce the diversity term in our proposed model to add diversity in the item list.
We particularly calculate the similarity f_ii' for pairs of items i, i' ∈ I to introduce diversity in the item list.
The diversity of the item list (i.e., the diversity term) is defined as the negative value of the summation of the items’ similarity degree f_ii' for overall adjacent items:


    D( x)= -∑_i ∈ I∑_i' ∈ I∑_j ∈ J∑_j' ∈ Jf_ii'd_jj'x_ijx_i'j' .

where d_jj' is the adjacent flag of the position j and j'; d_jj'=1 is for the adjacent positions; and d_jj'=0 is for the non-adjacent positions.
The value of function D( x) decreases because the high-similarity items are adjacent.
We solve a multi-objective optimization problem based on two values: the sales of individual products S( x) and diversity of the item list D( x).
Let w be a parameter used to determine the penalty for listing items with high similarity. This problem is formulated as follows:


    [                                                                                ∑_i ∈ I∑_j ∈ J s_ijx_ij - w ∑_i ∈ I∑_i' ∈ I∑_j ∈ J∑_j' ∈ Jf_ii'd_jj'x_ijx_i'j'                                                       ∑_i ∈ Ix_ij = 1,  j ∈ J,                                                       ∑_j ∈ Jx_ij = 1,  i ∈ I,                                                    x_ij∈{0,1},  i ∈ I,  j ∈ J. ]

Two methods for calculating the similarity f_ii' are introduced: explicit and implicit expressions.
In the explicit expression, semantic features, such as product category and average price, can be quantified using the distances between the feature vectors that can be calculated. The smaller the feature distance, the higher the similarity f_ii' between the items.
In the implicit expression, the higher the co-browsing number of an item (i.e., the number of times that the items were viewed in the same session by the same user on the website), the higher the similarity f_ii'.

The advantage of the first approach is that the interpretation of the result is straightforward; nevertheless, it suffers from a disadvantage in that appropriate semantic features must be created and quantified.
In contrast, the advantage of the second approach is that it involves easy calculations and can consider various information reflecting customer behavior; however, its disadvantage is that semantic interpretation might be difficult.
Nevertheless, we employed herein the latter method of counting the co-browsing number to calculate similarity f_ii'.
Section <ref> describes the detailed calculation for this.
As previously specified, the optimization problem (<ref>) is a QAP.

The QAP is well-known as an NP-hard problem that cannot be solved in polynomial time <cit.>.



 §.§ Formulating the Item Listing Optimization Problem as a QUBO

We utilize QA to solve our optimization problem.
The optimization problem must be expressed in the form of a QUBO to use QA as a solver. QUBO is given as follows <cit.>:


    [              x^T Q  x x∈{0,1}^N, ]

where Q ∈ℝ^N × N.
Thus, our optimization problem can be transformed into a QUBO by employing a penalty function for violating constraints and adding this penalty function to the objective function:

    [                                                                                 -∑_i ∈ I∑_j ∈ J s_ijx_ij + w ∑_i ∈ I∑_i' ∈ I∑_j ∈ J∑_j' ∈ Jf_ii'd_jj'x_ijx_i'j'                   +M ( ∑_i ∈ I(∑_j ∈ Jx_ij-1)^2  + ∑_j ∈ J( ∑_i ∈ Ix_ij-1 )^2 )                                                     x_ij∈{0,1},  i ∈ I,  j ∈ J. ]

where M is a parameter used to prevent the violation of the constraint conditions. 
This is ensured by setting an appropriate value for M.
In theory, M should take an extremely large value.
However, we cannot set M to such a large value because of the limitations of the current version of the quantum annealer used (i.e., D-Wave 2000Q).
Thus, for simplicity, we present M to the size of the largest element of the absolute value of Q in (<ref>).



 §.§ Decomposition Methods for Item Listing Problems

 is a software tool released by D-Wave Systems that enables solving a QUBO larger than one that can be processed using D-Wave 2000Q <cit.>.  is essentially a decomposing solver that divides a large problem into smaller parts, which can then be solved by D-Wave 2000Q.
Thus, when a large QUBO is inputted,  divides the problem and sends each part of the problem independently to D-Wave 2000Q for calculation to obtain partial solutions.
This process is repeated by selecting different parts of the problem using the tabu search until solution improvement stops.
See <cit.> for the detailed algorithm of .
Furthermore,  selects the subset of variables in the order of the energy impact of each variable for division of a problem.
However, in some cases, no feasible solution can be obtained when the target variables are extracted, regardless of the structure of the original problem.

Therefore, we focus herein on the structure of the assignment problem and propose a method to extract problems with feasible solutions.
Particularly in the case of an assignment problem, one condition involves each item being necessarily assigned to one position and another condition, in which each position is necessarily assigned to one item.
Therefore, while dividing the problem, we have to select variables with candidate combinations of items and positions that are already assigned.
Figure <ref> shows an example of the decomposition.


The original problem can be decomposed as follows if the number of items in the original problem is N_ org and the number of items solved by a partial problem is N_ sub:



  * Let 𝒩_s be the set of N_ sub items extracted from N_ org items.

  * Let 𝒫_s be the set of positions of items of 𝒩_s.

  * Let 𝒩_s ×𝒫_s be the variables of the decomposed problem.

This procedure involves _N_ org^2C_N_ sub^2 candidates for variable combinations in the selection of the subproblems; however, the number of solution candidates can be reduced to _N_ orgC_N_ sub exploiting the structure of the item list.

In practice, it is most important to determine the order in which items are listed in the upper positions of the item list because they are the items that are browsed most often.
Therefore, it is effective to solve the entire list as an integer programming problem as in Problem (<ref>) first, then only resolve the particularly important upper positions of the list using the QAP (<ref>).




§ RESULTS AND DISCUSSION



 §.§ Experimental Setup

For our experiments, we used the actual access log data of the online hotel reservation site Jalan [< https://www.jalan.net/en/japan_hotels_ryokan/>].
On this e-commerce website, a hotel list is created daily based on each area in Japan and the number of guests, including adults and children, that the hotels can accommodate in their rooms.
The access log includes the date and the time the customer accessed the item list screen, position of each item when the item list screen was accessed, and information on the hotel at which the customer made a reservation.
We estimated the sales s_ij and similarity f_ii' for the top 10 accessed areas on the hotel reservation website using the access log for the past six months.
Various methods can be considered for estimating s_ij.
The similarity f_ii' was estimated using the log of the co-browsed items in the same session of the customer.
Furthermore, s_ij and f_ii' were normalized such that their average was 0, and the standard deviation was 1 for each item list.

We conducted two experiments in this study:


  * evaluating the effect of the diversity term: comparison of solutions when the diversity control parameter is changed for the QAP (<ref>);

  * performance evaluation of problem decomposition: comparison of the objective values when the structure of the item list is considered for the  problem decomposition.

As previously specified, we used D-Wave 2000Q (DW_2000Q_VFYC_2) for our experiments.
Coupler strengths mapping logical to physical couplers with two physical couplers connecting each pair of logical qubits were set as 3.0.
Table <ref> lists the values set for the parameters of D-Wave 2000Q and .
, , and  represent the number of requests for problems, time per annealing, and number of times the main loop of the algorithm is repeated with no change in the optimal value before stopping, respectively. <cit.> 





 §.§ Effect of the Diversity Term

D-Wave 2000Q has less than 2048 qubits because the qubits typically have defects. In addition, as previously specified, the connection between the physical qubits is sparse and limited on the chimera graph, in which D-Wave 2000Q has been based on. Thus, we can consider the problem with eight items per subproblem because complete graph embedding can be applied to arbitrary problem graphs with less than 64 logical variables.
Therefore, we first compare the sales term S( x) and the diversity term D( x) obtained by solving the QAP (<ref>) by changing the diversity control parameter w in the problem when the number of items is 6 and 8.
The obtained solution will be the same as that obtained via solving (<ref>) if w is 0.

In Figures <ref> and <ref>, the horizontal axes represent the value of w in (<ref>); the right hand side in the figures indicates that the larger the similarity between the similar items listed together, the larger the penalty.
In addition, the vertical axes in Figures <ref> and <ref> represent the values of S( x) and D( x), respectively. The average value of all solutions is approximately 0 because s_ij and f_ii' are normalized for each area.
Figures <ref> and <ref> depict plots of S( x) and D( x), respectively, for three typical areas X, Y, and Z from among the 10 areas for the problem with eight items.
Figures <ref> and <ref> are plotted by aggregating the typical values for the 10 areas, wherein the marker and the bar represent the average and standard error values, respectively.



Figure <ref> shows that the average of the sales term S( x) is higher than 0 for any w. In addition, S( x) gradually decreases as w increases.
In other words, a trade-off relationship exists between increasing the sales term S( x) and ensuring that similar items are not placed adjacent to each other in the item list.
Furthermore, the decrease in the sales term value slows down as w increases.



Figure <ref> shows that when diversity is not considered for item listing, the diversity term D( x) is lower than the average value of 0 and increases with the increasing w.
Furthermore, the increase in the diversity term D( x) gradually slows down as w increases.
Note that the behavior of S( x) and D( x) based on w is the same whether or not the number of items is 6 or 8.

Figure <ref> shows the positional relationships between the top eight hotels and their type in area X of Figures <ref> and <ref>.
The top eight hotels in area X particularly include four city hotels and four budget hotels, which can also be classified as in the north or south region of area X; hence, it was chosen as an example.



Figure <ref> represents the actual hotel lists in area X obtained by solving (<ref>) for different values of w.
The top five hotels in the northern area are consecutively listed when w is 0. The city hotels are listed in succession among the top four; thus, the list is biased.
In contrast, only the top two hotels are the same in terms of both area and hotel type when w is between 0.1 and 0.6, that is, when diversity is considered.
Furthermore, no similarity exists between the consecutively listed hotels both in terms of area and hotel type when w is larger than 0.7.
Thus, despite not using semantic information (e.g., area and hotel type) in the actual calculations, a semantically diverse item list was created using the number of co-browsers as the similarity measure.





 §.§ Performance Evaluation of Problem Decomposition

We compared the method of extracting partial problems by  and our method of extracting partial problems considering the problem structure of the item list.
Our experiments were performed by comparing the objective values of (<ref>) when solving the problem of 12, 16, 20 and 24 items using each method.

Table <ref> lists the average of the objective values for the 10 areas obtained by solving the problem for different numbers of items and the gap between the objective values of the proposed and original  methods. 
The objective of the problem was to minimize the objective value, indicating that the smaller the objective value, the better the performance.
The gap between the objective values of the proposed and original  methods increased as the number of items increased. 
The effectiveness of the proposed method also increased.
We conducted a Wilcoxon signed-rank test on the difference between the objective values of the original  and proposed methods (Bonferroni correction was performed on the number of items).
Consequently, the null hypothesis in this case was rejected at a significance level of 5% (p-value = 0.020).



In terms of application, our proposed method is not only limited to the item list optimization problem, but can also be widely applied to other problems involving similar constraints, such as the assignment problem (<ref>).
For example, our method can be applied to the traveling salesman problem, which typically includes two constants A, B along with the following penalty terms:

    A ∑_i ∈ I( ∑_j ∈ Jx_ij-1) ^2  + B ∑_j ∈ J( ∑_i ∈ Ix_ij-1) ^2.





§ CONCLUSION

This study proposed a method of creating item lists on an e-commerce website as a QAP considering item sales and diversity.
We converted the QAP to a QUBO such that it can be directly solved by the quantum annealer, D-Wave 2000Q.
Direct manipulation to solve the resulting QUBO was not possible in the case with a large number of items because of the limited number of qubits available in the current version of the quantum annealer and the restriction on specifying connections between the qubits.
Therefore, we proposed a decomposition technique exploiting the structure of the problem.
The original large problem was divided into several subproblems, which can eventually be solved by D-Wave 2000Q individually.
Our experiments using actual real-world data demonstrated the efficiency of our proposed approach.
A remarkable observation made from the experimental results was that the output item list changed based on the diversity control parameter.
Our formulation led to the antiferromagnetic Ising model with a random field.
The resulting lists were “aligned" along the random field when the diversity control parameter was small.
In contrast, increasing the diversity control parameter eliminated the order in the item list and introduced diversity.

However, note that our research has some limitations.
The item list created by our method will not work well when there are insufficient data to calculate the sales s_ij and similarity f_ii'; because of which, no reliable estimate values existed. 
Having access logs when items were placed at various positions was important in obtaining good estimate values for s_ij and f_ii'.
The determination of the diversity consideration parameter w was also a limitation. 
w depends on the scale of diversity a customer of an e-commerce service wants; thus, it should be adjusted by changing w and monitoring performance, which is cumbersome to implementation in real-world scenarios.

As the number of qubits available in the quantum annealer increases in the future, our method for division of a large-sized problem into smaller subproblems will become more useful.
Our experiments clearly showed that our method performed better than  when the size of the problems became large.
The structure of our problem is similar to the traveling salesman problem and the scheduling problem in that it takes the form of a QUBO with two quadratic functions based on two constraints (i.e., one for distance or time and the other for list locations or tasks).

Thus, our method has a wide range of applications involving optimization problems that can be solved via QUBO formulation, such as those using the QA, D-Wave 2000Q, and other types of QUBO solvers.
Our results indicate that not only the evolution of hardware devices, but also the development of better software based on the structure of problems are essential for future QA applications.



§ CONFLICT OF INTEREST STATEMENT

The authors declare that the research was conducted without any commercial or financial association that could be construed as potential conflicts of interest.



§ AUTHOR CONTRIBUTIONS

NN contributed to the conception and design of the study, performed all the experiments, and wrote the first draft of the manuscript. KT implemented the program of our problem decomposition method. MO and MJM  contributed to manuscript revision. All authors discussed this study, and then reviewed and approved the final version of the manuscript.



§ ACKNOWLEDGMENTS

The authors thank Recruit Lifestyle Co., Ltd. and Recruit Communications Co., Ltd. for its support in this exploratory research project.



§ DATA AVAILABILITY STATEMENT

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.


32
abdel2018comprehensive
Abdel-Basset, M., Manogaran, G., Rashad, H., and Zaied, A. N. H. (2018).
A comprehensive review of quadratic assignment problem: variants, hybrids and applications.
Journal of Ambient Intelligence and Humanized Computing,
1–24.

adomavicius2011maximizing
Adomavicius, G. and Kwon, Y. (2011).
Maximizing aggregate recommendation diversity: A graph-theoretic approach.
In Proceedings of the 1st International Workshop on Novelty and Diversity in Recommender Systems (DiveRS 2011) 
(pp. 3–10).

anstreicher2003recent
Anstreicher, K. M. (2003).
Recent advances in the solution of quadratic assignment problems.
Mathematical Programming 97,
27–42.

antikacioglu2017post
Antikacioglu, A. and Ravi, R. (2017).
Post processing recommender systems for diversity.
In Proceedings of the 23rd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (pp. 707–716).

Arai2018nn
Arai, S., Ohzeki, M., and Tanaka, K. (2018).
Deep neural network detects quantum phase transition.
Journal of the Physical Society of Japan 87,
033001.

Dwave2010b
Berkley, A. J., Johnson, M. W., Bunyk, P., Harris, R., Johansson, J., Lanting, T., Ladizinsky, E., Tolkacheva, E., Amin, M.H.S., and Rose, G. (2010).
A scalable readout system for a superconducting adiabatic quantum optimization system.
Superconductor Science and Technology 23,
105014.

booth2017partitioning
Booth, M., Reinhardt, S. P., and Roy, A. (2017).
Partitioning optimization problems for hybrid classical/quantum execution.
<https://www.dwavesys.com/resources/publications>.

boothby2016fast
Boothby, T., King, A. D., and Roy, A. (2016).
Fast clique minor generation in chimera qubit connectivity graphs.
Quantum Information Processing 15, 
495–508.

Dwave2014
Bunyk, P., Hoskinson, E. M., Johnson, M. W., Tolkacheva, E., Altomare, F., Berkley, A. J., Harris, R.,  Hilton, J. P., Lanting, T., Przybysz, A. J., and Whittaker, J. (2014).
Architectural considerations in the design of a superconducting quantum annealing processor.
IEEE Transactions on Applied Superconductivity 24, 
1–10.

Crawford2016
Crawford, D., Levit, A., Ghadermarzy, N., Oberoi, J. S., and Ronagh, P. (2016).
Reinforcement learning using quantum boltzmann machines.
arXiv preprint,
arXiv:1612.05695.

dwave2017solver
D-Wave Systems Inc. (2017).
D-wave solver properties and parameters reference.
<https://www.dwavesys.com/resources/publications>.

Dwave2010c
Harris, R., Johnson, M. W., Lanting, T., Berkley, A. J., Johansson, J., Bunyk, P., Tolkacheva, E., Ladizinsky, E., Ladizinsky, N., Oh, T., Cioata, F., Perminov, I., Spear, P., Enderud, C., Rich, C., Uchaikin, S., Thom, M. C., Chapple, E. M., Wang, J., Wilson, B., Amin, M. H. S., Dickson, N., Karimi, K., Macready, B., Truncik, C. J. S., and Rose, G. (2010).
Experimental investigation of an eight-qubit unit cell in a superconducting optimization processor.
Physical Review B 82,
024511.

Henderson2018
Henderson, M., Novak, J., and Cook, T. (2018).
Leveraging adiabatic quantum computation for election forecasting.
arXiv preprint,
arXiv:1802.00069.

Hernandez2017
Hernandez, M. and Aramon, M. (2017).
Enhancing quantum annealing performance for the molecular similarity problem.
Quantum Information Processing 16, 
133.

Dwave2010a
Johnson, M. W., Bunyk, P., Maibaum, F., Tolkacheva, E., Berkley, A. J., Chapple, E. M., Harris, R., Johansson, J., Lanting, T., Perminov, I., Ladizinsky, E., Oh, T., and Rose, G. (2010).
A scalable control system for a superconducting adiabatic quantum optimization processor.
Superconductor Science and Technology 23, 
065004.

Kadowaki1998
Kadowaki, T. and Nishimori, H. (1998).
Quantum annealing in the transverse Ising model.
Physical Review E 58, 
5355.

Khoshaman2018
Khoshaman, A., Vinci, W., Denis, B., Andriyash, E., and Amin, M. H. (2018).
Quantum variational autoencoder.
Quantum Science and Technology 4, 
014001.

Richard2018
Li, R. Y., Felice, R. D., Rohs, R., and Lidar, D. A. (2018).
Quantum annealing versus classical machine learning applied to a simplified computational biology problem.
npj Quantum Information 4, 
14.

long2014relevance
Long, B. and Chang, Y. (2014).
Relevance ranking for vertical search engines.
Newnes.


Lucas2014
Lucas, A. (2014).
Ising formulations of many NP problems.
Frontiers in Physics 2, 
5.

Neukart2017
Neukart, F., Compostella, G., Seidel, C., Von Dollen, D., Yarkoni, S., and Parney, B. (2017).
Traffic flow optimization using a quantum annealer.
Frontiers in ICT 4, 
29.

Neukart2018
Neukart, F., Dollen, D. V., Seidel, C., and Compostella, G. (2018).
Quantum-enhanced reinforcement learning for finite-episode games with discrete state spaces.
Frontiers in Physics 5,
71.

Ohzeki2019
Ohzeki, M., Miki, A., Miyama, M. J., and Terabe, M. (2019).
Control of automated guided vehicles without collision by quantum annealer and digital devices.
arXiv preprint,
arXiv:1812.01532.

Ohzeki2018
Ohzeki, M., Okada, S., Terabe, M., and Taguchi, S. (2018).
Optimization of neural networks via finite-value quantum fluctuations.
Scientific Reports 8, 
9950.

Ohzeki2018NOLTA
Ohzeki, M., Takahashi, C., Okada, S., Terabe, M., Taguchi, S., and Tanaka, K. (2018).
Quantum annealing: next-generation computation and how to implement it when information is missing.
Nonlinear Theory and Its Applications, IEICE 9, 
392–405.



Okada2019
Okada, S., Ohzeki, M., Terabe, M., and Taguchi, S. (2019).
Improving solutions by embedding larger subproblems in a d-wave quantum annealer.
Scientific Reports 9, 
2098.

Perdomo2012
Perdomo-Ortiz, A., Dickson, N., Drew-Brook, M., Rose, G., and Aspuru-Guzik, A. (2012).
Finding low-energy conformations of lattice protein models by quantum annealing.
Scientific Reports 2, 
571.


Rosenberg2016
Rosenberg, G., Haghnegahdar, P., Goddard, P., Carr, P., Wu, K., and Prado, M. L. (2016).
Solving the optimal trading trajectory problem using a quantum annealer.
IEEE Journal of Selected Topics in Signal Processing 10,
1053–1060.

Takahashi2018
Takahashi, C., Ohzeki, M., Okada, S., Terabe, M., Taguchi, S., and Tanaka, K. (2018).
Statistical-mechanical analysis of compressed sensing for hamiltonian estimation of ising spin glass.
Journal of the Physical Society of Japan 87, 
074001.


vazirani2013approximation
Vazirani, V. V. (2013).
Approximation algorithms.
Springer Science & Business Media.


Venturelli2015
Venturelli, D., Dominic, J. M., and Rojo., G. (2015).
Quantum annealing implementation of job-shop scheduling.
arXiv preprint,
arXiv:1506.08479.


wang2016beyond
Wang, Y., Yin, D., Jie, L., Wang, P., Yamada, M., Chang, Y., and Mei, Q. (2016).
Beyond ranking: Optimizing whole-page presentation.
In Proceedings of the Ninth ACM International Conference on Web Search and Data Mining (pp. 103–112).



