Next Article in Journal
Enhanced Intrusion Detection Systems Performance with UNSW-NB15 Data Analysis
Next Article in Special Issue
An Attention-Based Method for the Minimum Vertex Cover Problem on Complex Networks
Previous Article in Journal
Program Code Generation with Generative AIs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Group Dynamics in Memory-Enhanced Ant Colonies: The Influence of Colony Division on a Maze Navigation Problem

Department of Mathematics and Computer Science, University of Catania, v.le Andrea Doria 6, 95125 Catania, Italy
*
Authors to whom correspondence should be addressed.
Algorithms 2024, 17(2), 63; https://doi.org/10.3390/a17020063
Submission received: 20 December 2023 / Revised: 24 January 2024 / Accepted: 27 January 2024 / Published: 1 February 2024
(This article belongs to the Special Issue Algorithms for Network Analysis: Theory and Practice)

Abstract

:
This paper introduces an agent-based model grounded in the ACO algorithm to investigate the impact of partitioning ant colonies on algorithmic performance. The exploration focuses on understanding the roles of group size and number within a multi-objective optimization context. The model consists of a colony of memory-enhanced ants (ME-ANTS) which, starting from a given position, must collaboratively discover the optimal path to the exit point within a grid network. The colony can be divided into groups of different sizes and its objectives are maximizing the number of ants that exit the grid while minimizing path costs. Three distinct analyses were conducted: an overall analysis assessing colony performance across different-sized groups, a group analysis examining the performance of each partitioned group, and a pheromone distribution analysis discerning correlations between temporal pheromone distribution and ant navigation. From the results, a dynamic correlation emerged between the degree of colony partitioning and solution quality within the ACO algorithm framework.

1. Introduction

Network analysis is typically understood as a collection of methodologies and heuristics that allow us to understand the relationships among entities within complex systems. Network analysis finds application in human fields such as sociology [1,2], economics [3,4,5], biology [6,7], medicine [8,9], and many others. For example, in sociology, network analysis can identify online communities [10,11] and social interactions [12] or measure social influence within a friendship network. In economics, we can analyze business relationships and understand how the interconnection between companies affects market stability. In biology, we can study interactions between proteins within cells to better understand biological mechanisms [13,14].
In this paper, we delve into a specific aspect of network analysis within the framework of an Ant Colony Optimization (ACO) algorithm, focusing on the dynamic relationships and interactions among individual ants and various groups within a colony. The colony functions as a complex network, where ants collaborate to find the optimal path between a food source and their nest, communicating through a simple yet indirect mechanism using chemical signals known as pheromones. This phenomenon serves as a prime example of a complex adaptive system (CAS), where the collective behavior transcends the actions of individual elements, showing the colony’s ability to thrive through cooperation [15]. In essence, it is not the solitary ant that discovers the optimal path, but rather the intricate interactions among groups of ants within the colony, facilitated by pheromones, that yield the overall outcome [16]. Our exploration delves into the interactions through pheromones and the structural organization of the colony. By analyzing how various forms of grouping, coupled with factors such as pheromone evaporation, influence intra-colony relations, we aim to comprehensively grasp the diverse outcomes generated by the algorithm in question. Consequently, our investigation extends beyond the realm of pheromone-driven interactions to encompass the substantial impact of colony subdivision on the overall dynamics.

1.1. The ACO Paradigm and Applications

Ant Colony Optimization (ACO) derives its inspiration from the foraging behavior observed in natural ant colonies, where decentralized communication through chemical signals known as pheromones plays a pivotal role in finding the shortest path between their nest and a source of food. Initially proposed in [17], ACO has demonstrated remarkable efficacy in addressing a multitude of complex combinatorial optimization problems, including scheduling and, particularly, routing problems [18,19,20,21,22,23]. More generally speaking, agent-based or swarm intelligence heuristics have been applied as well to coloring challenges [24,25,26], robot path planning for inaccessible areas [27,28], transportation complexities and traffic problems [29,30], and feature selection across diverse fields [31]. Beyond optimization, ACO demonstrates utility in modeling crowd dynamics. In [32], the authors translated ACO rules into conceptual frameworks rooted in social dynamics. They proposed an agent-based model to study collective crowd behaviors. Furthermore, ref. [33] provides a sensitivity analysis.
The notion of dividing ant colonies into groups has received particular attention because of its potential to enhance optimization performance. In [34], an adaptive ACO was implemented, where the colony was separated into three groups: ordinary ants, using heuristic information; abnormal ants, seeking low pheromone paths; and random ants, exploring randomly. This collaborative approach improved the optimization performance, and was able, for instance, to effectively address the traveling salesman problem.
Extending this exploration into multi-objective optimization problems, ref. [35] emphasized the performance benefits of employing multiple ant groups. In their study, a grouping approach was applied to enhance the efficacy of the ACO paradigm when dealing with multi-objective optimization problems.
In a related effort, ref. [36] proposed a multi-objective ACO algorithm for community detection in complex networks. By partitioning the ant colony into groups and coordinating their exploration, the algorithm achieved a more efficient search ability for solutions. Further corroborating these findings, Ref. [37] introduced a group-based ACO algorithm for bi-objective optimization, where ants were organized into groups based on parameters related to the objective function, resulting in an enhanced overall algorithm performance. Recent advancements include the Grouping Ant Colony Optimization (GACO) algorithm presented by [38]. GACO specifically addresses the multi-UAV reconnaissance task allocation problem by dividing the ant colony into distinct subgroups. This division proves crucial as different UAVs focus on distinct targets, ultimately leading to superior optimality in the results when compared to existing methods reported in the literature. Lastly, we stress the possibility of future applications of the GACO algorithm to networks modeling the mobility of groups of people, providing, for instance, personalized recommendations of places of interest [39] and possibly analyzing the impact of restrictions to mobility due to extreme circumstances [40].
Building upon the advancements in multi-objective optimization techniques, some recent studies have introduced innovative algorithms to address specific problems in different domains. For instance, in [30] the authors propose a modified Competitive Mechanism Multi-Objective Particle Swarm Optimization (MCMOPSO) algorithm for optimizing signalized traffic problems. This approach incorporates a dynamic leader selection mechanism, guiding particles toward elite solutions through angle comparisons. Similarly, in [41] a single deep reinforcement learning model is employed to tackle the multi-objective traveling salesman problem (MOTSP). The method makes use of an encoder–decoder framework with a routing encoder, efficiently extracting and aggregating information to generate approximate Pareto-optimal solutions in parallel.

1.2. Objective of Our Research

The exploration of dividing ant colonies into groups as a strategy to optimize Ant Colony Optimization (ACO) raises some compelling questions: what is the optimal degree of partitioning to achieve the best solution to the problem at hand? Are there any potential correlations between the solution quality and variables within the ACO framework? To address these questions, an agent-based model, rooted in the ACO paradigm but where the agents (ants) have memory, was employed to tackle a maze navigation problem. Drawing on insights from existing research, the present study investigates the impact of partitioning the colony into groups on the ACO algorithm’s performance, particularly examining the roles of group size and number in the context of a multi-objective optimization problem. The decision to adopt an agent-based approach is motivated by the necessity to explore the network dynamics inherent in the collaborative behaviors of the ant colony. In related studies [42], the same approach was used to explore the network dynamics of both collaborative and competitive ant colonies. The flexibility of this model proves crucial in capturing the complex and dynamic interactions within the colony, where each agent, representing an individual ant, interacts with the environment and each other. These interactions, facilitated by pheromones, lead to the emergence of collective behaviors, offering a deeper understanding of how the network evolves and adapts to changing conditions while exploring diverse scenarios.
In conclusion, the present study focuses on a dynamic scenario where a colony of memory-enhanced ants, starting from a designated position, collaboratively attempts to discover the optimal path to the exit within a grid network. The colony is divided into groups of varying sizes, with the overarching goal of maximizing the number of ants successfully exiting the maze while minimizing path costs.
The paper is structured as follows: Section 2 clarifies the model, offering a comprehensive understanding of it. In Section 3, the experimental setup and results of the experiments are provided. Section 4 critically discusses the obtained results. Finally, in Section 5 we present the concluding remarks.

2. The Model

We used the NetLogo software [43] to implement our model. NetLogo is an agent-based programming language as well as an integrated development environment (IDE) that provides a comprehensive environment, from the implementation of the simulation model to the creation of the graphical interface supporting the simulation. The graphical interface typically includes various elements such as switches, sliders, choosers, and inputs, chosen according to the user’s needs for model development. These elements allow a real-time interaction to observe changes in the dynamics without directly modifying the code. Interface elements can correspond to either simple model variables or code blocks and are arbitrary, as they may not be present for the model’s executability, as long as the associated variables or code are properly defined in the source code. The graphical interface must include two elements: one designed to invoke the code initializing the simulation environment (typically named “setup”), and another that calls the code initiating the dynamics (commonly labeled as “start”).
The environment is a weighted network defined as a graph G = ( V , E , w ) , where V is the set of vertices, E is the set of edges, and w : V × V R + is the weight function which assigns a positive cost to each edge of the graph. The weight measures the difficulty of crossing a specific edge. The starting point is a randomly selected node on one side (e.g., left side), while the exit point is another randomly chosen node on the opposite side (e.g., right side). The graph is undirected, therefore, all the edges can be crossed in both directions.
Let A i = { j V : ( i , j ) E } represent the set of vertices adjacent to vertex i and let π k ( t ) = ( π 1 , π 2 , , π t ) denote the set of vertices visited by the ant k at a specific time t, where ( π i , π i + 1 ) E for i = 1 , , t 1 . The probability p i j k ( t ) which measures the likelihood that ant k , placed on a vertex i , chooses as destination one of its neighbor vertices j at time t is defined as
p i j k ( t ) = τ i j ( t ) α · η i j β l J i k τ i l ( t ) α · η i l β if j J i k 0 otherwise ,
where J i k = A i { π t k } represents all the possible displacements of ant k from vertex i, τ i j ( t ) is the pheromone intensity on edge ( i , j ) , and η i j ( t ) is the desirability of edge ( i , j ) at a given time t. The parameters α and β determine, respectively, the importance of pheromone intensity and the importance of the desirability of an edge. Each link undergoes asynchronous updates through two types of rules based on the number of ticks T (the tick is the unit of time in the simulation software). The first rule is a local updating rule, which adjusts the pheromone levels after each movement of the ant as follows:
τ i j ( t + 1 ) = τ i j ( t ) + K .
Here, K denotes the pheromone left by each ant after traversing an edge ( i , j ) , and τ i j ( t ) represents the pheromone amount on the link at time t.
The second rule is a global updating rule, which updates the pheromone amount on all network links every T ticks as follows:
τ i j ( t + 1 ) = ( 1 ρ ) · τ i j ( t ) ,
where τ i j ( t ) is the trace amount on edge ( i , j ) at time t, and ρ is the evaporation decay parameter. At any given time t, the desirability η i j ( t ) measures the attractiveness of an edge ( i , j ) . This desirability is associated with the information I i j ( t ) about the weight of the edge. It is shared by each ant on the endpoint of the crossed edge ( i , j ) and initialized to 0, as in Algorithm 1 line 3.

2.1. Memory-Enhanced Ants and Knowledge Acquisition

In our model, ants employ a navigation strategy that combines the use of pheromones and the calculated desirability of edges to make informed decisions about their movement within the environment. Pheromones act as a communication device among ants, directing them towards paths with higher concentrations of pheromones. Conversely, desirability is internally determined based on acquired information I i j ( t ) about edge weights, achieved through two mechanisms.
The first memory mechanism employed by ants involves traversing edges and storing the weights in their memory as prior knowledge  w ¯ p k . As ants explore their environment, they accumulate information about the weights of the edges they traverse. The prior knowledge represents the average weight of the visited edges, computed by summing the weights along a generic path ( π ( t ) ) and dividing by the number of edges in that path (n). This process, outlined in Equation (4), enables ants to establish a historical understanding of the environment. The prior knowledge mechanism allows ants to form preferences for certain paths based on their past experiences. In essence, this memory mechanism reflects the ants’ ability to internalize and use information about previously encountered paths, shaping their decision-making in navigating the environment.
w ¯ p k = 1 n i = 1 n w ( π i , π i + 1 ) ;
The second memory mechanism involves the acquisition of information about edge weights at the nearest endpoint, leading to the formation of their local knowledge  w ¯ l k . Ants discern information ( I i j ( t ) ) about the weights of edges in their immediate vicinity, considering the number of neighboring edges (m) at the position where this information is present (Figure 1).
The local knowledge is then determined as the average weight of these neighboring edges, as expressed in Equation (5). This mechanism allows ants to adapt their decision-making based on real-time information acquired from their immediate surroundings. By considering the specific context at the nearest endpoint, ants enhance their awareness of the current environment and adjust their desirability calculations accordingly. The combination of prior knowledge and local knowledge mechanisms enables ants to strike a balance between past experience and current information, facilitating efficient navigation and adaptation to their dynamic surroundings,
w ¯ l k = 1 m i = 1 m w ( π i , π i + 1 ) .
The overall desirability value η i j ( t ) is then determined using a decision tree (Equation (6)), considering the availability of information I i j ( t ) and the type of knowledge:
η i j ( t ) = 1 w i j if I i j 0 and T 0 1 w ¯ if I i j = 0 and T 0 1 if I i j = 0 and T = 0
The global knowledge ( w ¯ ) is the mean of the prior and local knowledge:
w ¯ = w ¯ p k + w ¯ l k 2 .
When an ant is on a node i, it calculates η i j ( t ) considering the edge weight, with higher desirability for lower weights. If information I i j ( t ) is available at either nearest endpoint, the ant inversely relates desirability to the edge weight. In scenarios where information I i j ( t ) about neighboring edges is absent, the ant estimates desirability as the inverse of its global knowledge. This estimation is derived by averaging prior and local knowledge. Figure 2b illustrates this.
In cases of no prior and local knowledge, as depicted in Figure 2a, the ant assigns an equal desirability value of 1 to all edges.
In essence, as ants traverse edges, their prior and local knowledge contributes to the calculation of the desirability values. These values, representing the attractiveness of different paths, guide ants in choosing routes that align with both the pheromone trail and the desirability criteria. The combination of pheromone-based communication and desirability-driven decision-making allows the ants to navigate efficiently in a complex environment, adapting to changing conditions and optimizing their collective exploration.
The total path cost from the initial to the destination point is determined by the following expression:
i = 1 t 1 w ( π i , π i + 1 ) ,
where π 1 and π t represent the starting and destination points, respectively.

2.2. Description of the ME-ANTS Algorithm

The ME-ANTS algorithm, whose pseudocode is given in Algorithm 1, takes as input a weighted graph G ( V , E , w ) , which represents the environment, with set of vertices V of size | V | = m , set of edges E, and edge weights w. It also takes as inputs several other parameters, namely, α and β , which determine the importance of pheromone intensity and edge desirability; K, which is the amount of pheromone released by each ant; ρ is the evaporation rate; T d and T m a x represent the exit interval between groups and maximum simulation time; A = { a 1 , a 2 , , a n } is the set of ants, partitioned into g groups S 1 , , S g of equal cardinality, with n being the total number of ants, g the number of groups, and Δ the collection of the g groups. The algorithm outputs two key values: num_exited, which denotes the number of ants successfully reaching the exit node, and path_cost, which is the average cost of the traversed paths.
Two matrices are used: P for pheromone levels and I as a Boolean matrix for information release. P is initialized with all entry values equal to 1, while I is initialized with all entry values equal to 0. I plays a crucial role in calculating edge ( v i , v j ) desirability ( η i j ) as ants traverse the edges and release information. Next, the following variables are initialized to manage performance metrics: num_exited counts successful exits, path_costs collects path costs, and each ant a k has a list ( L w k ) associated with it, tracking its traversed edges.
The while loop governs the progression of the simulation, which continues until the current simulation time t until it reaches the maximum T m a x . Within this loop, an if statement checks whether the current time t is a multiple of the total number of vertices | V | . If true, a process is launched involving a designated group S h in Δ and the global updating rule (as defined by Equation (2)). In the network, each ant a k in the selected group S h moves through the network, allowing for the simultaneous exploration of different nodes or edges. The Ant-Movement-Async function runs in parallel for each ant in the network and defines how each ant explores the network.
As the ant traverses the graph, starting from initial node v 1 , it dynamically selects adjacent nodes using the transition probability equation (Equation (1)), considering both the pheromone matrix P and the Boolean matrix I. Once an ant moves between nodes, the visited edge is recorded in the path variable π . Concurrently, pheromone levels on traversed edges are updated, and information is released. A proportional wait time on each node, determined by the edge weight w i j , is incorporated to simulate realistic movement dynamics. If an ant successfully reaches the exit node v m , the count of exited ants (num_exited) is increased, and the cost of the traversed path π is calculated using the defined path cost function (Equation (8)).
The return statement concludes the algorithm, and it provides the count of successful exits (num_exited) and the mean of path costs (mean(path_costs)).
Algorithm 1: ME-ANTS Pseudocode
Algorithms 17 00063 i001

3. Experiments and Results

In the analyses, a graph with | V | = 225 and | E | = 501 was employed. In this graph, each node was connected to a maximum of eight neighbors, and edge weights were real numbers uniformly chosen from the range ( 0 , 1 ] . The simulations involved distributing n = 200 ants across g groups (see Algorithm 1 line 7), with g taking values from the set { 1 , 2 , 5 , 10 , 20 , 50 } and n g representing the number of ants in the group. The specific set considered was a subset of all possible divisors of n, that was selected based on practical considerations and to cover a range of scenarios. Smaller group sizes ( g = 10 , 20 , 50 ) provided insight into the performance of small groups, while larger group sizes ( g = 1 , 2 , 5 ) simulated scenarios where collaboration involved more ants. This approach struck a balance between granularity analysis and computational efficiency, ensuring a thorough exploration of various group sizes while maintaining practical feasibility in terms of computational resources. Exploration by groups started after a fixed time T e = | V | following the previous group, and an overall time limit T m a x was established for the entire colony to reach the exit. This limit was defined as T m a x = c × | V | , with c set to 150. Individual ants had specific time windows to reach the exit, determined by T i = T m a x ( T e × ( i 1 ) ) . Pheromone reduction occurred at a fixed degradation interval of T d = | V | , and the global updating rule took place every T d ticks with evaporation rates of ρ = 0.01 and ρ = 0.001 . The initial pheromone value was τ i j ( t = 0 ) = 1.0 , as expressed in Algorithm 1 line 2, and the parameter determining the amount of pheromone released by ants after crossing an edge was K = 0.1 (as in Algorithm 1 line 32). To assess the impact of group size, 10 independent simulations were conducted for each value of g. Three different kinds of analysis were performed: an overall analysis to evaluate the performance of the colony, a group analysis for every group configuration to evaluate the performance of the single groups, and a pheromone distribution analysis to identify correlations between the temporal diffusion of pheromones and the navigation behavior of ants. All variables and parameters of interest in our model are listed in Table 1 for clarity.

3.1. Overall Analysis

In the overall analysis phase, the costs were normalized using the success rate (SR). The SR is calculated as the ratio between the number of ants that successfully exited and the total number of ants involved in the exploration. This normalization provides a fair and informative perspective, facilitating an accurate assessment of the performance relative to the overall success of ants in completing the exploration of the virtual environment. The tables below present the obtained results based on two evaluation metrics: the success rate at the top of Table 2 and the cost of the discovered paths at the bottom. These values represent the averages across 10 simulations. The values in bold represent the best results for each evaluation metric.
Examining the success rate in Table 2, it is observed that at ρ = 0.01 optimization occurs with 10 groups, while at ρ = 0.001 optimization shifts to 20 groups. Regarding the path cost at the bottom of Table 2, it is noted that, in general, it decreases with an increase in the number of groups. However, the optimal value is at 10 groups for ρ = 0.01 and 20 groups for ρ = 0.001 .

3.2. Group Analysis

In the group analysis phase, adjustments were made to the path costs considering each group’s success rate. The success rate, denoted as S R g , is determined by the ratio of ants that successfully reached the exit to the total number of ants within the group.
When the colony is composed of a single group, as in Table 3, the success rates at both evaporation rates are comparable ( 0.849 at ρ = 0.01 and 0.848 at ρ = 0.001 ), indicating consistent performance. The normalized cost considering the success rate allows for a more accurate evaluation. For ρ = 0.01 the cost is 1293.53 , while at ρ = 0.001 the corresponding value is 1311.197 . These results suggest that despite minor variations the colony maintains relatively consistent performance in terms of success rates and cost across different evaporation rates.
Let us analyze now the results when we have more than one group. In what follows, the index of the group will denote its turn in starting the exploration. Thus, group 1 starts first, then group 2, and so on.
When the colony is divided into two groups, as in Table 4, for ρ = 0.01 group 1, which is the one that starts first, has the highest success rate of 0.856 and a lower cost of 1270.571 . In comparison, the second group, group 2, exhibits a slightly lower success rate of 0.854 and a marginally higher cost of 1281.639 . Shifting to ρ = 0.001 , group 1 maintains a consistently high success rate of 0.860 . However, this comes with a higher cost of 1312.964 . On the other hand, group 2 shows a notable improvement, achieving the highest success rate in this context at 0.876 , accompanied by the lowest cost of 1220.156 .
When the colony is divided into five distinct groups ( g = 5 ), as in Table 5, under ρ = 0.01 group 1 shows a remarkable success rate of 0.9125 accompanied by a cost of 1201.690 . The subsequent groups (3, 4, and 5) show signs of performance degradation, as evidenced by varying success rates and costs. Transitioning to ρ = 0.001 , group 1 maintains the best performance across both metrics, with a success rate of 0.9100 and a cost of 1210.888 . However, similarly to the ρ = 0.01 scenario, groups 2, 3, 4, and 5 exhibit varying degrees of success rates and costs, denoting a slight general degradation in performance as we move across the groups.
When the colony is divided into ten groups ( g = 10 ), as in Table 6, under ρ = 0.01 group 1 distinguishes itself with an impressive success rate of 0.925 , showing a high level of efficiency in reaching the exit. Following closely, group 2 achieves a slightly lower success rate of 0.850 but maintains competitive values for the cost. Group 3 stands out with the lowest cost of 1178.144 , emphasizing its cost-effectiveness despite a moderate success rate. Transitioning to ρ = 0.001 , group 1 continues to lead with the highest success rate of 0.905 . Group 2 shows improvement, achieving a success rate of 0.890 with the lowest cost of 1155.843 . Despite these strengths, there is a slight trend of performance degradation as we progress through the groups, as one can see in both ρ scenarios.
When the colony is divided into g = 20 groups, as in Table 7, and ρ = 0.01 , group 1 emerges as the top performer, with a notable success rate of 0.92 accompanied by a competitive cost of 1171.355 . Group 2 closely follows, achieving a success rate of 0.90 and showing competitive values for the cost. Noteworthy is group 6, standing out with the lowest cost of 1115.392 , emphasizing its cost-effectiveness. Transitioning to ρ = 0.001 , group 1 maintains a high success rate of 0.92 with a slightly increased cost of 1226.367 . Group 5 exhibits the lowest cost of 1059.165 . Conversely, group 4 boasts the highest success rate of 0.96 coupled with a competitive value for the path cost. Additionally, there is a noticeable trend of performance degradation, especially under ρ = 0.01 .
Finally, when the colony is divided into g = 50 groups, as in Table 8, under ρ = 0.01 group 1 stands out, with a perfect success rate of 1.000 along with a competitive cost of 1162.1375 . Subsequent groups closely follow with relatively high success rates and competitive values for cost and exit time. Notably, group 11 excels with the lowest path cost of 985.3046 . However, there is a moderate degradation in performance across successive groups. Transitioning to ρ = 0.001 , group 2 maintains a perfect success rate of 1.000 , coupled with a path cost of 1111.2629 . Meanwhile, group 6 emerges as the top performer, showing a perfect success rate of 1.000 and the lowest cost of 938.1332 . This highlights the group’s remarkable ability to efficiently navigate and exit the environment. Overall, for both values of ρ , a diverse range of success rates, costs, and exit times is observed across the groups. Some groups excel in achieving high success rates, while others prioritize minimizing costs or achieving faster exit times, providing valuable insights into optimizing the management of larger ant colonies.

3.3. Pheromone Distribution Analysis

We now study the correlation between the temporal distribution of pheromones on the edges and the number of ants navigating the network during the simulation. This involves the systematic computation of the mean and standard deviation of pheromones on the edges, as well as counting how many ants are still in the network.
The mean and standard deviation of pheromones over time provide insights into how ants distribute the pheromone along the edges. The average quantity of pheromone deposited by ants is directly proportional to the number of ants in the network. As the number of ants increases, the quantity of pheromone deposited along the edges also increases. Conversely, when some ants exit the network, the pheromone is deposited in a lower quantity. The standard deviation of the pheromone, on the other hand, measures the degree of heterogeneity of the pheromone on the edges. A low value indicates a higher uniformity of the pheromone, leading to greater difficulty for the ants to exit the network since the edges have more or less the same level of pheromone. Conversely, high values denote a higher degree of heterogeneity of the pheromone and this information can be used by ants to exit the network. Finally, the number of ants that are present in the network at each instant of time allows us to better visualize the groups of ants that enter the network at regular intervals and their persistence in the network. This process was carried out for each group configuration and both values of ρ . The resulting data are presented through three plots discussed before, which are obtained by selecting the best experiment that maximizes the number of ants that exit the network, considering a total of 10 experiments performed for each configuration.
In Figure 3, the outcomes for g = 1 and n g = 200 are illustrated. In the initial phase of the simulation, there is an increase in the mean of pheromones, coinciding with active exploration by a substantial number of ants in the environment. Subsequently, a decline occurs, more pronounced, with an evaporation rate of 0.01 and a slower stabilization at 0.001 . This pattern underscores the deep connection between mean pheromone levels and the presence of ants in the network. As ants exit the network, there is a corresponding decrease in the pheromone quantity on edges, particularly clear at a higher evaporation rate. This decrease is attributed to the reduced ability of the remaining ants to counteract the ongoing pheromone evaporation. The second plot mirrors these dynamics, illustrating changes in the standard deviation. Also, in this scenario there is an initial rise in the value during the simulation’s early phase, signifying an increase in the heterogeneity of pheromones as ants release them at the edges. As the ants gradually exit the network, the standard deviation decreases when ρ = 0.01 and stabilizes when ρ = 0.001 . This is consistent with the idea that a smaller ant population in the network results in a diminished ability to counteract the evaporation rate, ultimately leading to a more uniform distribution of pheromone levels across the edges. The third plot depicts the evolving number of ants in the network over time, starting at 200, corresponding to the initial exploration group in the g = 1 configuration. As time progresses, ants exit with varied trends influenced by the ρ value. Notably, these trends undergo an inversion around 14 , 000 ticks.
In Figure 4, the plots for g = 2 and n g = 100 exhibit a trend similar to the previous configuration. Both mean pheromones and standard deviation show an initial increase during the simulation. Subsequently, a decline is observed, but only when ρ = 0.01 . Under ρ = 0.001 , both mean and standard deviation continue to rise over time. This discrepancy can be attributed to the more efficient exploratory activity of ants in this setup, resulting in fewer ants being present in the network at any given moment to counteract ongoing pheromone evaporation, particularly with the higher evaporation rate of ρ = 0.01 . The third plot illustrates the evolving number of ants in the network over time, starting at 100 and peaking at 200, indicating the entry of groups 1 and 2. As time progresses, ants exit more rapidly under ρ = 0.001 and slower under ρ = 0.01 .
The plots for g = 5 and n g = 40 in Figure 5 follow a similar trend. Both the mean and the standard deviation exhibit an initial increase during the simulation, followed by a decline, but only when ρ = 0.01 . Under ρ = 0.001 , both metrics reach stabilization. As ants gradually exit the environment, their ability to counteract pheromone evaporation is more effective with a low evaporation rate ( ρ = 0.001 ), resulting in a higher degree of pheromone heterogeneity. Conversely, when the evaporation rate is higher ( ρ = 0.01 ), the remaining ants in the network struggle to counteract pheromone evaporation, leading to a decrease in the mean and heterogeneity. Concerning the number of ants in the environment, it initiates at 40, with the fluctuations in the plot representing the intermittent entry of the five groups. Over time, ants exit at varying rates depending on the value of ρ .
In Figure 6 and Figure 7, the plots for the configurations g = 10 and N g = 20 , as well as g = 20 and n g = 10 , exhibit a trend similar to previous configurations. For both configurations, both the mean and the standard deviation show an initial increase during the simulation, followed by a decline, but only under ρ = 0.01 . Under ρ = 0.001 , both metrics in both cases reach a stabilization. In other words, as ants gradually exit the environment their ability to counteract pheromone evaporation is more effective with a low evaporation rate ( ρ = 0.001 ), resulting in a prolonged period of higher pheromone heterogeneity. Conversely, when the evaporation rate is higher ( ρ = 0.01 ) the remaining ants in the network struggle to counteract pheromone evaporation, leading to a decrease in mean and standard deviation. Regarding the number of ants in the environment, they initiate at 20 and 10, respectively, with the fluctuations in the plot representing the intermittent entry of the other groups (10 and 20, respectively). For g = 10 , it is clear that over time ants exit more rapidly when ρ = 0.01 than when ρ = 0.001 , suggesting that the ants increasingly benefit from the remaining pheromones on the edges as time unfolds. In contrast, for g = 20 this pattern holds true but only until around 11,000 ticks, beyond which the two trends reverse.
In Figure 8, for the configuration g = 50 and n g = 4 , both the mean and standard deviation show a continuous increase during the simulation when ρ = 0.001 . However, with ρ = 0.01 , the trend in the standard deviation, after an initial phase of increase, starts to decrease, followed by stabilization and a slight subsequent increase. This is in contrast to the mean, which experiences a more pronounced decrease. Furthermore, in this configuration a high degree of heterogeneity is sustained for a more extended period. These behaviors are closely linked to the evolving number of ants in the network, as observed in the third plot. Initially (for 0 < T < 5000 ), groups of four ants engage in exploration, swiftly exiting with no noticeable peaks in the curve. Concurrently, other groups are launched, but their rapid exit is not reflected in the curve. In the subsequent time interval ( 5000 < T < 10,000), there is an increase in the number of ants, marked by fluctuations in the curves. Notably, higher values are observed under ρ = 0.01 , while lower values are seen at ρ = 0.001 . Around 11,000 ticks, the number of ants gradually starts to decrease.

4. Discussion

The observed trends across the two metrics in the overall analysis are deeply connected to the interplay between the number of groups and the evaporation rate. In terms of success rate, the dynamics shift at different evaporation rates. At a higher rate ( ρ = 0.01 ) the optimal performance is achieved with 10 groups, suggesting that a more consolidated effort among a moderate number of groups is favored for efficient exploration. Conversely, at a lower evaporation rate ( ρ = 0.001 ) the extended lifespan of pheromones promotes sustained collaboration across 20 groups, leading to improved overall efficiency. The trend in path cost reflects a delicate balance between the collaboration benefits and potential drawbacks of increased group size. At a higher evaporation rate ( ρ = 0.01 ) the optimal path cost is achieved with 10 groups, denoting an optimal level of cooperation. Conversely, at a lower evaporation rate ( ρ = 0.001 ), sustained collaboration across 20 groups becomes beneficial due to the extended period of information retention. The slower evaporation rate facilitates distributed collaboration, enhancing exploration efficiency and reducing overall path costs. In this scenario, the persistence of pheromones enables distributed ant groups to communicate better and share information.
In the group analysis, operating as a single group, the colony displayed a remarkable level of consistency, showing comparable success rates and costs at different evaporation rates ( ρ values). This consistency suggests that the collective behavior of the colony is inherently robust, adapting seamlessly to varying environmental conditions. Upon division into two groups, a clear dynamic emerged. Group 1 outperformed group 2 when the evaporation rate was higher and vice versa. However, the colony showed an overall consistent performance, hinting at the robustness of its collective decision-making processes. Expanding the analysis to five groups revealed a standout performance by group 1, achieving the highest success rate along with a competitive cost. However, the subsequent groups exhibited signs of performance degradation, highlighting the complexity of managing larger ant colonies with diverse group compositions. Similar trends persisted when the colony was divided into ten groups, with specific groups distinguishing themselves in terms of success rates and costs. Notably, there was a discernible trend of performance degradation progress through the groups. This suggests that as the number of groups increases, the coordination challenges become more pronounced, potentially affecting the overall effectiveness of the colony. Further dividing the colony into twenty and fifty groups exposed a more intricate landscape of performance variations. While certain groups consistently excelled, demonstrating high success rates and cost-effectiveness, others displayed a noticeable decline in performance. Overall, it is noteworthy that the optimization process does not reach its best with the final group of ants in each configuration. In other words, the group achieving the optimal values in terms of success rate or path cost is typically one of the initial groups. Subsequent groups tend to exhibit a decline in performance after one of these early groups attains optimal values in both evaluation metrics.
Finally, in the pheromone distribution analysis, in general, the plots across all configurations exhibit comparable trends. Initially, there is a noticeable increase in mean pheromone levels, corresponding to active exploration and a significant presence of ants. However, a subsequent decline becomes apparent, particularly significant when there is a higher evaporation rate ( ρ = 0.01 ). This decrease can be attributed to the reduced capacity of the remaining ants to counteract the ongoing evaporation of pheromones. In contrast, lower evaporation rates ( ρ = 0.001 ) result in more extended periods of elevated pheromone levels. The standard deviation mirrors this pattern, suggesting heightened heterogeneity during exploration and a more uniform distribution as ants exit. The number of ants in the network undergoes fluctuations over time, reaching a peak corresponding to the entry of the last group into the environment. As various groups enter and exit, the network experiences dynamic changes. Generally, a lower evaporation rate leads to a swifter and more efficient exit of ants from the network. These findings shed light on why, in group analysis, the groups that perform exceptionally well are typically among the first rather than the last in each configuration. What emerges is that once a particular group within the considered configuration discovers or achieves the optimal value in terms of success rate and cost, subsequent groups exhibit a decline in their performance. This phenomenon could be attributed to the non-constant presence of groups and ants within the environment over time. Groups initiate exploration at regular intervals, but as time progresses ants exit the environment. Consequently, the activity of pheromone evaporation becomes increasingly significant because the remaining ants within the environment struggle to counteract this process. Consequently, there is yet again an increase in the homogeneity of pheromone distribution across the edges. As a result, each path becomes similar for the ants, impeding their ability to optimize path choices.

5. Conclusions

This paper introduces an agent-based model based on the ACO paradigm where ants are memory-enhanced to explore how partitioning a colony into groups influences the algorithm’s performance. The study specifically delves into the roles of group size and number in the context of a network analysis. Within the model, a collective of ants, starting from a specified location, collaboratively works to find the optimal path leading to the exit within a grid network. The overarching objectives of the colony include maximizing the number of ants successfully leaving and minimizing path costs. Three different kinds of analysis were conducted: an overall analysis, to analyze the performance of the whole colony concerning the two evaluation metrics when divided into groups of different sizes; a group analysis, to investigate the performance of each group into which the colony is partitioned; and a pheromone distributions analysis, to discern any correlation between the temporal distribution of pheromones on the edges and the number of ants navigating the network during the simulation.
In the overall analysis, a nuanced relationship emerged between the number of groups and the algorithm’s performance. At a higher evaporation rate, the optimal configuration involved 10 groups, obtaining the best success rate and path costs. It follows that under conditions of accelerated pheromone decay a moderate level of colony partitioning enhances effective collaboration. Conversely, at a lower evaporation rate, a superior performance was observed with 20 groups, emphasizing the role of prolonged pheromone lifespan in facilitating distributed collaboration and reducing overall costs. This finding underscores the significance of pheromone persistence, particularly in scenarios where environmental conditions permit longer-lasting pheromones.
The group analysis delved deeper into the performance of individual groups. Notably, the colony operating as a single group demonstrated consistent adaptability across different evaporation rates. However, as the number of groups increased, a complex landscape of performance rates emerged. Optimal values in success rate or path cost were typically achieved by one of the initial groups, with subsequent groups showing a decline in performance. This highlights the coordination challenges that escalate with the number of groups and their interactions, emphasizing the need for careful consideration in partitioning colonies, especially in larger-scale scenarios, to mitigate challenges associated with managing diverse groups effectively.
The pheromone distribution analysis provided additional insights into the correlation between pheromone dynamics, ant presence, and group behavior. As the number of ants increased, mean pheromone levels collectively rose, indicating active exploration. However, a subsequent decline, particularly noticeable at higher evaporation rates, emphasized the challenge of sustaining the pheromone quantity over time. The standard deviation analysis added depth to this understanding, revealing that a more uniform distribution, as indicated by lower values, could make it harder for ants to differentiate paths and exit efficiently. Conversely, higher standard deviation values suggested greater heterogeneity, enabling ants to effectively use pheromone signals for their exploration. The temporal dynamics of group entry and exit further highlighted the transient nature of group presence. Peaks in mean pheromone levels aligned with group entry, followed by decreases indicating exits. This temporal pattern underscores the dynamic nature of group interactions and their impact on pheromone distribution.
In summary, our study reveals a correlation between the partitioning of ant colonies and the quality of solutions within an ACO framework, shedding light on the complex interplay between individual ant behaviors and collective dynamics. The adaptability of ant colonies, discerned through the optimal degree of partitioning influenced by pheromone lifespan and group entry into the environment, underscores their responsiveness to diverse environmental conditions. The importance of the evaporation rate of pheromones emerges as a key factor in comprehending group dynamics during maze navigation. When pheromones persist for a longer time, dividing the colony into more groups becomes advantageous. This is because it allows for a more distributed and efficient exploration of the maze. The surplus of pheromones in the network, if not properly managed, could hinder exploration by leading ants in unproductive directions. Conversely, the heightened performance of initial groups highlights the critical role of the early exploration stages. The delicate balance between pheromone presence and evaporation rate fosters effective communication, enabling these groups to efficiently identify and reinforce optimal paths. Yet, as the simulation progresses and ants exit the network, subsequent groups experience performance declines due to a diminished ability to counteract pheromone evaporation. The observed variability in optimal configurations underscores the need for adaptive strategies to optimize ACO performance in evolving environmental conditions. The demonstrated adaptability of the Memory-Enhanced Ant Colony Optimization (ME-ACO) framework holds practical implications, particularly in fields such as logistics, supply chain management, and energy distribution. Dynamically adjusting group configurations could optimize transportation routes in response to changing traffic conditions and seasonal demands. Similarly, in smart grids and energy management, adaptability could enhance efficiency by responding to fluctuating demands and renewable energy availability. Expanding beyond optimization, the ME-ACO framework can be used in modeling and understanding collective behaviors in diverse scenarios. Its flexibility might be crucial in capturing the complex and dynamic interactions among social agents as they interact with the environment and each other. Utilizing an indirect communication mechanism, these interactions give rise to the emergence of collective behaviors, providing a valuable perspective on social dynamics across various contexts. Future extensions of the presented framework may include a sensitivity analysis of parameters to explore model dynamics under different configurations, additional experiments with more complex environments, possibly incorporating benchmark networks from the literature, and the consideration of a dynamic version of the framework where group number and size adapt dynamically to environmental conditions.

Author Contributions

All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the project Future Artificial Intelligence Research (FAIR)-PNRR MUR Cod. PE0000013—CUP: E63C22001940006.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Plotnikova, M.; Ulceluse, M. Inequality as a driver of migration: A social network analysis. Popul. Space Place 2022, 28, e2497. [Google Scholar] [CrossRef]
  2. McAndrew, D. The structural analysis of criminal networks. In The Social Psychology of Crime; Routledge: London, UK, 2021; pp. 51–94. [Google Scholar]
  3. Miyakoshi, T.; Shimada, J.; Li, K.W. A network analysis on country and financial center attractiveness: Evidence from Asian economies, 2001–2018. Int. Rev. Econ. Financ. 2023, 87, 418–432. [Google Scholar] [CrossRef]
  4. Engel, J.; Nardo, M.; Rancan, M. Network Analysis for Economics and Finance: An Application to Firm Ownership. In Data Science for Economics and Finance: Methodologies and Applications; Springer International Publishing: Cham, Switherland, 2021; pp. 331–355. [Google Scholar] [CrossRef]
  5. Gönçer-Demiral, D.; İnce-Yenilmez, M. Network analysis of international export pattern. Soc. Netw. Anal. Min. 2022, 12, 156. [Google Scholar] [CrossRef]
  6. Matchado, M.S.; Lauber, M.; Reitmeier, S.; Kacprowski, T.; Baumbach, J.; Haller, D.; List, M. Network analysis methods for studying microbial communities: A mini review. Comput. Struct. Biotechnol. J. 2021, 19, 2687–2698. [Google Scholar] [CrossRef] [PubMed]
  7. Desquilles, L.; Musso, O. Metabolic Networks: Weighted Gene Correlation Network Analysis. In Metabolic Reprogramming: Methods and Protocols; Springer: Berlin/Heidelberg, Germany, 2023; pp. 317–325. [Google Scholar] [CrossRef]
  8. Jiang, D.; Armour, C.R.; Hu, C.; Mei, M.; Tian, C.; Sharpton, T.J.; Jiang, Y. Microbiome multi-omics network analysis: Statistical considerations, limitations, and opportunities. Front. Genet. 2019, 10, 995. [Google Scholar] [CrossRef]
  9. Kosvyra, A.; Ntzioni, E.; Chouvarda, I. Network analysis with biological data of cancer patients: A scoping review. J. Biomed. Inform. 2021, 120, 103873. [Google Scholar] [CrossRef] [PubMed]
  10. Jo, W.; Lee, J.; Park, J.; Kim, Y. Online information exchange and anxiety spread in the early stage of the novel coronavirus (COVID-19) outbreak in South Korea: Structural topic model and network analysis. J. Med. Internet Res. 2020, 22, e19455. [Google Scholar] [CrossRef]
  11. Urman, A.; Katz, S. What they do in the shadows: Examining the far-right networks on Telegram. Information, Commun. Soc. 2022, 25, 904–923. [Google Scholar] [CrossRef]
  12. Cutello, V.; Fargetta, G.; Pavone, M.; Scollo, R.A. Optimization algorithms for detection of social interactions. Algorithms 2020, 13, 139. [Google Scholar] [CrossRef]
  13. Wu, Z.; Liao, Q.; Liu, B. A comprehensive review and evaluation of computational methods for identifying protein complexes from protein–protein interaction networks. Briefings Bioinform. 2020, 21, 1531–1548. [Google Scholar] [CrossRef]
  14. Wimalagunasekara, S.S.; Weeraman, J.W.; Tirimanne, S.; Fernando, P.C. Protein-protein interaction (PPI) network analysis reveals important hub proteins and sub-network modules for root development in rice (Oryza sativa). J. Genet. Eng. Biotechnol. 2023, 21, 1–15. [Google Scholar] [CrossRef]
  15. Bonabeau, E. Social Insect Colonies as Complex Adaptive Systems. Ecosystems 1998, 1, 437–443. [Google Scholar] [CrossRef]
  16. Carmichael, T.; Hadžikadić, M. The Fundamentals of Complex Adaptive Systems; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar] [CrossRef]
  17. Dorigo, M.; Gambardella, L. Ant colony system: A cooperative learning approach to the traveling salesman problem. IEEE Trans. Evol. Comput. 1997, 1, 53–66. [Google Scholar] [CrossRef]
  18. Deng, W.; Xu, J.; Zhao, H. An Improved Ant Colony Optimization Algorithm Based on Hybrid Strategies for Scheduling Problem. IEEE Access 2019, 7, 20281–20292. [Google Scholar] [CrossRef]
  19. Wang, C.; Zhao, B.; Peng, W.; Wu, C.; Gong, Z. Routing Algorithm Based on Ant Colony Optimization for DTN Congestion Control. In Proceedings of the 2012 15th International Conference on Network-Based Information Systems, Melbourne, VIC, Australia, 26–28 September 2012; pp. 715–720. [Google Scholar] [CrossRef]
  20. Lin, T.L.; Chen, Y.S.; Chang, H.Y. Performance Evaluations of an Ant Colony Optimization Routing Algorithm for Wireless Sensor Networks. In Proceedings of the 2014 Tenth International Conference on Intelligent Information Hiding and Multimedia Signal Processing, Kitakyushu, Japan, 27–29 August 2014; pp. 690–693. [Google Scholar] [CrossRef]
  21. Samà, M.; D’Ariano, A.; Pacciarelli, D.; Pellegrini, P.; Rodriguez, J. Ant colony optimization for train routing selection: Operational vs tactical application. In Proceedings of the 2017 5th IEEE International Conference on Models and Technologies for Intelligent Transportation Systems (MT-ITS), Naples, Italy, 26–28 June 2017; pp. 297–302. [Google Scholar] [CrossRef]
  22. Sun, Y.; Dong, W.; Chen, Y. An Improved Routing Algorithm Based on Ant Colony Optimization in Wireless Sensor Networks. IEEE Commun. Lett. 2017, 21, 1317–1320. [Google Scholar] [CrossRef]
  23. Jia, Y.H.; Mei, Y.; Zhang, M. A Bilevel Ant Colony Optimization Algorithm for Capacitated Electric Vehicle Routing Problem. IEEE Trans. Cybern. 2021, 52, 10855–10868. [Google Scholar] [CrossRef]
  24. Consoli, P.; Collerà, A.; Pavone, M. Swarm Intelligence heuristics for Graph Coloring Problem. In Proceedings of the 2013 IEEE Congress on Evolutionary Computation, Cancun, Mexico, 20–23 June 2013; pp. 1909–1916. [Google Scholar] [CrossRef]
  25. Fidanova, S.; Pop, P. An improved hybrid ant-local search algorithm for the partition graph coloring problem. J. Comput. Appl. Math. 2016, 293, 55–61. [Google Scholar] [CrossRef]
  26. Mostafaie, T.; Modarres Khiyabani, F.; Navimipour, N.J. A systematic study on meta-heuristic approaches for solving the graph coloring problem. Comput. Oper. Res. 2020, 120, 104850. [Google Scholar] [CrossRef]
  27. Akka, K.; Khaber, F. Mobile robot path planning using an improved ant colony optimization. Int. J. Adv. Robot. Syst. 2018, 15, 1729881418774673. [Google Scholar] [CrossRef]
  28. Zhang, D.; You, X.; Liu, S.; Pan, H. Dynamic Multi-Role Adaptive Collaborative Ant Colony Optimization for Robot Path Planning. IEEE Access 2020, 8, 129958–129974. [Google Scholar] [CrossRef]
  29. Jovanovic, R.; Tuba, M.; Voß, S. An efficient ant colony optimization algorithm for the blocks relocation problem. Eur. J. Oper. Res. 2019, 274, 78–90. [Google Scholar] [CrossRef]
  30. Man-Chung Yuen, S.C.N.; Leung, M.F. A Competitive Mechanism Multi-Objective Particle Swarm Optimization Algorithm and Its Application to Signalized Traffic Problem. Cybern. Syst. 2021, 52, 73–104. [Google Scholar] [CrossRef]
  31. Peng, H.; Ying, C.; Tan, S.; Hu, B.; Sun, Z. An Improved Feature Selection Algorithm Based on Ant Colony Optimization. IEEE Access 2018, 6, 69203–69209. [Google Scholar] [CrossRef]
  32. Crespi, C.; Fargetta, G.; Pavone, M.; Scollo, R.A. An Agent-Based Model to Investigate Different Behaviours in a Crowd Simulation. In Proceedings of the Bioinspired Optimization Methods and Their Applications: 10th International Conference, BIOMA 2022, Maribor, Slovenia, 17–18 November 2022; Springer: Berlin/Heidelberg, Germany, 2022; pp. 1–14. [Google Scholar] [CrossRef]
  33. Crespi, C.; Scollo, R.A.; Fargetta, G.; Pavone, M. A sensitivity analysis of parameters in an agent-based model for crowd simulations. Appl. Soft Comput. 2023, 146, 110684. [Google Scholar] [CrossRef]
  34. Ping, G.; Chunbo, X.; Yi, C.; Jing, L.; Yanqing, L. Adaptive ant colony optimization algorithm. In Proceedings of the 2014 International Conference on Mechatronics and Control (ICMC), Jinzhou, China, 3–5 July 2014; pp. 95–98. [Google Scholar] [CrossRef]
  35. Ning, J.; Zhang, C.; Sun, P.; Feng, Y. Comparative study of ant colony algorithms for multi-objective optimization. Information 2018, 10, 11. [Google Scholar] [CrossRef]
  36. Mu, C.; Zhang, J.; Liu, Y.; Qu, R.; Huang, T. Multi-objective ant colony optimization algorithm based on decomposition for community detection in complex networks. Soft Comput. 2019, 23, 12683–12709. [Google Scholar] [CrossRef]
  37. Chen, W.N.; Tan, D.Z.; Yang, Q.; Gu, T.; Zhang, J. Ant colony optimization for the control of pollutant spreading on social networks. IEEE Trans. Cybern. 2019, 50, 4053–4065. [Google Scholar] [CrossRef]
  38. Gao, S.; Wu, J.; Ai, J. Multi-UAV reconnaissance task allocation for heterogeneous targets using grouping ant colony optimization algorithm. Soft Comput. 2021, 25, 7155–7167. [Google Scholar] [CrossRef]
  39. Cavallaro, C.; Verga, G.; Tramontana, E.; Muscato, O. Multi-Agent Architecture for Point of Interest Detection and Recommendation. In Proceedings of the CEUR Workshop Proceedings; 2019; Volume 2404, pp. 98–104. Available online: https://ceur-ws.org/Vol-2404/paper15.pdf (accessed on 22 January 2024).
  40. Cavallaro, C.; Bujari, A.; Foschini, L.; Di Modica, G.; Bellavista, P. Measuring the impact of COVID-19 restrictions on mobility: A real case study from italy. J. Commun. Netw. 2021, 23, 340–349. [Google Scholar] [CrossRef]
  41. Wang, Z.; Yao, S.; Li, G.; Zhang, Q. Multiobjective Combinatorial Optimization Using a Single Deep Reinforcement Learning Model. IEEE Trans. Cybern. 2023. [Google Scholar] [CrossRef]
  42. Crespi, C.; Scollo, R.A.; Pavone, M. Effects of Different Dynamics in an Ant Colony Optimization Algorithm. In Proceedings of the 2020 7th International Conference on Soft Computing Machine Intelligence (ISCMI2020), Stockholm, Sweden, 14–15 November 2020; pp. 8–11. [Google Scholar] [CrossRef]
  43. Wilensky, U. NetLogo; Center for Connected Learning and Computer-Based Modeling, Northwestern University: Evanston, IL, USA, 1999. [Google Scholar]
Figure 1. Ant k positioned on node i perceives information I i j ( t ) (green dot) regarding the weight of edge ( i , s ) only when it is located at the nearest endpoint (a); otherwise, it remains unaware of it (b).
Figure 1. Ant k positioned on node i perceives information I i j ( t ) (green dot) regarding the weight of edge ( i , s ) only when it is located at the nearest endpoint (a); otherwise, it remains unaware of it (b).
Algorithms 17 00063 g001
Figure 2. In (a), an ant k situated at node i assesses the desirability of adjacent edges as equal to 1 in the absence of both prior ( w p k = 0 ) and local knowledge ( w l k = 0 ). In (b), when ant k reaches node i from node i 1 , it evaluates adjacent edges. If information I i j ( t ) on edge weights is either missing or not visible, the desirability is determined as the inverse of global knowledge. When the information I i j ( t ) is visible, the desirability is assessed as the inverse of the weight of the link.
Figure 2. In (a), an ant k situated at node i assesses the desirability of adjacent edges as equal to 1 in the absence of both prior ( w p k = 0 ) and local knowledge ( w l k = 0 ). In (b), when ant k reaches node i from node i 1 , it evaluates adjacent edges. If information I i j ( t ) on edge weights is either missing or not visible, the desirability is determined as the inverse of global knowledge. When the information I i j ( t ) is visible, the desirability is assessed as the inverse of the weight of the link.
Algorithms 17 00063 g002
Figure 3. Plots obtained with one group ( g = 1 ) and two hundred ants per group ( n g = 200 ).
Figure 3. Plots obtained with one group ( g = 1 ) and two hundred ants per group ( n g = 200 ).
Algorithms 17 00063 g003
Figure 4. Plots obtained with two groups ( g = 2 ) and one hundred ants per group ( n g = 100 ).
Figure 4. Plots obtained with two groups ( g = 2 ) and one hundred ants per group ( n g = 100 ).
Algorithms 17 00063 g004
Figure 5. Plots obtained with five groups ( g = 5 ) and forty ants per group ( n g = 40 ).
Figure 5. Plots obtained with five groups ( g = 5 ) and forty ants per group ( n g = 40 ).
Algorithms 17 00063 g005
Figure 6. Plots obtained with ten groups ( g = 10 ) and twenty ants per group ( n g = 20 ).
Figure 6. Plots obtained with ten groups ( g = 10 ) and twenty ants per group ( n g = 20 ).
Algorithms 17 00063 g006
Figure 7. Plots obtained with twenty groups ( g = 20 ) and ten ants per group ( n g = 10 ).
Figure 7. Plots obtained with twenty groups ( g = 20 ) and ten ants per group ( n g = 10 ).
Algorithms 17 00063 g007
Figure 8. Plots obtained with fifty groups ( g = 50 ) and four ants per group ( n g = 4 ).
Figure 8. Plots obtained with fifty groups ( g = 50 ) and four ants per group ( n g = 4 ).
Algorithms 17 00063 g008
Table 1. Description of variables and parameters used in the proposed model.
Table 1. Description of variables and parameters used in the proposed model.
VariableDescription
w ( i , j ) weight of an edge
p i j k ( t ) transition probability of the ants
τ i j ( t ) pheromone intensity on the edge
η i j ( t ) desirability of an edge
I i j ( t ) information about the weight of an edge
α importance of pheromone intensity
β importance of desirability
Kquantity of pheromone intensity left by ants
ρ evaporation rate
nnumber of ants
gnumber of groups
n g number of ants per group
w ¯ p k prior knowledge of an ant
w ¯ l k local knowledge of an ant
w ¯ global knowledge of an ant
T e exit interval between groups
T i maximum time to reach the exit for ith group
T m a x maximum time of simulations
Table 2. Results obtained based on two evaluation metrics: success rate (on the top) and path cost (on the bottom).
Table 2. Results obtained based on two evaluation metrics: success rate (on the top) and path cost (on the bottom).
Success Rate
g ρ = 0.01 ρ = 0.001
10.8490.848
20.8550.868
50.8620.8765
100.86850.8665
200.81250.8875
500.7950.8645
Path Cost
g ρ = 0.01 ρ = 0.001
11293.531311.20
21276.091265.81
51249.521251.05
101242.631264.48
201339.241201.36
501289.451227.27
Table 3. Results obtained with one group ( g = 1 ) .
Table 3. Results obtained with one group ( g = 1 ) .
ρ = 0.01 ρ = 0.001
GroupSuccess RateCostSuccess RateCost
10.8491293.530.8481311.197
Table 4. Results obtained with two groups ( g = 2 ) .
Table 4. Results obtained with two groups ( g = 2 ) .
ρ = 0.01 ρ = 0.001
GroupSuccess RateCostSuccess RateCost
10.8561270.5710.8601312.964
20.8541281.6390.8761220.156
Table 5. Results obtained with five groups ( g = 5 ) .
Table 5. Results obtained with five groups ( g = 5 ) .
ρ = 0.01 ρ = 0.001
GroupSuccess RateCostSuccess RateCost
10.91251201.6900.91001210.888
20.87001206.3680.87751235.365
30.83501266.7260.86501301.793
40.87001243.5340.85751268.863
50.82251336.4510.87251240.825
Table 6. Results obtained with ten groups ( g = 10 ) .
Table 6. Results obtained with ten groups ( g = 10 ) .
ρ = 0.01 ρ = 0.001
GroupSuccess RateCostSuccess RateCost
10.9251251.0040.9051227.731
20.8501288.6510.8901155.843
30.8701178.1440.8801217.516
40.8951214.6120.8601287.937
50.8851180.4990.8701299.610
60.9051282.3150.8801217.064
70.8301259.4520.8451306.233
80.8501198.3900.8401311.303
90.8301310.1300.8551303.442
100.8451254.9350.8601331.005
Table 7. Results obtained with twenty groups ( g = 20 ) .
Table 7. Results obtained with twenty groups ( g = 20 ) .
ρ = 0.01 ρ = 0.001
GroupSuccess RateCostSuccess RateCost
10.921171.3550.921226.367
20.901233.4590.911118.648
30.881255.3810.951083.556
40.861294.7750.961080.811
50.831263.1560.901059.165
60.891115.3920.901166.186
70.881246.5940.891075.165
80.771394.6490.891217.577
90.851350.0210.921310.244
100.831334.2830.911169.145
110.771439.8800.881286.262
120.801439.4350.911191.761
130.901129.5090.871176.769
140.771504.2200.831423.822
150.751419.9980.841276.752
160.691720.0040.811240.255
170.771385.7820.891247.258
180.721397.6290.871226.935
190.751353.4570.861240.497
200.721571.2630.841258.229
Table 8. Results obtained with fifty groups ( g = 50 ) .
Table 8. Results obtained with fifty groups ( g = 50 ) .
ρ = 0.01 ρ = 0.001
GroupSuccess RateCostSuccess RateCost
11.0001162.13750.9501135.9887
20.9251158.70231.0001111.2629
30.9751163.79950.9751131.4784
40.9251206.36860.9001380.4290
50.9501048.18860.9001298.5476
60.8751135.87861.000938.1332
70.9251233.23640.9251068.3007
80.9251049.11580.9751144.3801
90.9251188.19060.9501105.7798
100.9251127.65580.9501077.9356
110.950985.30460.9751173.0694
120.8501138.14770.9251045.7094
130.8501327.33940.9001314.4445
140.9001183.35870.9251185.8384
150.8501346.72650.9501078.2593
160.8001243.68180.8751234.9214
170.8001344.95760.9501236.5639
180.8251319.19440.8501407.0209
190.9001128.43600.8751211.4137
200.8751343.88260.9001131.3349
210.6751666.87710.9001112.5776
220.7751218.70890.9251141.0186
230.7751496.62440.9001312.2767
240.8251388.54350.9001088.1684
250.6751435.58380.8751215.7378
260.8001136.25910.9501107.3724
270.6501282.38040.8501325.3072
280.8001288.66050.9001056.2103
290.8001151.41920.8751104.8473
300.7751276.22880.8751242.3154
310.72513729.1380.7501361.5867
320.8501198.27540.8501176.4607
330.7251352.96150.8251108.4895
340.7751269.99170.8251483.7002
350.7251479.83740.8751050.3742
360.7001409.37470.8001428.6403
370.7001245.40930.7251462.6775
380.8001120.86370.7751484.8115
390.7251501.10830.7001359.2875
400.7001484.29700.8251152.9818
410.7501403.84160.8751118.3799
420.7501351.12520.8001334.9016
430.6251568.72990.7001496.3216
440.7501352.69760.7001432.3804
450.7001427.53470.8751116.7630
460.5501753.95320.7001416.9457
470.7251300.91520.6751488.6332
480.6251443.72470.77515,315.673
490.6001653.46990.87513,732.195
500.7501155.25760.7001360.1492
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cavallaro, C.; Crespi, C.; Cutello, V.; Pavone, M.; Zito, F. Group Dynamics in Memory-Enhanced Ant Colonies: The Influence of Colony Division on a Maze Navigation Problem. Algorithms 2024, 17, 63. https://doi.org/10.3390/a17020063

AMA Style

Cavallaro C, Crespi C, Cutello V, Pavone M, Zito F. Group Dynamics in Memory-Enhanced Ant Colonies: The Influence of Colony Division on a Maze Navigation Problem. Algorithms. 2024; 17(2):63. https://doi.org/10.3390/a17020063

Chicago/Turabian Style

Cavallaro, Claudia, Carolina Crespi, Vincenzo Cutello, Mario Pavone, and Francesco Zito. 2024. "Group Dynamics in Memory-Enhanced Ant Colonies: The Influence of Colony Division on a Maze Navigation Problem" Algorithms 17, no. 2: 63. https://doi.org/10.3390/a17020063

APA Style

Cavallaro, C., Crespi, C., Cutello, V., Pavone, M., & Zito, F. (2024). Group Dynamics in Memory-Enhanced Ant Colonies: The Influence of Colony Division on a Maze Navigation Problem. Algorithms, 17(2), 63. https://doi.org/10.3390/a17020063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop