Heuristics for the Weighted K-Chinese/Rural Postman Problem with a Hint of Fixed Costs with Applications to Urban Snow Removal Kaj Holmberg Lith-MAT-R--2015/13--SE

✬ ✩

Department of Mathematics

Heuristics for the Weighted k-Chinese/rural Postman Problem with a Hint of Fixed Costs with Applications to Urban Snow Removal Kaj Holmberg LiTH-MAT-R--2015/13--SE

✫ ✪ Department of Mathematics Linköping University S-58183 Linköping, Sweden. Heuristics for the weighted k-Chinese/rural postman problem with a hint of fixed costs with applications to urban snow removal

Kaj Holmberg

Department of Mathematics Linköping Institute of Technology SE-581 83 Linköping, Sweden

October 22, 2015

Abstract: We describe a weighted version of the k-Chinese or k-rural postman problem that occurs in the context of snow removal. The problem concerns the questions of which vehicle shall do each task and how the vehicles shall travel between tasks. We also consider diﬀerent numbers of vehicles, in view of a ﬁxed cost for each vehicle. We describe and discuss heuristic solution approaches, based on usable substructures, such as Chinese/rural postman problems, meta-heuristics, k-means clustering and local search improvements by moving cycles. The methods have been implemented and tested on real life examples.

1 Introduction

Snow removal can be an important and costly issue, in certain places at certain times. The yearly amounts of snow has recently varied much in several countries. Some years there is a lot, other years hardly anything. Hence the demand of snow removal varies a lot. Snow removal can be quite difficult and expensive, if there is very much snow. Winters without snow make it difficult to motivate keeping a large fleet of snow removal vehicles available. Complaints from the public about snow removal are common.

A conclusion is that one cannot expect to use the plans from last year. New situations require new plans, which require planning tools.

The situation (in Sweden) is as follows. A city is divided into areas, and for each area, a contractor is hired. As the contractors work independently, such an area gives natural borders for the optimization. We therefore consider one such area, and focus on the question how to allocate the streets to the vehicles available.

1 In the paper Holmberg (2014), we give a thorough and detailed description of the snow removal problem, which is very difficult. Here we deal with a simplified problem. The treatment of a street includes many details, such as several sweeps, depending on the width of the street, additional clearing at crossings, and even considerations of the fact that it might take a longer time turning a vehicle around than continuing straight ahead. In this paper, we omit all these details, and approximate the work needed to clear a street by one operation. Furthermore, we consider only normal roads, not pedestrian roads and cycling roads, which may require different vehicles.

We assume that the vehicles are identical, and the main task is to decide which links each vehicle should take care of. We don’t use depots, since we don’t know where the vehicles come from. The areas to be treated are often parts of larger areas, so sometimes the vehicles may come from others tasks, and when they are done, the will go to yet other tasks, and we don’t know where these tasks are located. Since there is no time information in our model, the vehicle can enter the tour at any node.

If there was only one vehicle to do all snow removal, the problem would be reduced to a Chinese postman problem, which simply is to ﬁnd a shortest round trip that covers all links in a graph.

Related problems are the following. The rural postman problem is to ﬁnd a shortest round trip that covers certain links. This corresponds to one vehicle clearing a subset of the links.

The k-Chinese Postman Problem is to ﬁnd k shortest round trips that covers all links, which corresponds to k vehicles clearing all links.

The k-Rural Postman Problem is to ﬁnd k shortest round trips that covers certain links, which corresponds to k vehicles clearing a subset of the links.

Our simpliﬁed snow removal planning problem turns out to be a version of the k-Chinese Postman Problem. More precisely, it is a weighted combination of the original k-Chinese Postman Problem, where the total distance is minimized, and the min-max k-Chinese Postman Problem, where the maximal distance traveled by a vehicle is minimized. The two objectives are here simply added together with certain weights. In practical terms, it is a compromise between minimizing the cost for all travel (and equipment) and minimizing the time until all snow is removed.

Furthermore it could well be case of k-rural Postman Problem, since there could be links that not need to be treated. There could be links where the snow already has been removed, or where the responsibility of snow removal lies elsewhere. Our methods will be able to solve such cases just as well.

We also wish to compare the result of using diﬀerent numbers of (identical) vehicles. Clearly, it is faster to let several vehicles work at the same time, but one has to pay for the vehicles. Each used vehicle has a ﬁxed cost (which does not depend on how much the vehicle is used, only on the fact that it is available).

A vehicle may also need to transport itself on a street, even when the street has been cleared, in order to get to an uncleared street. This can often be done faster than

2 clearing. We will assume that the ﬁrst time a link is used, it is cleared, and the following passes will only be transportation, which is quicker.

However, as motivated more later, this has limited influence on the optimization, since the actual work (snow removal) to be done is constant. If all links have to be cleared, the time needed for this is the total length of all links multiplied with the unit snow removal time. The total time for the whole operation will be this constant plus the time needed for additional transportation, so the optimization can only minimize the transportation time (when not clearing). Therefore the difference in speed will only occur between the snow removal term (which is constant) and the transportation term (which is minimized). Thus this difference does not affect which solution is best, as all feasible solutions have the same constant base. It is only when calculating the final end time this difference must be taken into account.

Summing up, the problem is the following. We have a number of streets in a city network that needs to be cleared from snow, and a limited number of identical vehicles that shall do the job. One main task is to allocate links to vehicles, i.e. decide which subarea each vehicle will take care of. In practice, a driver is allocated to an area, and is then often free to plan the tour by himself. Obviously it could be useful for the driver to have a suggested tour, but it is not certain that this tour will be used. Therefore, the area to be covered by each vehicle is a more important result than how the tours exactly look like. However, we need to ﬁnd the tours in order to evaluate an allocation.

There are several possible objectives for this problem. We could minimize the time when everything is ready, i.e. minimize the maximal time of a vehicle. This leads to the min-max k-Chinese postman problem. We could also minimize total cost for the clearing, i.e. minimize the sum of the times for the vehicles. This leads to the original k-Chinese postman problem.

Another possibility is to minimize the transportation distance, when the vehicles are not clearing. This is more or less the same as minimizing the pollution caused by the vehicles, so this objective is environmentally friendly. It is also in principle the same goal as minimizing the total distance traveled. The best objective function is, we believe, a weighted combination of those above.

We need data for this problem. The street network and all link characteristics can nowa- days be found in diﬀerent data bases. In this paper we use data from OpenStreetMap. In Holmberg (2015b), we give more details about what data there is and how to extract the data to a useful form.

Times for vehicles are known approximately, as an average speed. This motivates the approximation of combining small tasks into larger.

The operators will in the future be required to equip each vehicle with a GPS. This gives data from which we may extract times for all tasks. So task times will be available with better accuracy. Analyzing the GPS tracks will require map matching, which is described in Holmberg (2015a).

3 2 Survey

The Chinese postman problem is a well-known problem where each link in a graph needs to be covered. It can be solved optimally and polynomially by a well-known method, see for example Eiselt, Gendreau, and Laporte (1995a). The rural postman is a similar problem, where however only a subset of the links need to be covered. This problem is NP-hard, but there exist eﬃcient heuristics, Christoﬁdes, Campos, Corberán, and Mota (1986), Eiselt, Gendreau, and Laporte (1995b), Cook, Schoenfeld, and Wainwright (1998), Ghiani, Laganá, and Musmanno (2006), Holmberg (2010).

The k-Chinese postman problem is treated in Pearn (1994), Ahr and Reinelt (2002), Osterhues and Mariak (2005), Ahr and Reinelt (2006), and the k-rural postman in Benavent, Corberán, Plana, and Sanchis (2009). Other papers on related arc routing problems are Kramberger and Zerovnik (2007), Benavent and Soler (1999), Corberan, Marti, Martinez, and Soler (2002), and Corberán and Prins (2010).

Winter road maintenance is treated in Campbell and Langevin (1995), Labelle, Langevin, and Campbell (2002), Trudel (2005) and Fu, Trudel, and Kim (2009). Rural snow removal is discussed in Razmara (2004) and Islam (2010). A comprehensive survey of problems around winter road maintenance is found in papers Perrier, Langevin, and Campbell (2006a), Perrier, Langevin, and Campbell (2006b), Perrier, Langevin, and Campbell (2007a), and Perrier, Langevin, and Campbell (2007b). In Perrier, Langevin, and Amaya (2008) a model for urban snow removal is given.

3 Modeling

An example of a model for the min-max k-rural postman problem is given in Benavent et al. (2009). Here we don’t plan to use MIP-codes for solving the problem, so we will not go into all the details of such a model. However, some parts of the model will help in the description of the heuristic methods.

We will use the following indices: k ∈ Q: the set of vehicles, l ∈ L: the set of links (streets), and i ∈ N: the set of nodes (points, crossings). We will use q = |Q|, n = |N| and m = |L|. The links are assumed to be undirected, i.e. that the combined operation discussed previously can be done in either direction. Link j has length lj, which is obtained from the digital maps.

The set of links that shall be cleared of snow are denoted LC ⊆ L. Often LC = L, i.e. all links in the network need to be cleared, but sometimes there are links that will be cleared by another operator, or by other types of vehicles, but still are available to use for transport of our vehicles. As will be clear later, for our method it will not matter much if LC = L or not.

The value of q, the number of vehicles, may be ﬁxed, but we also consider the case when it can be varied. We will not make a model where this parameter is a variable, but will rather do a parametric analysis, i.e. solve the problem for diﬀerent values of q

4 and compare the results, as there are only a few interesting values.

We also need the following data.

tj : the time needed for a vehicle to clear link j.

dj : the time needed for driving a vehicle along link j (while not clearing).

A A We will consider the average speed to be v m/s, so it will take dj = lj/v seconds to drive the whole street j. The time to clean a road is estimated in average to be p times the length of the road divided by the speed. A value of p = 2 is reasonable for an ordinary street. This time includes two or three sweeps (middle, right side, left side) as well as cleaning of turning space and/or crossroads. All of this is seen as one operation that starts in one end of the street and ﬁnishes in the other end. For narrow streets or bicycle paths, p = 1 is more appropriate.

A It then takes tj = pdj = plj/v seconds to clean road j. We will assume that the cost for driving a vehicle is proportional to the time it is driven, so letting cD be the cost D D A for driving one second, the cost for cleaning link j is cj = c tj = c plj/v .

T OT The total time needed for clearing is t = Pj∈LC tj, which, as mentioned earlier, is a constant. As a matter of fact, the optimization will only concern the rest of the solution, which is traveling around with the vehicles without removing snow. That means that the optimization will only deal with the times d. Therefore, the diﬀerence between t and d (however big it is) is unimportant when it comes to which solution is best. However, in order to calculate correct times and costs of solutions, we need to take the diﬀerence into account.

We will also use the following, when considering the number of vehicles.

f : the ﬁxed cost for using a vehicle.

The main variables are the following.

C xjk = 1 if vehicle k treats link j, for all j ∈ L , k ∈ Q.

F zjk the number of times vehicle k is driven forwards on link j, for all j, k.

B zjk the number of times vehicle k is driven backwards on link j, for all j, k.

We may also use the following variables as an eﬃcient representation of an allocation.

κj : the vehicle (index) that is assigned to task j.

Obviously this means that xjκj = 1. We will also use Sk as the set of links allocated to vehicle k, i.e. Sk = {j : κj = k} = {j : xjk = 1}.

The main constraints state that each link must be treated by one vehicle.

C X xjk = 1 for all j ∈ L (1) k

5 The vehicles have to drive on the streets that are treated.

F B C zjk + zjk ≥ zjk for all j ∈ L , k ∈ Q (2)

We also need node equilibrium constraints, including zF and zB, but we don’t spell them out here. Slightly more complicated (and numerous) are the constraints that make the tours connected, but we won’t go into the details of this either. All in all we require that z forms connected round trips in the network, which is denoted by

z ∈ R. (3)

Note that constraints 3 must allow subtours, i.e. repeated visits to nodes.

Furthermore we need some variables and constraints for times.

vR : the time when all vehicles are ready.

U vk : the total time vehicle k is used.

C vk : the total time vehicle k is used for clearing snow.

L vk : the total time vehicle k is used for transportation.

These times are deﬁned as follows.

C vk = X tjxjk (4) j

L F B vk = X dj(zjk + zjk − xjk) (5) j

U C L vk = vk + vk (6)

R U v ≥ vk for all k (7)

Some possible objective functions are the following.

Minimize ﬁnal end time:

min vR

This leads to the min-max k-Chinese/rural postman problem.

Minimize sum of total times for vehicles:

U min X vk k

This leads to the original k-Chinese/rural postman problem.

Minimize sum of transport times:

6 L min X vk k

As there is a ﬁxed proportion between time and cost (cD), the last two objective functions above could just as well be regarded as minimizing the cost. They are often accompanied by a constraint on vR (the snow must be removed within a given time frame).

In this paper, we will use an objective function that is a weighted sum of the maximal time for a vehicle and the total cost for the whole operation.

A R C D U min w v + w (c X vk + fq) k

where wA and wC are weights reﬂecting the importance of time relative to cost. If wA is large, it is important to ﬁnish early, i.e. to minimize the total time until all links are cleared, while if wC is large, it is more important to minimize the total costs for the whole operation.

When we solve the problem, fq will be constant, but we may solve the problem for a couple of diﬀerent values of q, and compare the results, and then this term will be important.

If we wish to include the number of vehicles explicitly as a variable in the model, we would need to use the variables

yk = 1 if vehicle k is used, 0 if not,

and the constraints

xjk ≤ yk for all k

X yk ≤ q k

0 ≤ q ≤ |Q|, integer.

This is possible, but in practice we are interested in only a few values of q, so it would not be worth while.

4 Lower bounds

We will later describe heuristics for ﬁnding feasible solutions. This yields upper bounds on the optimal objective function value. It can then be quite useful to obtain some lower bounds on the optimal objective function value, by solving a relaxation of the problem.

We can then compare the upper and lower bounds, in order to estimate the possible potential for improvement. If we are very lucky, the lower bound may be close enough

7 to the upper bound, so that we may conclude that the solution is close enough to the optimum.

A possible use for a relaxation is to construct a branch-and-bound method around it. By branching, the lower bound will get closer to the upper bound, and ﬁnally indicate optimality. A possible binary branching is to ﬁx a certain xjk to 1 in one branch and to 0 in the other. Another possibility is to create q branches for a certain link j, where xjk = 1 in branch k.

The simplest lower bound is obtained by dividing the total work that needs to be done by the number vehicles. Here all transportation is omitted and the need for forming connected cycles is ignored.

tT OT l0 = wA + wC (cDtT OT + fq) q

A slightly better lower bound is to ﬁnd the cost of one Chinese postman tour, vP , i.e. the cost for one vehicle to do all, and divide by the number of vehicles.

1 t 1 l1 = wA + wC (cDt + fq) q

E 1 v where t = + tT OT , since the extra work, vE = vP − tT OT , is only transportation. p

The time for obtaining l0 is negligible, while the time needed to obtain l1 is somewhat longer, but still much shorter than what is needed for solving the problem, i.e. ﬁnding a good feasible solution.

5 Solution methods

We will focus on heuristic solution methods, such as constructive heuristics, local search and metaheuristics. We propose heuristic solution approaches with the following parts:

1. Initial allocation, dividing the tasks among the vehicles.

2. Find a feasible solution, i.e. a route for each vehicle.

3. Improve the routes by local search.

4. Improve the routes by metaheuristics.

We recall that routing a vehicle on a given set of undirected tasks is a rural postman problem, and eﬃcient heuristics for this problem is presented in Holmberg (2010).

A practical aspect omitted in the problem formulation is that it may take a longer time to turn a vehicle around than to continue around the block. However, we can do a slight eﬀort to avoid turning as follows.

8 When solving a rural postman problem, additional arcs/edges are added, in order to form a graph where an Euler tour exists. The procedure ends by ﬁnding an Euler tour in the extended graph. In this stage, the Euler tour is seldom unique. All Euler tours are equally good, so we may choose any one of them. Then it is possible to avoid a U-turn as much as possible, by simply avoiding going back to the recently visited node if possible.

A solution is thus evaluated as follows. The time for each vehicle tour is calculated, vU . (This also gives the cost as cDvU .) The total solution time is then equal to the maximal time over the vehicles. The total cost is the sum of the costs for the vehicle tours, plus the ﬁxed costs. The total objective function value is then calculated as a weighted combination of the maximal time and the total cost.

5.1 Local search and metaheuristics

A local search method may be used to try to improve the allocation of vehicles to tasks. The neighborhood is deﬁned by the following changes.

• Move a task from one vehicle to another.

• Swap tasks between two vehicles.

This can also be used in metaheuristics, such as simulated annealing or tabu search. Even though immediate improvement is not obtained, better solutions might be obtained after a number of changes.

The changes can easily be made by explicitly working with κj. A move means changing one value, from κj = k1 to κj = k2. A swap means changing two, but can also be seen as changing places of two elements in the vector. (This is interesting only if these two elements do not have the same value.)

When ﬁnding the routes for the vehicles, we solve a rural postman problem for each vehicle, and one question is the order in which the vehicles are handled. This is unin- teresting, if the vehicles do not aﬀect each other (in this phase). However, one possible heuristic improvement is the following.

When the tour of the ﬁrst vehicle is found, one usually ﬁnds that the vehicle passes L links that are not allocated to be cleared by that vehicle. (This happens when vk > 0, F B for links where zjk + zjk > 0 and xjk = 0.)

Then it is possible to “disobey” the allocation, and let that vehicle clear all links it passes. These links are then done, and later vehicles need not pass them. It is not certain that this is a good idea, since it could mean that the first vehicle does a lot of work while the following do less and less. (Here the difference between c and d matters. It is probably a better idea if this difference is small.) It could however be worth to try.

Consider the example in ﬁgure 1. First we give the allocation for a red and a blue vehicle and then the cycle of the red vehicle. Then we change the allocation of the blue

9 Figure 1: Changing allocation. links that the red vehicle passes to red. After this, the blue cycle only need to cover the middle link. Assuming length 1 of each link, the red cycle has length 6, and the blue cycle length 2, which sums up to 8. In the last picture, we show how long the blue tour would be if the allocation was not changed. Then the blue cycle would get length 6, so the sum of lengths would be 12.

If this is done, the order of the vehicles obviously matters. Therefore we consider changing the order of the vehicles. We let oi be the vehicle that is treated as number i in the order. Changes that can be done are swapping two elements of o, or shifting o one or more steps to the right or the left.

Some important choices are the following.

1. How to generate the initial allocation.

2. How to change the order and allocation between the iterations.

3. What kind of improvements of the solutions to do.

5.2 Manual interaction

Since the problem addressed in this paper previously was solved completely manually, it is not far-fetched to consider manual interaction with our proposed methods. Our goal here is to ﬁnd out how much this might help.

Since the street networks are not exceedingly large, it is not impossible to do the initial allocation by hand. We will later describe how one can, in our tool Vineopt, select a set of links and designate that set as allocation for a vehicle. This needs to be done q times, and for rather small areas, the number of vehicles is rather limited.

Another possibility is to modify an existing allocation. If the links in the areas for different vehicles are shown in different colors on a map, it is easier for a human to find promising changes of the allocation. The idea is to move links from vehicles that do much to vehicles that do less.

Informal testing in a project in an advanced course at Linköping university gave the result that manual improvements gave signiﬁcant improvements, see section 6.3.

10 5.3 k-means clustering

Given the coordinates of the nodes N, the k-means clustering problem means ﬁnding k centers and allocating each node to a center, so that the sum of distances from each node to its allocated center is minimal. The problem is NP-hard, but eﬃcient heuristics exist, see for example Lloyd (1982).

This can be used for our problem with k = q, i.e. to create one cluster for each vehicle. This gives a partitioning of the nodes to the vehicles. Since we are more interested in the links, we must do the following. If both the starting and ending node of a link belongs to the same cluster, the link is allocated to that cluster. However, if the starting node and the ending node of a link belongs to diﬀerent clusters, we must decide to which cluster of the two the link should belong. In our computational test, we have tried some simple heuristic ways, namely to allocate the link to the cluster that has the highest index of these two, or to the lowest index, or randomly to one of the two clusters.

We have also tried the following. If link j has its two end nodes in diﬀerent clusters, sj ∈ A and ej ∈ B, we calculate the “cost” of moving each node to the other cluster, as the distance from each cluster center to the other endpoint of the link. If this distance is large, the vehicle has to travel far from its center to include the link. Obviously this is no exact measure at all, but it may give some indication. We then choose the allocation that minimizes this measure.

Let d(A, sj ) denote the distance between center of cluster A and node sj, and d(B, ej) the distance between center of cluster B and node ej. If node ej is moved to cluster A, we get the distance d(A, sj)+lj, and if node sj is moved to cluster B, we get the distance d(B, ej )+ lj. Obviously, if d(A, sj) < d(B, ej ), then the “cost” is least if we choose the ﬁrst alternative. Thus, if d(A, sj ) < d(B, ej ), then we move node ej to cluster A, and if d(A, sj ) > d(B, ej ), then we move node sj to cluster B.

In the rest of the paper, we will use “vehicle” and “cluster” interchangeably, both indexed with k.

One might think that is will never be better to let the different vehicles go into each others clusters. However, in figure 2 we give a red and a blue cycle, with lengths on the paths. Swapping the red and blue vehicles on the shorter paths as shown to the right in the figure, makes the longest cycle, the red one, shorter, so it improves the solution. Thus we see that it may be beneficial to let vehicles cross paths.

5.4 Details of the method

At the heart of the method we have the code for solving a rural postman problem, Holmberg (2010). The set of required links for vehicle k is simply those links that have κj = k.

Two important vectors are o, ord, which is the order in which the vehicles are considered, and κ, vehdo, which states which vehicle should clean each link.

11 7 7 10 10 4 4 3 3

Figure 2: Changing allocation.

The method contains the following parts.

1. Read the network.

2. Construct a starting order.

3. Construct an initial allocation of links to vehicles.

4. Solve the rural postman problem for each vehicle.

5. Optional: Make a vehicle clean all uncleaned links passed by it, even if they were allocated to another vehicle.

6. Evaluate the objective function.

7. Optional: Improve the solution by moving cycles from one tour to another one.

8. Optional: Change the allocation to the current solution.

9. Change the order for the vehicles.

10. Change the allocation of links to vehicles.

11. Repeat steps 4 - 10.

12. Save the allocation of links to vehicles.

The total cost, which is the sum of the costs for each vehicle’s tour plus the ﬁxed costs for vehicle, is computed at each iteration. Moreover the time for the solution is computed, which is the maximal over the vehicles’ times. The objective function value is then computed as a weighted sum of the cost and the time.

Simulated annealing is used for moving to a new allocation (for κ), i.e. better points are always accepted and worse points are accepted with a certain probability.

Steps 3, 9 and 10 can be done in various ways, as described later. We also have to choose which of the optional steps 5, 7 and 8 to do. Below we give a list of the user- controlled parameters, with abbreviations, and possible values. (SA stands for simulated annealing.)

12 Parameter Abbr Description maxiter MI Number of main iterations markdone MD Let all passed links be cleared by the vehicle, (0,1) dorevsol RS Change the allocation following the solution, (0,1) startvar SV Variant for starting allocation, (0-16) movevar MV Variant for changing allocation, (0-7) moalloc MA Number of steps for vehdo to be moved chordvar ChV Variant for changing order, (0-5) moord MO Number of steps for ord to be moved impvar IV Variant for the improvement heuristic, (0-3) cymvar CyV Variant for moving cycles, (0-4) cymimp CI Move cycles only if if gives improvement, (0,1) tcymvar TC Variant for two cycle move, (0,1) fixcore FC Use fixed core, (0,1) fixcoredelta FD Fixed core proportion, δF C unfixext UF Use vehicle area reduction, (0,1) ufixdelta UD Vehicle area reduction proportion, δAR starttemp ST Starting temperature in SA intiter II Number of iterations with fixed temperature in SA rcold RC Cooling factor in SA

In the following sections, we will describe these parameters in more detail.

5.4.1 Mark links done

The method first finds a tour for the first vehicle, which passes all links that are allocated for the vehicle. If markdone=1, all links passed by the vehicle are marked as done, i.e. these links do not need to be treated by other vehicles. This may mean that the next vehicle has less work to do. If we change the order in which tours for vehicles are decided, a different solution could be generated.

Therefore the vector ord can be changed to get diﬀerent solutions. There are several diﬀerent approaches to (randomly) change ord, see below.

Step 8 in the algorithm is only interesting if step 5 is used. It decides if the changed allocation should be saved for later iterations, or if the old allocation (which was not followed) should be kept.

5.4.2 Starting allocation

It can be very important which starting allocation is used. One possible use of the method is to use only the starting allocation, without subsequent tries to improve it.

We work with the vector κ, vehdo. It is controlled by the parameter startvar.

13 startvar Action 0 Start from the current allocation (assuming there is one). 1 Create a random allocation. 2 Read a link set from Vineopt and use as a starting allocation. 3 Let vehicle 1 do everything. 4 Sort in x for starting node. 5 Sort in x for ending node. 6 Sort in x for middle point on links. 7 Sort in y for starting node. 8 Sort in y for ending node. 9 Sort in y for middle point on links. 10 Sort by distance from the origin to starting node. 11 Sort by distance from the origin to ending node. 12 Sort by distance from the origin to middle point on links. 13 Use k-means, with boundary links allocated to the highest number. 14 Use k-means, with boundary links allocated to the lowest number. 15 Use k-means, with boundary links randomly allocated. 16 Use k-means, with boundary links allocated to minimize distance to center.

Using startvar = 0, we can repeatedly try to improve the current solution, while re- setting the simulated annealing temperature between the runs. The case startvar = 2 enables starting from a previously constructed allocation. This can be used iteratively, to continue with a previous allocation that has been improved by other means (in Vi- neopt). For startvar = 4 - 12, the links are sorted in various ways, and the allocation is made for one vehicle at a time, picking the closest links in order. It can be seen as a kind of (crude) geographical clustering. For startvar = 13 - 16, the nodes are clustered in q groups (with the code kmeans), thereby inducing an allocation of the links for the vehicles. The diﬀerence is how boundary links, i.e. links with the end nodes in diﬀerent clusters, are treated.

An interesting question is if one should try to ﬁnd the best possible starting solution, or if a worse starting solution may enable the metaheuristics to ﬁnd even better solutions.

5.4.3 Change allocation

The changes of the allocation of links to vehicles is done in the vector κ, controlled by the parameters movevar and moalloc. The changes are (with one exception) random, intended for a simulated annealing framework, and not aiming towards local improvement.

14 movevar Action 0 No change. 1 Swap two random links. 2 Change moalloc random links. 3 Move a link from a vehicle that does most to a vehicle that does least. 4 Move a link between two random vehicles that have both nodes in common. 5 Swap a link between two random vehicles that have both nodes in common. 6 Move a link, when one node is in common. 7 Swap a link, when one node is in common.

For movevar = 3, the number of links allocated to each vehicle is calculated, and a link is moved from the vehicle that has the highest number of links to the vehicle that has the least number of links. For movevar = 4 - 7, two random vehicles are chosen, and a change is done if the corresponding tours have one (6, 7) or two (4, 5) nodes in common. The change is either to move a link (4, 6) or to swap the two links (5, 7). In the case of two nodes in common, the previous tours are still feasible, but maybe one of them can be shortened, as the vehicle is not required to clear the link.

5.4.4 Change order of vehicles

To change the order of vehicles, we may swap elements in the vector o (ord), or shift it a number of steps to the left or to the right. In detail, we have the following possibilities, controlled by parameters chordvar and moord.

chordvar Action 0 No change of ord. 1 Swap two random elements in the vector ord. 2 Create a completely new order. 3 Shift the vector ord moord steps to the right. 4 Shift the vector ord moord steps to the left. 5 Swap the most and least used vehicles.

For chordvar = 5, we sum up the number of links allocated to each vehicle, and switch the positions of the vehicle that has the most links allocated to it and the vehicles that has the least. If we are using markdone = 1, this often means swapping places between the ﬁrst and the last vehicle, and will make a big diﬀerence for the next tours.

5.4.5 Improve solutions

In figure 3, we first give two tours, one short red (dashed) and one longer blue. Both vehicles pass the middle node, so one may move the cycle shown in green (dotted) in the middle figure from blue to red, without the need for any vehicle to move to any other node. The total distance is unchanged by this, but the longer cycle becomes shorter, and the shorter cycle becomes longer, so the maximal time is decreased.

Here we try to improve the solution by ﬁnding a cycle in a tour, checking if the cycle contains a node in another tour, and then moving the cycle to that tour. Note that this

15 Figure 3: Moving cycles. does not change the total cost for the solution, as precisely the same subtours are used. However, moving a subtour from one vehicle to another may decrease the ﬁnal time. It is controlled by the parameter impvar.

impvar Action 0 Do not move cycles. 1 Move cycles, one random try. 2 Move cycles, terminate when there is no improvement. 3 Move cycles, try all combinations.

Clearly impvar = 1 is very quick, but usually fails, while impvar = 3 takes quite a long time.

5.4.6 Cycle moving

This procedure is invoked from Improve solution, and performs the actual movement of the cycle. It is controlled by the parameter cymvar.

cymvar Action 0 No change. 1 Try to move a random cycle between two random vehicles. 2 Try to move a cycle from a long tour to a short one. 3 Try to move a cycle between two given vehicles. 4 Try to move a cycle from an expensive tour to a cheap one.

The choice cymvar = 3 is only used for impvar = 3. If cymvar = 4 succeeds, it usually gives improvement.

When two cycles are chosen, one try is done if the parameter tcymvar = 0. If tcymvar = 1, all possible moves are tried.

5.4.7 Fixed core

As the evaluation of all possible changes in a local search or metaheuristic may take a signiﬁcant time, it is interesting to try to reduce this time by reducing the number of changes to consider. One may argue that a link very close to the center of a cluster

16 should belong to this cluster, and not be moved to another. Therefore we may use a modiﬁcation called “ﬁxed core”, which simply means that links that lie within a certain distance of the center of the cluster are not moved, i.e. the metaheuristic skips these changes.

F C First let djk = max(d(Ck,sj), d(Ck, ej)), where Ck is the center of cluster k and d(Ck, i) F C is the distance between Ck and node i. We use djk as the distance between link j and FCmin F C the center of cluster k. Now we find the minimal of those, dj = mink djk , over the 1 FCmin2 F C clusters, yielding k as argument, and the second minimal, dj = mink:k6=k1 djk , FCmin FCmin2 and calculate the quotient between them, δj = dj /dj . This is compared to F C F C a given value, δ (called “fixcoredelta”). All links j with δj < δ are considered to be the fixed core, and we do not consider changing them. In other words, a link is in the fixed core if the closest cluster center is much closer then the second closest, i.e. if FCmin F C FCmin2 F C dj < δ dj . Setting δ = 0.5 would mean that the distance to the second closest should be more than twice the distance to the closest.

This can be used only when the starting allocation is based on clustering, and the calculations of the distances are only done once, directly after the starting allocation has been decided.

5.4.8 Vehicle area reduction

Another possibility, called “vehicle area reduction”, aims at decreasing the size of the rural postman problems that need to be solved. Recall that a rural postman problem is solved for each vehicle, with required edges those that are allocated to the vehicle. Usually this is done in the whole network. However, if the allocated links lie in a very restricted part of the whole network, it is a fair guess that the vehicle should not wander oﬀ into other distant parts of the network. Therefore we omit all links that lie further away from the allocated links than a certain limit. This way the rural postman problems become smaller, and are quicker to solve.

Considering cluster k, we ﬁnd the minimal distance between links in that cluster, Sk, and each link j not in that cluster.

AR djk = min (min(d(sj1 ,sj), d(ej1 ,sj), d(sj1 , ej), d(ej1 , ej))) for all j 6∈ Sk. j1∈Sk

This is compared to the maximal distance within the cluster, where we also take single links into account.

ARmax dk = max (max(d(sj1 ,sj2 ), d(ej1 ,sj2 ), d(sj1 , ej2 ), d(ej1 , ej2 ), lj1 , lj2 )). j1∈Sk,j2∈Sk

AR AR AR If djk > δ dk , link j is considered to lie far away from cluster k, and link j is not included in the rural postman problem for vehicle k. Values for δAR could be around 2, which means that the minimal distance between the link and a link in the set is twice as large as the maximal distance between two nodes within the set.

These calculations are based on the current clusters, and need to be redone when the

17 Figure 4: Node selection in Vineopt allocation is changed. It is done just before the rural postman problem is solved. AR AR AR Unfortunately it is somewhat time consuming, and if no links with djk > δ dk are found, it is a waste of time.

6 Computational details

6.1 Vineopt

We use the graphic network optimization tool Vineopt, written by Kaj Holmberg in Tcl/Tk, see www.vineopt.com. It has the ability to show networks and allow changes in them, as well as solve various optimization problems in graphs, in this case most notably the Chinese postman problem (optimally) and the rural postman problem (heuristi- cally). (The optimization codes are written in C.)

Vineopt can visualize solutions and link sets (in diﬀerent colors), over a map background obtained from OpenStreetMap. It also allows manual creation and changes of the link sets. Here the link sets represent the allocations.

We can use the codes within Vineopt for solving the Chinese postman problem and rural postman problem. The allocation of links to vehicles can in Vineopt be done by hand and/or with the help of heuristics.

18 6.2 Input data

The street networks have been downloaded from OpenStreetMap and converted to Vi- neopt format (as rural postman problems, undirected graphs with arc costs and possibility to deﬁne “required” links). The arc costs are equal to the length of road segments measured in meters. See Holmberg (2015b) for details.

In our computational tests, we consider the average speed to be vA = 7.2 km/h, which is 2 m/s. The time to clean a road is estimated in average to be twice the length of the road divided by the speed, i.e. we use p = 2. With these assumptions, the operation takes exactly one second per meter, i.e. it takes lj seconds to clean street j if it has length lj. We also give one second the cost of one, so times will be equal to costs. We assume that the ﬁxed cost is 200 for each vehicle. Therefore the total cost for a solution is the sum of the costs for each vehicle plus number of vehicles times 200. The total objective function value is here the sum of the maximal time and the total cost.

6.3 Preliminary student work

The basis for this work has been used in a project in an advanced course in Opti- mization at Linköping University. The students used a preliminary code, written in Matlab/Octave, called snowplan.

They used a small network, colonia, (student accommodations) with 24 nodes and 36 links, two slightly larger networks, skanninge-n (Skänninge norra), with 62 nodes and 80 links and atvid-s (Åtvidaberg södra), with 68 nodes and 95 links, as well as two larger, studentryd (part of Ryd where students live), with 387 nodes and 519 links and vadstena (a small town with a large castle), with 581 node and 778 links.

The following was done by the students. The small problem colonia was ﬁrst solved for one vehicle. Then a solution for two vehicles was found by doing link selection/allocation manually in Vineopt. The same was done for three vehicles. Then the best of these solutions was found, ﬁrst under the assumption that time and cost are of equal importance, i.e. wA = wC , and then under the assumption that time is twice as important as cost, i.e. wA = 2wC .

Then the same was done with snowplan, without manual intervention.

After this, the same was done for the slightly larger problem skanninge-n (Skänninge norra).

For the problem atvid-s (Åtvidaberg södra), solutions for 1, 2, 3 and 4 vehicles were found using snowplan. The best number of vehicles was then found for the two weight scenarios above.

Finally the goal was to ﬁnd a very good solution for 6 vehicles for the largest problem vadstena, with equal weights for cost and time. First snowplan was run, and the times for each vehicle noted, especially which vehicle took the longest time. Then this solution was improved in Vineopt, by moving arcs from one set to another. After this,

19 the allocation was imported into snowplan, and the program was run again.

The students were encouraged to do this a couple of times, by the announcement of a competition, to see which group produced the best solution. The results were quite varied, and some students did more manual work than others. Unfortunately many groups let snowplan do too few iterations, due to a limited practical experience of metaheuristics. Therefore the manual improvements were signiﬁcant.

6.4 New implementation

A new implementation has been done in Python, using Numpy and Scipy (providing the code kmeans). The tests were run on an Acer Aspire X3 X3995 3.4GHz, running Linux, Fedora Core 22. The machine has four CPUs, but only one was used in the test runs. For visual aid, Tkinter and networkx were used. For solving the rural postman problems, we use a heuristic method described in Holmberg (2010), implemented in C.

7 Computational tests and evaluation

7.1 Test problems

The methods have been tested on several diﬀerent problem instances, based on Open- StreetMap-data, obtained from the Internet. See Holmberg (2015b) for details of the data extraction procedure. The networks cover smaller cities around Linköping, and also diﬀerent parts of Linköping.

A certain ﬁlter was applied to the OSM-data, leaving only streets (highways) with the following label: “motorway”, “trunk”, “primary”, “secondary” and “tertiary”, but leaving out “road”, “cycleway”, “living_street”, “pedestrian”, “footway”, “path” and “residental”. These instances are thus based on real life networks, with real life distances, but with certain links removed, and nodes with degree two eliminated by simply adding the two adjacent links into one, summing up the costs. Therefore the curvature of the link is not visible, but the distance is correct.

We have, for practical reasons, divided the set into three parts, depending on the size. The ﬁrst part contains networks with less than 200 nodes. They are rather small cities or areas, and give rather easy problems. The second part contains networks with between 200 and 400 nodes. These problems are a bit more time consuming to solve. The third part contains the most diﬃcult problems, networks with more than 400 nodes.

In table 1, we give the number of nodes and links, the total distance of all the links (i.e. the total distance to be cleared from snow), the average node degree, and the number of inhabitants for the cities represented in the test set. For the examples based on parts of Linköping, no number of inhabitants can be speciﬁed. (The whole of Linköping has around 140 000 inhabitants.)

20 After this, we give pictures of some of the networks.

Name Nodes Links Total dist Deg What? Inhab atvid-s 68 95 2028 1.40 Åtvidaberg, south - brokind 179 199 2664 1.11 Brokind 502 colonia 24 36 131 1.50 Colonia, student homes - ekangen 149 184 2629 1.23 Ekängen 2037 mantorp-c 122 145 2243 1.19 Mantorp, central - mantorp 185 226 4152 1.22 Mantorp 3671 rimforsa-o 94 109 1474 1.16 Rimforsa, east - skanninge-c 122 152 2292 1.25 Skänninge, central - skanninge-n 62 80 974 1.29 Skänninge, north - skanninge 127 158 2322 1.24 Skänninge 3140 studentryd-e 76 104 543 1.37 Ryd E - sturefors 183 212 3283 1.16 Sturefors 2229 valla-liu 160 213 1183 1.33 Valla, university part - vikingstad 178 213 3390 1.20 Vikingstad 2096 atvid 220 290 5940 1.32 Åtvidaberg 6859 borensberg 211 257 4994 1.22 Borensberg 2886 kisa 311 359 4541 1.15 Kisa 3687 liu 291 394 2713 1.35 Linköping University - rimforsa 233 262 2837 1.12 Rimforsa 2238 ryd1 382 508 2200 1.33 Ryd - studentryd-int 235 316 1494 1.34 Ryd, student part - studentryd 387 519 2920 1.34 Ryd, student part - linghem 836 1129 6624 1.35 Linghem 518 ljungsbro 812 1024 14013 1.26 Ljungsbro 6620 malmslatt 489 626 6542 1.28 Malmslätt 5214 mjolby 419 524 11894 1.25 Mjölby 12245 ryd-m 1107 1495 10910 1.35 Ryd - ryd 952 1291 7903 1.36 Ryd - ryd2 987 1337 9174 1.35 Ryd - soderkoping 666 910 10998 1.37 Söderköping 6992 vadstena 581 778 7785 1.34 Vadstena 5613 valla 527 707 5583 1.34 Valla -

Table 1: City instances, sorted in three groups.

21 VINEOPT

Figure 5: Åtvidaberg (atvid), 224 nodes, 294 links.

VINEOPT

Figure 6: Brokind, 179 nodes, 199 links.

22 VINEOPT

Figure 7: Kisa, 311 nodes, 359 links.

VINEOPT

Figure 8: Linghem, 836 nodes, 1129 links.

23 VINEOPT

Figure 9: Ljungsbro, 812 nodes, 1024 links.

VINEOPT

Figure 10: Malmslätt (malmslatt), 489 nodes, 626 links.

24 VINEOPT

Figure 11: Mantorp central (mantorp-c), 122 nodes, 145 links.

VINEOPT

Figure 12: Skänninge (skanninge), 127 nodes, 158 links.

25 VINEOPT

Figure 13: Söderköping (soderkoping), 666 nodes, 910 links.

VINEOPT

Figure 14: Vadstena, 581 nodes, 778 links.

26 VINEOPT

Figure 15: Valla, 527 nodes, 707 links.

VINEOPT

Figure 16: Vikingstad, 178 nodes, 213 links.

27 Figure 17: Mantorp, G30, q = 4.

7.2 Preliminary tests

First we made some preliminary tests, in order to ﬁnd promising values of the parameters.

In figures 17 and 18, we give pictures of allocations, where different vehicles are represented by different colors. One may see parts of the allocations that do not seem to be optimal. In both these examples, the yellow vehicle seems to be driving around too much. However, one cannot be certain of this. If, for example, the yellow vehicle is ready much before the others, it makes sense for it to help with parts that might delay the others. (The pictures are only examples, and not certain to be good solutions.)

In the ﬁgure 19, we show a preliminary run with Snowplan. In the graph at the bottom of the picture, the blue line denotes the actual objective function value for each iteration, the green cross denotes the objective function value after the improvement procedure, and the dashed red line denotes the best objective function value obtained up to that iteration. The light blue horizontal line denotes the lower bound.

In the top of the picture, one can see the parameters that are possible to change, and their current values. At the top left, the total time for each vehicle is shown. There one can identify the vehicles that takes the longest and the shortest time, and may consider moving links between them. In this example, one may consider moving links from vehicles 1 and 2 to vehicles 3 and/or 4.

28 Figure 18: Åtvidaberg, A30, q = 4.

In the middle green part, the best objective function value is shown, together with some more information, such as at which iteration the best value was found, and the lower bound.

In the ﬁgures 20 - 22, we show some more runs, mostly to show the curves of objective function values. One can see that there is no monotone improvement. However, this depends on which starting solution was used.

In the ﬁgure 23, we show how the network and allocation are brought into Vineopt. Figure 24 shows zooming in, which allows single links to be changed in Vineopt. The resulting allocation can then be saved and used as starting solution in Snowplan.

In tables 31 - 32 in the Appendix, we present preliminary tests, made in order to ﬁnd interesting values of the parameters. The best results are shown ordered in table 33, and a group of faster settings in table 34.

These preliminary computational tests yield several interesting settings. One should remember that the method contains randomness, so all results should be regarded as approximate. In any case, in these tests, the best setting seems to be the following: markdone 1, dorevsol 1, startvar 13, movevar 3, impvar 3, chordvar 5, cymvar 4, cymimp 0, tcymvar 0, ﬁxcore 0, intiter 7, starttemp 1000, rcold 0.3. In words: Mark passed streets done, revise allocation after this, start with k-means, put border links in the highest number set, move a link from a vehicle that does most to one that does least, try all combinations of moving cycles, swap the most and least used vehicles, try to move a cycle from an expensive tour to a cheap one, move a cycle even if it does not give

29 Figure 19: Snowplan: Mantorp, prel tests, q = 4.

30 Figure 20: Snowplan: Åtvidaberg, G100, q = 4.

Figure 21: Snowplan: Åtvidaberg, G100, q = 4.

31 Figure 22: Snowplan: Åtvidaberg, G100, q = 4.

Figure 23: 1045: Vadstena l4, Vineopt, q = 6.

32 Figure 24: 1046: Vadstena l4, Vineopt, zoom, q = 6. improvement, try once for each pair of vehicles, do not use a ﬁxed core. In simulated annealing, do 7 interior iterations (with the same temperature), use start temperature 1000 och cooling factor 0.3.

We do not claim that these settings are the best, but only consider them a promising possibility. If this tool is to be used on a certain city, one would obviously do additional parameter tuning.

We have chosen a number of interesting settings to try more. They are shown in table 2.

We will also use the notation Ak,Bk etc, which denotes setting A with k iterations. For the setting A1, B1 etc, the method becomes a one-shot heuristic, which however may be rather good. If the starting solution is very good, the succeeding changes may not give much improvement, and the question is how many iterations one should spend on trying to improve the solution. In some cases, one may spend a fair amount of time for quite small improvements or none at all.

We also made some preliminary tests with fixed core and vehicle area reduction on the instance atvid, with q = 4, n = 220 and m = 290. Common settings are in principle A above, but with start variant 15. Preliminary conclusions are that the effect of fixed core is very small, and that vehicle area reduction doesn’t pay off. A disadvantage is that the check of which links to include in the rural postman problem must be redone after each change of the allocation, i.e. each iteration. The gain in time of a smaller rural postman problem is much less the the increase of time needed for the additional calculations.

33 Name MD RS SV IV ChV CyV TC II RC A 1 113 3 5 4 070.3 B 1 1 13 3 5 4 0 10 0.5 C 1 114 3 5 4 070.3 D 1 1 15 3 5 4 0 10 0.5 E 1 116 3 5 4 070.3 F 1 113 3 5 4 170.3 G 1 013 3 5 4 070.3 H 1 1 73 5 4 070.3 I 0 013 3 5 4 070.3 J 1 110 3 5 4 070.3 K 1 1 10 5 0 070.3 L 1 1 10 2 0 070.3

Table 2: Settings to try. Additionally MI=30, MV=3, CI=0, FC=0, ST=1000 for all.

Name FC FD UF UD Obj Time Q 0 - 0 - 10063.39 7.85 R 1 0.5 0 - 10553.78 7.77 S 0 - 1 1.5 10463.08 15.11 T 1 0.5 1 1.5 10667.33 14.55

Table 3: Tests with ﬁxed core and vehicle area reduction on atvid.

7.3 Tests with manual improvements

We have tried manual improvement on the example Vadstena, with q = 6. After the ﬁrst run of Snowplan, we had the objective function value 12152. We brought the allocation into Vineopt, and move some links from the vehicles that did most, to a vehicle that did much less. When we read this allocation into Snowplan, we got the value 12123. Repeating this several times, gave the values 12002, 12015, 11994, and 11961. We see that it gives signiﬁcant improvements. However, the improvements are not monotone, since we don’t know exactly how the vehicles will travel after the change. Finally we let Snowplan do several more iterations with the best solution, which improved it to 11934. We also have the lower bound l1 = 11791, so this solution is quite close to the optimal one. The relative error is around 1%.

7.4 Tests of one-step methods

In the tables 35 - 46 in the Appendix, we give the results for the settings with only one iteration, followed by one round of improvements. Here we use four vehicles. In the tables 4 and 6, we give the objective function values of the diﬀerent settings. In the tables 5 and 7, we give the times. Here we ﬁnd that the times are very reasonable.

In table 8, we give the minimal (best) and maximal (worst) objective function values for each problem, and also give which setting gave this. We do the same with the solution times, and ﬁnally give the lower bounds.

34 Name A1 B1 C1 D1 E1 F1 atvid-s 4292 4202 4098 4043 4002 4097 brokind 5761 5766 5865 5833 6060 6088 colonia 1021 1000 1012 1001 992 1014 ekangen 5077 5712 5653 5736 5567 4999 mantorp-c 4511 4556 4558 4757 4576 4611 mantorp 7598 7869 7974 7502 7685 7685 rimforsa-o 3485 3427 3388 3546 3373 3469 skanninge-c 4492 4673 4649 4642 4668 4688 skanninge-n 2433 2446 2451 2438 2362 2458 skanninge 4712 4810 4486 4772 4610 4834 studentryd-e 1690 1628 1639 1645 1588 1649 sturefors 6477 6633 6449 6576 6449 6495 valla-liu 2718 2676 2607 2686 2864 2676 vikingstad 7180 7178 7456 7647 7723 8918 atvid 10381 10501 11209 10581 10727 10622 borensberg 11423 9851 11940 10360 10761 10389 kisa 9133 9454 9650 9058 9089 8720 liu 5161 5568 5429 5629 5338 5205 rimforsa 6298 6010 6016 6099 5880 6028 ryd1 4110 4085 4352 4094 4073 4380 studentryd-int 3083 2995 2991 3034 3024 3126 studentryd 5439 5396 5545 5446 5498 5438 linghem 11731 11219 11997 11423 11259 11320 ljungsbro 23866 24235 23864 24506 23571 25426 malmslatt 10867 10712 10529 10973 10538 11030 mjolby 21177 21117 21225 23912 21185 21345 ryd-m 17561 17558 16718 17158 17732 17242 ryd 12612 12889 12341 12576 12273 12477 ryd2 14359 14400 14948 14250 14280 15438 soderkoping 18017 18440 18950 18818 17928 18680 vadstena 12548 12603 12799 12504 12531 12843 valla 11094 10217 9928 9733 9709 9713

Table 4: Results, objectives, with settings A1-F1.

35 Name A1 B1 C1 D1 E1 F1 atvid-s 0.06 0.06 0.06 0.07 0.06 0.07 brokind 0.20 0.20 0.20 0.20 0.20 0.19 colonia 0.04 0.04 0.03 0.03 0.04 0.04 ekangen 0.15 0.20 0.16 0.15 0.15 0.15 mantorp-c 0.11 0.11 0.10 0.11 0.11 0.10 mantorp 0.22 0.24 0.25 0.21 0.22 0.23 rimforsa-o 0.08 0.08 0.08 0.08 0.07 0.08 skanninge-c 0.12 0.13 0.11 0.11 0.10 0.12 skanninge-n 0.07 0.06 0.05 0.05 0.05 0.06 skanninge 0.11 0.11 0.12 0.11 0.11 0.12 studentryd-e 0.07 0.07 0.06 0.06 0.06 0.07 sturefors 0.21 0.27 0.21 0.22 0.21 0.21 valla-liu 0.18 0.17 0.17 0.21 0.17 0.17 vikingstad 0.22 0.21 0.21 0.22 0.22 0.21 atvid 0.33 0.33 0.34 0.32 0.34 0.32 borensberg 0.30 0.28 0.30 0.28 0.34 0.30 kisa 0.87 0.66 0.75 0.70 0.71 0.70 liu 0.63 0.60 0.63 0.60 0.61 0.63 rimforsa 0.36 0.35 0.42 0.35 0.37 0.36 ryd1 1.47 1.45 1.19 1.30 1.28 0.94 studentryd-int 0.43 0.43 0.39 0.38 0.41 0.39 studentryd 1.28 1.37 1.33 1.48 1.36 1.41 linghem 8.74 9.05 8.96 8.80 9.32 8.61 ljungsbro 8.24 8.25 8.24 8.28 8.39 8.16 malmslatt 2.16 2.30 2.11 2.23 2.37 2.11 mjolby 1.38 1.54 1.44 1.53 1.54 1.46 ryd-m 23.15 23.18 23.34 23.28 22.95 23.29 ryd 15.45 14.66 14.58 14.80 15.78 14.87 ryd2 16.52 16.61 16.62 16.61 16.17 16.51 soderkoping 5.03 5.07 4.69 4.61 4.80 5.05 vadstena 4.01 3.49 3.68 3.49 3.74 3.65 valla 2.65 2.50 2.52 2.56 2.55 2.78

Table 5: Results, times, with settings A1-F1.

36 Name G1 H1 I1 J1 K1 L1 atvid-s 4029 4020 4084 4069 4876 5263 brokind 5647 6567 5962 6531 9954 9803 colonia 1007 1011 1002 1024 1071 1065 ekangen 5170 5202 5141 5513 7998 8099 mantorp-c 4692 5018 4642 4906 7403 7247 mantorp 7925 8095 7573 7933 11861 12261 rimforsa-o 3466 3384 3362 3380 4972 5133 skanninge-c 4540 4872 4863 4887 6535 6372 skanninge-n 2438 2475 2476 2501 3083 3366 skanninge 4764 4898 4484 4930 6782 6875 studentryd-e 1694 1641 1666 1619 1936 1873 sturefors 6466 6908 6629 6958 10338 9429 valla-liu 2754 2949 2835 2983 3600 3426 vikingstad 7768 7915 7512 8250 11215 11352 atvid 10355 10650 10620 11035 15662 14978 borensberg 10358 10137 10922 10154 13794 13834 kisa 9491 9550 8942 9530 14403 14554 liu 5184 5274 5516 5298 6838 6845 rimforsa 5842 6747 5898 6612 9653 9807 ryd1 4147 4050 4106 4037 5457 5554 studentryd-int 3085 3158 3172 3151 4118 4206 studentryd 5254 5425 5388 5419 7301 7021 linghem 11764 11212 11564 11178 16745 16872 ljungsbro 24840 23960 23937 24001 36856 36365 malmslatt 11030 11358 11080 11373 16443 15939 mjolby 21076 23044 21775 23045 31275 32309 ryd-m 17744 17340 17491 17201 23525 23316 ryd 12833 12380 12316 12399 17273 17174 ryd2 14952 14774 14293 14984 19164 19687 soderkoping 19067 18081 18624 17976 24275 24289 vadstena 12543 12565 12795 12570 17762 17801 valla 10056 9764 9982 9797 13165 13054

Table 6: Results, objectives, with settings G1-L1.

37 Name G1 H1 I1 J1 K1 L1 atvid-s 0.05 0.06 0.06 0.06 0.04 0.05 brokind 0.25 0.21 0.17 0.21 0.18 0.19 colonia 0.04 0.03 0.03 0.04 0.03 0.03 ekangen 0.15 0.15 0.12 0.15 0.12 0.13 mantorp-c 0.11 0.11 0.09 0.11 0.08 0.08 mantorp 0.24 0.24 0.18 0.22 0.23 0.18 rimforsa-o 0.09 0.07 0.06 0.08 0.06 0.06 skanninge-c 0.11 0.10 0.08 0.11 0.09 0.08 skanninge-n 0.05 0.06 0.05 0.05 0.04 0.04 skanninge 0.12 0.11 0.10 0.11 0.09 0.09 studentryd-e 0.06 0.07 0.06 0.06 0.05 0.05 sturefors 0.29 0.22 0.18 0.22 0.18 0.18 valla-liu 0.17 0.19 0.15 0.19 0.14 0.13 vikingstad 0.22 0.22 0.18 0.25 0.18 0.18 atvid 0.32 0.34 0.30 0.36 0.33 0.30 borensberg 0.28 0.30 0.25 0.31 0.27 0.25 kisa 0.66 0.73 0.71 0.76 0.72 0.73 liu 0.72 0.61 0.54 0.66 0.67 0.72 rimforsa 0.32 0.39 0.29 0.37 0.32 0.37 ryd1 1.17 1.36 1.28 1.27 1.49 1.46 studentryd-int 0.36 0.40 0.31 0.37 0.39 0.59 studentryd 1.30 1.33 1.33 1.33 1.63 1.42 linghem 8.93 8.97 8.17 8.67 16.14 15.91 ljungsbro 8.20 8.31 7.58 8.44 13.52 14.09 malmslatt 2.18 2.21 2.01 2.16 2.79 2.86 mjolby 1.37 1.43 1.20 1.48 1.73 1.56 ryd-m 22.87 22.93 21.75 23.07 45.44 46.06 ryd 14.17 14.86 13.78 14.90 28.06 27.65 ryd2 16.25 16.44 15.51 16.52 30.60 31.26 soderkoping 4.54 5.02 4.38 4.84 7.79 7.28 vadstena 3.57 3.48 3.25 3.82 4.89 5.24 valla 2.45 2.55 2.27 2.55 3.45 3.37

Table 7: Results, times, with settings G1-L1.

38 Obj Time Lowb Name Min Max Min Max l0 l1 atvid-s E1 4002 L1 5263 K1 0.04 D1 0.07 3335 3690 brokind G1 5647 K1 9954 I1 0.17 G1 0.25 4130 5225 colonia E1 992 K1 1071 C1 0.03 A1 0.04 964 981 ekangen F1 4999 L1 8099 I1 0.12 B1 0.20 4086 4707 mantorp-c A1 4511 K1 7403 K1 0.08 A1 0.11 3604 4174 mantorp D1 7502 L1 12261 I1 0.18 C1 0.25 5990 6948 rimforsa-o I1 3362 L1 5133 I1 0.06 G1 0.09 2643 3105 skanninge-c A1 4492 K1 6535 I1 0.08 B1 0.13 3665 4259 skanninge-n E1 2362 L1 3366 K1 0.04 A1 0.07 2017 2248 skanninge I1 4484 L1 6875 K1 0.09 C1 0.12 3703 4300 studentryd-e E1 1588 K1 1936 K1 0.05 A1 0.07 1479 1559 sturefors C1 6449 K1 10338 I1 0.18 G1 0.29 4904 5826 valla-liu C1 2607 K1 3600 L1 0.13 D1 0.21 2279 2533 vikingstad B1 7178 L1 11352 I1 0.18 J1 0.25 5038 6341 atvid G1 10355 K1 15662 I1 0.30 J1 0.36 8225 9657 borensberg B1 9851 L1 13834 I1 0.25 E1 0.34 7042 8785 kisa F1 8720 L1 14554 B1 0.66 A1 0.87 6476 7981 liu A1 5161 L1 6845 I1 0.54 G1 0.72 4191 4931 rimforsa G1 5842 L1 9807 I1 0.29 C1 0.42 4346 5266 ryd1 J1 4037 L1 5554 F1 0.94 K1 1.49 3550 3926 studentryd-int C1 2991 L1 4206 I1 0.31 L1 0.59 2668 2926 studentryd G1 5254 K1 7301 A1 1.28 K1 1.63 4450 5074 linghem J1 11178 L1 16872 I1 8.17 K1 16.14 9080 10562 ljungsbro E1 23571 K1 36856 I1 7.58 L1 14.09 18316 22207 malmslatt C1 10529 K1 16443 I1 2.01 L1 2.86 8978 10187 mjolby G1 21076 L1 32309 I1 1.20 K1 1.73 15668 18851 ryd-m C1 16718 K1 23525 I1 21.75 L1 46.06 14437 16159 ryd E1 12273 K1 17273 I1 13.78 K1 28.06 10678 11910 ryd2 D1 14250 L1 19687 I1 15.51 L1 31.26 12268 13707 soderkoping E1 17928 L1 24289 I1 4.38 K1 7.79 14548 16913 vadstena D1 12504 L1 17801 I1 3.25 L1 5.24 10532 12147 valla E1 9709 K1 13165 I1 2.27 K1 3.45 7779 9039

Table 8: Results, best and worst objectives and times, settings A1-L1, q = 4.

39 We ﬁnd that the settings K1 and L1 give bad solutions, while good solutions are obtained by several other settings. Some solutions are actually close to the lower bound, but most are not very close.

Setting E1 gives the best solution for 8 instances, followed by G1 and C1, both giving the best for 5 instances.

40 7.5 Tests with 30 iterations

In the tables 47 - 58 in the Appendix, we present results of our main test with the chosen settings, A - L, using 30 iterations. In tables 9, 10, 11 and 12, we do a comparison for each instance.

Name A B C D E F atvid-s 4014 3842 3859 3921 3982 3862 brokind 5916 5761 5810 5701 5810 5748 colonia 1001 1000 1000 995 1009 1000 ekangen 5607 5299 5595 5214 5327 4995 mantorp-c 4557 4662 4481 4590 4672 4509 mantorp 7353 7520 7477 7650 7523 7777 rimforsa-o 3412 3337 3373 3389 3477 3506 skanninge-c 4375 4484 4458 4545 4676 4647 skanninge-n 2373 2370 2409 2376 2366 2372 skanninge 4470 4628 4476 4521 4562 4450 studentryd-e 1608 1604 1623 1623 1621 1599 sturefors 6988 6686 6482 6418 6744 6466 valla-liu 2692 2665 2712 2737 2845 2669 vikingstad 7574 7234 7360 7534 7547 7365 atvid 10444 10489 10594 10623 10595 10501 borensberg 9451 9481 10073 9894 10257 10190 kisa 9146 9519 8903 8962 9495 9459 liu 5236 5320 5249 5343 5323 5150 rimforsa 5875 5724 5981 5813 5872 5792 ryd1 4094 4028 3993 4038 4050 4046 studentryd-int 3102 3113 2995 3160 3087 3051 studentryd 5242 5304 5256 5243 5345 5336 linghem 11955 11336 11668 11542 11423 11236 ljungsbro 23290 23795 24009 24230 24403 23835 malmslatt 10658 10476 10521 10830 10411 10508 mjolby 21041 21524 21063 21446 21229 22548 ryd-m 16994 16920 16971 16996 16868 17226 ryd 12137 12240 12213 12539 12340 12683 ryd2 14876 14184 14087 14242 13948 14238 soderkoping 17566 18020 18621 17821 17659 19055 vadstena 12416 12582 12474 12492 12530 12343 valla 9897 9721 9593 9730 10022 9691

Table 9: Results, objectives, with settings A-F.

In table 13, we give the best and the worst of the results. Checking how many times each setting gave the best objective function value for q = 4 and 30 iterations, we get the following. A: 8, G: 6, B: 4, C: 4, E: 4, F: 4, D: 2, H: 1, I: 0, J: 0. The randomness makes it hard to pick one “best” setting, but here A seems to be a good choice, since it gives the best solution for 8 problems.

41 Name A B C D E F atvid-s 1.35 1.42 1.31 1.36 1.37 1.38 brokind 4.45 4.40 4.52 4.51 4.77 4.49 colonia 0.77 0.80 0.85 0.80 0.80 0.78 ekangen 3.56 3.72 3.45 3.24 3.41 3.16 mantorp-c 2.43 2.69 2.41 2.46 2.51 2.32 mantorp 4.83 5.52 5.75 5.18 4.88 5.06 rimforsa-o 1.77 1.88 1.76 1.77 1.76 1.93 skanninge-c 2.43 2.76 2.51 2.57 2.55 2.74 skanninge-n 1.18 1.24 1.33 1.36 1.32 1.32 skanninge 2.65 2.96 2.77 3.02 2.53 2.71 studentryd-e 1.46 2.27 1.49 1.54 1.48 1.57 sturefors 4.85 5.10 5.29 4.55 5.43 5.36 valla-liu 3.86 3.87 3.98 4.06 4.15 3.98 vikingstad 4.84 5.69 5.00 4.71 4.60 4.88 atvid 8.05 8.92 8.58 7.98 8.03 7.12 borensberg 6.49 7.11 7.23 6.48 8.54 7.26 kisa 18.09 15.88 17.17 17.34 16.52 15.56 liu 16.81 14.71 18.70 15.87 16.46 15.75 rimforsa 9.05 8.37 9.77 9.02 8.75 8.06 ryd1 30.48 25.58 27.74 30.46 30.17 27.94 studentryd-int 10.03 8.93 9.57 11.16 9.06 8.39 studentryd 30.41 29.56 33.61 30.25 31.41 28.50 linghem 182.20 182.50 179.77 179.49 190.59 190.22 ljungsbro 161.62 165.63 173.87 161.00 163.33 164.44 malmslatt 49.55 46.92 50.32 48.66 46.41 47.04 mjolby 34.11 38.65 37.27 37.98 35.47 34.24 ryd-m 420.20 422.05 416.06 417.10 419.66 416.16 ryd 277.52 270.36 275.99 270.01 269.26 273.43 ryd2 298.73 300.00 297.11 299.23 297.34 297.81 soderkoping 95.51 98.58 105.69 98.41 101.25 98.66 vadstena 69.41 74.77 74.51 79.45 76.18 72.93 valla 58.29 58.94 58.90 58.57 58.02 56.87

Table 10: Results, times, with settings A-F.

42 Name G H I J K L atvid-s 3957 3930 4290 4012 4881 4657 brokind 5859 6567 5878 6531 9134 9082 colonia 989 1005 1006 1020 1042 1034 ekangen 4933 5144 5585 5397 6770 7265 mantorp-c 4473 4904 4561 4855 6760 7093 mantorp 7555 7883 7432 7829 11477 11332 rimforsa-o 3441 3344 3435 3380 4686 4704 skanninge-c 4712 4867 4768 4828 6026 6095 skanninge-n 2390 2415 2458 2439 2791 2840 skanninge 4477 4866 4593 4708 6098 5923 studentryd-e 1621 1628 1619 1616 1720 1732 sturefors 6613 6895 6574 6958 8674 9398 valla-liu 2657 2829 2782 2830 3323 3286 vikingstad 7604 7877 8905 7921 10646 10238 atvid 10443 10600 10785 10801 13963 13890 borensberg 9901 10115 10675 9995 12660 12124 kisa 9213 9498 9807 9447 13040 12769 liu 5168 5274 5622 5227 6250 6353 rimforsa 5894 6592 6020 6612 8794 8837 ryd1 4030 4027 4190 4010 4753 4767 studentryd-int 3058 3158 3025 3099 3739 3697 studentryd 5295 5363 5493 5400 6287 6307 linghem 11611 11132 11748 11161 14713 14804 ljungsbro 23584 23960 24786 23878 32700 32543 malmslatt 10593 11306 11031 11250 13933 14130 mjolby 23003 22623 21558 22750 26573 28322 ryd-m 17202 17008 17715 17194 20340 20164 ryd 12367 12380 12582 12336 14486 14262 ryd2 14575 14774 14367 14958 16957 16699 soderkoping 18220 17840 18436 17802 22537 22223 vadstena 12365 12560 12692 12570 15545 15493 valla 9517 9764 9677 9612 11486 11371

Table 11: Results, objectives, with settings G-L.

43 Name G H I J K L atvid-s 1.29 1.34 1.40 1.73 1.00 1.02 brokind 4.31 4.59 4.66 5.12 3.02 3.15 colonia 0.77 0.77 0.79 0.96 0.69 0.77 ekangen 3.20 3.48 3.92 3.76 2.30 2.32 mantorp-c 2.33 3.68 2.43 2.64 1.79 1.67 mantorp 4.69 5.13 4.72 4.97 3.36 3.47 rimforsa-o 1.67 1.73 1.81 1.83 1.28 1.29 skanninge-c 2.35 2.41 3.39 2.47 1.68 1.71 skanninge-n 1.38 1.21 1.33 1.25 0.90 0.88 skanninge 2.52 2.50 2.59 2.88 1.79 1.82 studentryd-e 1.52 3.08 1.69 2.65 1.14 1.03 sturefors 4.11 4.92 5.69 4.92 3.39 3.24 valla-liu 3.90 4.13 4.15 4.56 2.68 2.59 vikingstad 4.85 5.12 5.01 5.21 3.21 3.63 atvid 7.59 8.45 7.77 7.57 5.40 5.40 borensberg 6.04 7.62 7.25 6.67 4.34 4.49 kisa 16.43 17.27 15.78 15.08 15.64 15.06 liu 16.36 19.08 14.78 13.81 14.91 15.69 rimforsa 8.00 8.35 7.89 14.07 6.18 6.84 ryd1 27.62 29.00 24.86 23.80 28.20 26.87 studentryd-int 9.32 8.94 9.20 9.61 7.15 7.22 studentryd 29.76 29.46 24.93 27.88 29.28 29.70 linghem 177.07 182.42 171.83 190.10 274.90 259.57 ljungsbro 166.58 165.24 154.99 185.94 229.68 224.47 malmslatt 45.78 48.02 41.85 45.16 50.67 47.65 mjolby 32.53 34.55 30.22 30.38 34.09 32.29 ryd-m 403.35 413.78 406.10 421.59 787.02 697.96 ryd 263.03 277.56 262.40 271.87 425.07 411.32 ryd2 286.28 302.82 290.70 308.63 522.63 465.33 soderkoping 95.05 100.43 98.56 101.08 131.28 129.88 vadstena 67.97 70.67 61.92 74.79 91.95 84.71 valla 55.48 54.24 49.57 58.56 65.81 63.72

Table 12: Results, times, with settings G-L.

44 Obj Time Lowb Name Min Max Min Max l0 l1 atvid-s B 3842 K 4881 K 1.00 J 1.73 3335 3690 brokind D 5701 K 9134 K 3.02 J 5.12 4130 5225 colonia G 989 K 1042 K 0.69 J 0.96 964 981 ekangen G 4933 L 7265 K 2.30 I 3.92 4086 4707 mantorp-c G 4473 L 7093 L 1.67 H 3.68 3604 4174 mantorp A 7353 K 11477 K 3.36 C 5.75 5990 6948 rimforsa-o B 3337 L 4704 K 1.28 F 1.93 2643 3105 skanninge-c A 4375 L 6095 K 1.68 I 3.39 3665 4259 skanninge-n E 2366 L 2840 L 0.88 G 1.38 2017 2248 skanninge F 4450 K 6098 K 1.79 D 3.02 3703 4300 studentryd-e F 1599 L 1732 L 1.03 H 3.08 1479 1559 sturefors D 6418 L 9398 L 3.24 I 5.69 4904 5826 valla-liu G 2657 K 3323 L 2.59 J 4.56 2279 2533 vikingstad B 7234 K 10646 K 3.21 B 5.69 5038 6341 atvid G 10443 K 13963 K 5.40 B 8.92 8225 9657 borensberg A 9451 K 12660 K 4.34 E 8.54 7042 8785 kisa C 8903 K 13040 L 15.06 A 18.09 6476 7981 liu F 5150 L 6353 J 13.81 H 19.08 4191 4931 rimforsa B 5724 L 8837 K 6.18 J 14.07 4346 5266 ryd1 C 3993 L 4767 J 23.80 A 30.48 3550 3926 studentryd-int C 2995 K 3739 K 7.15 D 11.16 2668 2926 studentryd A 5242 L 6307 I 24.93 C 33.61 4450 5074 linghem H 11132 L 14804 I 171.83 K 274.90 9080 10562 ljungsbro A 23290 K 32700 I 154.99 K 229.68 18316 22207 malmslatt E 10411 L 14130 I 41.85 K 50.67 8978 10187 mjolby A 21041 L 28322 I 30.22 B 38.65 15668 18851 ryd-m E 16868 K 20340 G 403.35 K 787.02 14437 16159 ryd A 12137 K 14486 I 262.40 K 425.07 10678 11910 ryd2 E 13948 K 16957 G 286.28 K 522.63 12268 13707 soderkoping A 17566 K 22537 G 95.05 K 131.28 14548 16913 vadstena F 12343 K 15545 I 61.92 K 91.95 10532 12147 valla G 9517 K 11486 I 49.57 K 65.81 7779 9039

Table 13: Results, best and worst objectives and times, settings A-L, q = 4.

45 Times are obviously longer here than with one iteration. However, it is still possible to use the method, from a practical perspective. Probably the planning does not have to be changed every day. If there are circumstances that require a new planning, it is still much faster to use this tool than to do it by hand. One may note that K is fast for small problems, but slow for large, and since in addition gives the worst solution for many problems, it is not the method of our choice.

In table 14, we investigate the setting A30 more thoroughly, by giving the results together with problem sizes, lower bounds and relative errors. The quality of the lower bounds may vary between instances.

Name n m q Obj Time l1 Err atvid-s 68 95 4 3828.73 1.49 3690.50 0.036 brokind 179 199 4 5766.56 4.24 5225.82 0.094 colonia 24 36 4 1001.79 0.75 981.29 0.020 ekangen 149 184 4 5502.08 3.41 4707.38 0.144 mantorp-c 122 145 4 4612.99 2.36 4174.39 0.095 mantorp 185 226 4 7592.03 4.86 6948.47 0.085 rimforsa-o 94 109 4 3381.52 1.75 3105.75 0.082 skanninge-c 122 152 4 4623.21 2.49 4259.22 0.079 skanninge-n 62 80 4 2324.27 1.31 2248.61 0.033 skanninge 127 158 4 4738.71 2.60 4300.15 0.093 studentryd-e 76 104 4 1602.75 1.45 1559.96 0.027 sturefors 183 212 4 6443.14 4.74 5826.71 0.096 valla-liu 160 213 4 2764.75 4.01 2533.03 0.084 vikingstad 178 213 4 8445.48 4.72 6341.68 0.249 atvid 220 290 4 10753.67 8.41 9657.40 0.102 borensberg 211 257 4 10337.00 6.30 8785.55 0.150 kisa 311 359 4 8765.30 16.63 7981.47 0.089 liu 291 394 4 5335.31 15.18 4931.24 0.076 rimforsa 233 262 4 5893.20 8.07 5266.41 0.106 ryd1 382 508 4 4040.36 29.37 3926.69 0.028 studentryd-int 235 316 4 3062.97 8.91 2926.23 0.045 studentryd 387 519 4 5298.02 28.83 5074.35 0.042 linghem 836 1129 4 11800.73 175.34 10562.68 0.105 ljungsbro 812 1024 4 23566.40 159.42 22207.49 0.058 malmslatt 489 626 4 10450.04 44.52 10187.53 0.025 mjolby 419 524 4 21320.07 31.40 18851.91 0.116 ryd-m 1107 1495 4 17133.13 439.34 16159.83 0.057 ryd 952 1291 4 12610.77 270.47 11910.95 0.055 ryd2 987 1337 4 14768.43 295.82 13707.77 0.072 soderkoping 666 910 4 18538.84 99.72 16913.25 0.088 vadstena 581 778 4 12389.38 74.34 12147.58 0.020 valla 527 707 4 9529.93 58.79 9039.05 0.052

Table 14: Results with setting A30, q = 4. Relative error to lower bound.

46 7.6 Diﬀerent number of iterations

In the tables 63 and 64 in the Appendix, we give results for methods A5 and A100, i.e. setting A with 5 or 100 iterations. In tables 15 and 16, we sum up the results for setting A with 1, 5, 30 and 100 iterations. The question here is the value of continued iterations. Since the methods contain randomness, we see cases where more iterations give a worse solution. That will obviously not happen within the same run.

Name A1 A5 A30 A100 atvid-s 4292 4005 4014 3971 brokind 5761 5804 5916 5585 colonia 1021 994 1001 990 ekangen 5077 5542 5607 5334 mantorp-c 4511 4718 4557 4382 mantorp 7598 7458 7353 7544 rimforsa-o 3485 3387 3412 3340 skanninge-c 4492 4638 4375 4616 skanninge-n 2433 2383 2373 2414 skanninge 4712 4651 4470 4671 studentryd-e 1690 1651 1608 1603 sturefors 6477 6495 6988 6428 valla-liu 2718 2676 2692 2624 vikingstad 7180 7561 7574 7287 atvid 10381 10506 10444 10179 borensberg 11423 10062 9451 10065 kisa 9133 9393 9146 9496 liu 5161 5242 5236 5197 rimforsa 6298 5908 5875 5700 ryd1 4110 4008 4094 4082 studentryd-int 3083 2993 3102 3062 studentryd 5439 5412 5242 5275 linghem 11731 11226 11955 - ljungsbro 23866 23525 23290 - malmslatt 10867 10492 10658 - mjolby 21177 21270 21041 - ryd-m 17561 16912 16994 - ryd 12612 12782 12137 - ryd2 14359 14186 14876 - soderkoping 18017 17935 17566 - vadstena 12548 12438 12416 - valla 11094 9648 9897 -

Table 15: Results, objectives, with setting A with 1, 5, 30 and 100 iterations.

Times grow almost linearly, as expected. As there can be several iterations without improvement of the objective function value, improvement of the objective function value is not very regular. Usually, doing 100 iterations gives better solutions than doing 30 or less, so it is mainly a question of how much time one wants to spend.

47 Name A1 A5 A30 A100 atvid-s 0.06 0.25 1.35 4.49 brokind 0.20 0.76 4.45 13.82 colonia 0.04 0.14 0.77 2.64 ekangen 0.15 0.64 3.56 11.19 mantorp-c 0.11 0.42 2.43 7.59 mantorp 0.22 0.95 4.83 18.71 rimforsa-o 0.08 0.32 1.77 6.24 skanninge-c 0.12 0.50 2.43 8.20 skanninge-n 0.07 0.22 1.18 3.94 skanninge 0.11 0.46 2.65 8.63 studentryd-e 0.07 0.26 1.46 4.76 sturefors 0.21 0.85 4.85 15.32 valla-liu 0.18 0.67 3.86 15.07 vikingstad 0.22 0.77 4.84 14.91 atvid 0.33 1.30 8.05 26.40 borensberg 0.30 1.19 6.49 22.56 kisa 0.87 3.46 18.09 65.91 liu 0.63 2.93 16.81 59.56 rimforsa 0.36 1.40 9.05 28.22 ryd1 1.47 5.33 30.48 100.28 studentryd-int 0.43 1.62 10.03 34.76 studentryd 1.28 5.27 30.41 106.38 linghem 8.74 32.94 182.20 - ljungsbro 8.24 31.05 161.62 - malmslatt 2.16 8.60 49.55 - mjolby 1.38 5.83 34.11 - ryd-m 23.15 77.47 420.20 - ryd 15.45 49.41 277.52 - ryd2 16.52 55.90 298.73 - soderkoping 5.03 18.30 95.51 - vadstena 4.01 13.21 69.41 - valla 2.65 9.93 58.29 -

Table 16: Results, times, with setting A with 1, 5, 30 and 100 iterations.

48 7.7 Tests with ﬁxed core and vehicle area reduction

Tests with ﬁxed core and vehicle area reduction, methods Q - T with 30 iterations, are reported in tables 59 - 62 in the Appendix, and summed up in tables 17, 18 and 19. Unfortunately these settings does not seem to give improvements.

Name Q R S T atvid-s 3848 3873 3881 3809 brokind 5723 5801 5963 5843 colonia 1010 1004 997 1003 ekangen 5362 5286 4972 5302 mantorp-c 4692 4553 4562 4754 mantorp 7558 7363 7511 7665 rimforsa-o 3471 3335 3393 3466 skanninge-c 4520 4576 4588 4543 skanninge-n 2417 2419 2402 2434 skanninge 4621 4998 4662 4471 studentryd-e 1594 1598 1583 1598 sturefors 7063 6807 7023 6960 valla-liu 2641 2709 2683 2808 vikingstad 7644 7698 7742 7726 atvid 10248 10336 10302 10785 borensberg 9608 9698 9650 10187 kisa 8856 9300 9274 9468 liu 5299 5280 5224 5340 rimforsa 5861 5897 5928 5800 ryd1 4159 4020 4159 4299 studentryd-int 3081 3151 3070 3037 studentryd 5344 5291 5235 5306 linghem 11636 12175 11795 11329 ljungsbro 23731 24535 24198 23980 malmslatt 10706 10637 10867 10636 mjolby 21181 20921 20850 20923 ryd-m 17115 16814 16544 17003 ryd 12637 12500 12532 12465 ryd2 14564 14435 14022 14117 soderkoping 18280 17995 18220 18161 vadstena 12411 12370 12536 12487 valla 9765 9966 9641 9785

Table 17: Results, objectives, with settings Q-T.

49 Name Q R S T atvid-s 1.37 1.30 2.34 2.36 brokind 4.32 4.10 8.11 7.75 colonia 0.80 0.72 1.02 0.98 ekangen 3.40 3.33 6.42 6.83 mantorp-c 2.46 2.41 4.76 4.51 mantorp 5.00 5.43 9.76 9.40 rimforsa-o 1.90 1.89 3.16 3.18 skanninge-c 2.60 2.43 4.77 4.97 skanninge-n 1.42 1.28 1.97 2.08 skanninge 2.69 2.32 4.95 5.03 studentryd-e 1.48 1.52 2.49 2.63 sturefors 4.99 4.73 8.63 8.77 valla-liu 4.07 4.01 8.73 7.97 vikingstad 4.72 4.80 9.14 8.77 atvid 7.90 7.90 15.35 15.04 borensberg 6.66 6.33 12.49 12.99 kisa 17.11 17.63 29.73 30.45 liu 15.51 16.57 29.74 29.93 rimforsa 8.06 9.22 13.88 17.71 ryd1 28.04 28.84 51.33 49.97 studentryd-int 9.08 8.92 17.82 18.68 studentryd 29.85 30.72 58.03 51.16 linghem 181.93 181.13 287.08 305.95 ljungsbro 161.45 163.81 251.54 244.40 malmslatt 48.84 49.56 81.69 81.18 mjolby 34.93 37.50 57.13 59.60 ryd-m 439.95 434.84 595.10 605.61 ryd 277.91 289.60 409.68 413.02 ryd2 311.75 302.58 445.76 448.99 soderkoping 100.78 101.29 164.85 166.28 vadstena 73.69 72.51 123.40 121.98 valla 58.39 58.80 102.45 99.26

Table 18: Results, times, with settings Q-T.

50 Obj Time Name Min Max Min Max atvid-s T 3809 S 3881 R 1.30 T 2.36 brokind Q 5723 S 5963 R 4.10 S 8.11 colonia S 997 Q 1010 R 0.72 S 1.02 ekangen S 4972 Q 5362 R 3.33 T 6.83 mantorp-c R 4553 T 4754 R 2.41 S 4.76 mantorp R 7363 T 7665 Q 5.00 S 9.76 rimforsa-o R 3335 Q 3471 R 1.89 T 3.18 skanninge-c Q 4520 S 4588 R 2.43 T 4.97 skanninge-n S 2402 T 2434 R 1.28 T 2.08 skanninge T 4471 R 4998 R 2.32 T 5.03 studentryd-e S 1583 R 1598 Q 1.48 T 2.63 sturefors R 6807 Q 7063 R 4.73 T 8.77 valla-liu Q 2641 T 2808 R 4.01 S 8.73 vikingstad Q 7644 S 7742 Q 4.72 S 9.14 atvid Q 10248 T 10785 Q 7.90 S 15.35 borensberg Q 9608 T 10187 R 6.33 T 12.99 kisa Q 8856 T 9468 Q 17.11 T 30.45 liu S 5224 T 5340 Q 15.51 T 29.93 rimforsa T 5800 S 5928 Q 8.06 T 17.71 ryd1 R 4020 T 4299 Q 28.04 S 51.33 studentryd-int T 3037 R 3151 R 8.92 T 18.68 studentryd S 5235 Q 5344 Q 29.85 S 58.03 linghem T 11329 R 12175 R 181.13 T 305.95 ljungsbro Q 23731 R 24535 Q 161.45 S 251.54 malmslatt T 10636 S 10867 Q 48.84 S 81.69 mjolby S 20850 Q 21181 Q 34.93 T 59.60 ryd-m S 16544 Q 17115 R 434.84 T 605.61 ryd T 12465 Q 12637 Q 277.91 T 413.02 ryd2 S 14022 Q 14564 R 302.58 T 448.99 soderkoping R 17995 Q 18280 Q 100.78 T 166.28 vadstena R 12370 S 12536 R 72.51 S 123.40 valla S 9641 R 9966 Q 58.39 S 102.45

Table 19: Results, best and worst objectives and times, settings Q-T.

51 7.8 Tests with ﬁve vehicles

Now present the 30 iteration tests with methods A - L and q = 5, in the same way as we previously did for q = 4.

Name A B C D E F atvid-s 3987 3967 4018 4155 3961 4156 brokind 5799 5777 5872 5847 5611 6114 colonia 1206 1200 1186 1202 1196 1187 ekangen 5096 5426 5467 5542 5311 5349 mantorp-c 4677 4805 4668 4688 4652 4819 mantorp 7705 7339 7381 7226 7672 7328 rimforsa-o 3381 3760 3480 3427 3576 3503 skanninge-c 4620 4608 4604 4526 4608 4620 skanninge-n 2561 2571 2629 2562 2600 2595 skanninge 4683 4737 4875 4868 4594 4655 studentryd-e 1818 1812 1798 1793 1791 1793 sturefors 6910 6724 6739 6870 6733 6634 valla-liu 2842 2836 2923 2852 2790 2831 vikingstad 7653 8244 7619 7529 7807 7928 atvid 10307 10802 9902 9919 10506 10270 borensberg 10159 10261 10244 10163 10343 9974 kisa 8885 9029 9140 9723 8920 9691 liu 5387 5343 5238 5426 5243 5326 rimforsa 5894 6130 6209 6147 5977 5744 ryd1 4083 4145 4137 4144 4124 4080 studentryd-int 3160 3174 3205 3197 3134 3179 studentryd 5342 5373 5396 5376 5258 5374 linghem 11166 11894 11846 11199 11046 11131 ljungsbro 23214 23435 23811 23716 22857 23029 malmslatt 10480 10604 10405 10566 10384 10704 mjolby 21525 21319 21406 21652 21330 21465 ryd-m 16663 16399 16628 16718 16249 16479 ryd 12243 12070 12401 12242 12124 12314 ryd2 14003 14247 14343 13826 13922 14130 soderkoping 17596 18856 17619 18233 17100 18411 vadstena 12246 12393 12095 12165 12355 12395 valla 9849 9574 9566 9767 9659 9633

Table 20: Results, objectives, with settings A-F.

In table 24, we give the best and the worst for q = 5.

52 Name A B C D E F atvid-s 1.90 1.95 1.89 1.96 2.10 2.20 brokind 7.29 7.80 8.55 6.61 8.53 7.20 colonia 1.03 0.98 0.93 1.06 1.09 1.16 ekangen 4.90 4.71 5.22 5.18 5.46 5.83 mantorp-c 3.43 3.44 3.47 3.54 3.67 4.17 mantorp 8.17 7.37 8.05 7.76 9.33 8.76 rimforsa-o 2.45 2.57 2.36 2.49 3.10 2.98 skanninge-c 4.26 3.40 3.49 3.56 4.34 5.31 skanninge-n 1.77 1.87 1.78 1.96 2.03 2.19 skanninge 3.91 3.80 4.10 3.99 4.53 4.33 studentryd-e 2.27 2.03 2.03 2.10 2.46 2.49 sturefors 8.54 7.54 7.49 8.04 7.73 9.11 valla-liu 6.12 6.24 6.11 6.52 5.97 6.72 vikingstad 7.61 8.52 8.94 7.56 7.70 9.44 atvid 11.66 13.42 13.26 11.96 12.16 12.52 borensberg 10.76 10.11 11.03 9.96 10.41 10.42 kisa 27.34 23.87 23.55 22.24 23.77 21.86 liu 21.91 23.33 22.61 21.59 26.22 20.69 rimforsa 11.79 12.44 12.71 11.89 13.46 11.56 ryd1 39.27 40.18 44.59 38.16 38.61 36.31 studentryd-int 13.64 14.44 15.04 12.36 13.79 16.53 studentryd 43.20 47.06 41.10 39.72 41.87 40.52 linghem 239.81 242.17 240.80 244.57 242.48 239.86 ljungsbro 216.77 216.53 216.44 215.32 231.03 220.74 malmslatt 64.59 66.09 63.79 60.84 63.48 63.31 mjolby 53.56 49.02 46.50 46.23 47.52 47.36 ryd-m 543.67 537.88 546.85 566.67 547.77 554.19 ryd 373.00 361.24 364.18 369.90 353.81 355.35 ryd2 394.04 394.64 395.50 393.49 386.62 411.53 soderkoping 131.95 134.87 146.26 132.26 130.04 136.10 vadstena 100.07 100.46 93.95 98.51 92.18 100.58 valla 81.10 78.77 75.75 79.92 74.75 78.47

Table 21: Results, times, with settings A-F.

53 Name G H I J K L atvid-s 4190 4334 4199 4380 5003 5117 brokind 6036 6883 6125 6794 9830 9743 colonia 1207 1205 1221 1202 1242 1212 ekangen 5342 5453 5739 5471 8320 7577 mantorp-c 4639 5141 4575 5067 7057 7188 mantorp 7298 8025 8278 8223 11688 12250 rimforsa-o 3554 3438 3468 3406 4946 4758 skanninge-c 4570 4818 4944 4807 6187 6518 skanninge-n 2616 2577 2674 2586 3134 3243 skanninge 4969 4744 4672 4894 6448 6573 studentryd-e 1782 1836 1838 1854 1991 1939 sturefors 6851 6878 6698 6992 9851 9632 valla-liu 2762 3125 2959 3151 3560 3521 vikingstad 7511 8147 7781 8281 11344 10952 atvid 10611 10664 10591 10697 15040 15102 borensberg 10103 10396 11271 10301 12365 12591 kisa 8708 9815 9379 9679 12604 13272 liu 5232 5396 5519 5452 6773 6616 rimforsa 6077 6851 6081 6879 9720 9683 ryd1 4160 4132 4104 4122 5156 4991 studentryd-int 3187 3316 3188 3291 4004 4038 studentryd 5320 5539 5432 5553 6831 6567 linghem 11212 11275 11928 11442 15964 15398 ljungsbro 23257 23664 23681 23781 33779 33471 malmslatt 10866 11415 10408 11166 15405 14713 mjolby 20916 22638 23925 22886 29317 28765 ryd-m 16696 16783 17008 16936 20644 20859 ryd 12185 12274 12801 12362 15049 14859 ryd2 14040 14627 14658 14611 17328 16852 soderkoping 17454 17271 18288 17515 21652 22675 vadstena 12295 12454 12940 12441 16426 16272 valla 9837 9811 9934 9822 12564 12063

Table 22: Results, objectives, with settings G-L.

54 Name G H I J K L atvid-s 1.84 2.19 1.93 1.97 1.14 1.25 brokind 6.91 8.68 6.41 6.90 3.61 4.08 colonia 1.05 1.10 1.05 1.00 0.83 0.85 ekangen 4.99 5.50 5.26 5.05 2.73 2.83 mantorp-c 3.36 4.24 3.25 3.66 1.93 2.24 mantorp 7.26 8.95 8.12 8.50 4.28 4.50 rimforsa-o 2.53 3.18 2.34 2.40 1.48 1.58 skanninge-c 3.94 4.04 4.04 3.86 1.99 2.34 skanninge-n 1.73 2.05 1.99 1.73 1.08 1.13 skanninge 4.09 4.17 3.75 5.56 2.34 2.37 studentryd-e 1.98 2.36 2.07 2.12 1.24 1.30 sturefors 7.57 7.98 8.18 7.54 3.93 7.11 valla-liu 6.01 7.72 6.44 6.57 3.14 3.45 vikingstad 7.57 8.94 7.10 8.19 3.69 6.53 atvid 14.16 14.20 14.74 13.16 8.48 7.04 borensberg 8.91 11.54 10.96 11.81 5.25 6.34 kisa 23.00 25.72 24.37 28.79 15.72 13.86 liu 22.97 24.64 22.95 23.74 15.96 16.54 rimforsa 11.85 14.56 12.61 14.30 7.31 8.42 ryd1 39.97 45.72 36.16 39.97 29.95 27.64 studentryd-int 13.87 15.89 13.70 15.14 7.83 8.40 studentryd 44.34 43.07 38.17 43.03 31.49 29.74 linghem 236.47 244.21 231.26 239.34 274.63 279.52 ljungsbro 210.47 242.81 204.00 223.22 236.16 235.81 malmslatt 63.85 65.61 55.93 66.44 54.88 51.85 mjolby 42.78 54.90 41.02 50.56 37.11 33.62 ryd-m 526.43 559.03 540.46 559.01 772.63 744.39 ryd 343.50 359.15 347.52 358.98 441.91 429.17 ryd2 384.70 398.83 378.31 393.22 500.45 496.10 soderkoping 142.92 135.22 128.92 135.78 133.75 139.24 vadstena 95.87 109.95 94.30 97.47 91.08 91.38 valla 78.91 81.16 70.91 79.98 68.02 65.19

Table 23: Results, times, with settings G-L.

55 Obj Time Lowb Name Min Max Min Max l0 l1 atvid-s E 3961 L 5117 K 1.14 F 2.20 3434 3774 brokind E 5611 K 9830 K 3.61 H 8.68 4197 5248 colonia C 1186 K 1242 K 0.83 F 1.16 1157 1174 ekangen A 5096 K 8320 K 2.73 F 5.83 4155 4751 mantorp-c I 4575 L 7188 K 1.93 H 4.24 3692 4239 mantorp D 7226 L 12250 K 4.28 E 9.33 5982 6902 rimforsa-o A 3381 K 4946 K 1.48 H 3.18 2769 3213 skanninge-c D 4526 L 6518 K 1.99 F 5.31 3750 4320 skanninge-n A 2561 L 3243 K 1.08 F 2.19 2168 2390 skanninge E 4594 L 6573 K 2.34 J 5.56 3786 4360 studentryd-e G 1782 K 1991 K 1.24 F 2.49 1651 1729 sturefors F 6634 K 9851 K 3.93 F 9.11 4940 5825 valla-liu G 2762 K 3560 K 3.14 H 7.72 2420 2663 vikingstad G 7511 K 11344 K 3.69 F 9.44 5068 6320 atvid C 9902 L 15102 L 7.04 I 14.74 8128 9503 borensberg F 9974 L 12591 K 5.25 J 11.81 6993 8666 kisa G 8708 L 13272 L 13.86 J 28.79 6449 7894 liu G 5232 K 6773 K 15.96 E 26.22 4256 4965 rimforsa F 5744 K 9720 K 7.31 H 14.56 4404 5287 ryd1 F 4080 K 5156 L 27.64 H 45.72 3640 4001 studentryd-int E 3134 L 4038 K 7.83 F 16.53 2793 3041 studentryd E 5258 K 6831 L 29.74 B 47.06 4504 5103 linghem E 11046 K 15964 I 231.26 L 279.52 8949 10372 ljungsbro E 22857 K 33779 I 204.00 H 242.81 17815 21551 malmslatt E 10384 K 15405 L 51.85 J 66.44 8851 10012 mjolby G 20916 K 29317 L 33.62 H 54.90 15273 18329 ryd-m E 16249 L 20859 G 526.43 K 772.63 14092 15745 ryd B 12070 K 15049 G 343.50 K 441.91 10483 11666 ryd2 D 13826 K 17328 I 378.31 K 500.45 12009 13391 soderkoping E 17100 L 22675 I 128.92 C 146.26 14198 16468 vadstena C 12095 K 16426 K 91.08 H 109.95 10342 11893 valla C 9566 K 12564 L 65.19 H 81.16 7700 8909

Table 24: Results, best and worst objectives and times, settings A-L.

56 7.9 Tests with six vehicles

Here we present the tests with methods A - J in 30 iteration for 6 vehicles. That many vehicles is only interesting for the larger networks.

Name A B C D E F linghem 11392 11102 11370 11125 11128 11132 ljungsbro 23934 23057 22485 23112 22677 22533 malmslatt 10296 10925 10582 10514 10354 10534 mjolby 21518 22089 21879 21727 21638 22123 ryd-m 16582 16699 16113 16281 16125 16396 ryd 12227 12387 12095 12425 11997 11943 ryd2 13929 13937 13697 13813 14134 13833 soderkoping 18309 17394 17342 18603 16820 17243 vadstena 12534 12107 12374 12218 12171 12389 valla 9367 9710 9968 9596 9368 9642

Table 25: Results, obj, with settings A-F, q = 6.

Name A B C D E F linghem 312.68 315.45 309.97 306.22 315.64 314.82 ljungsbro 275.13 273.34 282.53 303.09 277.73 279.56 malmslatt 81.40 86.22 84.32 85.99 92.05 83.41 mjolby 57.80 73.29 63.29 62.68 63.39 61.84 ryd-m 683.31 695.27 744.05 692.72 679.59 709.62 ryd 459.98 451.41 446.19 456.49 452.54 452.58 ryd2 501.52 503.46 490.39 494.95 500.29 495.19 soderkoping 178.67 179.65 176.52 177.44 197.14 178.38 vadstena 130.53 132.13 128.81 133.65 130.72 134.48 valla 104.00 105.20 107.07 105.19 105.57 104.15

Table 26: Results, times, with settings A-F, q = 6.

57 Name G H I J linghem 11057 11419 11705 11563 ljungsbro 22894 23396 24238 23203 malmslatt 10755 10776 10886 11350 mjolby 21290 23731 22178 23658 ryd-m 16482 17061 16947 16797 ryd 12557 12341 11935 12347 ryd2 14251 14638 14184 14681 soderkoping 17756 17389 17872 17539 vadstena 12313 12433 12413 12493 valla 9758 9934 10084 9814

Table 27: Results, obj, with settings G-J, q = 6.

Name G H I J linghem 304.92 317.15 299.70 317.10 ljungsbro 275.34 282.91 273.85 280.31 malmslatt 82.24 85.39 77.31 87.93 mjolby 61.99 66.96 59.40 74.90 ryd-m 671.67 712.64 667.29 697.00 ryd 443.92 493.46 440.63 456.35 ryd2 516.21 500.31 488.72 515.84 soderkoping 173.91 184.41 170.62 202.87 vadstena 129.14 131.32 140.17 130.95 valla 102.21 108.92 98.98 109.23

Table 28: Results, times, with settings G-J, q = 6.

Obj Time Lowb Name Min Max Min Max l0 l1 linghem G 11057 I 11705 I 299.70 H 317.15 8928 10311 ljungsbro C 22485 I 24238 B 273.34 D 303.09 17548 21180 malmslatt A 10296 J 11350 I 77.31 E 92.05 8833 9961 mjolby G 21290 H 23731 A 57.80 J 74.90 15077 18048 ryd-m C 16113 H 17061 I 667.29 C 744.05 13928 15535 ryd I 11935 G 12557 I 440.63 H 493.46 10420 11570 ryd2 C 13697 J 14681 I 488.72 G 516.21 11903 13247 soderkoping E 16820 D 18603 I 170.62 J 202.87 14031 16239 vadstena B 12107 A 12534 C 128.81 I 140.17 10283 11791 valla A 9367 I 10084 I 98.98 J 109.23 7714 8889

Table 29: Results, best and worst objectives and times, settings A-J, q = 6.

58 7.10 Diﬀerent numbers of vehicles

We have solved a couple of problems for diﬀerent numbers of vehicles, from one up to eight, and then added a ﬁxed cost for each vehicle. The results are reported in table 30. We used settings A.

mantorp atvid vadstena Veh Obj Time Obj Time Obj Time 1 10037 0.11 14372 0.20 18356 2.72 2 8161 2.22 11412 3.49 14203 43.08 3 7541 3.23 11026 4.98 12956 55.85 4 7233 4.77 10303 7.47 12475 74.48 5 7660 7.31 10404 13.42 12237 95.96 6 7364 11.61 10798 16.34 12306 139.74 7 7730 15.63 10567 25.07 12248 164.33 8 7812 20.14 10327 28.74 12130 211.49

Table 30: Diﬀerent numbers of vehicles.

The time for clearing obviously decreases when the number of vehicles increase, but the vehicle cost increases. One can then ﬁnd a minimum of the total cost.

For Mantorp we find that the total cost is minimized for 4 vehicles, while 6 also is good. Again we recall that the methods use randomness, so one should not put too much trust in small differences. For Åtvidaberg, 4, 5 or 8 vehicles seem to be best. For Vadstena, which is a somewhat larger city, the cost is minimal for 8 vehicles. (This is just an illustration of how our results can be used to dimension the fleet for a certain area.)

We also see clearly that the solution time increases with the number of vehicles.

8 Conclusions and future work

The weighted k-Chinese/rural postman problem appears when planning urban snow removal. We have developed heuristics for the problem, using Chinese/rural postman problems, metaheuristics, k-means clustering and local search improvements by moving cycles. The methods give fairly good solutions to real life problems. Solutions times are reasonable and allows real life usage of the code. Manual improvements are possible with the graphical software.

By adding a ﬁxed cost for each vehicle, we can compare using diﬀerent numbers of vehicles.

Future research includes improved local search by moving paths, not only cycles, and improved metaheuristics. We will also use the results of this simpliﬁed problem for a more detailed routing.

Acknowledgment This work has beneﬁted from feedback from students taking the course TAOP61 at Linköping University.

59 References

Ahr, D., and Reinelt, G. (2002), “New heuristics and lower bounds for the min-max k-Chinese postman problem”, in: Proceedings of the 10th Annual European Sym- posium on Algorithms, ESA ’02, 64–74 London, UK, UK. Springer-Verlag.

Ahr, D., and Reinelt, G. (2006), “A tabu search algorithm for the min-max k-Chinese postman problem”, Comput. Oper. Res. 33/12, 3403–3422.

Benavent, E., Corberán, A., Plana, I., and Sanchis, J. M. (2009), “Min-max k-vehicles windy rural postman problem”, Networks 54/4, 216–226.

Benavent, E., and Soler, D. (1999), “The directed rural postman problem with turn penalties”, Transportation Science 33, 408–418.

Campbell, J. F., and Langevin, A. (1995), “The snow disposal assignment problem”, The Journal of the Operational Research Society 46, 919–929.

Christoﬁdes, N., Campos, V., Corberán, A., and Mota, E. (1986), “An algorithm for the rural postman problem on a directed graph”, Mathematical Programming Study 26, 155–166.

Cook, C., Schoenfeld, D. A., and Wainwright, R. L. (1998), “Finding rural postman tours”, in: Proceedings of the 1998 ACM Symposium on Applied Computing, 318–326. ACM Press.

Corberan, A., Marti, R., Martinez, E., and Soler, D. (2002), “The rural postman problem on mixed graphs with turn penalties”, Computers and Operations Research 29, 887–903.

Corberán, A., and Prins, C. (2010), “Recent results on arc routing problems: An anno- tated bibliography”, Networks 56, 50–69.

Eiselt, H. A., Gendreau, M., and Laporte, G. (1995a), “Arc routing problems, part I: The Chinese postman problem”, Operations Research 43, 231–242.

Eiselt, H. A., Gendreau, M., and Laporte, G. (1995b), “Arc routing problems, part II: The rural postman problem”, Operations Research 43, 399–414.

Fu, L., Trudel, M., and Kim, V. (2009), “Optimizing winter road maintenance operations under real-time information”, European Journal of Operational Research 196, 332– 341.

Ghiani, G., Laganá, D., and Musmanno, R. (2006), “A constructive heuristic for the undirected rural postman problem”, Computers & Operations Research 33, 3450– 3457.

Holmberg, K. (2010), “Heuristics for the rural postman problem”, Computers & Opera- tions Research 37, 981–990.

Holmberg, K. (2014), “Urban snow removal: Modeling and relaxations”, Research report LiTH-MAT-R–2014/08–SE, Department of Mathematics, Linköping Institute of Technology, Sweden.

60 Holmberg, K. (2015a), “Map matching by optimization”, Research report LiTH-MAT- R–2015/01–SE, Linköping Intitute of Technology, Sweden.

Holmberg, K. (2015b), “On using open street maps for optimization”, Research report, Linköping Intitute of Technology.

Islam, S. (2010), “Application framework for snow removal routing problem”, Master’s thesis, Linköping University. LIU-IDA/LITH-EX-A–10/012–SE.

Kramberger, T., and Zerovnik, J. (2007), “Priority constrained Chinese postman problem”, Logistics and sustainable transport 1, 33–46.

Labelle, A., Langevin, A., and Campbell, J. (2002), “Sector design for snow removal and disposal in urban areas”, Socio-Economic Planning Sciences 36, 183–202.

Lloyd, S. (1982), “Least squares quantization in PCM”, Information Theory, IEEE Transactions on 28/2, 129–137.

Osterhues, A., and Mariak, F. (2005), “On variants of the k-Chinese postman problem”, Diskussionsbeiträge des fachgebiets operations research und wirtschaftsin- formatik 30, Universität Dortmund.

Pearn, W. L. (1994), “Solvable cases of the k-person Chinese postman problem”, Oper- ations Research Letters 16, 241–244.

Perrier, N., Langevin, A., and Amaya, C.-A. (2008), “Vehicle routing for urban snow plowing operations”, Transportation Science 42, 44–56.

Perrier, N., Langevin, A., and Campbell, J. F. (2006a), “A survey of models and algorithms for winter road maintenance: Part I: System design for spreading and plowing”, Computers and Operations Research 33, 209–238.

Perrier, N., Langevin, A., and Campbell, J. F. (2006b), “A survey of models and algorithms for winter road maintenance: Part II: System design for snow disposal”, Computers and Operations Research 33, 239–262.

Perrier, N., Langevin, A., and Campbell, J. F. (2007a), “A survey of models and algorithms for winter road maintenance: Part III: Vehicle routing and depot location for spreading”, Computers and Operations Research 34, 211 – 257.

Perrier, N., Langevin, A., and Campbell, J. F. (2007b), “A survey of models and algorithms for winter road maintenance: Part IV: Vehicle routing and ﬂeet sizing for plowing and snow disposal”, Computers and Operations Research 34, 258 – 294.

Razmara, G. (2004), Snow Removal Routing Problems - Theory and Applications. PhD dissertation, Linköping University, Sweden. Linköping Studies in Science and Technology. Dissertation no. 888.

Trudel, M. R. (2005), “A mathematical model for winter maintenance operations man- agement”, Master’s thesis, University of Waterloo.

61 Appendix

8.1 Preliminary tests

In tables 31 and 32, we present preliminary tests, made in order to ﬁnd interesting values of the parameters. The best results are shown ordered in 33, and a group of faster settings in in 34.

MD RS SV MV IV ChV CyV CI TC II ST RC Obj time 0 0 1 0 0 0 0 0 0 7 1000 0.3 7929 1.55 0 0 1 0 0 0 0 0 0 7 1000 0.3 7929 1.81 0 0 1 0 0 0 0 0 0 7 1000 0.3 7380 2.49 0 0 1 0 0 0 0 0 0 7 1000 0.3 7947 2.56 1 0 1 0 0 0 0 0 0 7 1000 0.3 6856 2.46 1 1 1 0 0 0 0 0 0 7 1000 0.3 6605 2.57 1 1 1 0 0 0 0 1 0 7 1000 0.3 7318 2.66 1 1 1 0 0 0 0 0 1 7 1000 0.3 6985 2.66 1 1 1 1 0 0 0 0 0 7 1000 0.3 6869 2.67 1 1 1 3 0 0 0 0 0 7 1000 0.3 6897 2.72 1 1 1 4 0 0 0 0 0 7 1000 0.3 7139 2.76 1 1 1 5 0 0 0 0 0 7 1000 0.3 7359 2.83 1 1 1 6 0 0 0 0 0 7 1000 0.3 6821 2.85 1 1 1 7 0 0 0 0 0 7 1000 0.3 7004 2.82 1 1 1 3 0 0 0 0 0 7 1000 0.3 6806 2.76 1 1 1 3 0 1 0 0 0 7 1000 0.3 6471 2.89 1 1 1 3 0 2 0 0 0 7 1000 0.3 6450 2.92 1 1 1 3 0 3 0 0 0 7 1000 0.3 6715 2.82 1 1 1 3 0 4 0 0 0 7 1000 0.3 7096 2.98 1 1 1 3 0 5 0 0 0 7 1000 0.3 6435 2.94 1 1 1 3 1 5 0 0 0 7 1000 0.3 5624 3.15 1 1 1 3 2 5 0 0 0 7 1000 0.3 5784 3.08 1 1 1 3 3 5 0 0 0 7 1000 0.3 5739 3.86 1 1 1 3 3 5 1 0 0 7 1000 0.3 5743 3.78 1 1 1 3 3 5 2 0 0 7 1000 0.3 5682 3.74 1 1 1 3 3 5 3 0 0 7 1000 0.3 5647 3.83 1 1 1 3 3 5 4 0 0 7 1000 0.3 5853 4.26 1 1 3 3 3 5 4 0 0 7 1000 0.3 5815 3.80 1 1 4 3 3 5 4 0 0 7 1000 0.3 5110 3.87

Table 31: Preliminary tests part 1 on instance mantorp-c, with n = 122, m = 145, q = 4, MI=30, MA=0, MO=0, FC=0, UF=0.

62 MD RS SV MV IV ChV CyV CI TC II ST RC Obj time 1 1 4 3 3 5 4 0 0 7 1000 0.3 5127 3.92 1 1 5 3 3 5 4 0 0 7 1000 0.3 5031 3.91 1 1 6 3 3 5 4 0 0 7 1000 0.3 5029 3.87 1 1 7 3 3 5 4 0 0 7 1000 0.3 4969 3.91 1 1 10 3 3 5 4 0 0 7 1000 0.3 4780 3.88 1 1 13 3 3 5 4 0 0 7 1000 0.3 4389 3.94 1 1 14 3 3 5 4 0 0 7 1000 0.3 4427 3.91 1 1 15 3 3 5 4 0 0 7 1000 0.3 4471 3.83 1 1 16 3 3 5 4 0 0 7 1000 0.3 4486 3.88 0 0 13 3 3 5 4 0 0 7 1000 0.3 4722 3.72 1 1 13 3 3 5 4 0 0 7 1000 0.3 4271 3.89 1 1 13 3 3 5 4 0 1 7 1000 0.3 4585 4.15 1 0 13 3 3 5 4 0 0 7 1000 0.3 4606 3.84 1 1 13 3 3 5 4 0 0 7 1000 0.3 4545 4.12 1 1 13 3 3 5 4 0 0 7 2000 0.3 4287 3.95 1 1 13 3 3 5 4 0 0 7 3000 0.3 4397 3.94 1 1 13 3 3 5 4 0 0 5 1000 0.3 4613 4.14 1 1 13 3 3 5 4 0 0 10 1000 0.3 4547 4.05 1 1 13 3 3 5 4 0 0 10 1000 0.5 4316 3.97 1 1 15 3 3 5 4 0 0 10 1000 0.5 4435 4.06 1 1 15 3 3 5 4 0 0 7 1000 0.2 4453 4.07 1 1 15 3 3 5 4 0 0 7 1000 0.2 4670 4.07

Table 32: Preliminary tests part 2 on instance mantorp-c, with n = 122, m = 145, q = 4, MI=30, MA=0, MO=0. FC=0, UF=0.

63 MD RS SV MV IV ChV CyV CI TC II ST RC Obj time 1 1 13 3 3 5 4 0 0 7 1000 0.3 4271 3.89 1 1 13 3 3 5 4 0 0 7 2000 0.3 4287 3.95 1 1 13 3 3 5 4 0 0 10 1000 0.5 4316 3.97 1 1 13 3 3 5 4 0 0 7 1000 0.3 4389 3.94 1 1 13 3 3 5 4 0 0 7 3000 0.3 4397 3.94 1 1 14 3 3 5 4 0 0 7 1000 0.3 4427 3.91 1 1 15 3 3 5 4 0 0 10 1000 0.5 4435 4.06 1 1 15 3 3 5 4 0 0 7 1000 0.2 4453 4.07 1 1 15 3 3 5 4 0 0 7 1000 0.3 4471 3.83 1 1 16 3 3 5 4 0 0 7 1000 0.3 4486 3.88 1 1 13 3 3 5 4 0 0 7 1000 0.3 4545 4.12 1 1 13 3 3 5 4 0 0 10 1000 0.3 4547 4.05 1 1 13 3 3 5 4 0 1 7 1000 0.3 4585 4.15 1 0 13 3 3 5 4 0 0 7 1000 0.3 4606 3.84 1 1 13 3 3 5 4 0 0 5 1000 0.3 4613 4.14 1 1 7 3 3 5 4 0 0 7 1000 0.3 4969 3.91 1 1 15 3 3 5 4 0 0 7 1000 0.2 4670 4.07 0 0 13 3 3 5 4 0 0 7 1000 0.3 4722 3.72 1 1 10 3 3 5 4 0 0 7 1000 0.3 4780 3.88 1 1 6 3 3 5 4 0 0 7 1000 0.3 5029 3.87 1 1 5 3 3 5 4 0 0 7 1000 0.3 5031 3.91 1 1 4 3 3 5 4 0 0 7 1000 0.3 5110 3.87 1 1 4 3 3 5 4 0 0 7 1000 0.3 5127 3.92 1 1 1 3 1 5 0 0 0 7 1000 0.3 5624 3.15 1 1 1 3 3 5 3 0 0 7 1000 0.3 5647 3.83 1 1 1 3 3 5 2 0 0 7 1000 0.3 5682 3.74 1 1 1 3 3 5 0 0 0 7 1000 0.3 5739 3.86 1 1 1 3 3 5 1 0 0 7 1000 0.3 5743 3.78 1 1 1 3 2 5 0 0 0 7 1000 0.3 5784 3.08 1 1 3 3 3 5 4 0 0 7 1000 0.3 5815 3.80 1 1 1 3 3 5 4 0 0 7 1000 0.3 5853 4.26

Table 33: Best: Preliminary tests on instance mantorp-c, sorted after objective function value.

64 MD RS SV MV IV ChV CyV CI TC II ST RC Obj time 1 1 1 3 0 5 0 0 0 7 1000 0.3 6435 2.94 1 1 1 3 0 2 0 0 0 7 1000 0.3 6450 2.92 1 1 1 3 0 1 0 0 0 7 1000 0.3 6471 2.89 1 1 1 0 0 0 0 0 0 7 1000 0.3 6605 2.57 1 1 1 3 0 3 0 0 0 7 1000 0.3 6715 2.82 1 1 1 3 0 0 0 0 0 7 1000 0.3 6806 2.76 1 1 1 6 0 0 0 0 0 7 1000 0.3 6821 2.85 1 0 1 0 0 0 0 0 0 7 1000 0.3 6856 2.46 1 1 1 1 0 0 0 0 0 7 1000 0.3 6869 2.67 1 1 1 3 0 0 0 0 0 7 1000 0.3 6897 2.72 1 1 1 0 0 0 0 0 1 7 1000 0.3 6985 2.66

Table 34: Faster: Preliminary tests on instance mantorp-c, sorted after objective function value.

65 8.2 Tests with one iteration

In the tables 35 - 46, we give the results for the settings using only one iteration.

Name n m q Obj Time atvid-s 68 95 4 4292.99 0.06 brokind 179 199 4 5761.22 0.20 colonia 24 36 4 1021.93 0.04 ekangen 149 184 4 5077.26 0.15 mantorp-c 122 145 4 4511.35 0.11 mantorp 185 226 4 7598.26 0.22 rimforsa-o 94 109 4 3485.69 0.08 skanninge-c 122 152 4 4492.03 0.12 skanninge-n 62 80 4 2433.58 0.07 skanninge 127 158 4 4712.01 0.11 studentryd-e 76 104 4 1690.63 0.07 sturefors 183 212 4 6477.29 0.21 valla-liu 160 213 4 2718.51 0.18 vikingstad 178 213 4 7180.76 0.22 atvid 220 290 4 10381.53 0.33 borensberg 211 257 4 11423.44 0.30 kisa 311 359 4 9133.56 0.87 liu 291 394 4 5161.28 0.63 rimforsa 233 262 4 6298.58 0.36 ryd1 382 508 4 4110.91 1.47 studentryd-int 235 316 4 3083.82 0.43 studentryd 387 519 4 5439.50 1.28 linghem 836 1129 4 11731.25 8.74 ljungsbro 812 1024 4 23866.51 8.24 malmslatt 489 626 4 10867.23 2.16 mjolby 419 524 4 21177.15 1.38 ryd-m 1107 1495 4 17561.68 23.15 ryd 952 1291 4 12612.44 15.45 ryd2 987 1337 4 14359.72 16.52 soderkoping 666 910 4 18017.16 5.03 vadstena 581 778 4 12548.05 4.01 valla 527 707 4 11094.55 2.65

Table 35: Results with setting A1.

66 Name n m q Obj Time atvid-s 68 95 4 4202.33 0.06 brokind 179 199 4 5766.56 0.20 colonia 24 36 4 1000.29 0.04 ekangen 149 184 4 5712.64 0.20 mantorp-c 122 145 4 4556.09 0.11 mantorp 185 226 4 7869.89 0.24 rimforsa-o 94 109 4 3427.18 0.08 skanninge-c 122 152 4 4673.22 0.13 skanninge-n 62 80 4 2446.20 0.06 skanninge 127 158 4 4810.49 0.11 studentryd-e 76 104 4 1628.94 0.07 sturefors 183 212 4 6633.19 0.27 valla-liu 160 213 4 2676.04 0.17 vikingstad 178 213 4 7178.84 0.21 atvid 220 290 4 10501.32 0.33 borensberg 211 257 4 9851.13 0.28 kisa 311 359 4 9454.17 0.66 liu 291 394 4 5568.01 0.60 rimforsa 233 262 4 6010.63 0.35 ryd1 382 508 4 4085.17 1.45 studentryd-int 235 316 4 2995.91 0.43 studentryd 387 519 4 5396.46 1.37 linghem 836 1129 4 11219.91 9.05 ljungsbro 812 1024 4 24235.39 8.25 malmslatt 489 626 4 10712.92 2.30 mjolby 419 524 4 21117.80 1.54 ryd-m 1107 1495 4 17558.14 23.18 ryd 952 1291 4 12889.70 14.66 ryd2 987 1337 4 14400.13 16.61 soderkoping 666 910 4 18440.67 5.07 vadstena 581 778 4 12603.38 3.49 valla 527 707 4 10217.07 2.50

Table 36: Results with setting B1.

67 Name n m q Obj Time atvid-s 68 95 4 4098.29 0.06 brokind 179 199 4 5865.04 0.20 colonia 24 36 4 1012.14 0.03 ekangen 149 184 4 5653.69 0.16 mantorp-c 122 145 4 4558.87 0.10 mantorp 185 226 4 7974.29 0.25 rimforsa-o 94 109 4 3388.56 0.08 skanninge-c 122 152 4 4649.71 0.11 skanninge-n 62 80 4 2451.72 0.05 skanninge 127 158 4 4486.88 0.12 studentryd-e 76 104 4 1639.85 0.06 sturefors 183 212 4 6449.09 0.21 valla-liu 160 213 4 2607.13 0.17 vikingstad 178 213 4 7456.45 0.21 atvid 220 290 4 11209.77 0.34 borensberg 211 257 4 11940.52 0.30 kisa 311 359 4 9650.40 0.75 liu 291 394 4 5429.01 0.63 rimforsa 233 262 4 6016.46 0.42 ryd1 382 508 4 4352.93 1.19 studentryd-int 235 316 4 2991.75 0.39 studentryd 387 519 4 5545.73 1.33 linghem 836 1129 4 11997.47 8.96 ljungsbro 812 1024 4 23864.45 8.24 malmslatt 489 626 4 10529.46 2.11 mjolby 419 524 4 21225.74 1.44 ryd-m 1107 1495 4 16718.41 23.34 ryd 952 1291 4 12341.18 14.58 ryd2 987 1337 4 14948.71 16.62 soderkoping 666 910 4 18950.13 4.69 vadstena 581 778 4 12799.22 3.68 valla 527 707 4 9928.93 2.52

Table 37: Results with setting C1.

68 Name n m q Obj Time atvid-s 68 95 4 4043.58 0.07 brokind 179 199 4 5833.49 0.20 colonia 24 36 4 1001.31 0.03 ekangen 149 184 4 5736.50 0.15 mantorp-c 122 145 4 4757.84 0.11 mantorp 185 226 4 7502.80 0.21 rimforsa-o 94 109 4 3546.42 0.08 skanninge-c 122 152 4 4642.49 0.11 skanninge-n 62 80 4 2438.07 0.05 skanninge 127 158 4 4772.83 0.11 studentryd-e 76 104 4 1645.31 0.06 sturefors 183 212 4 6576.28 0.22 valla-liu 160 213 4 2686.43 0.21 vikingstad 178 213 4 7647.02 0.22 atvid 220 290 4 10581.15 0.32 borensberg 211 257 4 10360.53 0.28 kisa 311 359 4 9058.64 0.70 liu 291 394 4 5629.90 0.60 rimforsa 233 262 4 6099.56 0.35 ryd1 382 508 4 4094.71 1.30 studentryd-int 235 316 4 3034.94 0.38 studentryd 387 519 4 5446.24 1.48 linghem 836 1129 4 11423.42 8.80 ljungsbro 812 1024 4 24506.00 8.28 malmslatt 489 626 4 10973.03 2.23 mjolby 419 524 4 23912.77 1.53 ryd-m 1107 1495 4 17158.90 23.28 ryd 952 1291 4 12576.72 14.80 ryd2 987 1337 4 14250.46 16.61 soderkoping 666 910 4 18818.98 4.61 vadstena 581 778 4 12504.13 3.49 valla 527 707 4 9733.04 2.56

Table 38: Results with setting D1.

69 Name n m q Obj Time atvid-s 68 95 4 4002.45 0.06 brokind 179 199 4 6060.45 0.20 colonia 24 36 4 992.83 0.04 ekangen 149 184 4 5567.42 0.15 mantorp-c 122 145 4 4576.28 0.11 mantorp 185 226 4 7685.16 0.22 rimforsa-o 94 109 4 3373.28 0.07 skanninge-c 122 152 4 4668.17 0.10 skanninge-n 62 80 4 2362.93 0.05 skanninge 127 158 4 4610.98 0.11 studentryd-e 76 104 4 1588.86 0.06 sturefors 183 212 4 6449.91 0.21 valla-liu 160 213 4 2864.05 0.17 vikingstad 178 213 4 7723.36 0.22 atvid 220 290 4 10727.33 0.34 borensberg 211 257 4 10761.30 0.34 kisa 311 359 4 9089.64 0.71 liu 291 394 4 5338.63 0.61 rimforsa 233 262 4 5880.67 0.37 ryd1 382 508 4 4073.47 1.28 studentryd-int 235 316 4 3024.74 0.41 studentryd 387 519 4 5498.14 1.36 linghem 836 1129 4 11259.86 9.32 ljungsbro 812 1024 4 23571.75 8.39 malmslatt 489 626 4 10538.78 2.37 mjolby 419 524 4 21185.18 1.54 ryd-m 1107 1495 4 17732.63 22.95 ryd 952 1291 4 12273.27 15.78 ryd2 987 1337 4 14280.31 16.17 soderkoping 666 910 4 17928.95 4.80 vadstena 581 778 4 12531.06 3.74 valla 527 707 4 9709.64 2.55

Table 39: Results with setting E1.

70 Name n m q Obj Time atvid-s 68 95 4 4097.01 0.07 brokind 179 199 4 6088.93 0.19 colonia 24 36 4 1014.46 0.04 ekangen 149 184 4 4999.19 0.15 mantorp-c 122 145 4 4611.48 0.10 mantorp 185 226 4 7685.53 0.23 rimforsa-o 94 109 4 3469.93 0.08 skanninge-c 122 152 4 4688.63 0.12 skanninge-n 62 80 4 2458.85 0.06 skanninge 127 158 4 4834.29 0.12 studentryd-e 76 104 4 1649.97 0.07 sturefors 183 212 4 6495.91 0.21 valla-liu 160 213 4 2676.04 0.17 vikingstad 178 213 4 8918.12 0.21 atvid 220 290 4 10622.87 0.32 borensberg 211 257 4 10389.60 0.30 kisa 311 359 4 8720.10 0.70 liu 291 394 4 5205.85 0.63 rimforsa 233 262 4 6028.30 0.36 ryd1 382 508 4 4380.72 0.94 studentryd-int 235 316 4 3126.50 0.39 studentryd 387 519 4 5438.38 1.41 linghem 836 1129 4 11320.12 8.61 ljungsbro 812 1024 4 25426.39 8.16 malmslatt 489 626 4 11030.91 2.11 mjolby 419 524 4 21345.96 1.46 ryd-m 1107 1495 4 17242.58 23.29 ryd 952 1291 4 12477.68 14.87 ryd2 987 1337 4 15438.62 16.51 soderkoping 666 910 4 18680.37 5.05 vadstena 581 778 4 12843.79 3.65 valla 527 707 4 9713.33 2.78

Table 40: Results with setting F1.

71 Name n m q Obj Time atvid-s 68 95 4 4029.70 0.05 brokind 179 199 4 5647.97 0.25 colonia 24 36 4 1007.38 0.04 ekangen 149 184 4 5170.89 0.15 mantorp-c 122 145 4 4692.08 0.11 mantorp 185 226 4 7925.07 0.24 rimforsa-o 94 109 4 3466.83 0.09 skanninge-c 122 152 4 4540.03 0.11 skanninge-n 62 80 4 2438.55 0.05 skanninge 127 158 4 4764.04 0.12 studentryd-e 76 104 4 1694.05 0.06 sturefors 183 212 4 6466.94 0.29 valla-liu 160 213 4 2754.42 0.17 vikingstad 178 213 4 7768.96 0.22 atvid 220 290 4 10355.88 0.32 borensberg 211 257 4 10358.61 0.28 kisa 311 359 4 9491.25 0.66 liu 291 394 4 5184.55 0.72 rimforsa 233 262 4 5842.28 0.32 ryd1 382 508 4 4147.45 1.17 studentryd-int 235 316 4 3085.42 0.36 studentryd 387 519 4 5254.86 1.30 linghem 836 1129 4 11764.30 8.93 ljungsbro 812 1024 4 24840.08 8.20 malmslatt 489 626 4 11030.40 2.18 mjolby 419 524 4 21076.14 1.37 ryd-m 1107 1495 4 17744.78 22.87 ryd 952 1291 4 12833.49 14.17 ryd2 987 1337 4 14952.63 16.25 soderkoping 666 910 4 19067.27 4.54 vadstena 581 778 4 12543.08 3.57 valla 527 707 4 10056.64 2.45

Table 41: Results with setting G1.

72 Name n m q Obj Time atvid-s 68 95 4 4020.54 0.06 brokind 179 199 4 6567.65 0.21 colonia 24 36 4 1011.12 0.03 ekangen 149 184 4 5202.10 0.15 mantorp-c 122 145 4 5018.77 0.11 mantorp 185 226 4 8095.74 0.24 rimforsa-o 94 109 4 3384.90 0.07 skanninge-c 122 152 4 4872.53 0.10 skanninge-n 62 80 4 2475.10 0.06 skanninge 127 158 4 4898.85 0.11 studentryd-e 76 104 4 1641.56 0.07 sturefors 183 212 4 6908.30 0.22 valla-liu 160 213 4 2949.16 0.19 vikingstad 178 213 4 7915.93 0.22 atvid 220 290 4 10650.21 0.34 borensberg 211 257 4 10137.26 0.30 kisa 311 359 4 9550.69 0.73 liu 291 394 4 5274.02 0.61 rimforsa 233 262 4 6747.74 0.39 ryd1 382 508 4 4050.60 1.36 studentryd-int 235 316 4 3158.73 0.40 studentryd 387 519 4 5425.97 1.33 linghem 836 1129 4 11212.04 8.97 ljungsbro 812 1024 4 23960.92 8.31 malmslatt 489 626 4 11358.14 2.21 mjolby 419 524 4 23044.18 1.43 ryd-m 1107 1495 4 17340.90 22.93 ryd 952 1291 4 12380.76 14.86 ryd2 987 1337 4 14774.93 16.44 soderkoping 666 910 4 18081.82 5.02 vadstena 581 778 4 12565.53 3.48 valla 527 707 4 9764.37 2.55

Table 42: Results with setting H1.

73 Name n m q Obj Time atvid-s 68 95 4 4084.17 0.06 brokind 179 199 4 5962.98 0.17 colonia 24 36 4 1002.25 0.03 ekangen 149 184 4 5141.96 0.12 mantorp-c 122 145 4 4642.37 0.09 mantorp 185 226 4 7573.32 0.18 rimforsa-o 94 109 4 3362.60 0.06 skanninge-c 122 152 4 4863.95 0.08 skanninge-n 62 80 4 2476.07 0.05 skanninge 127 158 4 4484.35 0.10 studentryd-e 76 104 4 1666.74 0.06 sturefors 183 212 4 6629.18 0.18 valla-liu 160 213 4 2835.91 0.15 vikingstad 178 213 4 7512.16 0.18 atvid 220 290 4 10620.93 0.30 borensberg 211 257 4 10922.50 0.25 kisa 311 359 4 8942.81 0.71 liu 291 394 4 5516.80 0.54 rimforsa 233 262 4 5898.01 0.29 ryd1 382 508 4 4106.73 1.28 studentryd-int 235 316 4 3172.94 0.31 studentryd 387 519 4 5388.61 1.33 linghem 836 1129 4 11564.77 8.17 ljungsbro 812 1024 4 23937.38 7.58 malmslatt 489 626 4 11080.24 2.01 mjolby 419 524 4 21775.31 1.20 ryd-m 1107 1495 4 17491.48 21.75 ryd 952 1291 4 12316.05 13.78 ryd2 987 1337 4 14293.34 15.51 soderkoping 666 910 4 18624.45 4.38 vadstena 581 778 4 12795.78 3.25 valla 527 707 4 9982.81 2.27

Table 43: Results with setting I1.

74 Name n m q Obj Time atvid-s 68 95 4 4069.84 0.06 brokind 179 199 4 6531.53 0.21 colonia 24 36 4 1024.96 0.04 ekangen 149 184 4 5513.36 0.15 mantorp-c 122 145 4 4906.11 0.11 mantorp 185 226 4 7933.01 0.22 rimforsa-o 94 109 4 3380.73 0.08 skanninge-c 122 152 4 4887.07 0.11 skanninge-n 62 80 4 2501.53 0.05 skanninge 127 158 4 4930.13 0.11 studentryd-e 76 104 4 1619.24 0.06 sturefors 183 212 4 6958.30 0.22 valla-liu 160 213 4 2983.92 0.19 vikingstad 178 213 4 8250.84 0.25 atvid 220 290 4 11035.99 0.36 borensberg 211 257 4 10154.99 0.31 kisa 311 359 4 9530.69 0.76 liu 291 394 4 5298.61 0.66 rimforsa 233 262 4 6612.25 0.37 ryd1 382 508 4 4037.38 1.27 studentryd-int 235 316 4 3151.55 0.37 studentryd 387 519 4 5419.07 1.33 linghem 836 1129 4 11178.15 8.67 ljungsbro 812 1024 4 24001.83 8.44 malmslatt 489 626 4 11373.68 2.16 mjolby 419 524 4 23045.48 1.48 ryd-m 1107 1495 4 17201.40 23.07 ryd 952 1291 4 12399.93 14.90 ryd2 987 1337 4 14984.91 16.52 soderkoping 666 910 4 17976.33 4.84 vadstena 581 778 4 12570.30 3.82 valla 527 707 4 9797.44 2.55

Table 44: Results with setting J1.

75 Name n m q Obj Time atvid-s 68 95 4 4876.64 0.04 brokind 179 199 4 9954.53 0.18 colonia 24 36 4 1071.43 0.03 ekangen 149 184 4 7998.76 0.12 mantorp-c 122 145 4 7403.59 0.08 mantorp 185 226 4 11861.96 0.23 rimforsa-o 94 109 4 4972.86 0.06 skanninge-c 122 152 4 6535.95 0.09 skanninge-n 62 80 4 3083.59 0.04 skanninge 127 158 4 6782.39 0.09 studentryd-e 76 104 4 1936.80 0.05 sturefors 183 212 4 10338.49 0.18 valla-liu 160 213 4 3600.53 0.14 vikingstad 178 213 4 11215.08 0.18 atvid 220 290 4 15662.36 0.33 borensberg 211 257 4 13794.94 0.27 kisa 311 359 4 14403.28 0.72 liu 291 394 4 6838.78 0.67 rimforsa 233 262 4 9653.41 0.32 ryd1 382 508 4 5457.20 1.49 studentryd-int 235 316 4 4118.59 0.39 studentryd 387 519 4 7301.20 1.63 linghem 836 1129 4 16745.74 16.14 ljungsbro 812 1024 4 36856.23 13.52 malmslatt 489 626 4 16443.23 2.79 mjolby 419 524 4 31275.01 1.73 ryd-m 1107 1495 4 23525.59 45.44 ryd 952 1291 4 17273.85 28.06 ryd2 987 1337 4 19164.47 30.60 soderkoping 666 910 4 24275.48 7.79 vadstena 581 778 4 17762.63 4.89 valla 527 707 4 13165.38 3.45

Table 45: Results with setting K1.

76 Name n m q Obj Time atvid-s 68 95 4 5263.07 0.05 brokind 179 199 4 9803.74 0.19 colonia 24 36 4 1065.17 0.03 ekangen 149 184 4 8099.05 0.13 mantorp-c 122 145 4 7247.68 0.08 mantorp 185 226 4 12261.55 0.18 rimforsa-o 94 109 4 5133.07 0.06 skanninge-c 122 152 4 6372.11 0.08 skanninge-n 62 80 4 3366.16 0.04 skanninge 127 158 4 6875.31 0.09 studentryd-e 76 104 4 1873.41 0.05 sturefors 183 212 4 9429.48 0.18 valla-liu 160 213 4 3426.59 0.13 vikingstad 178 213 4 11352.93 0.18 atvid 220 290 4 14978.87 0.30 borensberg 211 257 4 13834.87 0.25 kisa 311 359 4 14554.64 0.73 liu 291 394 4 6845.29 0.72 rimforsa 233 262 4 9807.20 0.37 ryd1 382 508 4 5554.93 1.46 studentryd-int 235 316 4 4206.16 0.59 studentryd 387 519 4 7021.03 1.42 linghem 836 1129 4 16872.48 15.91 ljungsbro 812 1024 4 36365.63 14.09 malmslatt 489 626 4 15939.10 2.86 mjolby 419 524 4 32309.30 1.56 ryd-m 1107 1495 4 23316.01 46.06 ryd 952 1291 4 17174.22 27.65 ryd2 987 1337 4 19687.52 31.26 soderkoping 666 910 4 24289.07 7.28 vadstena 581 778 4 17801.28 5.24 valla 527 707 4 13054.52 3.37

Table 46: Results with setting L1.

77 8.3 Tests with 30 iterations

In tables 47 - 62, we present results of runs with 30 iterations.

Name n m q Obj Time atvid-s 68 95 4 4014.35 1.35 brokind 179 199 4 5916.24 4.45 colonia 24 36 4 1001.01 0.77 ekangen 149 184 4 5607.90 3.56 mantorp-c 122 145 4 4557.82 2.43 mantorp 185 226 4 7353.81 4.83 rimforsa-o 94 109 4 3412.41 1.77 skanninge-c 122 152 4 4375.77 2.43 skanninge-n 62 80 4 2373.09 1.18 skanninge 127 158 4 4470.06 2.65 studentryd-e 76 104 4 1608.00 1.46 sturefors 183 212 4 6988.09 4.85 valla-liu 160 213 4 2692.72 3.86 vikingstad 178 213 4 7574.55 4.84 atvid 220 290 4 10444.18 8.05 borensberg 211 257 4 9451.83 6.49 kisa 311 359 4 9146.86 18.09 liu 291 394 4 5236.79 16.81 rimforsa 233 262 4 5875.31 9.05 ryd1 382 508 4 4094.66 30.48 studentryd-int 235 316 4 3102.48 10.03 studentryd 387 519 4 5242.03 30.41 linghem 836 1129 4 11955.03 182.20 ljungsbro 812 1024 4 23290.93 161.62 malmslatt 489 626 4 10658.54 49.55 mjolby 419 524 4 21041.42 34.11 ryd-m 1107 1495 4 16994.80 420.20 ryd 952 1291 4 12137.43 277.52 ryd2 987 1337 4 14876.40 298.73 soderkoping 666 910 4 17566.32 95.51 vadstena 581 778 4 12416.79 69.41 valla 527 707 4 9897.03 58.29

Table 47: Results with setting A30.

78 Name n m q Obj Time atvid-s 68 95 4 3842.55 1.42 brokind 179 199 4 5761.55 4.40 colonia 24 36 4 1000.49 0.80 ekangen 149 184 4 5299.37 3.72 mantorp-c 122 145 4 4662.50 2.69 mantorp 185 226 4 7520.47 5.52 rimforsa-o 94 109 4 3337.12 1.88 skanninge-c 122 152 4 4484.26 2.76 skanninge-n 62 80 4 2370.41 1.24 skanninge 127 158 4 4628.87 2.96 studentryd-e 76 104 4 1604.67 2.27 sturefors 183 212 4 6686.49 5.10 valla-liu 160 213 4 2665.54 3.87 vikingstad 178 213 4 7234.23 5.69 atvid 220 290 4 10489.46 8.92 borensberg 211 257 4 9481.08 7.11 kisa 311 359 4 9519.21 15.88 liu 291 394 4 5320.16 14.71 rimforsa 233 262 4 5724.08 8.37 ryd1 382 508 4 4028.04 25.58 studentryd-int 235 316 4 3113.80 8.93 studentryd 387 519 4 5304.90 29.56 linghem 836 1129 4 11336.45 182.50 ljungsbro 812 1024 4 23795.30 165.63 malmslatt 489 626 4 10476.14 46.92 mjolby 419 524 4 21524.04 38.65 ryd-m 1107 1495 4 16920.26 422.05 ryd 952 1291 4 12240.06 270.36 ryd2 987 1337 4 14184.12 300.00 soderkoping 666 910 4 18020.13 98.58 vadstena 581 778 4 12582.37 74.77 valla 527 707 4 9721.53 58.94

Table 48: Results with setting B30.

79 Name n m q Obj Time atvid-s 68 95 4 3859.21 1.31 brokind 179 199 4 5810.22 4.52 colonia 24 36 4 1000.97 0.85 ekangen 149 184 4 5595.41 3.45 mantorp-c 122 145 4 4481.26 2.41 mantorp 185 226 4 7477.57 5.75 rimforsa-o 94 109 4 3373.19 1.76 skanninge-c 122 152 4 4458.42 2.51 skanninge-n 62 80 4 2409.60 1.33 skanninge 127 158 4 4476.60 2.77 studentryd-e 76 104 4 1623.29 1.49 sturefors 183 212 4 6482.83 5.29 valla-liu 160 213 4 2712.90 3.98 vikingstad 178 213 4 7360.37 5.00 atvid 220 290 4 10594.96 8.58 borensberg 211 257 4 10073.31 7.23 kisa 311 359 4 8903.99 17.17 liu 291 394 4 5249.11 18.70 rimforsa 233 262 4 5981.56 9.77 ryd1 382 508 4 3993.24 27.74 studentryd-int 235 316 4 2995.91 9.57 studentryd 387 519 4 5256.90 33.61 linghem 836 1129 4 11668.18 179.77 ljungsbro 812 1024 4 24009.75 173.87 malmslatt 489 626 4 10521.85 50.32 mjolby 419 524 4 21063.56 37.27 ryd-m 1107 1495 4 16971.89 416.06 ryd 952 1291 4 12213.37 275.99 ryd2 987 1337 4 14087.49 297.11 soderkoping 666 910 4 18621.36 105.69 vadstena 581 778 4 12474.87 74.51 valla 527 707 4 9593.06 58.90

Table 49: Results with setting C30.

80 Name n m q Obj Time atvid-s 68 95 4 3921.47 1.36 brokind 179 199 4 5701.96 4.51 colonia 24 36 4 995.01 0.80 ekangen 149 184 4 5214.58 3.24 mantorp-c 122 145 4 4590.56 2.46 mantorp 185 226 4 7650.24 5.18 rimforsa-o 94 109 4 3389.12 1.77 skanninge-c 122 152 4 4545.44 2.57 skanninge-n 62 80 4 2376.28 1.36 skanninge 127 158 4 4521.35 3.02 studentryd-e 76 104 4 1623.73 1.54 sturefors 183 212 4 6418.14 4.55 valla-liu 160 213 4 2737.34 4.06 vikingstad 178 213 4 7534.31 4.71 atvid 220 290 4 10623.02 7.98 borensberg 211 257 4 9894.11 6.48 kisa 311 359 4 8962.49 17.34 liu 291 394 4 5343.10 15.87 rimforsa 233 262 4 5813.88 9.02 ryd1 382 508 4 4038.03 30.46 studentryd-int 235 316 4 3160.14 11.16 studentryd 387 519 4 5243.48 30.25 linghem 836 1129 4 11542.20 179.49 ljungsbro 812 1024 4 24230.80 161.00 malmslatt 489 626 4 10830.51 48.66 mjolby 419 524 4 21446.68 37.98 ryd-m 1107 1495 4 16996.78 417.10 ryd 952 1291 4 12539.28 270.01 ryd2 987 1337 4 14242.56 299.23 soderkoping 666 910 4 17821.89 98.41 vadstena 581 778 4 12492.97 79.45 valla 527 707 4 9730.34 58.57

Table 50: Results with setting D30.

81 Name n m q Obj Time atvid-s 68 95 4 3982.22 1.37 brokind 179 199 4 5810.68 4.77 colonia 24 36 4 1009.46 0.80 ekangen 149 184 4 5327.88 3.41 mantorp-c 122 145 4 4672.97 2.51 mantorp 185 226 4 7523.94 4.88 rimforsa-o 94 109 4 3477.86 1.76 skanninge-c 122 152 4 4676.32 2.55 skanninge-n 62 80 4 2366.70 1.32 skanninge 127 158 4 4562.91 2.53 studentryd-e 76 104 4 1621.75 1.48 sturefors 183 212 4 6744.84 5.43 valla-liu 160 213 4 2845.84 4.15 vikingstad 178 213 4 7547.22 4.60 atvid 220 290 4 10595.47 8.03 borensberg 211 257 4 10257.59 8.54 kisa 311 359 4 9495.02 16.52 liu 291 394 4 5323.21 16.46 rimforsa 233 262 4 5872.40 8.75 ryd1 382 508 4 4050.00 30.17 studentryd-int 235 316 4 3087.24 9.06 studentryd 387 519 4 5345.46 31.41 linghem 836 1129 4 11423.91 190.59 ljungsbro 812 1024 4 24403.26 163.33 malmslatt 489 626 4 10411.43 46.41 mjolby 419 524 4 21229.63 35.47 ryd-m 1107 1495 4 16868.44 419.66 ryd 952 1291 4 12340.27 269.26 ryd2 987 1337 4 13948.04 297.34 soderkoping 666 910 4 17659.23 101.25 vadstena 581 778 4 12530.68 76.18 valla 527 707 4 10022.07 58.02

Table 51: Results with setting E30.

82 Name n m q Obj Time atvid-s 68 95 4 3862.02 1.38 brokind 179 199 4 5748.66 4.49 colonia 24 36 4 1000.34 0.78 ekangen 149 184 4 4995.67 3.16 mantorp-c 122 145 4 4509.49 2.32 mantorp 185 226 4 7777.73 5.06 rimforsa-o 94 109 4 3506.40 1.93 skanninge-c 122 152 4 4647.16 2.74 skanninge-n 62 80 4 2372.18 1.32 skanninge 127 158 4 4450.56 2.71 studentryd-e 76 104 4 1599.65 1.57 sturefors 183 212 4 6466.94 5.36 valla-liu 160 213 4 2669.97 3.98 vikingstad 178 213 4 7365.45 4.88 atvid 220 290 4 10501.47 7.12 borensberg 211 257 4 10190.50 7.26 kisa 311 359 4 9459.82 15.56 liu 291 394 4 5150.66 15.75 rimforsa 233 262 4 5792.07 8.06 ryd1 382 508 4 4046.42 27.94 studentryd-int 235 316 4 3051.13 8.39 studentryd 387 519 4 5336.90 28.50 linghem 836 1129 4 11236.89 190.22 ljungsbro 812 1024 4 23835.25 164.44 malmslatt 489 626 4 10508.90 47.04 mjolby 419 524 4 22548.45 34.24 ryd-m 1107 1495 4 17226.02 416.16 ryd 952 1291 4 12683.28 273.43 ryd2 987 1337 4 14238.05 297.81 soderkoping 666 910 4 19055.82 98.66 vadstena 581 778 4 12343.26 72.93 valla 527 707 4 9691.83 56.87

Table 52: Results with setting F30.

83 Name n m q Obj Time atvid-s 68 95 4 3957.36 1.29 brokind 179 199 4 5859.65 4.31 colonia 24 36 4 989.41 0.77 ekangen 149 184 4 4933.69 3.20 mantorp-c 122 145 4 4473.70 2.33 mantorp 185 226 4 7555.98 4.69 rimforsa-o 94 109 4 3441.06 1.67 skanninge-c 122 152 4 4712.75 2.35 skanninge-n 62 80 4 2390.49 1.38 skanninge 127 158 4 4477.75 2.52 studentryd-e 76 104 4 1621.92 1.52 sturefors 183 212 4 6613.57 4.11 valla-liu 160 213 4 2657.70 3.90 vikingstad 178 213 4 7604.28 4.85 atvid 220 290 4 10443.68 7.59 borensberg 211 257 4 9901.56 6.04 kisa 311 359 4 9213.12 16.43 liu 291 394 4 5168.66 16.36 rimforsa 233 262 4 5894.86 8.00 ryd1 382 508 4 4030.41 27.62 studentryd-int 235 316 4 3058.81 9.32 studentryd 387 519 4 5295.67 29.76 linghem 836 1129 4 11611.55 177.07 ljungsbro 812 1024 4 23584.37 166.58 malmslatt 489 626 4 10593.92 45.78 mjolby 419 524 4 23003.33 32.53 ryd-m 1107 1495 4 17202.70 403.35 ryd 952 1291 4 12367.75 263.03 ryd2 987 1337 4 14575.19 286.28 soderkoping 666 910 4 18220.98 95.05 vadstena 581 778 4 12365.77 67.97 valla 527 707 4 9517.01 55.48

Table 53: Results with setting G30.

84 Name n m q Obj Time atvid-s 68 95 4 3930.96 1.34 brokind 179 199 4 6567.65 4.59 colonia 24 36 4 1005.92 0.77 ekangen 149 184 4 5144.47 3.48 mantorp-c 122 145 4 4904.36 3.68 mantorp 185 226 4 7883.04 5.13 rimforsa-o 94 109 4 3344.20 1.73 skanninge-c 122 152 4 4867.49 2.41 skanninge-n 62 80 4 2415.41 1.21 skanninge 127 158 4 4866.04 2.50 studentryd-e 76 104 4 1628.39 3.08 sturefors 183 212 4 6895.41 4.92 valla-liu 160 213 4 2829.86 4.13 vikingstad 178 213 4 7877.91 5.12 atvid 220 290 4 10600.44 8.45 borensberg 211 257 4 10115.10 7.62 kisa 311 359 4 9498.92 17.27 liu 291 394 4 5274.02 19.08 rimforsa 233 262 4 6592.28 8.35 ryd1 382 508 4 4027.60 29.00 studentryd-int 235 316 4 3158.73 8.94 studentryd 387 519 4 5363.96 29.46 linghem 836 1129 4 11132.73 182.42 ljungsbro 812 1024 4 23960.92 165.24 malmslatt 489 626 4 11306.74 48.02 mjolby 419 524 4 22623.82 34.55 ryd-m 1107 1495 4 17008.88 413.78 ryd 952 1291 4 12380.76 277.56 ryd2 987 1337 4 14774.93 302.82 soderkoping 666 910 4 17840.78 100.43 vadstena 581 778 4 12560.22 70.67 valla 527 707 4 9764.37 54.24

Table 54: Results with setting H30.

85 Name n m q Obj Time atvid-s 68 95 4 4290.40 1.40 brokind 179 199 4 5878.69 4.66 colonia 24 36 4 1006.73 0.79 ekangen 149 184 4 5585.79 3.92 mantorp-c 122 145 4 4561.04 2.43 mantorp 185 226 4 7432.72 4.72 rimforsa-o 94 109 4 3435.82 1.81 skanninge-c 122 152 4 4768.31 3.39 skanninge-n 62 80 4 2458.03 1.33 skanninge 127 158 4 4593.10 2.59 studentryd-e 76 104 4 1619.44 1.69 sturefors 183 212 4 6574.11 5.69 valla-liu 160 213 4 2782.96 4.15 vikingstad 178 213 4 8905.74 5.01 atvid 220 290 4 10785.35 7.77 borensberg 211 257 4 10675.62 7.25 kisa 311 359 4 9807.60 15.78 liu 291 394 4 5622.63 14.78 rimforsa 233 262 4 6020.91 7.89 ryd1 382 508 4 4190.90 24.86 studentryd-int 235 316 4 3025.56 9.20 studentryd 387 519 4 5493.14 24.93 linghem 836 1129 4 11748.01 171.83 ljungsbro 812 1024 4 24786.03 154.99 malmslatt 489 626 4 11031.11 41.85 mjolby 419 524 4 21558.72 30.22 ryd-m 1107 1495 4 17715.57 406.10 ryd 952 1291 4 12582.26 262.40 ryd2 987 1337 4 14367.13 290.70 soderkoping 666 910 4 18436.87 98.56 vadstena 581 778 4 12692.96 61.92 valla 527 707 4 9677.31 49.57

Table 55: Results with setting I30.

86 Name n m q Obj Time atvid-s 68 95 4 4012.52 1.73 brokind 179 199 4 6531.53 5.12 colonia 24 36 4 1020.65 0.96 ekangen 149 184 4 5397.55 3.76 mantorp-c 122 145 4 4855.36 2.64 mantorp 185 226 4 7829.94 4.97 rimforsa-o 94 109 4 3380.73 1.83 skanninge-c 122 152 4 4828.63 2.47 skanninge-n 62 80 4 2439.63 1.25 skanninge 127 158 4 4708.60 2.88 studentryd-e 76 104 4 1616.23 2.65 sturefors 183 212 4 6958.30 4.92 valla-liu 160 213 4 2830.16 4.56 vikingstad 178 213 4 7921.57 5.21 atvid 220 290 4 10801.84 7.57 borensberg 211 257 4 9995.68 6.67 kisa 311 359 4 9447.31 15.08 liu 291 394 4 5227.43 13.81 rimforsa 233 262 4 6612.25 14.07 ryd1 382 508 4 4010.59 23.80 studentryd-int 235 316 4 3099.69 9.61 studentryd 387 519 4 5400.81 27.88 linghem 836 1129 4 11161.93 190.10 ljungsbro 812 1024 4 23878.90 185.94 malmslatt 489 626 4 11250.31 45.16 mjolby 419 524 4 22750.39 30.38 ryd-m 1107 1495 4 17194.45 421.59 ryd 952 1291 4 12336.72 271.87 ryd2 987 1337 4 14958.87 308.63 soderkoping 666 910 4 17802.88 101.08 vadstena 581 778 4 12570.30 74.79 valla 527 707 4 9612.47 58.56

Table 56: Results with setting J30.

87 Name n m q Obj Time atvid-s 68 95 4 4881.41 1.00 brokind 179 199 4 9134.04 3.02 colonia 24 36 4 1042.11 0.69 ekangen 149 184 4 6770.16 2.30 mantorp-c 122 145 4 6760.74 1.79 mantorp 185 226 4 11477.70 3.36 rimforsa-o 94 109 4 4686.30 1.28 skanninge-c 122 152 4 6026.85 1.68 skanninge-n 62 80 4 2791.65 0.90 skanninge 127 158 4 6098.02 1.79 studentryd-e 76 104 4 1720.43 1.14 sturefors 183 212 4 8674.44 3.39 valla-liu 160 213 4 3323.40 2.68 vikingstad 178 213 4 10646.28 3.21 atvid 220 290 4 13963.04 5.40 borensberg 211 257 4 12660.79 4.34 kisa 311 359 4 13040.81 15.64 liu 291 394 4 6250.21 14.91 rimforsa 233 262 4 8794.24 6.18 ryd1 382 508 4 4753.12 28.20 studentryd-int 235 316 4 3739.98 7.15 studentryd 387 519 4 6287.69 29.28 linghem 836 1129 4 14713.11 274.90 ljungsbro 812 1024 4 32700.15 229.68 malmslatt 489 626 4 13933.76 50.67 mjolby 419 524 4 26573.34 34.09 ryd-m 1107 1495 4 20340.54 787.02 ryd 952 1291 4 14486.85 425.07 ryd2 987 1337 4 16957.66 522.63 soderkoping 666 910 4 22537.79 131.28 vadstena 581 778 4 15545.60 91.95 valla 527 707 4 11486.48 65.81

Table 57: Results with setting K30.

88 Name n m q Obj Time atvid-s 68 95 4 4657.57 1.02 brokind 179 199 4 9082.59 3.15 colonia 24 36 4 1034.15 0.77 ekangen 149 184 4 7265.83 2.32 mantorp-c 122 145 4 7093.67 1.67 mantorp 185 226 4 11332.70 3.47 rimforsa-o 94 109 4 4704.52 1.29 skanninge-c 122 152 4 6095.91 1.71 skanninge-n 62 80 4 2840.06 0.88 skanninge 127 158 4 5923.98 1.82 studentryd-e 76 104 4 1732.61 1.03 sturefors 183 212 4 9398.53 3.24 valla-liu 160 213 4 3286.13 2.59 vikingstad 178 213 4 10238.10 3.63 atvid 220 290 4 13890.50 5.40 borensberg 211 257 4 12124.70 4.49 kisa 311 359 4 12769.55 15.06 liu 291 394 4 6353.01 15.69 rimforsa 233 262 4 8837.82 6.84 ryd1 382 508 4 4767.51 26.87 studentryd-int 235 316 4 3697.74 7.22 studentryd 387 519 4 6307.04 29.70 linghem 836 1129 4 14804.78 259.57 ljungsbro 812 1024 4 32543.04 224.47 malmslatt 489 626 4 14130.35 47.65 mjolby 419 524 4 28322.90 32.29 ryd-m 1107 1495 4 20164.37 697.96 ryd 952 1291 4 14262.77 411.32 ryd2 987 1337 4 16699.76 465.33 soderkoping 666 910 4 22223.33 129.88 vadstena 581 778 4 15493.38 84.71 valla 527 707 4 11371.01 63.72

Table 58: Results with setting L30.

89 Name n m q Obj Time atvid-s 68 95 4 3848.21 1.37 brokind 179 199 4 5723.70 4.32 colonia 24 36 4 1010.20 0.80 ekangen 149 184 4 5362.64 3.40 mantorp-c 122 145 4 4692.24 2.46 mantorp 185 226 4 7558.31 5.00 rimforsa-o 94 109 4 3471.46 1.90 skanninge-c 122 152 4 4520.06 2.60 skanninge-n 62 80 4 2417.24 1.42 skanninge 127 158 4 4621.02 2.69 studentryd-e 76 104 4 1594.14 1.48 sturefors 183 212 4 7063.86 4.99 valla-liu 160 213 4 2641.81 4.07 vikingstad 178 213 4 7644.38 4.72 atvid 220 290 4 10248.13 7.90 borensberg 211 257 4 9608.79 6.66 kisa 311 359 4 8856.03 17.11 liu 291 394 4 5299.67 15.51 rimforsa 233 262 4 5861.77 8.06 ryd1 382 508 4 4159.98 28.04 studentryd-int 235 316 4 3081.31 9.08 studentryd 387 519 4 5344.33 29.85 linghem 836 1129 4 11636.02 181.93 ljungsbro 812 1024 4 23731.08 161.45 malmslatt 489 626 4 10706.37 48.84 mjolby 419 524 4 21181.79 34.93 ryd-m 1107 1495 4 17115.10 439.95 ryd 952 1291 4 12637.77 277.91 ryd2 987 1337 4 14564.45 311.75 soderkoping 666 910 4 18280.58 100.78 vadstena 581 778 4 12411.59 73.69 valla 527 707 4 9765.44 58.39

Table 59: Results with setting Q30.

90 Name n m q Obj Time atvid-s 68 95 4 3873.55 1.30 brokind 179 199 4 5801.88 4.10 colonia 24 36 4 1004.94 0.72 ekangen 149 184 4 5286.14 3.33 mantorp-c 122 145 4 4553.66 2.41 mantorp 185 226 4 7363.08 5.43 rimforsa-o 94 109 4 3335.25 1.89 skanninge-c 122 152 4 4576.26 2.43 skanninge-n 62 80 4 2419.27 1.28 skanninge 127 158 4 4998.31 2.32 studentryd-e 76 104 4 1598.55 1.52 sturefors 183 212 4 6807.50 4.73 valla-liu 160 213 4 2709.89 4.01 vikingstad 178 213 4 7698.10 4.80 atvid 220 290 4 10336.61 7.90 borensberg 211 257 4 9698.73 6.33 kisa 311 359 4 9300.94 17.63 liu 291 394 4 5280.46 16.57 rimforsa 233 262 4 5897.43 9.22 ryd1 382 508 4 4020.74 28.84 studentryd-int 235 316 4 3151.08 8.92 studentryd 387 519 4 5291.06 30.72 linghem 836 1129 4 12175.71 181.13 ljungsbro 812 1024 4 24535.63 163.81 malmslatt 489 626 4 10637.03 49.56 mjolby 419 524 4 20921.04 37.50 ryd-m 1107 1495 4 16814.76 434.84 ryd 952 1291 4 12500.55 289.60 ryd2 987 1337 4 14435.55 302.58 soderkoping 666 910 4 17995.44 101.29 vadstena 581 778 4 12370.34 72.51 valla 527 707 4 9966.27 58.80

Table 60: Results with setting R30.

91 Name n m q Obj Time atvid-s 68 95 4 3881.50 2.34 brokind 179 199 4 5963.77 8.11 colonia 24 36 4 997.27 1.02 ekangen 149 184 4 4972.71 6.42 mantorp-c 122 145 4 4562.24 4.76 mantorp 185 226 4 7511.88 9.76 rimforsa-o 94 109 4 3393.50 3.16 skanninge-c 122 152 4 4588.23 4.77 skanninge-n 62 80 4 2402.70 1.97 skanninge 127 158 4 4662.69 4.95 studentryd-e 76 104 4 1583.59 2.49 sturefors 183 212 4 7023.50 8.63 valla-liu 160 213 4 2683.97 8.73 vikingstad 178 213 4 7742.96 9.14 atvid 220 290 4 10302.02 15.35 borensberg 211 257 4 9650.71 12.49 kisa 311 359 4 9274.66 29.73 liu 291 394 4 5224.62 29.74 rimforsa 233 262 4 5928.61 13.88 ryd1 382 508 4 4159.71 51.33 studentryd-int 235 316 4 3070.78 17.82 studentryd 387 519 4 5235.90 58.03 linghem 836 1129 4 11795.56 287.08 ljungsbro 812 1024 4 24198.84 251.54 malmslatt 489 626 4 10867.44 81.69 mjolby 419 524 4 20850.33 57.13 ryd-m 1107 1495 4 16544.36 595.10 ryd 952 1291 4 12532.73 409.68 ryd2 987 1337 4 14022.28 445.76 soderkoping 666 910 4 18220.62 164.85 vadstena 581 778 4 12536.96 123.40 valla 527 707 4 9641.52 102.45

Table 61: Results with setting S30.

92 Name n m q Obj Time atvid-s 68 95 4 3809.50 2.36 brokind 179 199 4 5843.82 7.75 colonia 24 36 4 1003.67 0.98 ekangen 149 184 4 5302.86 6.83 mantorp-c 122 145 4 4754.22 4.51 mantorp 185 226 4 7665.46 9.40 rimforsa-o 94 109 4 3466.07 3.18 skanninge-c 122 152 4 4543.93 4.97 skanninge-n 62 80 4 2434.36 2.08 skanninge 127 158 4 4471.80 5.03 studentryd-e 76 104 4 1598.17 2.63 sturefors 183 212 4 6960.59 8.77 valla-liu 160 213 4 2808.27 7.97 vikingstad 178 213 4 7726.33 8.77 atvid 220 290 4 10785.91 15.04 borensberg 211 257 4 10187.33 12.99 kisa 311 359 4 9468.65 30.45 liu 291 394 4 5340.83 29.93 rimforsa 233 262 4 5800.12 17.71 ryd1 382 508 4 4299.86 49.97 studentryd-int 235 316 4 3037.44 18.68 studentryd 387 519 4 5306.95 51.16 linghem 836 1129 4 11329.21 305.95 ljungsbro 812 1024 4 23980.01 244.40 malmslatt 489 626 4 10636.01 81.18 mjolby 419 524 4 20923.71 59.60 ryd-m 1107 1495 4 17003.85 605.61 ryd 952 1291 4 12465.97 413.02 ryd2 987 1337 4 14117.40 448.99 soderkoping 666 910 4 18161.23 166.28 vadstena 581 778 4 12487.50 121.98 valla 527 707 4 9785.71 99.26

Table 62: Results with setting T30.

93 8.4 More tests with setting A

Name n m q Obj Time atvid-s 68 95 4 4005.00 0.25 brokind 179 199 4 5804.49 0.76 colonia 24 36 4 994.05 0.14 ekangen 149 184 4 5542.43 0.64 mantorp-c 122 145 4 4718.49 0.42 mantorp 185 226 4 7458.55 0.95 rimforsa-o 94 109 4 3387.66 0.32 skanninge-c 122 152 4 4638.72 0.50 skanninge-n 62 80 4 2383.14 0.22 skanninge 127 158 4 4651.02 0.46 studentryd-e 76 104 4 1651.93 0.26 sturefors 183 212 4 6495.91 0.85 valla-liu 160 213 4 2676.04 0.67 vikingstad 178 213 4 7561.85 0.77 atvid 220 290 4 10506.32 1.30 borensberg 211 257 4 10062.74 1.19 kisa 311 359 4 9393.75 3.46 liu 291 394 4 5242.79 2.93 rimforsa 233 262 4 5908.57 1.40 ryd1 382 508 4 4008.33 5.33 studentryd-int 235 316 4 2993.41 1.62 studentryd 387 519 4 5412.68 5.27 linghem 836 1129 4 11226.23 32.94 ljungsbro 812 1024 4 23525.80 31.05 malmslatt 489 626 4 10492.20 8.60 mjolby 419 524 4 21270.19 5.83 ryd-m 1107 1495 4 16912.03 77.47 ryd 952 1291 4 12782.49 49.41 ryd2 987 1337 4 14186.95 55.90 soderkoping 666 910 4 17935.21 18.30 vadstena 581 778 4 12438.91 13.21 valla 527 707 4 9648.21 9.93

Table 63: Results with setting A5.

94 Name n m q Obj Time atvid-s 68 95 4 3971.17 4.49 brokind 179 199 4 5585.59 13.82 colonia 24 36 4 990.21 2.64 ekangen 149 184 4 5334.80 11.19 mantorp-c 122 145 4 4382.86 7.59 mantorp 185 226 4 7544.77 18.71 rimforsa-o 94 109 4 3340.55 6.24 skanninge-c 122 152 4 4616.46 8.20 skanninge-n 62 80 4 2414.27 3.94 skanninge 127 158 4 4671.74 8.63 studentryd-e 76 104 4 1603.34 4.76 sturefors 183 212 4 6428.52 15.32 valla-liu 160 213 4 2624.69 15.07 vikingstad 178 213 4 7287.21 14.91 atvid 220 290 4 10179.22 26.40 borensberg 211 257 4 10065.06 22.56 kisa 311 359 4 9496.54 65.91 liu 291 394 4 5197.29 59.56 rimforsa 233 262 4 5700.51 28.22 ryd1 382 508 4 4082.52 100.28 studentryd-int 235 316 4 3062.97 34.76 studentryd 387 519 4 5275.56 106.38

Table 64: Results with setting A100.

95 8.5 Tests with 30 iterations and 5 vehicles

In the following tables we give the results for methods A - L with 30 iterations with q = 5.

Name n m q Obj Time atvid-s 68 95 5 3987.88 1.90 brokind 179 199 5 5799.34 7.29 colonia 24 36 5 1206.92 1.03 ekangen 149 184 5 5096.89 4.90 mantorp-c 122 145 5 4677.59 3.43 mantorp 185 226 5 7705.51 8.17 rimforsa-o 94 109 5 3381.00 2.45 skanninge-c 122 152 5 4620.91 4.26 skanninge-n 62 80 5 2561.69 1.77 skanninge 127 158 5 4683.98 3.91 studentryd-e 76 104 5 1818.81 2.27 sturefors 183 212 5 6910.18 8.54 valla-liu 160 213 5 2842.18 6.12 vikingstad 178 213 5 7653.05 7.61 atvid 220 290 5 10307.55 11.66 borensberg 211 257 5 10159.21 10.76 kisa 311 359 5 8885.59 27.34 liu 291 394 5 5387.26 21.91 rimforsa 233 262 5 5894.07 11.79 ryd1 382 508 5 4083.34 39.27 studentryd-int 235 316 5 3160.53 13.64 studentryd 387 519 5 5342.14 43.20 linghem 836 1129 5 11166.46 239.81 ljungsbro 812 1024 5 23214.98 216.77 malmslatt 489 626 5 10480.68 64.59 mjolby 419 524 5 21525.15 53.56 ryd-m 1107 1495 5 16663.00 543.67 ryd 952 1291 5 12243.81 373.00 ryd2 987 1337 5 14003.50 394.04 soderkoping 666 910 5 17596.59 131.95 vadstena 581 778 5 12246.02 100.07 valla 527 707 5 9849.49 81.10

Table 65: Results with setting A.

96 Name n m q Obj Time atvid-s 68 95 5 3967.34 1.95 brokind 179 199 5 5777.49 7.80 colonia 24 36 5 1200.88 0.98 ekangen 149 184 5 5426.71 4.71 mantorp-c 122 145 5 4805.23 3.44 mantorp 185 226 5 7339.66 7.37 rimforsa-o 94 109 5 3760.25 2.57 skanninge-c 122 152 5 4608.45 3.40 skanninge-n 62 80 5 2571.16 1.87 skanninge 127 158 5 4737.19 3.80 studentryd-e 76 104 5 1812.18 2.03 sturefors 183 212 5 6724.02 7.54 valla-liu 160 213 5 2836.53 6.24 vikingstad 178 213 5 8244.67 8.52 atvid 220 290 5 10802.35 13.42 borensberg 211 257 5 10261.07 10.11 kisa 311 359 5 9029.33 23.87 liu 291 394 5 5343.30 23.33 rimforsa 233 262 5 6130.88 12.44 ryd1 382 508 5 4145.66 40.18 studentryd-int 235 316 5 3174.40 14.44 studentryd 387 519 5 5373.19 47.06 linghem 836 1129 5 11894.24 242.17 ljungsbro 812 1024 5 23435.67 216.53 malmslatt 489 626 5 10604.59 66.09 mjolby 419 524 5 21319.18 49.02 ryd-m 1107 1495 5 16399.97 537.88 ryd 952 1291 5 12070.30 361.24 ryd2 987 1337 5 14247.75 394.64 soderkoping 666 910 5 18856.25 134.87 vadstena 581 778 5 12393.61 100.46 valla 527 707 5 9574.62 78.77

Table 66: Results with setting B.

97 Name n m q Obj Time atvid-s 68 95 5 4018.42 1.89 brokind 179 199 5 5872.64 8.55 colonia 24 36 5 1186.71 0.93 ekangen 149 184 5 5467.64 5.22 mantorp-c 122 145 5 4668.25 3.47 mantorp 185 226 5 7381.57 8.05 rimforsa-o 94 109 5 3480.41 2.36 skanninge-c 122 152 5 4604.30 3.49 skanninge-n 62 80 5 2629.85 1.78 skanninge 127 158 5 4875.59 4.10 studentryd-e 76 104 5 1798.16 2.03 sturefors 183 212 5 6739.48 7.49 valla-liu 160 213 5 2923.14 6.11 vikingstad 178 213 5 7619.68 8.94 atvid 220 290 5 9902.72 13.26 borensberg 211 257 5 10244.35 11.03 kisa 311 359 5 9140.79 23.55 liu 291 394 5 5238.76 22.61 rimforsa 233 262 5 6209.94 12.71 ryd1 382 508 5 4137.34 44.59 studentryd-int 235 316 5 3205.26 15.04 studentryd 387 519 5 5396.94 41.10 linghem 836 1129 5 11846.70 240.80 ljungsbro 812 1024 5 23811.71 216.44 malmslatt 489 626 5 10405.54 63.79 mjolby 419 524 5 21406.98 46.50 ryd-m 1107 1495 5 16628.76 546.85 ryd 952 1291 5 12401.26 364.18 ryd2 987 1337 5 14343.00 395.50 soderkoping 666 910 5 17619.56 146.26 vadstena 581 778 5 12095.35 93.95 valla 527 707 5 9566.57 75.75

Table 67: Results with setting C.

98 Name n m q Obj Time atvid-s 68 95 5 4155.64 1.96 brokind 179 199 5 5847.63 6.61 colonia 24 36 5 1202.97 1.06 ekangen 149 184 5 5542.87 5.18 mantorp-c 122 145 5 4688.94 3.54 mantorp 185 226 5 7226.81 7.76 rimforsa-o 94 109 5 3427.61 2.49 skanninge-c 122 152 5 4526.29 3.56 skanninge-n 62 80 5 2562.10 1.96 skanninge 127 158 5 4868.95 3.99 studentryd-e 76 104 5 1793.84 2.10 sturefors 183 212 5 6870.01 8.04 valla-liu 160 213 5 2852.37 6.52 vikingstad 178 213 5 7529.92 7.56 atvid 220 290 5 9919.69 11.96 borensberg 211 257 5 10163.24 9.96 kisa 311 359 5 9723.47 22.24 liu 291 394 5 5426.24 21.59 rimforsa 233 262 5 6147.06 11.89 ryd1 382 508 5 4144.10 38.16 studentryd-int 235 316 5 3197.19 12.36 studentryd 387 519 5 5376.62 39.72 linghem 836 1129 5 11199.33 244.57 ljungsbro 812 1024 5 23716.11 215.32 malmslatt 489 626 5 10566.58 60.84 mjolby 419 524 5 21652.35 46.23 ryd-m 1107 1495 5 16718.72 566.67 ryd 952 1291 5 12242.10 369.90 ryd2 987 1337 5 13826.04 393.49 soderkoping 666 910 5 18233.48 132.26 vadstena 581 778 5 12165.64 98.51 valla 527 707 5 9767.27 79.92

Table 68: Results with setting D.

99 Name n m q Obj Time atvid-s 68 95 5 3961.92 2.10 brokind 179 199 5 5611.07 8.53 colonia 24 36 5 1196.38 1.09 ekangen 149 184 5 5311.27 5.46 mantorp-c 122 145 5 4652.91 3.67 mantorp 185 226 5 7672.56 9.33 rimforsa-o 94 109 5 3576.94 3.10 skanninge-c 122 152 5 4608.32 4.34 skanninge-n 62 80 5 2600.08 2.03 skanninge 127 158 5 4594.41 4.53 studentryd-e 76 104 5 1791.13 2.46 sturefors 183 212 5 6733.57 7.73 valla-liu 160 213 5 2790.34 5.97 vikingstad 178 213 5 7807.42 7.70 atvid 220 290 5 10506.48 12.16 borensberg 211 257 5 10343.40 10.41 kisa 311 359 5 8920.98 23.77 liu 291 394 5 5243.24 26.22 rimforsa 233 262 5 5977.58 13.46 ryd1 382 508 5 4124.13 38.61 studentryd-int 235 316 5 3134.73 13.79 studentryd 387 519 5 5258.73 41.87 linghem 836 1129 5 11046.62 242.48 ljungsbro 812 1024 5 22857.41 231.03 malmslatt 489 626 5 10384.51 63.48 mjolby 419 524 5 21330.33 47.52 ryd-m 1107 1495 5 16249.87 547.77 ryd 952 1291 5 12124.06 353.81 ryd2 987 1337 5 13922.75 386.62 soderkoping 666 910 5 17100.48 130.04 vadstena 581 778 5 12355.20 92.18 valla 527 707 5 9659.27 74.75

Table 69: Results with setting E.

100 Name n m q Obj Time atvid-s 68 95 5 4156.19 2.20 brokind 179 199 5 6114.44 7.20 colonia 24 36 5 1187.74 1.16 ekangen 149 184 5 5349.28 5.83 mantorp-c 122 145 5 4819.10 4.17 mantorp 185 226 5 7328.24 8.76 rimforsa-o 94 109 5 3503.40 2.98 skanninge-c 122 152 5 4620.38 5.31 skanninge-n 62 80 5 2595.01 2.19 skanninge 127 158 5 4655.15 4.33 studentryd-e 76 104 5 1793.79 2.49 sturefors 183 212 5 6634.18 9.11 valla-liu 160 213 5 2831.35 6.72 vikingstad 178 213 5 7928.99 9.44 atvid 220 290 5 10270.10 12.52 borensberg 211 257 5 9974.57 10.42 kisa 311 359 5 9691.25 21.86 liu 291 394 5 5326.13 20.69 rimforsa 233 262 5 5744.86 11.56 ryd1 382 508 5 4080.61 36.31 studentryd-int 235 316 5 3179.72 16.53 studentryd 387 519 5 5374.59 40.52 linghem 836 1129 5 11131.78 239.86 ljungsbro 812 1024 5 23029.67 220.74 malmslatt 489 626 5 10704.15 63.31 mjolby 419 524 5 21465.38 47.36 ryd-m 1107 1495 5 16479.26 554.19 ryd 952 1291 5 12314.00 355.35 ryd2 987 1337 5 14130.20 411.53 soderkoping 666 910 5 18411.90 136.10 vadstena 581 778 5 12395.77 100.58 valla 527 707 5 9633.83 78.47

Table 70: Results with setting F.

101 Name n m q Obj Time atvid-s 68 95 5 4190.59 1.84 brokind 179 199 5 6036.10 6.91 colonia 24 36 5 1207.01 1.05 ekangen 149 184 5 5342.70 4.99 mantorp-c 122 145 5 4639.32 3.36 mantorp 185 226 5 7298.48 7.26 rimforsa-o 94 109 5 3554.81 2.53 skanninge-c 122 152 5 4570.41 3.94 skanninge-n 62 80 5 2616.16 1.73 skanninge 127 158 5 4969.37 4.09 studentryd-e 76 104 5 1782.10 1.98 sturefors 183 212 5 6851.67 7.57 valla-liu 160 213 5 2762.16 6.01 vikingstad 178 213 5 7511.70 7.57 atvid 220 290 5 10611.22 14.16 borensberg 211 257 5 10103.54 8.91 kisa 311 359 5 8708.15 23.00 liu 291 394 5 5232.30 22.97 rimforsa 233 262 5 6077.21 11.85 ryd1 382 508 5 4160.04 39.97 studentryd-int 235 316 5 3187.97 13.87 studentryd 387 519 5 5320.78 44.34 linghem 836 1129 5 11212.82 236.47 ljungsbro 812 1024 5 23257.03 210.47 malmslatt 489 626 5 10866.99 63.85 mjolby 419 524 5 20916.35 42.78 ryd-m 1107 1495 5 16696.79 526.43 ryd 952 1291 5 12185.99 343.50 ryd2 987 1337 5 14040.83 384.70 soderkoping 666 910 5 17454.80 142.92 vadstena 581 778 5 12295.61 95.87 valla 527 707 5 9837.86 78.91

Table 71: Results with setting G.

102 Name n m q Obj Time atvid-s 68 95 5 4334.23 2.19 brokind 179 199 5 6883.37 8.68 colonia 24 36 5 1205.78 1.10 ekangen 149 184 5 5453.42 5.50 mantorp-c 122 145 5 5141.02 4.24 mantorp 185 226 5 8025.43 8.95 rimforsa-o 94 109 5 3438.31 3.18 skanninge-c 122 152 5 4818.61 4.04 skanninge-n 62 80 5 2577.62 2.05 skanninge 127 158 5 4744.95 4.17 studentryd-e 76 104 5 1836.05 2.36 sturefors 183 212 5 6878.04 7.98 valla-liu 160 213 5 3125.10 7.72 vikingstad 178 213 5 8147.27 8.94 atvid 220 290 5 10664.48 14.20 borensberg 211 257 5 10396.56 11.54 kisa 311 359 5 9815.42 25.72 liu 291 394 5 5396.69 24.64 rimforsa 233 262 5 6851.96 14.56 ryd1 382 508 5 4132.63 45.72 studentryd-int 235 316 5 3316.70 15.89 studentryd 387 519 5 5539.42 43.07 linghem 836 1129 5 11275.81 244.21 ljungsbro 812 1024 5 23664.12 242.81 malmslatt 489 626 5 11415.54 65.61 mjolby 419 524 5 22638.72 54.90 ryd-m 1107 1495 5 16783.62 559.03 ryd 952 1291 5 12274.16 359.15 ryd2 987 1337 5 14627.38 398.83 soderkoping 666 910 5 17271.14 135.22 vadstena 581 778 5 12454.10 109.95 valla 527 707 5 9811.03 81.16

Table 72: Results with setting H.

103 Name n m q Obj Time atvid-s 68 95 5 4199.93 1.93 brokind 179 199 5 6125.43 6.41 colonia 24 36 5 1221.98 1.05 ekangen 149 184 5 5739.72 5.26 mantorp-c 122 145 5 4575.36 3.25 mantorp 185 226 5 8278.86 8.12 rimforsa-o 94 109 5 3468.92 2.34 skanninge-c 122 152 5 4944.06 4.04 skanninge-n 62 80 5 2674.04 1.99 skanninge 127 158 5 4672.25 3.75 studentryd-e 76 104 5 1838.66 2.07 sturefors 183 212 5 6698.96 8.18 valla-liu 160 213 5 2959.45 6.44 vikingstad 178 213 5 7781.38 7.10 atvid 220 290 5 10591.06 14.74 borensberg 211 257 5 11271.80 10.96 kisa 311 359 5 9379.75 24.37 liu 291 394 5 5519.10 22.95 rimforsa 233 262 5 6081.54 12.61 ryd1 382 508 5 4104.09 36.16 studentryd-int 235 316 5 3188.09 13.70 studentryd 387 519 5 5432.97 38.17 linghem 836 1129 5 11928.61 231.26 ljungsbro 812 1024 5 23681.34 204.00 malmslatt 489 626 5 10408.68 55.93 mjolby 419 524 5 23925.11 41.02 ryd-m 1107 1495 5 17008.45 540.46 ryd 952 1291 5 12801.69 347.52 ryd2 987 1337 5 14658.91 378.31 soderkoping 666 910 5 18288.09 128.92 vadstena 581 778 5 12940.92 94.30 valla 527 707 5 9934.63 70.91

Table 73: Results with setting I.

104 Name n m q Obj Time atvid-s 68 95 5 4380.76 1.97 brokind 179 199 5 6794.68 6.90 colonia 24 36 5 1202.17 1.00 ekangen 149 184 5 5471.02 5.05 mantorp-c 122 145 5 5067.42 3.66 mantorp 185 226 5 8223.74 8.50 rimforsa-o 94 109 5 3406.33 2.40 skanninge-c 122 152 5 4807.61 3.86 skanninge-n 62 80 5 2586.49 1.73 skanninge 127 158 5 4894.21 5.56 studentryd-e 76 104 5 1854.81 2.12 sturefors 183 212 5 6992.01 7.54 valla-liu 160 213 5 3151.67 6.57 vikingstad 178 213 5 8281.00 8.19 atvid 220 290 5 10697.34 13.16 borensberg 211 257 5 10301.98 11.81 kisa 311 359 5 9679.93 28.79 liu 291 394 5 5452.27 23.74 rimforsa 233 262 5 6879.24 14.30 ryd1 382 508 5 4122.35 39.97 studentryd-int 235 316 5 3291.25 15.14 studentryd 387 519 5 5553.05 43.03 linghem 836 1129 5 11442.34 239.34 ljungsbro 812 1024 5 23781.65 223.22 malmslatt 489 626 5 11166.23 66.44 mjolby 419 524 5 22886.68 50.56 ryd-m 1107 1495 5 16936.12 559.01 ryd 952 1291 5 12362.26 358.98 ryd2 987 1337 5 14611.03 393.22 soderkoping 666 910 5 17515.91 135.78 vadstena 581 778 5 12441.52 97.47 valla 527 707 5 9822.06 79.98

Table 74: Results with setting J.

105 Name n m q Obj Time atvid-s 68 95 5 5003.78 1.14 brokind 179 199 5 9830.67 3.61 colonia 24 36 5 1242.82 0.83 ekangen 149 184 5 8320.46 2.73 mantorp-c 122 145 5 7057.41 1.93 mantorp 185 226 5 11688.71 4.28 rimforsa-o 94 109 5 4946.19 1.48 skanninge-c 122 152 5 6187.37 1.99 skanninge-n 62 80 5 3134.36 1.08 skanninge 127 158 5 6448.76 2.34 studentryd-e 76 104 5 1991.01 1.24 sturefors 183 212 5 9851.81 3.93 valla-liu 160 213 5 3560.91 3.14 vikingstad 178 213 5 11344.06 3.69 atvid 220 290 5 15040.49 8.48 borensberg 211 257 5 12365.87 5.25 kisa 311 359 5 12604.19 15.72 liu 291 394 5 6773.62 15.96 rimforsa 233 262 5 9720.68 7.31 ryd1 382 508 5 5156.94 29.95 studentryd-int 235 316 5 4004.03 7.83 studentryd 387 519 5 6831.48 31.49 linghem 836 1129 5 15964.74 274.63 ljungsbro 812 1024 5 33779.64 236.16 malmslatt 489 626 5 15405.53 54.88 mjolby 419 524 5 29317.41 37.11 ryd-m 1107 1495 5 20644.68 772.63 ryd 952 1291 5 15049.67 441.91 ryd2 987 1337 5 17328.92 500.45 soderkoping 666 910 5 21652.53 133.75 vadstena 581 778 5 16426.24 91.08 valla 527 707 5 12564.23 68.02

Table 75: Results with setting K.

106 Name n m q Obj Time atvid-s 68 95 5 5117.14 1.25 brokind 179 199 5 9743.97 4.08 colonia 24 36 5 1212.12 0.85 ekangen 149 184 5 7577.26 2.83 mantorp-c 122 145 5 7188.69 2.24 mantorp 185 226 5 12250.25 4.50 rimforsa-o 94 109 5 4758.60 1.58 skanninge-c 122 152 5 6518.97 2.34 skanninge-n 62 80 5 3243.13 1.13 skanninge 127 158 5 6573.83 2.37 studentryd-e 76 104 5 1939.13 1.30 sturefors 183 212 5 9632.97 7.11 valla-liu 160 213 5 3521.31 3.45 vikingstad 178 213 5 10952.31 6.53 atvid 220 290 5 15102.13 7.04 borensberg 211 257 5 12591.97 6.34 kisa 311 359 5 13272.01 13.86 liu 291 394 5 6616.56 16.54 rimforsa 233 262 5 9683.27 8.42 ryd1 382 508 5 4991.79 27.64 studentryd-int 235 316 5 4038.76 8.40 studentryd 387 519 5 6567.92 29.74 linghem 836 1129 5 15398.13 279.52 ljungsbro 812 1024 5 33471.03 235.81 malmslatt 489 626 5 14713.82 51.85 mjolby 419 524 5 28765.31 33.62 ryd-m 1107 1495 5 20859.98 744.39 ryd 952 1291 5 14859.84 429.17 ryd2 987 1337 5 16852.39 496.10 soderkoping 666 910 5 22675.43 139.24 vadstena 581 778 5 16272.33 91.38 valla 527 707 5 12063.28 65.19

Table 76: Results with setting L.

107 8.6 Tests with 30 iterations and 6 vehicles

In the following tables, we present the tests with 30 iterations for methods A - J, for 6 vehicles. This is only done for the larger networks.

Name n m q Obj Time linghem 836 1129 6 11392.13 312.68 ljungsbro 812 1024 6 23934.07 275.13 malmslatt 489 626 6 10296.34 81.40 mjolby 419 524 6 21518.60 57.80 ryd-m 1107 1495 6 16582.52 683.31 ryd 952 1291 6 12227.31 459.98 ryd2 987 1337 6 13929.70 501.52 soderkoping 666 910 6 18309.84 178.67 vadstena 581 778 6 12534.69 130.53 valla 527 707 6 9367.20 104.00

Table 77: Results with setting A.

Name n m q Obj Time linghem 836 1129 6 11102.31 315.45 ljungsbro 812 1024 6 23057.87 273.34 malmslatt 489 626 6 10925.02 86.22 mjolby 419 524 6 22089.57 73.29 ryd-m 1107 1495 6 16699.06 695.27 ryd 952 1291 6 12387.43 451.41 ryd2 987 1337 6 13937.20 503.46 soderkoping 666 910 6 17394.15 179.65 vadstena 581 778 6 12107.50 132.13 valla 527 707 6 9710.82 105.20

Table 78: Results with setting B.

108 Name n m q Obj Time linghem 836 1129 6 11370.23 309.97 ljungsbro 812 1024 6 22485.16 282.53 malmslatt 489 626 6 10582.96 84.32 mjolby 419 524 6 21879.88 63.29 ryd-m 1107 1495 6 16113.94 744.05 ryd 952 1291 6 12095.26 446.19 ryd2 987 1337 6 13697.39 490.39 soderkoping 666 910 6 17342.14 176.52 vadstena 581 778 6 12374.84 128.81 valla 527 707 6 9968.37 107.07

Table 79: Results with setting C.

Name n m q Obj Time linghem 836 1129 6 11125.35 306.22 ljungsbro 812 1024 6 23112.16 303.09 malmslatt 489 626 6 10514.00 85.99 mjolby 419 524 6 21727.01 62.68 ryd-m 1107 1495 6 16281.24 692.72 ryd 952 1291 6 12425.14 456.49 ryd2 987 1337 6 13813.80 494.95 soderkoping 666 910 6 18603.38 177.44 vadstena 581 778 6 12218.11 133.65 valla 527 707 6 9596.32 105.19

Table 80: Results with setting D.

Name n m q Obj Time linghem 836 1129 6 11128.36 315.64 ljungsbro 812 1024 6 22677.11 277.73 malmslatt 489 626 6 10354.32 92.05 mjolby 419 524 6 21638.14 63.39 ryd-m 1107 1495 6 16125.77 679.59 ryd 952 1291 6 11997.35 452.54 ryd2 987 1337 6 14134.14 500.29 soderkoping 666 910 6 16820.99 197.14 vadstena 581 778 6 12171.60 130.72 valla 527 707 6 9368.08 105.57

Table 81: Results with setting E.

109 Name n m q Obj Time linghem 836 1129 6 11132.51 314.82 ljungsbro 812 1024 6 22533.56 279.56 malmslatt 489 626 6 10534.22 83.41 mjolby 419 524 6 22123.02 61.84 ryd-m 1107 1495 6 16396.77 709.62 ryd 952 1291 6 11943.41 452.58 ryd2 987 1337 6 13833.37 495.19 soderkoping 666 910 6 17243.02 178.38 vadstena 581 778 6 12389.54 134.48 valla 527 707 6 9642.39 104.15

Table 82: Results with setting F.

Name n m q Obj Time linghem 836 1129 6 11057.88 304.92 ljungsbro 812 1024 6 22894.70 275.34 malmslatt 489 626 6 10755.72 82.24 mjolby 419 524 6 21290.91 61.99 ryd-m 1107 1495 6 16482.94 671.67 ryd 952 1291 6 12557.74 443.92 ryd2 987 1337 6 14251.51 516.21 soderkoping 666 910 6 17756.65 173.91 vadstena 581 778 6 12313.33 129.14 valla 527 707 6 9758.34 102.21

Table 83: Results with setting G.

Name n m q Obj Time linghem 836 1129 6 11419.59 317.15 ljungsbro 812 1024 6 23396.21 282.91 malmslatt 489 626 6 10776.61 85.39 mjolby 419 524 6 23731.64 66.96 ryd-m 1107 1495 6 17061.52 712.64 ryd 952 1291 6 12341.53 493.46 ryd2 987 1337 6 14638.89 500.31 soderkoping 666 910 6 17389.80 184.41 vadstena 581 778 6 12433.01 131.32 valla 527 707 6 9934.98 108.92

Table 84: Results with setting H.

110 Name n m q Obj Time linghem 836 1129 6 11705.06 299.70 ljungsbro 812 1024 6 24238.56 273.85 malmslatt 489 626 6 10886.59 77.31 mjolby 419 524 6 22178.55 59.40 ryd-m 1107 1495 6 16947.61 667.29 ryd 952 1291 6 11935.03 440.63 ryd2 987 1337 6 14184.68 488.72 soderkoping 666 910 6 17872.35 170.62 vadstena 581 778 6 12413.86 140.17 valla 527 707 6 10084.40 98.98

Table 85: Results with setting I.

Name n m q Obj Time linghem 836 1129 6 11563.84 317.10 ljungsbro 812 1024 6 23203.57 280.31 malmslatt 489 626 6 11350.88 87.93 mjolby 419 524 6 23658.92 74.90 ryd-m 1107 1495 6 16797.32 697.00 ryd 952 1291 6 12347.11 456.35 ryd2 987 1337 6 14681.74 515.84 soderkoping 666 910 6 17539.79 202.87 vadstena 581 778 6 12493.50 130.95 valla 527 707 6 9814.37 109.23

Table 86: Results with setting J.

111