Journal of Applied Mathematics

Volume 2014 (2014), Article ID 818479, 8 pages

http://dx.doi.org/10.1155/2014/818479

Research Article

A Novel Algorithm of Quantum Random Walk in Server Traffic Control and Task Scheduling

Dong Yumin¹ and Xiao Shufen²

¹Network Center, Qingdao Technological University, Qingdao 266033, China

²College of Automobile and Transportation, Qingdao Technological University, Qingdao 266033, China

Received 18 January 2014; Accepted 19 February 2014; Published 7 April 2014

Academic Editor: Feng Gao

Copyright © 2014 Dong Yumin and Xiao Shufen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A quantum random walk optimization model and algorithm in network cluster server traffic control and task scheduling is proposed. In order to solve the problem of server load balancing, we research and discuss the distribution theory of energy field in quantum mechanics and apply it to data clustering. We introduce the method of random walk and illuminate what the quantum random walk is. Here, we mainly research the standard model of one-dimensional quantum random walk. For the data clustering problem of high dimensional space, we can decompose one -dimensional quantum random walk into one-dimensional quantum random walk. In the end of the paper, we compare the quantum random walk optimization method with GA (genetic algorithm), ACO (ant colony optimization), and SAA (simulated annealing algorithm). In the same time, we prove its validity and rationality by the experiment of analog and simulation.

1. Introduction

The server cluster technology may be connecting multiple independent servers, and, in the same time, it must provide services as a whole by a cluster. In the server cluster, how to solve the problem of server traffic control and task scheduling is very important.

In order to reduce the access time, optimize the overall performance and achieve parallel program in a high efficiency; the task request must be allocated to each on the server. So, load balancing mechanism is the core of cluster technologies.

In the literature [1], it expands the analogies employed on the development of quantum evolutionary algorithms by putting forward quantum-inspired Hadamard walks, called QHW. In order to solve combinatorial optimization problems, a quantum evolutionary algorithm, abbreviated HQEA, is proposed. From the results of the experiments carried out on the knapsack problem, HQEA performs noticeably better than a conventional genetic algorithm, in terms of convergence speed and accuracy. The literature [2] explores how a spectral technique suggested through coined quantum walks can be used to differentiate between graphs which are cospectral as for standard matrix representations. This algorithm runs in polynomial time; it can differentiate many graphs for there is no subexponential time algorithm which is proven to be able to differentiate between them.

By the literature [3], they propose a quantum algorithm to evaluate formulas by an extended gate set, including two- and three-bit binary gates. This algorithm is more optimal on read, once formulas for that each gate’s inputs are balanced in a certain sense. It describes a very compact triaxial instrument in the literature [4]. The triaxial instrument is based on a rhombic dodecahedral geometry that can accommodate three nonplanar ring light paths with orthogonal sensing axes. Component count can be substantially reduced by a discharge of layout to use a single cathode and two anodes running all three axes in balanced plasma currents. Two Monte Carlo-based approaches to assess parameter uncertainty with complex hydrologic models are considered in the literature [5]. The importance sampling has been carried out in the generalized likelihood uncertainty estimation framework by Beven and Binley. The metropolis algorithm is different from importance sampling which uses a random walk that adapts to the real probability distribution describing parameter uncertainty.

Because existing search protocols for unstructured peer-to-peer systems to create huge burden on communications or cause long response time and result in unreliable performance. In the literature [6], in order to discover service providers, it reports that an important function of a peer-to-peer system is a distributed message relaying. They present an incentive mechanism which not only relieves the free-riding problem but also accomplishes good system efficiency in message relaying for peer discovery. The passed along message propagation process is promised rewards in the mechanism.

In the literature [7], it analyzes the discrete-time quantum walk by separating the quantum evolution equation into Markovian and its interference terms. Because of this separation, it is possible to show analytically which quadratic increase in the variation of the position of quantum walker with time is a direct aftermath of the coherence of the quantum evolution. As expected, the variation is shown to increase linearly with time, if the evolution is decoherent, as in the classical case. Moreover, it shows that the system has an evolving operator analogous to which of a resonant quantum kicked rotor. At the same time, the rotator can be described by evolution of the quantum walker.

Quantum random walks on a graph, which is analogous to classical stochastic walk, form the basis for many of the recent quantum algorithms that promise to obviously outperform existing classical random walk algorithms. A number of studies have been done on the many applications of quantum random walk to some important computing problems. There are two kinds of quantum random walk algorithms: continuous-time and discrete-time. It is reported that a quantum arithmetic is defined by a sequence of the operations that runs on an actual model of quantum computation in the literature [8]. It proposes quantum circuit designs for both kinds of random walk algorithms which operate on various graphs. It considers two important problems to which random walk arithmetic are applicable: the triangle finding problem and binary tree problem. Because of it a few research works that are related to quantum random walk circuit design on graphs exist; the circuit designs they present here are the first of their kind. At the same time, they also provide an estimate of the quantum cost of the circuits of quantum systems. And it is based on the number of execution cycles and quantum operations.

In the literature [9], the natural random walk causing Brownian motion occurs to be always biased in a very delicate way: emphasizing some possibilities by only approximative maximal uncertainty principle. It introduces a new method of stochastic model, and they use the merely maximizing entropy choice of transition probabilities.

Berry and Wang show numerically that a discrete time quantum random walk of two irrelevant particles is able to differentiate some nonisomorphic powerfully regular graphs from the same household in the literature [10]. They analytically show how it is possible for the walks to differentiate such graphs, while the continuous time quantum walks of two irrelevant particles cannot.

It is reported that the quantum walks are quantum mechanical theory analog of random walks in the literature [11]. By traversing the edges of a graph, a quantum “walker” progresses between initial and final states. They present a hybrid model for general quantum computing in which a quantum walker gets discrete steps of continuous evolution.

Effective server traffic control can extend the “capacity” of the server, and the task scheduling can improve system throughput. In early research methods of it, such as Min-Min algorithm, Max-Min algorithm, genetic algorithm (GA), round robin (RR), simulated annealing algorithm (SAA), dynamic feedback algorithm (DFA), and ant colony optimization algorithm (ACO). These arithmetics have some improvements in different degree on the task scheduling.

But these algorithms have this or that problem, such as local premature problem and divergence problem.

In order to overcome the instability above the algorithms, the quantum random walk algorithm is proposed, and it is proved better than above GA, ACO, and SAA by simulation experiments.

2. Quantum Random Walk

2.1. Random Walk

Random walk is a mathematical method to study the formation of trajectory by a random sequence of continuous; it is not only a means to study mathematics but also a basic tool in the natural sciences. Any stage of the random walk behavior is not limited to previous history of migration; the process is also called Markov process. Random walk can be simply described as follows.

Suppose in a straight line, there is a moving particle; it is at the origin to move left or right one unit of distance; the probability is and , respectively, each time the particle in accordance with the probability to move a unit distance to the left or right. Here we assume that the probability of the particle is equal to the left or right; that is, ; random variable can be used to represent the probability ; its value is as follows:

If the particles every moment in a straight line position constitute an independent identical distribution of random variables sequence, denoted is a sequence of independent and identically distributed variables on meet , . is its first terms and is called a random walk. After the particles move steps, the probability of it being found in position is

Among them, take only integer, ; other cases were 0.

Although the classical random walk has a broad application, but compared with the quantum random walk, it feels much ashamed of its inferiority.

Quantum random walk is a quantum computing model proposed in recent years; scholars have also become increasingly interested in research.

2.2. Quantum Random Walk

For discrete quantum random walk, the system added an extra degree of freedom; some literature defines it as chirality that can build an adaptation to global local unitary process. This walk is also called quantum Hadamard walk. The only possible remained unchanged in the global process of unitary transformation is the only mobile operator between adjacent lattice points to the left or the right.

Significantly different between quantum Hadamard random walk and classical quantum is the interference, the diffusion rate of quantum walk square magnitude faster than classical square. Due to the existence of quantum superposition states, in quantum random walk, position of the particle from the probability distribution may be seen; particles may be in several locations simultaneously with different probability. Quantum random walk process is accomplished by a unitary matrix transformation [12].

2.2.1. One-Dimensional Quantum Random Walk

The random walk model sets up corresponding quantum algorithms and does quantum information processing. People commonly used coined quantum walks; it corresponds to Hilbert space which can be expressed as follows:

where is the random walk of grid space; it corresponds to a classic case of -degree regular graphs , which is a coin flip operator space (coin space). The total unitary evolution matrix is by the two independent parts, namely, flipping a coin and conditional replacement

The first step of the quantum walk is to perform a rotation operation in coin space, equivalent to the classical random walk in a coin toss, through this operation to get a coin superposition state [13]. Then, the replacement operator makes the particles by a coin to decide an edge vertex adjacent to move to the next. Starting from the initial state , repeat the walk after steps and obtain the probability distribution of each vertex as follows:

Quantum random walks of a variety of ways common are one-dimensional linear walk, ring walk, hypercube walk, and so forth. For one-dimensional linear walk, , replacement operator applied to the base is expressed as

Starting from the initial state , continuous action ; after each step, the distribution of every point is as follows:

Not only is the distribution situation different with the classical random walk, but one-dimensional linear walk is also higher and faster than the classical random walk in the diffusion velocity.

2.2.2. -Dimensional Hypercube Quantum Random Walk

For the -dimensional hypercube quantum random walk, it has vertices each vertex can be marked by an binary string,

. Each vertex is degrees; therefore, coin space , using to mark the coin space basis vectors, ; it indicates the direction of the next step; each corresponds to a -dimensional vector ,

, the number of bit is 1 and that of the other bits is 0. The replacement operator walk on the hypercube is applied to the base vectors and can be expressed as

It represents two points marked by quantum bits; only when they have only one bit is not the same when they are connected directly via a side (e.g., 001101 and 011101 are communicating). There are many ways to select hypercube walking coin flip operators; however, the following form is usually taken to maintain certain symmetry:

The operator of such forms has the characteristics that in all directions are a permutation invariant; it retains hypercube replacement invariance. Such a form is a commonly used Grover diffusion operator selection

Among them, are equally weighted superposition states in all directions. Grover operator is one of the permutation invariant operators, which is the farthest with unit transformation ; It will effectively mix all of the any given initial state into the superposition of them. The total random walk evolution operator can be expressed as

Define the Hamming distance, starting from any point to another point of the minimum number of edges experience (i.e., the required number of steps), with . Marking each vertex string the number of “1” is called Hamming weights, for example, Hamming weight 010 to 1. Starting from , any Hamming weighing the same point total can reach at the same number of steps. When the coins are formed to the symmetry type such as (10), it has the same probability starting from this point to reaching the point with the same weight. This allows us to put all Hamming weight of the same point “accumulation” to a point, thereby reducing the symmetry of the random walk; the walk on the hypercube is becoming walk in a straight line. It is noteworthy that this walk difference on the straight is not unbiased; it differs from the previously discussed one-dimensional linear walk. The number of total vertices after walking on the variable linear is ; is the hypercube dimension.

By putting a hypercube walk into walk on straight line, many problems can be resolved to simplify and get results. Moor and Russell found that when , random walk is a balanced distribution. Kempe through the research on the hitting time found hypercube quantum random walk to reach vertical angles of time relative to the classic case which is an exponential acceleration; this shows that the quantum random walk has the potential to make quantum algorithm acceleration.

3. The Model of Server Traffic Control and Task Scheduling

In cluster services, the task scheduling can be described as follows: tasks need to be allocated to nodes (these nodes are the servers) with different handling capacity; the goal is finding an optimization schedule to minimize the total completion time. The system model is shown as follows.

We suppose there are nodes (or servers) and tasks. Every task should be assigned to only one node. We use denoting the nodes (or servers), in this paper, where denotes one of the nodes (or servers); expresses the current load, where expresses the current load of node . For instance, means that the node (server) has a current load of 0; in other words, the node is idle. Here tasks are expressed by, where is one of the tasks. A matrix is built between servers and tasks: , where is one of the elements. So, there are two states as follows:

where, , .

We use to express the time of processing on one task; in other words, the time of task processed on node . The processing time is denoted as follows:

where, , .

It is not difficult to see that is also an matrix.

Here, we define the optimal state occurring with these conditions: (a) the total system has a relative short time of processing; (b) the throughput of system is relatively larger in unite time. We can describe this state using the following equations:

where is the new task, is the current total load at the node, is the length of ready queue at node , is the average processing time at node , are constant, is a function which can show the ability of node processing. The system is on the optimal running state, when the capacity of processing tasks (or task scheduling) reaches the maximal matching at one node.

4. The Method of Task Scheduling Based on Quantum Random Walk (QRW) Clustering Algorithm

In the paper, we mainly research the standard model of one-dimensional quantum random walk. For the data clustering problem of high dimensional space, we can decompose one -dimensional quantum random walk into one-dimensional quantum random walk.

4.1. Clustering Algorithm Based on One-Dimensional Quantum Random Walk, Referred to as Quantum Random Walk Clustering Algorithm (QRWA)

Step 1. Assume an unlabeled data set , where each data point with m features.

Step 2. Each data point in the data set can be considered as a particle that transfers in the space according to the probability.

Step 3. Establish clustering algorithm based on the one-dimensional quantum random walk.

The clustering algorithm uses a distributed control strategy, that is, each data point of the data set only affected by its neighbor within the neighborhood. The neighbor of data points available -nearest neighbor method or the method of default scope of to determine and use indicates the set of neighbors of a data point in time.

Step 4. Calculate the probability for each data point transfer to all neighbors in the neighborhood ; the formula is as follows:

Among them, and , respectively, represent the degree of current and initial data points; similarly, and , respectively, represent the current and initial data point distances between and .

Step 5. Find the maximum transition probability and the neighbors of greatest probability of metastasis , ,  .

4.2. The Server Traffic Control Clustering Method of Quantum Random Walk

As previously mentioned, -dimensional quantum random walk is decomposed into one-dimensional quantum random walk. For each dimension, data points can only move a step left or right, or . Therefore, the maximum transition probability is mapped to the interval , so, the probability of transfer in the opposite direction is . When , and probability of metastasis in both directions is equal, available aforementioned Hadamard transforms as a coin matrix. However, in normal conditions, , therefore, the coin matrix is used in the algorithm is

It is easy to verify that the matrix is unitary matrix meeting the reversibility requirements of quantum mechanics.

Since the quantum random walk clustering algorithm will use consecutive transformation for initial state, thus, every time it changes left and right transfer step and to the original 1/2, that is, and , Then, conditional operator will press type structure:

As is known, in quantum mechanics, each one of superposition states can be seen as a position of particle and indicates the probability of finding the particle at this location. If repeatedly used transforms for initial state, then the resulting superposition state will contain more items; this increases possible appearing position of the particle, and this is not present in the classical random walk. It is these possible positions that increase the searching range of solution space and provide an opportunity for better results. To calculate the probability of multiple locations of particles and their appearance in the corresponding locations, a unitary operation is sufficient because of the quantum parallelism, but in the classical world, you need multiple operations to complete; it also reflects an aspect of quantum computing to accelerate the classical computing.

If the initial state of the particle is , then, after applying times transform , get superposition state which is

From (19) particles can be found not only with probability and at the same time appear in the probability of and appearing in another new position could be . At this point, if projection measurement of the superposition state , it will collapse to one of these three positions according to the probability; then, the component of particle in the dimension is updated with the following formula:

As data points are random walk in the space, its position and its nearest neighbor are constantly changing with time. Therefore, in the process of walking, the distance of the data points and the degree of it need to be recalculated. Repeat the entire process above until the sum of moving length of all the particles is less than some preset threshold . At this time, some separating section of the natural emergence in the space can be observed; each section corresponds a separate cluster.

5. Analog and Simulation Experiments

5.1. The Experimental Environment

In order to compare quantum random walk clustering algorithm (QRWA), genetic algorithm (GA), ant colony optimization (ACO), and simulated annealing algorithm (SAA), we select six servers as nodes. In the experiments, we select the number of task from 0 to 2500 (or 3000). We compare the results of these schemes by Matlab. The correlation parameters of selected servers for experiments are in Table 1.

Table 1: Parameters of selected servers.

SM
Model CPU
MS
NA
OS

SUN SPARC Enterprise T5120
UltraSPARC T2

1.2 G
32.0 G
2 * 1000 M
Solaris

IBM System p5520Q
POWER5+

2.66 G
32.0 G
4 * 1000 M
Linux

HP rx4640
intel Itanium 2.

1.5 G 
32.0 G
2 * 1000 M
Unix

In Table 1, OS represents operation system; NA represents network adapter; MS represents memory size; SM represents specifications and models.

The topological structure of the network servers is as shown in Figure 1.

Figure 1: The topological structure of the network servers.

5.2. Results

Figure 2 shows the system flow control rate of QRWA is better than GA, SAA, and ACO. And the more the task quantity is, the closer the flow control rate is. The task quantity is from 0 to 2500.

Figure 2: The flow control rate of SUN SPARC Enterprise T5120 in GA, SAA, ACO, and QRWA.

Figure 3 shows the server traffic of GA, SAA, ACO, and QRWA. That is to say, the QRWA is bigger than GA, SAA, and ACO.

Figure 3: The server traffic of SUN SPARC Enterprise T5120 in GA, SAA, ACO, and QRWA.

The Figure 5 shows that the throughput of QRWA is bigger than ACO, GA, and SAA.

The Figure 6 shows that the throughput of QRWA is bigger than ACO, GA, and SAA.

Figures 4, 5, and 6 show the performances of QRWA are better than GA, SAA, and ACO.

Figure 4: In SUN SPARC Enterprise T5120, the relationship of the flow control rate, task quantity and time. And the QRWA is better than GA, SAA, and ACO.

Figure 5: The throughput of HP rx4640 in ACO, GA, SAA, and QRWA during one day.

Figure 6: The throughput of IBM System p5520Q in ACO, GA, SAA, and QRWA during one week.

From the results, it is clear that quantum random walk algorithm (QRWA) is better in server traffic control and task scheduling than genetic algorithm (GA), simulated annealing algorithm (SAA), and ant colony optimization (ACO). QRWA is more effective in task scheduling.

6. Conclusions

The paper gives a quantum random walk model and algorithm on server traffic control and task scheduling. We mainly research the standard model of one-dimensional quantum random walk. For the data clustering problem of high dimensional space, we can decompose one -dimensional quantum random walk into one-dimensional quantum random walk.

The simulation results demonstrate the effectiveness and superiority of QRWA.

The model and algorithm increases the throughput and efficiency of the system, and it had some merits than traditional model and arithmetics.

We will research the two directions in the future. The first one is the effects of noise on the scheme and model; the second one is the method of how to apply in the field of intelligence.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by the National Natural Science Foundation of China (61173056).

References

S. Koppaka and A. R. Hota, “Superior exploration-exploitation balance with quantum-inspired Hadamard walks,” in Proceedings of the 12th Annual Genetic and Evolutionary Computation Conference (GECCO '10), pp. 2093–2094, Companion Publication, July 2010.

D. Emms, S. Severini, R. C. Wilson, and E. R. Hancock, “Coined quantum walks lift the cospectrality of graphs and trees,” Pattern Recognition, vol. 42, no. 9, pp. 1988–2002, 2009.

B. W. Reichardt and R. Špalek, “Span-Program-based quantum algorithm for evaluating formulas,” in Proceedings of the 40th Annual ACM Symposium on Theory of Computing (STOC '08), pp. 103–112, May 2008.

G. J. Martin, S. C. Gillespie, and C. H. Volk, “Litton 11 cm triaxial zero-lock gyro,” in Proceedings of the IEEE Position Location and Navigation Symposium (PLANS '96), pp. 45–55, April 1996.

G. Kuczera and E. Parent, “Monte Carlo assessment of parameter uncertainty in conceptual catchment models: the Metropolis algorithm,” Journal of Hydrology, vol. 211, no. 1-4, pp. 69–85, 1998.

C. Li, B. Yu, and K. Sycara, “An incentive mechanism for message relaying in unstructured peer-to-peer systems,” Electronic Commerce Research and Applications, vol. 8, no. 6, pp. 315–326, 2009.

A. Romanelli, A. C. S. Schifino, R. Siri, G. Abal, A. Auyuanet, and R. Donangelo, “Quantum random walk on the line as a Markovian process,” Physica A: Statistical Mechanics and its Applications, vol. 338, no. 3-4, pp. 395–405, 2004.

A. Chakrabarti, C. Lin, and N. K. Jha, “Design of quantum circuits for random walk algorithms,” in Proceedings of the IEEE Computer Society Annual Symposium on VLSI (ISVLSI '12), pp. 135–140, 2012.

J. Duda, “From maximal entropy random walk to quantum thermodynamics,” Journal of Physics: Conference Series, vol. 361, no. 1, Article ID 012039, 2012.

K. Rudinger, J. K. Gamble, E. Bach, M. Friesen, R. Joynt, and S. N. Coppersmith, “Comparing algorithms for graph isomorphism using discrete-and continuous-time quantum random walks,” Journal of Computational and Theoretical Nanoscience, vol. 10, no. 7, pp. 1653–1661, 2013.

M. S. Underwood and D. L. Feder, “Universal quantum computation by discontinuous quantum walk,” Physical Review A: Atomic, Molecular, and Optical Physics, vol. 82, no. 4, Article ID 042304, 2010.

L. Jun, Investigations on Quantum Random Walk Search Algorithm, 2006.

N. Shenvi, J. Kempe, and K. B. Whaley, “Quantum random-walk search algorithm,” Physical Review A: Atomic, Molecular, and Optical Physics, vol. 67, no. 5, Article ID 052307, 2003.