Natural systems are a source of inspiration for computer algorithms designed to solve optimisation problems. Yet most ‘nature-inspired’ algorithms take only superficial inspiration from biology, and little is known about how real biological systems solve difficult problems. Moreover, ant algorithms, neural networks and similar methods are usually applied to static problems, whereas most biological systems have evolved to perform under dynamically changing conditions. We used the Towers of Hanoi puzzle to test whether Argentine ants can solve a potentially difficult optimisation problem. We also tested whether the ants can adapt to dynamic changes in the problem. We mapped all possible solutions to the Towers of Hanoi on a single graph and converted this into a maze for the ants to solve. We show that the ants are capable of solving the Towers of Hanoi, and are able to adapt when sections of the maze are blocked off and new sections installed. The presence of exploration pheromone increased the efficiency of the resulting network and increased the ants' ability to adapt to changing conditions. Contrary to previous studies, our study shows that mass-recruiting ant species such as the Argentine ant can forage effectively in a dynamic environment. Our results also suggest that novel optimisation algorithms can benefit from stronger biological mimicry.

Routing telephone calls through busy networks while minimising connection time, constructing complex machinery while keeping costs and build-time low, and finding the most efficient set-down and pick-up routes for delivery vehicles are all examples of combinatorial optimisation problems (Gonzalez and Sahni, 1976; Alon and Srinivasan, 1997; Bell and McMullen, 2004). These problems become increasingly difficult to solve as the size of the problem increases. In many cases, the difficulty of the problem increases exponentially with the number of components in the system. In the worst case these problems are ‘NP-hard’ (‘non-deterministic polynomial-time hard’), so that there is no known algorithm for solving general instances of the problem within a reasonable time frame. As a result, programmers must exploit ‘heuristic’ algorithms that find near-optimal solutions in a reasonable time.

Natural systems have proved to be a rich source of inspiration for computer scientists in designing optimisation algorithms (Bonabeau et al., 2000; Vassiliadis and Dounias, 2009). Many biological systems have been refined through millions of years of natural selection to efficiently exploit the ephemeral, often fiercely contested and spatially isolated resources of their environment. Such systems often hinge upon the construction of efficient transportation networks connecting the resources. Such networks have been studied at all levels of biological organisation: colonies of ants (Dussutour et al., 2004; Buhl et al., 2009) and termites (Perna et al., 2008), fungal mycelia (Bebber et al., 2007), acellular slime mould (Nakagaki et al., 2000; Nakagaki et al., 2001; Nakagaki et al., 2004; Latty and Beekman, 2009; Tero et al., 2010) and even the vasculature of plants and animals (Banavar et al., 1999; West et al., 1999). These complex systems are constructed in the absence of any central control by many individual autonomous components possessing only local information. It is this ‘swarm intelligence’ (Bonabeau et al., 1999) that has been widely used in a variety of computing techniques.

Probably the best-known nature-inspired algorithm used for NP-hard problems is Ant Colony Optimisation (ACO) (Dorigo and Stützle, 2004). ACO was inspired by the foraging behaviour of trail-laying ants (Dorigo et al., 1996). Many ant species construct foraging networks by laying pheromone trails towards food sources (Hölldobler and Wilson, 1990). Scouts deposit pheromones linking a newly discovered food source to the nest. The presence of the pheromone recruits nestmates to follow the trail, laying their own pheromone in turn and leading to a process of amplification. In addition, the trail pheromone is volatile and gradually evaporates, thus requiring constant reinforcement. As a result, longer trails are less competitive and eventually only shorter trails will be selected (Goss et al., 1989; Beckers et al., 1990). The ACO method for solving shortest path problems deploys virtual ‘ants’ to explore all possible routes, depositing pheromone on each edge they travel. The amount of pheromone is inversely proportional to the length of the tour, so that shorter tours receive the most pheromone. Evaporation of pheromone after each tour and many reiterations of the process means eventually only the shortest tour remains, or one close to it (Dorigo et al., 1996). The ant-based system has since been widely adopted for use in other combinatorial optimisation problems (Liang and Smith, 2004; Akay and Toksari, 2009; Zhang et al., 2009).

Although ACO and other nature-inspired algorithms often prove competitive in solving specific optimisation problems, apart from the original studies on binary choice shortest path problems there has been surprisingly little work on how real ants solve combinatorial optimisation problems. Notable exceptions are the studies of Vittori et al. (Vittori et al., 2006) and Garnier et al. (Garnier et al., 2009), which examined how ants used the branching angles in shortest-path mazes to decide which path to take. They showed that Argentine ants use the polarity of their trail network to determine whether the path leads towards or away from the nest. Unfed foragers have a high tendency to make a U-turn when they encounter an asymmetrical bifurcation, whereas fed foragers will tend to make a U-turn when they confront a symmetrical bifurcation (Gerbier et al., 2008).

However, these studies did not provide the ants with a difficult optimisation problem. Moreover, little work has focused on how ants solve dynamic problems, that is, those that change over time. The problem of finding the shortest path to food in situations in which the food changes position or paths become blocked and reopened has been faced by ant colonies for millions of years. Understanding how the ants solve such dynamic problems could prove a fruitful source of new inspiration for ACO and other computing methods.

The current convention remains that many mass-recruiting ant species are unable to adapt to dynamic changes, such as exploiting a newly discovered superior food source when an inferior one is currently being used. The apparent inability of the ants to adapt to changing conditions is supported by laboratory experiments (Goss et al., 1989; Beckers et al., 1992; Traniello et al., 1995) and mathematical models (Goss et al., 1989; Nicolis and Deneubourg, 1999; Camazine et al., 2001). However, more recent studies also suggest that trail-laying ants are better able to adjust their foraging efforts than previously assumed (Dussutour et al., 2009a; Dussutour et al., 2009b). In particular, Dussutour et al. (Dussutour et al., 2009b) argue that allowing the ants to first explore an environment before food is placed within it improves their ability to react to changes in the position of the food.

Here, we use the ‘Towers of Hanoi’ to examine how ants construct a solution to a dynamic optimisation problem. In the Towers of Hanoi, N discs of different diameters are stacked upon one of three pegs in order of decreasing diameter from the base. The problem is solved optimally when the stack of discs is transferred to one of the remaining two pegs using the smallest number of possible moves while obeying the following rules: only one disc can be moved at a time; only the top disc in a stack can be moved; and a disc can never be stacked upon another of smaller diameter. The optimal solution to this version of the problem is a simple recursive algorithm requiring exactly 2N-1 moves (Buneman and Levy, 1980). However, with the addition of extra discs and a single extra peg, the problem becomes NP-hard (Korf and Felner, 2007). We used the three-disc, three-rod version of the Towers of Hanoi problem to test the optimisation abilities of Argentine ants (Linepithema humile Mayr). We chose this version of the problem because the set of possible moves and configurations forms an undirected graph that can easily be transformed into a maze with a clear entry point (the nest) and goal (the food) when linked to its mirror image [see Fig. 1A, adapted from Romik (Romik, 2006), and 1B]. We could thus translate the Towers of Hanoi problem into a shortest-path problem. Our maze contains 32,768 possible paths from one end to the other without turning an ant back to directly face the origin of travel. The two shortest paths are along the outside edges of the maze (Fig. 1B). We tested the ants' ability to adapt to changing conditions by obstructing the most-used shortest path(s) once a solution had been constructed by the ants. As the presence of ‘exploration’ pheromone (pheromone deposited while the ants explore their environment in the absence of food) might help a colony to better forage dynamically (Dussutour et al., 2009b), we conducted trials both in the presence and absence of exploration pheromone. Lastly, we tested the effect of colony size, as the number of ants is known to affect the ability of a colony to form pheromone trails (Beekman et al., 2001).

Fig. 1.

(A) Undirected graph of three-disc, three-rod Towers of Hanoi [adapted from Romik (Romik, 2006)]. The numbers 1, 2 and 3 represent discs of increasing size and their positions on the three rods a, b and c. The nodes signify distributions of discs, whereas the branches represent moves. (B) Towers of Hanoi maze. The dashed lines depict the minimal length paths from nest to food source.

Fig. 1.

(A) Undirected graph of three-disc, three-rod Towers of Hanoi [adapted from Romik (Romik, 2006)]. The numbers 1, 2 and 3 represent discs of increasing size and their positions on the three rods a, b and c. The nodes signify distributions of discs, whereas the branches represent moves. (B) Towers of Hanoi maze. The dashed lines depict the minimal length paths from nest to food source.

General procedures

We used Argentine ants collected from the grounds of the University of Sydney, New South Wales, Australia. We housed colonies in plastic boxes with sides coated with Fluon to prevent escape. Each box contained up to 20 cardboard-wrapped test tubes partially filled with water and a cotton wool plug. We fed the colonies with a vitamin-enriched diet (Bhatkar and Whitcomb, 1970) twice per week, supplemented with mealworms.

All experiments were performed in a temperature-controlled room with diffuse overhead lighting and with colonies of different ages. Experiments were run from May 2009 until March 2010. The mazes were uncovered to prevent pheromone saturation. At the end of each experiment, the mazes were sprayed with ethanol (90%) and washed in bleach (80%) to remove any trace of pheromone. The experimental setup was surrounded by white walls to avoid any disturbance and to prevent the ants from using visual cues to navigate.

Experimental procedure

Towers of Hanoi

We starved 30 colonies containing 500 workers and one queen, and 30 colonies containing 1000 workers and two queens for 3 days before each experiment. Each colony contained varied numbers of brood and males. To determine the effect of the presence of exploration pheromone on the ability of the ants to solve the Towers of Hanoi, half of the 500 and 1000 ant colonies (15 each) were given access to the experimental arena for at least two hours prior to the start of a trial. These were our pre-exposure trials.

We used a Perspex maze 1 m in length consisting of 64 identical hexagons aligned in a rhomboid array (Fig. 3). The underside of the maze and the stilts it was mounted upon were coated in Fluon to ensure the ants could only walk on the upper side of the maze. Each bifurcation within the maze was made symmetrical when coming from the nest, and asymmetrical when returning from the food source (Vittori et al., 2006; Gerbier et al., 2008; Garnier et al., 2009). To convert the maze into the undirected graph of the Towers of Hanoi problem, several paths were blocked off with a perpendicular wall of Fluon-coated plastic (Fig. 1B).

Ants were given access to the maze via a wooden skewer that was placed in the nest and rested against the inner side of the first bifurcation. We placed a food source (wax insect feeder containing 2 mol l–1 sucrose solution) at the distal end of the maze at the beginning of each experiment and recorded the construction of trails as described below for 1 h. After 1 h, we blocked off the connecting bridge that was most used by the ants using Fluon-coated plastic walls. This configuration is hereafter referred to as ‘Towers of Hanoi’ and is presented in Fig. 2A. We again recorded the construction of trails for 1 h, this time to determine how capable the ants were of constructing an alternative solution after the maze had been changed. This was performed for all 30 trials.

We alternated the order in which we tested the different treatments (with and without pre-exposure, 500 and 1000 ants). We also alternated which side of the maze contained the nest and the food source.

Variations on the Towers of Hanoi

To determine whether Argentine ants have an intrinsic bias when constructing solutions through a maze, we also performed experiments using an open 64-hexagon maze (referred to as ‘open maze’ hereafter, Fig. 3). Twelve colonies of 1000 workers and two queens were starved for 3 days without access to the maze. After 3 days, we placed food at the distal end of the maze and recorded the construction of trails for each of the 12 colonies for 1 h.

As ants in an open maze appeared to prefer constructing trails along the outside edges of the maze (which corresponds to the optimal solution to the Towers of Hanoi maze), we performed a variation on the Towers of Hanoi maze. As before, we allowed the ants to construct a solution for 1 h using the Towers of Hanoi maze. After 1 h, we blocked off the two connecting bridges and opened a new bridge (central bridge) in the middle of the maze. Thus, in this version of the maze, the ants needed to construct a solution that did not follow the outside edges. This version is hereafter referred to as ‘Modified Towers of Hanoi’ (Fig. 2B). We used 24 colonies of 1000 workers and two queens. All colonies were starved for 3 days prior to the experiments, with half of the colonies (12) having access to the maze to lay exploration pheromone. These were our pre-exposure trials.

Behaviour of individual ants

Because we found an effect of exploration pheromone on the networks that the ants constructed, we next compared the behaviour of individual ants in the presence and absence of exploration pheromone. We starved approximately 100 ants for 3 days in a separate queenless colony. We then placed a single worker from this colony on either the proximal or distal end of an open maze using a wooden skewer. Because we needed a smaller field of view in order to perform our analyses of individual behaviour, we used a truncated maze; the area available to the ant did not exceed 42 cm along the outside edge of the maze. We collected data as described below for 8 min, after which time the ant was removed and another placed on the opposite end of the maze to ensure the pheromone laid by the first individual did not interfere with the behaviour of the second. This was repeated until 15 ants had been tested. The entire maze was washed with ethanol (90%) and bleach (80%) to remove any pheromone after the testing of each pair of ants. To obtain a maze coated in exploration pheromone, the maze was offered to a colony of 500 workers and one queen for at least 2 h over the starvation period, and just prior to testing. The same colony was used to mark the mazes in these pre-exposure trials. We then tested 15 ants as described above. In the treatment with exploration pheromone, we did not place consecutive ants on opposite ends of the maze, but used a new maze after every 3 ants tested. This was done to prevent excessive evaporation of the pheromone over repeated experiments, as the mean lifetime of L. humile pheromone has been shown to be approximately 30 min (Deneubourg et al., 1990).

Fig. 2.

(A) Towers of Hanoi maze after dynamic change. In this case the upper connecting bridge has been blocked. Arena dimensions: a, 90 deg; b, 135 deg; c, 135 deg; d, 3 cm; e, 4 cm; f, 1.5 cm; g, 35 cm; h, 1 m. (B) Modified Towers of Hanoi maze. Both connecting bridges have been blocked, and a central bridge has been inserted in the centre of the maze. In both A and B, the dashed lines represent the new minimal length path(s). For B, there are four possible minimal paths, representing the possible combinations of upper and lower minimal paths in each half of the maze. The shaded paths in B represent example paths with corrective moves, which have been circled.

Fig. 2.

(A) Towers of Hanoi maze after dynamic change. In this case the upper connecting bridge has been blocked. Arena dimensions: a, 90 deg; b, 135 deg; c, 135 deg; d, 3 cm; e, 4 cm; f, 1.5 cm; g, 35 cm; h, 1 m. (B) Modified Towers of Hanoi maze. Both connecting bridges have been blocked, and a central bridge has been inserted in the centre of the maze. In both A and B, the dashed lines represent the new minimal length path(s). For B, there are four possible minimal paths, representing the possible combinations of upper and lower minimal paths in each half of the maze. The shaded paths in B represent example paths with corrective moves, which have been circled.

Fig. 3.

Open maze analysis. An observed trail accumulated a value dependent on the proximity of its constituent path segments relative to the outside edges of the maze. The dashed lines provide examples of two paths at opposite ends of the spectrum. The long-dashed trail travels through the centre of the maze and thus accumulates the lowest edge proximity score; the short-dashed trail traverses the outer edge of the open maze, and hence accumulates the largest score. The average edge proximity score represents the calculated proximity score for an ant making random left and right choices when travelling from one side of the maze to the other. Note that all possible paths continuing from proximal to distal end are of the same physical length.

Fig. 3.

Open maze analysis. An observed trail accumulated a value dependent on the proximity of its constituent path segments relative to the outside edges of the maze. The dashed lines provide examples of two paths at opposite ends of the spectrum. The long-dashed trail travels through the centre of the maze and thus accumulates the lowest edge proximity score; the short-dashed trail traverses the outer edge of the open maze, and hence accumulates the largest score. The average edge proximity score represents the calculated proximity score for an ant making random left and right choices when travelling from one side of the maze to the other. Note that all possible paths continuing from proximal to distal end are of the same physical length.

Data collection and analysis

Towers of Hanoi

We collected still images of the entire maze every 10 s with a Logitech 9000 webcam mounted 1 m above the maze. Images were collected for 1 h after adding the food source. After blocking the connecting bridge used by the ants we collected another hour of data. If both connecting bridges were used, we blocked one.

In order to easily visualise the trail networks, we condensed the collected images into 10 min overlays (each consisting of 60 images). We defined a trail as any continuous (i.e. no break in the dark area exceeded six pixels) procession of dark pixelation between the nest and food source. Pilot studies showed that this criterion allows us to distinguish between established trails and individual ants displayed in successive images in the overlay.

At each time point of 10 min, we calculated the proportion of all colonies with definable trails (refer to the n values in Figs 4, 6) that had at least one minimal solution. We also determined the total length of the trail network for each period of 10 min by summing the lengths (cm) of all segments containing a trail between nest and food source. Dead-end branches of trails that connected the nest and food source were included in the analysis and the dead-end segments counted twice to reflect that ants must turn back to continue on to the nest/food source. The number of times the ants constructed trails along the minimal path(s) of each maze was also counted. The final analysis involved counting the number of ‘corrective moves’ employed in the trail network. These are any moves between nodes in a direction perpendicular to the long axis of the maze. For example, moving between the node ‘cbc’ and the node ‘abc’ in the bottom right of Fig. 1A involves one corrective move, moving between ‘aac‘ and ‘ccc’ involves three corrective moves, and so on. Further examples are given in Fig. 2B. Any sub-optimal solution to the problem will thus contain at least one corrective move, the number of which will scale directly with the divergence from the optimal solution. We calculated the mean total network length (cm), the mean number of minimal paths and the mean number of corrective moves across all replicates for each 10 min period. Each metric was graphed over the 12 time points of the experiments (six time points before one trail was blocked and six afterwards). The first and second hour of the experiments were treated separately in the analyses. We used the Shapiro–Wilk test for normality, which showed none of the data to be sufficiently normal. Hence, we applied a Poisson transformation to the data [x′=√(x+0.5)] and then performed three-way ANOVA with repeated measures. The factors were colony size (500/1000), pre-exposure (with/without), and time as the repeated measure. We also compared the partial eta squared (η2p) values to compare the size of effect of each factor (colony size and pre-exposure).

Variations on the Towers of Hanoi

For the open maze experiments, trail segments were given a value based on their proximity to the edge of the maze, as shown in Fig. 3, providing a method with which to quantify any edge-biased behaviour exhibited by the ants. For comparison with data, we identified the path that gave the lowest possible edge proximity total, i.e. through the centre of the graph (proximity score=95), and that giving the highest possible edge proximity total, i.e. along the edges (proximity score=155) (Figs 3, 5). We also calculated the average edge proximity of the path taken by ants choosing left and right paths at random as they move through the maze. Ants choosing paths at random as they leave the nest will initially be biased toward the centre of the open maze. In particular, at the point of maximum maze width (i.e. where there are nine parallel segments), the probability an ant arrives on a segment i from the top is given by:
(1)
However, in the second half of the maze there is a bias towards ants moving to the edges of the maze, as once an ant has arrived on the edge it no longer faces any choices. The probability that an ant on a segment i from the top (assuming for now that i<5) is on a segment with proximity score 5 after L choices is:
(2)
which increases with L. We combined these first and second half calculations to find the probabilities that each edge in the graph is used, and then found the average proximity score of a randomly moving ant. This gave a proximity score of 113.6.
Fig. 4.

Results for Towers of Hanoi. (A) Proportion of colonies with at least one minimal path. (B) Total network length. ‘Minimal length’ indicates the length of the shortest path through the maze (108 cm). (C) Number of minimal paths used. (D) Number of corrective moves in the network. The vertical dashed line indicates the time of dynamic change. Shown in B–D are means ± s.e.m. The N values are the number of colonies that constructed trails at each time point. At the beginning of the experiment not all colonies had constructed trails, whereas at the end of the experiment some colonies ceased foraging due to satiation.

Fig. 4.

Results for Towers of Hanoi. (A) Proportion of colonies with at least one minimal path. (B) Total network length. ‘Minimal length’ indicates the length of the shortest path through the maze (108 cm). (C) Number of minimal paths used. (D) Number of corrective moves in the network. The vertical dashed line indicates the time of dynamic change. Shown in B–D are means ± s.e.m. The N values are the number of colonies that constructed trails at each time point. At the beginning of the experiment not all colonies had constructed trails, whereas at the end of the experiment some colonies ceased foraging due to satiation.

For the Modified Towers of Hanoi, we performed two-way ANOVA on Poisson-transformed data with factors colony size (500/1000) and time as the repeated measure. In addition, in analysing the number of minimal path solutions, all observed combinations of minimal paths were scored. For example, where one minimal path linked the nest and the central bridge, and two minimal paths linked the central bridge and the food source, two minimal paths were scored, as this was the number of minimal paths available to any individual.

Behaviour of individual ants

We recorded video at 10 frames per second with a Logitech 9000 webcam placed 30 cm above the section of the maze that the ant was allowed to explore. We omitted the first minute of data – just after the focal ant had been introduced to the maze – from the analysis. The remaining 7 min of video were analysed using the open source machine-vision program ‘Ctrax’ (Branson et al., 2009). For each frame, the body axis and movement direction of the ant was used to calculate its orientation, and to generate 4200 data points for velocity, acceleration and turning speed for each ant. We then determined an average of each measure for each ant, and calculated a mean for each treatment group of 15 replicates. These means were compared via one-way ANOVA.

Towers of Hanoi

All ant colonies tested succeeded in solving the Towers of Hanoi in the sense that they all constructed a trail between the nest and the food source. More importantly, 93.3% of the colonies with definable trails in all treatments succeeded in finding a shortest path solution by the end of the first hour, except for the 500, no pre-exposure treatment, of which only 43% found a minimal path (Fig. 4A). When one of the minimal paths was then blocked, 90 to 93.3% of the colonies with exposure were able to find the remaining minimal solution, whereas 73.3 to 78.6% of the colonies without pre-exposure achieved this.

Fig. 5.

Results of open maze analysis. The dashed lines represent the maximum (corresponding with following the edge as shown in Fig. 3), average, and minimum (when ants construct a trail though the centre in the maze as illustrated in Fig. 3) possible edge values. Shown are means ± s.e.m. The results show a strong intrinsic edge following behaviour, with the mean edge values at all time points remaining near the upper bound.

Fig. 5.

Results of open maze analysis. The dashed lines represent the maximum (corresponding with following the edge as shown in Fig. 3), average, and minimum (when ants construct a trail though the centre in the maze as illustrated in Fig. 3) possible edge values. Shown are means ± s.e.m. The results show a strong intrinsic edge following behaviour, with the mean edge values at all time points remaining near the upper bound.

When exploration pheromone was present, definable trails were observed in the first 10 min (Fig. 4B). This was not the case when exploration pheromone was absent (Fig. 4B). Despite the network being longer during the first 10 to 20 min, the final network length after 60 min was shorter in the presence of exploration pheromone than in its absence (F(1,56)=13.718, P<0.001, η2p=0.197). Increasing colony size resulted in a longer total network length both in the presence and in the absence of exploration pheromone (F(1,56)=28.484, P<0.001, η2p=0.337). The larger η2p value for colony size indicates that colony size had a greater effect on total network length than exploration pheromone.

During the first hour, the mean number of minimal paths used (Fig. 4C) was significantly lower in colonies of 500 ants with no pre-exposure than in all other treatments. The presence of exploration pheromone had a stronger effect on the number of minimal paths (F(1,56)=29.060, P<0.001, η2p=0.342) than did colony size (F(56,1)=21.758, P<0.001, η2p=0.280). Pre-exposure had an effect on the quality of the solution constructed, with pre-exposed treatments showing far fewer corrective moves than the others (F(1,56)=28.965, P<0.001, η2p=0.341) (Fig. 4D), whereas there was no significant effect of colony size (F(1,56)=1.703, P=0.197, η2p=0.030).

After blocking the trail at the connecting bridge (after 60 min), total network length initially either increased slightly or remained stable. This is because the high level of traffic at this stage enabled the ants to build new, definable trails to reconnect the nest and food source over this 10 min period. Ten to 20 minutes after this dynamic change, we saw a gradual pruning of the network (Fig. 4B). Pre-exposure had a larger effect on total network length in the second hour (F(1,56)=26.674, P<0.001, η2p=0.323) than did colony size (F(1,56)=13.200, P=0.001, η2p=0.191). At the end of the second hour, all treatments reached a plateau between a mean of 0.6 and 1 minimal path (Fig. 4C), with no significant difference between the presence and absence of exploration pheromone (F(1,56)=1.524, P=0.222, η2p=0.027) or between colony sizes (F(1,56)=1.524, P=0.222, η2p=0.027).

All treatments showed an increase in corrective moves after the blockage, followed by a rapid decrease to very near to the number of corrective moves made before the change (Fig. 4D). Pre-exposure reduced the number of corrective moves required (F(1,56)=63.578, P<0.001, η2p=0.532), as did colony size (F(1,56)=8.219, P=0.006, η2p=0.006), although this latter effect was far weaker than the former. Interestingly, colonies of 500 ants without pre-exposure required less corrective moves than did unexposed colonies of 1000 ants. This effect might be related to the trails of the 1000 ant colonies being less focused on minimal paths directly before the obstruction.

Variations on Towers of Hanoi

Ants showed a strong intrinsic bias for forming trails along the outer edges of the open maze. Mean edge proximity score remained close to maximal at all time points (Fig. 5). The ants thus preferentially built trails near or along the outside edges of the open maze.

As in the original Towers of Hanoi experiment, all colonies tested succeeded in producing at least one trail between nest and food source in the Modified Towers of Hanoi. At the end of the first hour, 100% of the pre-exposed and 83.3% of the no exposure colonies used at least one shortest path (Fig. 6A). They were also able to respond to the blockage at the end of the first hour. Over the first 10 min after both connecting bridges were blocked, 25% of the pre-exposure colonies began to use a minimal path, whereas this was seen for 8.3% of the no pre-exposure colonies. By the end of the second hour, 100% of the pre-exposure colonies and 85.7% of the no pre-exposure colonies were using a minimal path.

Unlike in the original Towers of Hanoi experiment, in the initial 10 min of the Modified Towers of Hanoi experiments definable trails were observed for the no exposure treatment. These differences are probably accounted for by the differences in activity levels due to colony age and time of year at which the experiments were conducted. Otherwise the first hour of the Modified Towers of Hanoi experiments showed a similar relationship between the total length of the network and exposure (Fig. 6B). Total network length was lowest in the pre-exposed treatment (F(1,22)=7.363, P=0.013). We found no significant differences in the second hour between the pre-exposure and no pre-exposure treatments (F(1,22)=1.478, P=0.237). By the end of the experiment, both treatments converged to a point at or near the minimal network length.

Both the pre-exposure and no pre-exposure treatments showed the same pattern in terms of number of minimal paths used (Fig. 6C), with no significant difference between treatments in either the first (F(1,22)=3.449, P=0.077) or second (F(1,22)=1.066, P=0.313) hour. On average, more than one minimal path was used over the first hour, and directly following the dynamic change this number dropped to near zero owing to the blocking of these paths. Over the second hour, both treatments converged to use a mean of just one of the four possible minimal paths.

No significant difference was found between treatments in the mean number of corrective moves (Fig. 6D) for either the first (F(1,22)=1.219, P=0.281) or second (F(1,22)=0.945, P=0.342) hour. The number of corrective moves remained constant and low over the first hour, and then increased immediately following the dynamic change. The mean number of corrective moves returned to the original low levels by the end of the experiment.

Behaviour of individual ants

Ants introduced onto a maze containing exploration pheromone moved significantly faster (Fig. 7A) than did those on a clean maze (10.61 mm s–1 compared with 8.61 mm s–1, respectively; F(1,28)=6.899, P=0.014). Similarly, in the presence of exploration pheromone, acceleration (42.48 mm s–2 compared with 33.9 mm s–2; F(1,28)=9.829, P=0.004; Fig. 7B) and turning speed (1.57 rad s–1 compared with 1.37 rad s–1; F(1,28)=7.643, P=0.010; Fig. 7C) were significantly higher than for ants exploring a clean maze.

Fig. 6.

Results for Modified Towers of Hanoi. (A) Proportion of colonies with a minimal path. (B) Total network length. (C) Number of minimal paths used. (D) Number of corrective moves in the network. The vertical dashed line indicates the time of dynamic change. Shown in B–D are means ± s.e.m.

Fig. 6.

Results for Modified Towers of Hanoi. (A) Proportion of colonies with a minimal path. (B) Total network length. (C) Number of minimal paths used. (D) Number of corrective moves in the network. The vertical dashed line indicates the time of dynamic change. Shown in B–D are means ± s.e.m.

Our results show that Argentine ants can find the optimal solution to dynamic shortest path problems. Moreover, the ants had little problem finding an efficient alternative solution when their main trail was obstructed. The open maze experiments showed that the ants displayed a strong bias towards the outside edges. This bias is consistent with ants attempting to explore the full 360 deg radius of the environment. The intrinsic edge bias raised the question of whether the ants really recognise the shortest path, or whether the trail simply converges to the edge of the maze. The Modified Towers of Hanoi was designed to test this, by blocking both connecting bridges and opening a new path through the centre of the maze, thus forcing the ants to overcome their edge preference. Clearly, when the shortest path does not follow the edge, the ants are still able to construct the optimal solution.

Although increasing colony size generally increased the total length of the network, it was exploration pheromone that had the most striking effect on solution efficiency. In almost all cases, the presence of exploration pheromone led to smaller, less convoluted networks, the use of closer to just one minimal path, and fewer corrective moves. Increased solution efficiency can be at least partly attributed to the significantly greater velocity, acceleration and turning speed of individual ants when exploration pheromone was present. A similar effect has been found previously in Lasius niger (Beckers et al., 1992), but never before in L. humile. Strong behavioural evidence for such an exploration pheromone has been found for a similar mass-recruiting ant species, the Big-headed ant Pheidole megacephala (Fabricius). The two-pheromone model built on this experimental evidence suggests that a short-lived foraging pheromone, coupled with a longer-lasting exploration pheromone allows the ants to quickly redistribute their foragers to exploit new resources, even in the presence of a strong foraging trail to where a food source has just been removed (Dussutour et al., 2009b). Gerbier et al. (Gerbier et al., 2008) also found that when exploration pheromone was present ants tended to move more directly towards the food source when unfed, and towards the nest when fed, by continuing to follow their initial direction of travel. Hence, in the presence of exploration pheromone ants appear to use a more direct path both to and from the food source.

The initial response of the ants to blocking their trail was to construct a new trail along the zig-zagging path perpendicular to the long axis of the maze. This can be seen in Fig. 8B, which was taken 30 min after the blockage. This perpendicular path constitutes a large number of corrective moves. In the Towers of Hanoi configuration, the perpendicular path is the shortest path between the blocked area and the open connecting bridge, and hence efficiently redirects traffic. Ultimately, as traffic increases on the remaining minimal path at the expense of the longer, perpendicular path, evaporation of the trail pheromone leads to the extinction of the longer path (Fig. 8C).

Fig. 7.

Boxplot results of the individual ant locomotion experiments. (A) Velocity. (B) Acceleration. (C) Turning speed. The upper and lower sides of the box represent the interquartile range and the dark line is the median. The bars extending from the box indicate the highest and lowest values, excluding outliers, which are presented as black dots. The velocity, acceleration and turning speed of individual workers was significantly greater in the presence than in absence of exploration pheromone.

Fig. 7.

Boxplot results of the individual ant locomotion experiments. (A) Velocity. (B) Acceleration. (C) Turning speed. The upper and lower sides of the box represent the interquartile range and the dark line is the median. The bars extending from the box indicate the highest and lowest values, excluding outliers, which are presented as black dots. The velocity, acceleration and turning speed of individual workers was significantly greater in the presence than in absence of exploration pheromone.

The response to the blockage violates a hypothesis of ‘simple’ trail-following behaviour in mass-recruiting ant species. On encountering the blockage, an individual ant is forced to complete a 180 deg turn to continue travel. At this point, both the asymmetry of the bifurcation (Gerbier et al., 2008; Garnier et al., 2009) and the strong, established pheromone trail (Goss et al., 1989; Beckers et al., 1990; Traniello et al., 1995; Nicolis and Deneubourg, 1999; Camazine et al., 2001) might be thought to impel her to retrace her steps directly back along the edge she has just travelled. This was not what we observed. Instead, the ants searched their local area around the now blocked path for a way around the blockage, despite the violation of up to seven bifurcation symmetry ‘preferences’ before the open connecting bridge was reached. This indicates that information provided by the bifurcation asymmetry is ignored by the ants in favour of information from some internal compass (Gerbier et al., 2008) that might provide the ants with a sense of directional intention. For individual foragers to return on a direct path to the nest, despite a tortuous outbound path, would require a compass, as well as an intrinsic odometer (Wittlinger et al., 2006). This ability has yet to be shown in a mass-recruiting species.

Fig. 8.

Example results for Towers of Hanoi (500 ants, no exploration pheromone) at different time points: (A) 60 min, (B) 90 min, (C) 120 min. Note that the usage of the perpendicular path (and hence number of corrective moves) increases after the dynamic change, and then returns to initial levels.

Fig. 8.

Example results for Towers of Hanoi (500 ants, no exploration pheromone) at different time points: (A) 60 min, (B) 90 min, (C) 120 min. Note that the usage of the perpendicular path (and hence number of corrective moves) increases after the dynamic change, and then returns to initial levels.

Fig. 9.

Example results for Modified Towers of Hanoi (with exploration pheromone) at different time points: (A) 60 min, (B) 90 min, (C) 120 min. At 90 min, the perpendicular path and the new minimal path show similar traffic levels, yet by 120 min the ants solely use the minimal path.

Fig. 9.

Example results for Modified Towers of Hanoi (with exploration pheromone) at different time points: (A) 60 min, (B) 90 min, (C) 120 min. At 90 min, the perpendicular path and the new minimal path show similar traffic levels, yet by 120 min the ants solely use the minimal path.

The bifurcation asymmetry, suggested by Garnier et al. (Garnier et al., 2009), could play a role in the creation of efficient paths in the modified Towers of Hanoi experiment. Unlike the standard Towers of Hanoi, the modified version does not have an established trail leading away from the new central bridge. Rather trails build up along the perpendicular paths on both sides of the blockage (Fig. 9B). However, rather than being attracted to these new pheromone trails and thus along the path perpendicular to their goals, the ants produce a new minimal path (Fig. 9C). This solution is consistent with the ants ignoring pheromone in favour of continuing in an established direction (i.e. bifurcation asymmetry).

Our study gives us an opportunity to think again what computer scientists and biologists can learn from one another. In many ways the field of Ant Colony Optimisation has moved away from its biological roots, because when applied to real problems the simple pheromone following ants do not solve difficult combinatorial problems (Dorigo and Stützle, 2004). Instead, various pieces of heuristic information are added to the ACO solution scheme in order to better match the solution to the problem in hand. We would raise three important points with respect to better linking experiment and algorithm development. Firstly, we have found that ants are capable of solving dynamic optimisation problems and suggest that it is such problems, rather than static NP-hard problems, that should be the focus of algorithm development. Secondly, our and several previous studies have identified exploration pheromone, laid out in advance of starting the problem-solving exercise, as an important component in reaching good solutions. Finally, as in the addition of heuristic information to ACO schemes, our results suggest that Argentine ants might have more advanced navigational abilities than they are often credited with (Aron et al., 1993). Although in general the link between experiments on ants and ACO is still rather weak, these types of observations show that experiments can continue to provide new inspiration for novel optimisation algorithms and that algorithms can inspire our investigation of biology.

Our experiment was inspired by a suggestion by Dave Broomhead and this article is dedicated to him on the occasion of his 60th birthday. We thank members of the Behaviour and Genetics of Social Insects lab (Sydney) for constructive comments on an earlier version of this manuscript and the anonymous referees for helpful comments and criticisms.

     
  • ACO

    ant colony optimisation

  •  
  • NP-hard

    non-deterministic polynomial-time hard

  •  
  • η2p

    partial eta squared

This work was funded by the Human Frontier Science Program.

Akay
D.
,
Toksari
M. D.
(
2009
).
Ant colony optimization approach for classification of occupational low back disorder risks
.
Hum. Factor. Ergon. Man.
19
,
1
-
14
.
Alon
N.
,
Srinivasan
A.
(
1997
).
Improved parallel approximation of a class of integer programming problems
.
Algorithmica
17
,
449
-
462
.
Aron
S.
,
Beckers
R.
,
Deneubourg
J. L.
,
Pasteels
J. M.
(
1993
).
Memory and chemical communication in the orientation of two mass-recruiting ant species
.
Insectes Soc.
40
,
369
-
380
.
Banavar
J. R.
,
Maritan
A.
,
Rinaldo
A.
(
1999
).
Size and form in efficient transportation networks
.
Nature
399
,
130
-
132
.
Bebber
D. P.
,
Hynes
J.
,
Darrah
P. R.
,
Boddy
L.
,
Fricker
M. D.
(
2007
).
Biological solutions to transport network design
.
Proc. R. Soc. B
274
,
2307
-
2315
.
Beckers
R.
,
Deneubourg
J. L.
,
Goss
S.
,
Pasteels
J. M.
(
1990
).
Collective decision-making through food recruitment
.
Insectes Soc.
37
,
258
-
267
.
Beckers
R.
,
Deneubourg
J. L.
,
Goss
S.
(
1992
).
Trail laying behavior during food recruitment in the ant Lasius niger
.
Insectes Soc.
39
,
59
-
72
.
Beekman
M.
,
Sumpter
D. J. T.
,
Ratnieks
F. L. W.
(
2001
).
Phase transition between disordered and ordered foraging in Pharaoh’s ants
.
Proc. Natl. Acad. Sci. USA
98
,
9703
-
9706
.
Bell
J. E.
,
McMullen
P. R.
(
2004
).
Ant colony optimization techniques for the vehicle routing problem
.
Adv. Eng. Inform.
18
,
41
-
48
.
Bhatkar
A.
,
Whitcomb
W. H.
(
1970
).
Artificial diet for rearing various species of ants
.
Fla. Entomol.
53
,
229
-
232
.
Bonabeau
E.
,
Dorigo
M.
,
Theraulaz
G.
(
1999
).
Swarm Intelligence: From Natural to Artificial Systems
.
New York
:
Oxford University Press
.
Bonabeau
E.
,
Dorigo
M.
,
Theraulaz
G.
(
2000
).
Inspiration for optimization from social insect behaviour
.
Nature
406
,
39
-
42
.
Branson
K.
,
Robie
A. A.
,
Bender
J.
,
Perona
P.
,
Dickinson
M. H.
(
2009
).
High-throughput ethomics in large groups of Drosophila
.
Nat. Methods
6
,
U451
-
U477
.
Buhl
J.
,
Hicks
K.
,
Miller
E.
,
Persey
S.
,
Alinvi
O.
,
Sumpter
D.
(
2009
).
Shape and efficiency of wood ant foraging networks
.
Behav. Ecol. Sociobiol.
63
,
451
-
460
.
Buneman
P.
,
Levy
L.
(
1980
).
The Towers of Hanoi problem
.
Inform. Process. Lett.
10
,
243
-
244
.
Camazine
S.
,
Deneubourg
J.-L.
,
Franks
N. R.
,
Sneyd
J.
,
Theraulaz
G.
,
Bonabeau
E.
(
2001
).
Self-Organization in Biological Systems
.
538
pp.
NJ
:
Princeton University Press
.
Deneubourg
J. L.
,
Aron
S.
,
Goss
S.
,
Pasteels
J. M.
(
1990
).
The self-organizing exploratory pattern of the Argentine ant
.
J. Insect Behav.
3
,
159
-
168
.
Dorigo
M.
,
Stützle
T.
(
2004
).
Ant Colony Optimization
.
Cambridge, MA
:
MIT Press
.
Dorigo
M.
,
Maniezzo
V.
,
Colorni
A.
(
1996
).
Ant system: optimization by a colony of cooperating agents
.
IEEE T. Syst. Man. Cy. B
26
,
29
-
41
.
Dussutour
A.
,
Fourcassie
V.
,
Helbing
D.
,
Deneubourg
J. L.
(
2004
).
Optimal traffic organization in ants under crowded conditions
.
Nature
428
,
70
-
73
.
Dussutour
A.
,
Beekman
M.
,
Nicolis
S. C.
,
Meyer
B.
(
2009a
).
Noise improves collective decision-making by ants in dynamic environments
.
Proc. R. Soc. B
276
,
4353
-
4361
.
Dussutour
A.
,
Nicolis
S. C.
,
Shephard
G.
,
Beekman
M.
,
Sumpter
D. J. T.
(
2009b
).
The role of multiple pheromones in food recruitment by ants
.
J. Exp. Biol.
212
,
2337
-
2348
.
Garnier
S.
,
Guerecheau
A.
,
Combe
M.
,
Fourcassie
V.
,
Theraulaz
G.
(
2009
).
Path selection and foraging efficiency in Argentine ant transport networks
.
Behav. Ecol. Sociobiol.
63
,
1167
-
1179
.
Gerbier
G.
,
Garnier
S.
,
Rieu
C.
,
Theraulaz
G.
,
Fourcassie
V.
(
2008
).
Are ants sensitive to the geometry of tunnel bifurcation?
Anim. Cogn.
11
,
637
-
642
.
Gonzalez
T.
,
Sahni
S.
(
1976
).
Open shop scheduling to minimize finish time
.
J. ACM
23
,
665
-
679
.
Goss
S.
,
Aron
S.
,
Deneubourg
J. L.
,
Pasteels
J. M.
(
1989
).
Self-organized shortcuts in the Argentine ant
.
Naturwissenschaften
76
,
579
-
581
.
Hölldobler
B.
,
Wilson
E. O.
(
1990
).
The Ants
.
Berlin
:
Springer-Verlag
.
Korf
R.
,
Felner
A.
(
2007
).
Recent progress in heuristic search: a case study of the four-peg Towers of Hanoi problem
.
Int. Joint Conf. Art. Intell.
2324
-
2329
.
Latty
T.
,
Beekman
M.
(
2009
).
Food quality affects search strategy in the acellular slime mould, Physarum polycephalum
.
Behav. Ecol.
20
,
1160
-
1167
.
Liang
Y. C.
,
Smith
A. E.
(
2004
).
An ant colony optimization algorithm for the redundancy allocation problem (RAP)
.
IEEE Trans. Rel.
53
,
417
-
423
.
Nakagaki
T.
,
Yamada
H.
,
Toth
A.
(
2000
).
Maze-solving by an amoeboid organism
.
Nature
407
,
470
.
Nakagaki
T.
,
Yamada
H.
,
Toth
A.
(
2001
).
Path finding by tube morphogenesis in an amoeboid organism
.
Biophys. Chem.
92
,
47
-
52
.
Nakagaki
T.
,
Yamada
H.
,
Hara
M.
(
2004
).
Smart network solutions in an amoeboid organism
.
Biophys. Chem.
107
,
1
-
5
.
Nicolis
S. C.
,
Deneubourg
J. L.
(
1999
).
Emerging patterns and food recruitment in ants: an analytical study
.
J. Theor. Biol.
198
,
575
-
592
.
Perna
A.
,
Valverde
S.
,
Gautrais
J.
,
Jost
C.
,
Sole
R.
,
Kuntz
P.
,
Theraulaz
G.
(
2008
).
Topological efficiency in three-dimensional gallery networks of termite nests
.
Physica A
387
,
6235
-
6244
.
Romik
D.
(
2006
).
Shortest paths in the Tower of Hanoi graph and finite automata
.
SIAM Discret. Math.
20
,
610
-
622
.
Tero
A.
,
Takagi
S.
,
Saigusa
T.
,
Ito
K.
,
Bebber
D. P.
,
Fricker
M. D.
,
Yumiki
K.
,
Kobayashi
R.
,
Nakagaki
T.
(
2010
).
Rules for biologically inspired adaptive network design
.
Science
327
,
439
-
442
.
Traniello
J. F. A.
,
Robson
S. K.
,
Carde
R. T.
,
Bell
W. J.
(
1995
).
Trail and territorial communication in social insects
.
Chem. Ecol. Insects
2
,
241
-
286
.
Vassiliadis
V.
,
Dounias
G.
(
2009
).
Nature-inspired intelligence: a review of selected methods and applications
.
Int. J. Artif. Intell. Tools
18
,
487
-
516
.
Vittori
K.
,
Talbot
G.
,
Gautrais
J.
,
Fourcassie
V.
,
Araujo
A. F. R.
,
Theraulaz
G.
(
2006
).
Path efficiency of ant foraging trails in an artificial network
.
J. Theor. Biol.
239
,
507
-
515
.
West
G. B.
,
Brown
J. H.
,
Enquist
B. J.
(
1999
).
The fourth dimension of life: fractal geometry and allometric scaling of organisms
.
Science
284
,
1677
-
1679
.
Wittlinger
M.
,
Wehner
R.
,
Wolf
H.
(
2006
).
The ant odometer: stepping on stilts and stumps
.
Science
312
,
1965
-
1967
.
Zhang
J.
,
Chung
H. S. H.
,
Lo
A. W. L.
,
Huang
T.
(
2009
).
Extended ant colony optimization algorithm for power electronic circuit design
.
IEEE Trans. Power Electron.
24
,
147
-
162
.