^{1}

^{1}

^{2}

This paper proposes a Hybridized Ant Colony Optimization (HACO) algorithm. It integrates the advantages of Ant System (AS) and Ant Colony System (ACS) of solving optimization problems. The main focus and core of the HACO algorithm are based on annexing the strengths of the AS, ACO and the Max-Min Ant System (MMAS) previously proposed by various researchers at one time or the order. In this paper, the HACO algorithm for solving optimization problems employs new Transition Probability relations with a Jump transition probability relation which indicates the point or path at which the desired optimum value has been met. Also, it brings to play a new pheromone updating rule and introduces the pheromone evaporation residue that calculates the amount of pheromone left after updating which serves as a guide to the successive ant traversing the path and diverse local search approaches. Regarding the computational efficiency of the HACO algorithm, we observe that the HACO algorithm can find very good solutions in a short time, as the algorithm has been tested on a number of combinatorial optimization problems and results shown to compare favourably with analytical results. This strength can be combined with other metaheuristic approaches in the future work to solve complex combinatorial optimization problems.

Ant System algorithms have been founded on the observation of real ant colonies by various researchers to which Marco Dorigo a key factor as it was found in his Ph.D. Thesis of 1992, where the concept was first brought to light. By living in colonies, ants’ social behavior is directed more to the survival of the colony than that of a single individual member of the colony. An interesting and significantly important behavior of ant colonies is their foraging behavior, and in particular, their ability to find the shortest route between their nest and a food source, realizing that they are almost blind. Without any leader that guides the ants to optimal trajectories, the ants manage to find these optimal trajectories with time, by interacting with their local environment.

The ants initially search for food in a random fashion such that, when they have found some, they return home depositing chemicals, called pheromones. According to [

It follows that when a number of paths are available from the nest to a food source, a colony of ants may be able to exploit the pheromone trail deposited by the individual members of the colony to discover the shortest path from the nest to the food source and back [

In general, Ant Colony Optimization (ACO) algorithms employ a finite size of artificial agents with defined characteristics which collectively search for good quality solutions to the problem under consideration. Starting from an initial state selected according to some case-dependent criteria, each ant builds a solution, which is similar to a chromosome in a genetic algorithm. While building its own solution, each ant collects information on its own performance and uses this information to modify the representation of the problem, as seen by the other ants [

The amount of pheromone deposited is made proportional to the goodness of the solution an artificial ant has developed (or is developing). Rapid drift of all the ants towards the same part of the search space is avoided by employing the stochastic component of the choice decision policy and the pheromone evaporation mechanism. To simulate pheromone evaporation, the pheromone persistence coefficient, ρ, is defined which enables greater exploration of the search space and minimizes the chance of premature convergence to suboptimal solutions.

A probabilistic decision policy is also used by the ants to direct their search towards the most interesting regions of the search space. The level of stochasticity in the policy and the strength of the updates in the pheromone trail determine the balance between the exploration of new points in the state space and the exploitation of accumulated knowledge [

This paper is structured as follows: Section 2 describes and briefly reviews variants Ant systems that have been developed over years by various researchers. While Section 3 is dedicated to the introduction and description in detail of the Hybridized Ant Colony Optimization (HACO) algorithm. Section 4 presents a step by step computational procedure of the HACO algorithm. Section 5 concludes the paper.

The concept of metaheuristic is a general purposes algorithmic framework that can be applied to different optimization problems with relatively few modifications to make them adaptable to a specific problem. In the view of [

A) How do the ants construct a solution?

B) How do the ants choose the next step to take?

C) How is the solution post-optimized?

D) How is the pheromone deposit updated?

Any variant of ACO algorithm that will be regathered as amenable will find answers to all the afore mentioned questions. There are numerous examples of Combinatorial Optimization problems and the variants being handled using Ant Colony Optimization algorithms. This is a type of algorithm that seeks to model the emergent behaviour observed in ant colonies and utilize this behaviour to solve problems. Before we roll out the HACO algorithm there is need to have an outline to various ACO developed by various researchers. This leads us to variants ACO algorithms that have been developed over years by various researchers and improvements that ACO algorithms have witnessed.

In the last three decades, there have been variants Ant Colony Optimizations due to modification and improvements in one respect or the other. It must be noted that, the several variants Ant Colony Optimization differing most importantly on when and how pheromones are updated. However, some of the variants are highlighted and discussed in what follows.

Ant Colony Optimization (ACO), called Ant System in [

Ant System is the first ACO algorithm proposed in literature. The AS is based on the behavior of real ants searching for food. Real ants communicate with each other using an aromatic essence, called pheromone, that they lay down on the path they traverse. By [

Its main characteristics are as follows: at each iteration, the pheromone values are updated by all the m ants that have built a solution in the iteration itself. At each step, an ant at node i chooses to go to a node j through the edge ( i , j ) and this arc is added to the solution. The repetition of this step stops when the ant has completed its tour. The pheromone, τ i j , associated with the edge joining cities, i and j, is updated using the relation:

τ i j ← ( 1 − ρ ) τ i j + ρ ∑ k = 1 m Δ τ i j k (1)

where ρ is the evaporation rate, m is the number of ants, and Δ τ i j k is the quantity of pheromone laid on edge ( i , j ) by ant k:

Δ τ i j k = { 1 / L k if ant k used ( i , j ) in its tour 0 otherwise (2)

where L k is the length of the tour constructed by ant k.

The ant colony system (ACS) was developed as in [

Once all the ants have built a solution, a pheromone trail is added to the arcs used by the ant that has found the best tour from the beginning of the trial. Thus, instead of allowing all the m ants to update the pheromones as in AS, only the ant that has found the best solution in ACS deposits pheromone on the arcs of the best tour. The local pheromone update is performed by all the ants after each construction step. Each ant applies update only to the last edge traversed as opined by [

τ i j = ( 1 − ρ ) τ i j + ρ τ 0 (3)

where ρ ∈ ( 0 , 1 ] is the pheromone decay coefficient, and τ 0 is the initial value of the pheromone.

The offline pheromone update, similarly to Max-Min Ant System (MMAS), is applied at the end of each iteration by the ant with the shortest tour, which can be either the iteration-best or the best-so-far. The update formula is:

τ i j ← { ( 1 − ρ ) τ i j + ρ Δ τ i j , if ( i , j ) ∈ t τ i j , otherwise (4)

In MMAS, not all pheromone levels are evaporated, as with the AS, but only those that also receive a pheromone deposit.

In MMAS, only one single ant is used to update the pheromone trails after every iteration. The algorithm achieves a strong exploitation of the search history by allowing only the best solutions to add pheromone during the pheromone trail update. Its characterizing elements are that only the best ant updates the pheromone trails and that the value of the pheromone is bounded by | τ min − τ max | . Also, the use of a rather simple mechanism for limiting the strengths of the pheromone trails effectively avoids premature convergence of the search.

The pheromone update is implemented as follows:

τ i j = [ ( 1 − ρ ) ⋅ τ i j + τ i j b e s t ] τ min τ max (5)

where τ max and τ min are respectively the upper and lower bounds imposed on the pheromone; the operator [ x ] b a is defined as:

[ x ] b a = { a if x > a b if x < b x otherwise (6)

and τ i j b e s t is defined by:

τ i j b e s t = { 1 / L b e s t if ( i , j ) belongs to the best tour 0 otherwise (7)

where L b e s t is the length of the tour of the best ant. This may be either the best tour found in the current iteration i.e. iteration-best, L i b or the best solution found since the start of the algorithm-best-so-far, L b s or a combination of both [

The Hybridized Ant Colony Optimization hinges on the strengths of the afore mentioned variants and improves on the drawbacks leading to the hybridized algorithm stressing on the Transition Probability, Jump Transition Probability, generating Randomized Numbers that guides in determining the path traversed by the ants, pheromone updating rule, pheromone evaporation residue and diverse local search approaches are proposed hence, evaluating the function values to which the best and worst function values is identified will be discussed in the next section.

The HACO algorithm is described as follows:

Step 1: Initialize, N , n , ρ , τ 0 . Set ξ = 0 . 0 0 0 0 1 and ζ = 2 .

where A a , a = 1 , 2 , ⋯ , N represent the number of Ants in the colony. { x j } represents a set of permissible discrete values of the Design Variables, x j , and j = 1 , 2 , ⋯ , n . The interval, α ≤ x j ≤ β or [ α , β ] , is divided into n equal parts such that, each part is assigned to each of the design variables, x j . The pheromone evaporation rate is assumed to be ρ = 0.5 . It is assumed that all the ants will deposit equal amount of pheromone at the initial stage hence, τ 0 = 1 .

Step 2: a) Compute the Transition Probability, T P i , j , using:

T P i , j = τ i ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1 (8)

where f b e s t is the highest function value, ∑ f r e q o f f b e s t is the sum of the frequency of the f b e s t where more than one exist and the Jump Transition Probability, P i , j j u m p is computed using:

T P i , j j u m p = τ i + 1 ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1 (9)

where the subscripts, i and j of T P i , j and T P i , j j u m p are the iterations and paths counters respectively. The cumulative transition probabilities’ values are computed and it is being associated with the paths traversed by the ant based on the value of the transition probability obtained using the expressions (8) and (9).

N.B: At i = 0 , the Jump Transition Probability, T P i , j j u m p , is not computed since each path has equiprobable chances of being traversed. Also, the ∑ f r e q o f f b e s t = 0 and τ i + 1 = 0 . The Jump Transition Probability, T P i , j j u m p , is aimed at speeding up the transition probability of the path traversed after a path have been identified as being best owning to the function value.

b) Compute the Cumulative Transition Probability, C T P i , j .

At iteration i = 0 , the Cumulative Transition Probability, C T P i , j along the path, j, is given by:

C T P i , 1 = T P i , 1 C T P i , 2 = C T P i , 1 + T P i , 1 C T P i , 3 = C T P i , 2 + T P i , 1 ⋮ C T P i , n = C T P i , n − 1 + T P i , 1 } (10)

At iteration i > 0 , the Cumulative Transition Probability, C T P i , j along the path, j, is given by:

C T P i , 1 = T P i , 1 C T P i , 2 = C T P i , 1 + T P i , 1 ⋮ C T P i , m = C T P i , m − 1 + T P i , 1 C T P i , m + 1 = C T P i , m + T P i , j j u m p C T P i , m + 2 = C T P i , m + 1 + T P i , 1 C T P i , m + 3 = C T P i , m + 2 + T P i , 1 ⋮ C T P i , n = C T P i , ( n − 1 ) + T P i , 1 } (11)

Step 3: a) Generate N random numbers.

The specific path chosen by a particular ant will then be determined using the Roulette Wheel selection process of randomization. From the closed interval [ 0 , 1 ] , N random numbers:

r 1 , r 2 , ⋯ , r N (12)

are generated such that, an ant, A a , is assigned a random number each. The assigned randomized number is then further used as a guide towards the path traversed by the ant.

b) Correlate each of the randomized numbers, r a .

Correlate each of the randomized numbers, r a , with a transition probability cumulative value, C T P i , j , such that the randomized number, r a , falls within the range of value of the C T P i , j along a particular path thus:

r a ≃ C T P i , j (13)

Step 4: a) Evaluate the objective function value, f ( x j ) .

The objective function value, f ( x j ) is evaluated with the value assigned to the design variable, x j , that corresponds to the cumulative transition probability’s paths, j, of the design variable path. i.e.

f a = C T P i , j ( j ) = f ( x j ) (14)

b) Determine the best and worst paths.

The best and worst paths among the N paths chosen by different ants are determined using the relation:

f b e s t = M a x ( f ( x 1 ) , f ( x 2 ) , ⋯ , f ( x N ) ) (15)

f w o r s t = M i n ( f ( x 1 ) , f ( x 2 ) , ⋯ , f ( x N ) ) (16)

Step 5: Test for convergence of the process.

The process is assumed to have converged if all the ants take the same path. i.e.

f ( x 1 ) = f ( x 2 ) = ⋯ = f ( x N ) or (17a)

f ( x N ) − f ( x N − 1 ) = f ( x N − 1 ) − f ( x N − 2 ) = ⋯ = f ( x 2 ) − f ( x 1 ) = ξ ≤ 0.00001 (17b)

If the convergence criterion has been satisfied, then stop. Else, it will be assumed that, all the ants return to the colony and start the search for a shorter path to the food source again.

Step 6: a) Calculate the Pheromone Residue, M:

M = ( ∑ i = 1 N τ i ) − ρ . (18)

where ρ , rate of pheromone evaporation, is a user defined parameter.

b) Evaluate the objective function value’s ratio

The range ratio of the best objective function value, f b e s t , to the worst objective function value, f w o r s t , in maximizing a function given by Δ τ j k , is evaluated using the relation:

Δ τ j k = ζ ( ∑ f r e q o f f b e s t ) ( f b e s t ) f w o r s t (19)

where the summation extends over all the best ants, A a (if multiple ants take the same best path). The subscripts, j, is the path preceding node and k is the path succeeding node.

c) Compute the pheromone deposit update:

To compute the pheromone deposits update set i = i + 1 and use the following relations:

τ i = ( 1 − ρ ) τ i − 1 (20)

τ i + 1 = τ i + Δ τ j k (21)

where τ i + 1 is the pheromone on the node link, ( j , k ) , after updating and τ i is the pheromone on the link ( j , k ) . This is before the update is computed, then, go to Step 2.

The HACO algorithm is now applied to solve the following test problems:

Problem 1: Find the Maximum of

f ( x ) = − 200 x 2 + 92000 x − 8400000 (22)

within the range 220 - 300 using HACO algorithm.

Step 1: Let A a , a = 1 , 2 , ⋯ , N represent the number of Ants in the colony.

Let { x j } represent a set of permissible discrete values of the Design Variable, x j , where j = 1 , 2 , ⋯ , n . The range, α ≤ x j ≤ β or [ α , β ] , is divided into n equal parts such that, a part each is assigned to each of the design variables, x j . Let ρ = 0.5 represent the pheromone evaporation rate. It is assumed that all the ants will deposit equal amount of pheromone at the initial stage hence, τ 1 = 1 .

Number of Ants = 4, n = 9 , are assumed within the range for x i as:

x 1 = 220 , x 2 = 230 , x 3 = 240 , x 4 = 250 , x 5 = 260 , x 6 = 270 , x 7 = 280 , x 8 = 290 , x 9 = 300

Step 2: Compute the Transition Probability, T P i , j , of selecting the arc or ray (or the discrete values) x j as:

T P i , j = τ i ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1 , j = 1 , 2 , ⋯ , n (23)

T P 1 , 1 = 1 ( 9 − 0 ) × 1 + 0 = 1 9 = 0.11111111111 (24)

The cumulative of the transition probability values are computed associated with different paths based on the transition probabilities values obtained from the expressions (8) and (9). The Cumulative Transition Probability, C T P i , j at iteration i and path j is given by:

C T P i , 1 = T P i , 1 C T P i , 2 = C T P i , 1 + T P i , 1 C T P i , 3 = C T P i , 2 + T P i , 1 ⋮ C T P i , n = C T P i , n − 1 + T P i , 1 }

Hence,

C T P 1 , 2 = 0.11111111111 + 0.11111111111 = 0.222222222222 ,

C T P 1 , 3 = 0.11111111111 + 0.22222222222 = 0.333333333333 ,

C T P 1 , 4 = 0.11111111111 + 0.33333333333 = 0.444444444444 ,

C T P 1 , 5 = 0.11111111111 + 0.44444444444 = 0.555555555556 ,

C T P 1 , 6 = 0.11111111111 + 0.55555555556 = 0.666666666667 ,

C T P 1 , 7 = 0.11111111111 + 0.66666666667 = 0.777777777778 ,

C T P 1 , 8 = 0.11111111111 + 0.77777777778 = 0.888888888889 ,

C T P 1 , 9 = 0.11111111111 + 0.88888888889 = 1.000000000000.

N. B: Since none of the paths have been identified as f b e s t , the jump transition probability needs not be computed at this stage.

Step 3: a) Generate N random numbers, r 1 , r 2 , ⋯ , r N , within the [ 0 , 1 ] . One random number for each ant thus:

r 1 = 0.3122 , r 2 = 0.8701 , r 3 = 0.4729 , r 4 = 0.6190 (25)

b) Correlate each of the randomized numbers, r a . with a transition probability cumulative value, C T P i , j , which is of the most approximate value or the cumulative transition probability value’s range to which the randomized numbers value, r a , fall:

r a ≅ C T P i , j . (26)

Ant 1: r 1 = 0.3122 ≈ C T P 1 , 3 = 0.333333333333

Ant 2: r 2 = 0.8701 ≈ C T P 1 , 8 = 0.888888888889

Ant 3: r 3 = 0.4729 ≈ C T P 1 , 5 = 0.555555555556

Ant 4: r 4 = 0.6190 ≈ C T P 1 , 6 = 0.666666666667

Step 4: a) Evaluate the objective function value, f ( x j ) .

The objective function values, f ( x j ) are evaluated using values assigned to the design variable values, x j , corresponding to the cumulative of the transition probability’s paths, j.

f a = C T P i , j at j = f ( x j ) (27)

Ant 1: f 1 = C T P 1 , 3 = 0.333333333333.

Then, x 3 = 240 . Hence, f ( x 1 ) = f ( 240 ) = 2160000 .

Ant 2: f 2 = C T P 1 , 8 = 0.888888888889.

Then, x 8 = 290 . Hence, f ( x 2 ) = f ( 290 ) = 1460000 .

Ant 3: f 3 = C T P 1 , 5 = 0.555555555556.

Then, x 5 = 260 . Hence, f ( x 3 ) = f ( 260 ) = 2000000 .

Ant 4: f 4 = C T P 1 , 6 = 0.666666666667.

Then, x 6 = 270 . Hence, f ( x 4 ) = f ( 270 ) = 1860000 .

b) Determine the best and worst paths.

The best and worst paths among the N paths chosen by different ants are determined using the relation:

f b e s t = M a x ( f ( x 1 ) , f ( x 2 ) , ⋯ , f ( x N ) ) (28)

f w o r s t = M i n ( f ( x 1 ) , f ( x 2 ) , ⋯ , f ( x N ) ) (29)

Hence, f b e s t = f ( x 3 ) = 216000 at x b e s t = x 3 = 240 and

f w o r s t = f ( x 8 ) = 146000 at x w o r s t = x 8 = 290 .

Step 5: Test for convergence of the process. The process is assumed to have converged if all the N ants take the same path. i.e.

f ( x 1 ) = f ( x 2 ) = ⋯ = f ( x N ) or (30)

f ( x N ) − f ( x N − 1 ) = f ( x N − 1 ) − f ( x N − 2 ) = ⋯ = f ( x 2 ) − f ( x 1 ) = ξ ≤ 0.00001 (31)

However, if convergence has not been achieved, it will be assumed that all the ants will have to return home and start again in search of food. Since neither of the convergence conditions is met, the process has to continue.

Step 6: a) Calculate the Pheromone Evaporation constant, M.

Since the convergence criterion has not been met, there is need to calculate the pheromone evaporation constant, M, using the relation:

M = ( ∑ i = 1 N τ i ) − ρ . (32)

M = 4 × 1 − 0.5 = 3.5

b) Evaluate the objective function value’s ratio:

Δ τ j k = ζ ( ∑ f r e q o f f b e s t ) ( f b e s t ) f w o r s t (33)

Δ τ 2 , 4 = 1 × 2 × 2160000 1460000 = 2.95890410959

where ζ = 2 , ρ = 0.5 .

c) Compute the Pheromone deposit update: Then, set the new iteration number as i = i + 1 , and update the pheromone deposits on different arcs (or discrete values of design variables) where τ i − 1 denotes the previous pheromone amount of the iteration after evaporation, which is taken as τ i can be computed by:

At i = i + 1 = 2 ,

τ i = ( 1 − ρ ) τ i − 1

τ 2 = ( 1 − 0.5 ) τ 1 = 0.5 × 1.0 = 0.5

τ i + 1 = τ i + Δ τ j k , (34)

where j = 1 , 2 , 4 , 5 , 6 , 7 , 8 , 9 . It must be noted that the path, j = 3 is omitted. This is because, that was the path at which f b e s t has been identified.

τ 3 = τ 2 + Δ τ 2 , 4 = 1.0 + 2.95890410959 = 3.95890410959

With this, we go back to step 2.

Step 2: Compute the Transition Probability, ( T P i , j ) and Transition Probability Jump, ( T P i , j j u m p ) of selecting the arc:

T P i , j = τ i ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1 and

T P i , j j u m p = τ i + 1 ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1

where f r e q o f f b e s t is the frequency of the f b e s t , τ 2 = 0.5 , τ 3 = 3.95890410959 and j = 1 , 2 , 4 , 5 , 6 , 7 , 8 , 9

This gives:

T P 2 , 1 = τ 2 ( N − ∑ f r e q o f f b e s t ) ( τ 2 ) + τ 3 = 0.5 ( 9 − 1 ) × 0.5 + 3.95890410959 = 0.06282271945

T P 2 , 3 j u m p = τ 3 ( N − ∑ f r e q o f f b e s t ) ( τ 2 ) + τ 3 = 3.95890410959 7.95890410959 = 0.49741824441

The Transition Probabilities would have been equiprobable but for the fact that, the transition probability has to jump at P 2 , 3 . Hence,

C T P 2 , 1 = 0.06282271945 ,

C T P 2 , 2 = 0.06282271945 + 0.06282271945 = 0.1256454389 ,

Note that, there is a jump at this point. Hence,

C T P 2 , 2 j u m p = 0.49741824441 + 0.12564543890 = 0.62306368330 ,

C T P 2 , 4 = 0.06282271945 + 0.62306368330 = 0.68588640276 ,

C T P 2 , 5 = 0.06282271945 + 0.68588640276 = 0.74870912221 ,

C T P 2 , 6 = 0.06282271945 + 0.74870912221 = 0.81153184166 ,

C T P 2 , 7 = 0.06282271945 + 0.81153184166 = 0.87435456111 ,

C T P 2 , 8 = 0.06282271945 + 0.87435456111 = 0.93717728056 ,

C T P 2 , 9 = 0.06282271945 + 0.93717728056 = 1.00000000001

Step 3: Generate N random numbers r 1 , r 2 , ⋯ , r N in the [ 0 , 1 ] , one for each ant

r 1 = 0.3688 , r 2 = 0.8577 , r 3 = 0.0776 , r 4 = 0.5791

a) Evaluate the objective function values corresponding to the complete paths

Ant 1 : f 1 = f ( x 1 ) = f ( 240 ) = 2160000

Ant 2 : f 2 = f ( x 2 ) = f ( 280 ) = 1680000

Ant 3 : f 3 = f ( x 3 ) = f ( 230 ) = 2180000

Ant 4 : f 4 = f ( x 4 ) = f ( 240 ) = 2160000

Step 4: Determine the best and worst paths among the N paths chosen by different ants

x b e s t = x 3 = 230 ; f b e s t = f 3 = 218000

x w o r s t = x 7 = 280 ; f w o r s t = f 2 = 168000

Step 5: At i = 3 , τ i + 1 = τ i + Δ τ j k ,

Δ τ j k = ζ ( ∑ f r e q o f f b e s t ) ( f b e s t ) f w o r s t = 2 × 2180000 1680000 = 2.59523809524 ,

where ζ = 2 , ρ = 0.5 and τ i can be computed as:

τ i = ( 1 − ρ ) τ i − 1 = ( 1 − 0.5 ) τ i − 1 = 0.5 × 0.5 = 0.25 ; j = 1 , 3 , 4 , 5 , 6 , 7 , 8 , 9

with the path 2 missed out. Hence,

τ 4 = τ 3 + Δ τ 1 , 3 = 3.95890410959 + 2.59523809524 = 6.55414220483

for j = 2 and τ 3 = 0.25

With this, we set i = 4 and go back to step 2.

Step 2: Compute the Transition Probability ( T P i j ) of selecting the arc

T P i , j = τ i ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1 and

T P i , j j u m p = τ i + 1 ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1 ,

where τ 3 = 0.25 ; j = 1 , 3 , 4 , 5 , 6 , 7 , 8 , 9 and τ 4 = 6.55414220483

This gives T P 3 , 1 = 0.25 8 × 0.25 + 6.55414220483 = 0.25 8.55414220483

T P 3 , 1 = 0.0292256072

T P 3 , 2 j u m p = 6.55414220483 8.55414220483 = 0.76619514241

C T P 3 , 3 = 0.82464635681 , C T P 3 , 4 = 0.85387196401 ,

C T P 3 , 5 = 0.88309757121 , C T P 3 , 6 = 0.91232317841 ,

C T P 3 , 7 = 0.94154878561 , C T P 3 , 8 = 0.97077439281 , C T P 3 , 9 = 1.00000000001

Step 3: a) Generate N random numbers r 1 , r 2 , ⋯ , r N in the [ 0 , 1 ] , one for each ant:

r 1 = 0.9251 , r 2 = 0.4928 , r 3 = 0.0632 , r 4 = 0.2064

b) Evaluate the objective function values corresponding to the complete paths:

Ant 1 : f 1 = f ( x 1 ) = f ( 280 ) = 1680000

Ant 2 : f 2 = f ( x 2 ) = f ( 230 ) = 2180000

Ant 3 : f 3 = f ( x 3 ) = f ( 230 ) = 2180000

Ant 4 : f 4 = f ( x 4 ) = f ( 230 ) = 2180000

Step 4: Determine the best and worst paths among the N paths chosen by different ants

x b e s t = x 1 , 2 = 230 ; f b e s t = f 2 = f 3 = f 4 = 2180000

x w o r s t = x 1 , 7 = 280 ; f w o r s t = f 1 = 1680000

Step 5: i = 4 ,

Δ τ 1 , 5 = ζ ( ∑ f r e q o f f b e s t ) ( f b e s t ) f w o r s t = 3 × 2 × 2180000 1680000 = 7.78571428571

where ζ = 2 , ρ = 0.5 , τ i and τ i + 1 can be computed as

τ i = ( 1 − ρ ) τ i − 1 = ( 1 − 0.5 ) τ 3 = 0.5 × 0.25 = 0.125 ; j = 1 , 5 , 6 , 7 , 8 , 9

τ i + 1 = τ i + Δ τ j k

τ 5 = τ 4 + Δ τ 1 , 5 = 6.55414220483 + 7.7857142857 = 4.33985649054

for j = 2 and τ 4 = 0.125

With this we back go to step 2.

Step 2: Compute the Cumulative transition probability ( C T P i j ) of selecting the arc:

T P i , j = τ i ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1 and

T P i , j j u m p = τ i + 1 ( N − ∑ f r e q o f f b e s t ) ( τ i ) + τ i + 1 ,

where τ i j = 0.125 ; j = 1 , 3 , 4 , 5 , 6 , 7 , 8 , 9 and τ 4 , 2 = 14.33985649054

This gives:

T P 4 , 1 = 0.125 6 × 0.125 + 14.33985649054 = 0.125 15.08985649054 = 0.00828371032

T P 4 , 2 j u m p = 14.33985649054 15.08985649054 = 0.95029773806

C T P 4 , 1 = 0.00828371032 ,

C T P 4 , 2 j u m p = 0.95858144838 , C T P 4 , 3 = 0.9668651587 ,

C T P 4 , 4 = 0.97514886902 , C T P 4 , 5 = 0.98343257934 ,

C T P 4 , 6 = 0.99171628966 , C T P 4 , 7 = 0.99999999998 ,

C T P 4 , 8 = 1.0082837103 , C T P 4 , 9 = 1.01656742062

Step 3: a) Generate N random numbers r 1 , r 2 , ⋯ , r N in the [ 0 , 1 ] , one for each ant:

r 1 = 0.5684 , r 2 = 0.8385 ,

r 3 = 0.4191 , r 4 = 0.250

b) Evaluate the objective function values corresponding to the complete paths:

Ant 1 : f 1 = f ( x 1 ) = f ( 230 ) = 2180000

Ant 2 : f 2 = f ( x 2 ) = f ( 230 ) = 2180000

Ant 3 : f 3 = f ( x 3 ) = f ( 230 ) = 2180000

Ant 4 : f 4 = f ( x 4 ) = f ( 230 ) = 2180000

It can be seen that all the ants now follow the same path, that is, Ant 1, Ant 2, Ant 3 and Ant 4 follow the same path. Hence, this is acclaimed as the best path with:

x b e s t = x 2 = 230 , f b e s t = f 1 = f 2 = f 3 = f 4 = 2180000

τ 4 , 2 = 14.33985649054 and f ( x * ) = 2180000 .

In this paper, we tried to annex the strengths of AS, ACS and MMAS, Ant Colony Optimization algorithms, to obtain solution of optimization problems. To this means, due to high computational complexity in solving these problems using the earliest ACO algorithms, we modify the relation for determining the Transition Probability, formulate the Jump Transition Probability and Pheromone Evaporation Residue relations and introduce the process of correlating the randomized numbers, r a , with the cumulative transition probability cumulative value, C T P i , j , into the algorithm to obtain the best solution to build a Hybridized ACO such that the randomized number, r a , falls within the range of value of the C T P i , j along a particular path. Implementation of the method and numerical results that show the accuracy of the method, speedy convergence and low occupation of the computer memory are also discussed.

The authors declare no conflicts of interest regarding the publication of this paper.

Adebayo, K.J., Aderibigbe, F.M. and Dele-Rotimi, A.O. (2019) On the Development of a Hybridized Ant Colony Optimization (HACO) Algorithm. American Journal of Computational Mathematics, 9, 358-372. https://doi.org/10.4236/ajcm.2019.94026