Recyclables Collection Route Balancing Problem with Heterogeneous Fleet

Nowadays, robust and efficient solid waste collection is crucial to motivate citizens to participate in the circular economy by sorting recyclable solid waste. Vocational vehicles, including garbage trucks, contribute significantly to CO2 emissions; therefore, it is strongly recommended, and in the European Union it is mandatory, to replace conventional-fuel-based garbage trucks with electric ones. For providing sustainable and energy-efficient solid waste collection with a heterogeneous fleet, in-depth mathematical computations are needed to support solving complex decision-making problems, including crew rostering and vehicle routing, because the distance and capacity of electric garbage trucks differ from conventional-fuel-based ones. However, the literature on solid waste collection using electric garbage trucks is still relatively scarce. The main contribution of this paper is developing an optimization problem for balancing travel distance assigned to each garbage truck of a heterogeneous fleet. The problem is based on specific requirements of the Municipal Solid Waste Management in Cracow, Poland, where the working time of routes is balanced and the total time of collection service can be minimized. For the problem, an MIP program was developed to generate optimal crew schedules, so that the hitherto network of segregated solid waste pickup nodes can be served using a heterogeneous fleet in which the share of electric garbage trucks is up to 30%. We study the impact of the changed composition of the fleet on modifications in crew rostering due to the shorter range of an electric vehicle compared to a conventional-fuel-based one.


Introduction
The selected solid waste collection process has been the subject of numerous research as it is an increasingly urgent and important issue faced by municipalities all over the world, because effective Solid Waste Management (SWM) is considered as an essential method for resource and energy recovery [1][2][3][4]. The existing substantial literature particularly concentrates on the Vehicle Routing Problem for garbage trucks performing the service of solid waste collection. Practical aspects of application OR-based methods (the most frequently VRPs) for improving decision-making in SWM were studied using cases of different municipalities, e.g., Rio de Janeiro, Brazil [5]; Austin, Texas, USA [6]; Malaga, Spain [7]; Mashhad, Iran [8]; Cracow, Poland [9][10][11][12]; La Palma, Spain [13]; Ipoh City, Malaysia [14]; Wageningen, the Netherlands [15]; and many cities in Portugal [16].
According to the The European Green Deal, the EU Member Countries are obliged to achieve climate neutrality by 2050 [17]; therefore, it is strongly recommended-and in the European Union, to some extent, even mandatory-to replace conventional-fuel-based vocational vehicles, including garbage trucks with zero-emission ones, e.g., electric vehicles [18,19]. In Poland, where electromobility still holds a small share of the transportation sector [20], the Act on Electromobility and Alternative Fuels obliges municipal service providers to guarantee that 30% of their fleets are comprised of zero-emission vehicles (ZEV) by 2028; the act foresees the following thresholds to reach: 5% by 2021, 10% in 2022, and 20% in 2025 [21]. Fleet management finds augmenting the share of electric vehicles (EV) challenging, due to the differences in the main characteristics of an electric truck compared to conventional ones [22][23][24], e.g., range, capacity, or noise level. In terms of municipal solid waste collection, electric garbage truck performance may result in rerouting (because of the shorter electric vehicle range) and rescheduling as solid waste collection at night can be allowed (because of electric vehicles' low noise level). With the growing share of electric garbage trucks in the fleet, the need for effective tools for routing a heterogeneous fleet also becomes more urgent. As solid waste management is getting increasingly complicated due to laws and regulations, there is a need for robust and efficient decision support systems in this area. Although the VRP for a heterogeneous fleet has been widely studied (e.g., by Coelho et al. [25], Schneider et al. [26], Taillard [27], Matthopoulos and Sofianopoulou [28], Li et al. [29], Berghida and Boukra [30], Kaewman and Akararungruangkul [31], Liu and Cao [32], Baldacci et al. [33], Toth and Vigo [34], Sousa et al. [35], Liong et al. [36]) the body of knowledge on the VRPs dedicated to solid waste management is still not impressive; to the authors' best knowledge, the garbage truck routing as a combinatorial optimization problem has neither been studied from the perspective of balancing the workload assigned to garbage trucks crew nor balancing the routes assigned to vehicles [37][38][39][40][41].
There are numerous approaches to organizing the municipal segregated solid waste collection system [42][43][44][45][46][47]. Roughly two approaches to selected recyclables collection can be distinguished. The first system is based on picking up segregated solid waste packed separately in dedicated bin bags directly from private pickup points, i.e., an individual real estate such as a private house, an apartment house, organization headquarters, companies, office buildings and other public facilities. This system is suitable for smaller towns or agglomeration districts dominated by detached residential houses. The second system uses a network of open access pickup points equipped with big recycling bins (BRB) dedicated to segregating different types of recyclable waste. BRB locations are not linked to any real estate and citizens have unlimited open access to them. Such a system works effectively in cities or districts which are built up with residential towers. In this paper, we concentrate only on the segregated solid waste collection based on BRBs in the form adopted in the municipality of Cracow, Poland. We formulate the Recyclables Collection Route Balancing Problem (RCRBP) and solve it using Mixed-Integer Programming (MIP) for an SWM with heterogeneous fleet. The main contribution of this paper is the adoption of an objective function representing the aim to balance the route kilometrage of the serving BRBs, i.e., balancing the workload on serving BRBs assigned to garbage collector teams. Such an approach is not typical for the VRPs, but is crucial for satisfying the specific requirements of the segregated solid waste collection system under investigation [48,49].
The remainder of this paper is as follows: In Section 2, we briefly review the recent advantages in Solid Waste Management, i.e., selective recyclable waste collection and sustainable energy-efficient fleet management including using green garbage trucks. We give special attention to optimization models dealing with routing problems with heterogeneous fleets, as these problems are crucial for decision making in Solid Waste Management (SWM) when reaching the mandatory share of ZEV within the fleet. In Section 3, we present the newly developed Recyclables Collection Route Balancing Problem. The problem fulfils the specific requirements of the BRB service performed by the Municipal Solid Waste Management in Cracow, Poland; therefore, it aims at balancing the workload assigned to each garbage truck instead of minimizing the total cost, which makes this approach unique among the VRPs. Next, in Section 4, we formulate an extension of the Recyclables Collection Route Balancing Problem, so that the requirements of using a heterogeneous fleet composed of conventional-fuel-based and electric garbage trucks can be met. Then, the results of computational experiments are reported and discussed. As the share of ZEV in the fleet serving solid waste collection is expected to grow and have a significant impact on SWM, in Section 5 we highlight some directions for future research.

Recent Advances in Decision Support Systems for Solid Waste Management
The SWM becomes an increasingly complex system and that makes it more expensive and requires elaborating logistics management using in-depth mathematical computation. Electric garbage trucks become increasingly popular in agglomerations all over the world as a cleaner alternative to conventional-fuel-based vehicles. As already mentioned, the European Union obliges member countries to use green garbage trucks, which significantly influences crew rostering and vehicle routing where the solid waste collection fleet is engaged. Such a situation results in further advances on decision support tools based on modeling processes in SWM systems and the VRPs with a heterogeneous fleet (i.e., mix VRP). A detailed review in this area is provided by Hannan et al. [50] and Sulemana et al. [51], which cover different types of optimization, objectives and constraints on solid waste collection, optimization methods and techniques, as well as research gaps, key issues, and challenges in this field.
On the table are both well-known and completely new issues. An example of the former is the trade-off between cost optimization and risk minimization in an integrated municipal SWM network, which is addressed by Yousefloo and Babazadeh [52] using a multi-objective mixed integer linear programming (MILP) model, solved using theconstraint method. The authors also included environmental aspects in the proposed model. Dolinina et al. [53] propose a mathematical model for rerouting in a solid waste collection system, while Puspita et al. [54] use an open capacitation VRP. An SWM is studied by Gdowska et al. [12] from another perspective. Selective solid waste collection is performed by a heterogeneous fleet and the kilometrage assigned to the collector crew must be balanced. Markov et al. [55] approach the recyclable waste collection problem as a vehicle routing problem with intermediate facilities, so that a heterogeneous fixed fleet and a flexible assignment of destination depots can be integrated. The problem is modeled as an MILP with several valid inequalities. Asefi et al. [56] investigate a VRP for a heterogeneous fleet with the aim of cost optimization by flexibly adjusting the fleet size. The problem is modeled with a bi-objective MILP model minimizing the transportation cost in the entire waste management system and the total deviation from the fair load allocation to transfer stations.
A new problem is introduced by Kang et al. [57], where modeling SWM is focused on the selective collection process of electric and electronics, because the amount of electric waste grows and becomes a challenge due to toxic chemicals and heavy metals in it. On the other hand, electric waste is considered as a valuable source of precious and base metals. The authors discuss the e-waste collection model for smart collection systems. Another innovation is proposed by Nowakowski et al. [58], where an algorithm and a productive model of the online system for scheduling solid waste collection are presented. The system includes parametric models of simulated annealing, tabu search, greedy, and bee colony optimization. An Internet platform enables comprehensive communication for people that request waste equipment for collection, registering of data, and solving the vehicle routing problem with time windows (VRPTW).
Erdinç et al. [59] present a Mixed-Integer Linear Programming (MILP)-based route optimization model as the first attempt explicitly for the waste-collection-oriented electric garbage truck routing process. The objective is to minimize total energy consumption, so that the solid waste collection process can be greener. Another mathematical model developed by Bányai et al. [60] combines routing, assignment, and scheduling problems for the solid waste collection process. The authors concentrated on optimizing the assignment of waste sources to garbage trucks, scheduling the waste collection through routing of each garbage truck to minimize the total operation cost, and increasing reliability, meanwhile considering comprehensive environmental indicators that have great impact on public health. A binary bat algorithm is used for solving the problem. Earlier, a simulation-based approach was used to model VRPs with a heterogeneous fleet partially composed with electric garbage trucks; Vonolfen et al. [61] presented simulation models for analyzing the scenarios of replacing conventional trucks with electric ones in the collection process of municipal glass waste.

Methodology and Materials
In this section, we present the Recyclables Collection Route Balancing Problem (RCRBP), which is the next one in the series of solid waste collection problems presented in the previous publications [11,12,62,63]. Here we consider a selective recyclables collection system which separates collection from open access BRBs and collection from recycling bins belonging directly to real estate (private pickup points); it means that to each recyclable collection route belong either BRB locations or pickup points located at real estates. The recyclables collection service is performed by one company for the entire agglomeration. The company uses a heterogeneous fleet of garbage trucks dedicated to collecting recyclables separately: (1) if a garbage truck is assigned to collect paper, it collects paper along the routes comprised of BRB locations as well as along the routes to which belong private pickup points, but (2) if a garbage truck is assigned to collect paper during the given go, it can be assigned to collect plastic on the next go, certainly after cleanup. During a working day, a garbage collector team can serve both a route in the BRB network and a route to which belong only private pickup points. In this paper, we focus exclusively on the emptying service dedicated to the network of open access BRBs. The BRB network is served with a heterogeneous fleet in which the share of electric garbage trucks constantly increases. We consider an SWM with a heterogeneous fleet of garbage trucks in which ZEVs can occupy any share. The objective is to find recyclables collection routes to be performed during a planning horizon (i.e., a certain number of planning periods), so that (1) every type of recyclable must be collected from every BRB location, (2) the working time of routes is balanced, and (3) the total time of collection service can be minimized.
We assume the demand for emptying service is deterministic, because it is estimated using historical data on filling in dynamics. The average fill level of BRBs was used also to set the frequency of emptying service during the planning horizon. The proposed approach does not consider the stochastic filling level of BRBs, because the monitoring of this index is out of the scope of the routing problem. In the SWM system on which our research is based, the filling level is not monitored by another unit and they set the size of BRBs in every location and the frequency of emptying them. For the routing problem, the size of BRBs and the frequency of service is a parameter, whose value cannot be changed. In our case, the frequency of emptying each BRB is once in a week (7 days), because the company found out that during a week the average filling level reaches ca. 60-80%. Moreover, garbage collector teams visually inspect the fill level every time and report it, so if BRBs in a given location are repeatedly filled below 50% or are overfilled, the management can analyze the data and undertake some actions, e.g., replace BRBs with smaller or bigger ones, relocate BRBs, increase the emptying frequency, etc. Finally, the analysis can show the need to reroute recyclables' collection.
In 2019, in the agglomeration of Cracow, Poland, BRBs were located in 639 locations. According to GIS data, the distance between BRB locations ranges from 1 km to 29 km, with the average equal to 7.3 km. In every location, the citizens could dump three types of recyclable material (r = 3): plastic and metal (yellow BRBs), paper and cardboard (blue BRBs), and glass (green BRBs). The recyclable collection service was provided everyday on weekdays (from Monday until Friday), so the schedule of segregated recyclable collection comprised of 5 days (i.e., five planning periods, t = 5). The service was performed by 12 garbage trucks (v = 12). Moreover, according to the SWM requests, each BRBs location can visit at most once a day, so that the citizens may not consider this service too arduous. There is a simplification adopted by decision makers to enable generating schedules manually without any dedicated decision support tool: each garbage truck has a fixed assignment to one type of recyclables as well as every garbage truck should serve only one route every day, so that balancing the working time of routes can be easier to obtain.
Using their knowledge and experience, the decision makers developed a Formula (1) for estimating the number of routes that guarantee collecting r types of recyclables with the fleet of v vehicles in t planning. Note that L is not proved to be the minimum number of balanced routes to serve the area. If decision makers are not supported with any optimization tool capable of minimizing the number of balanced routes of recyclables collection from BRBs, they can take it for granted that L is a reasonable number of feasible routes to be performed in a relatively short time and the working time of each route is similar. Formula (1) should rather be considered as a useful tool for making a quick decision on the number of routes to design and the number of garbage trucks to engage in serving BRBs.
In the SWM system under investigation, r = 3, t = 5, and v = 12, hence L = 20, which means that in the network of 639 BRB locations we should find 20 balanced routes for the selected recyclables collection. In Table 1, we present an exemplary schedule of the selected recyclable collection prepared by decision makers using manual tools. The fleet of 12 vehicles (i.e., 12 garbage collectors teams) is divided into three groups and each group is assigned to one recyclable: Teams 1-4 to plastic, Teams 5-8 to paper, and Teams 9-12 to glass. Roman numbers (I-XX) represent 20 routes which cover 639 BRB locations. Each route is used three times a week: for collecting separately plastic, paper, and glass. Routes are assigned to teams and vehicles with respect to the total pickup demand of the route, vehicle capacity, vehicle range, etc.
Note that in the schedule presented in Table 1 each garbage collector team (i.e., each vehicle) has a fixed assignment to one recyclable. After careful analysis of fleet capacity, fleet range, distances between BRB locations, and historical data on pickup demand that the number of vehicles engaged in serving BRBs should be reasonably low (and not necessarily minimum), but must guarantee proper service at every BRB location. Moreover, a route assigned to a garbage truck must be performed in one go. This is why Teams 1-12 every day serve only one route each, because for the rest of the working day the collectors teams and vehicles are engaged in performing other activities. The heuristic method of assignment and routing presented above is handy and provides feasible solutions with acceptable effort and computation time. However, the obtained solutions cannot be considered optimal, because the method assumes using every available collectors team and vehicle, even if a smaller number of teams was sufficient to perform the service. Moreover, the method does not guarantee indisputable balancing the working time assigned to routes, because the routing is separated from the assignment problem. For these reasons, we would like to propose an OR-based hierarchical approach combining clustering BRB locations and balancing the working time of garbage collectors teams engaged in serving BRBs.

Recyclables Collection Route Balancing Problem
In our research, we take the SWM system described in Section 3 as the reference for formulating the optimization problem called the Recyclables Collection Route Balancing Problem (RCRBP). The objective is to find routes for segregated solid waste collection to be performed in planning horizon T . Recyclables are represented by set R, so the decision maker may decide on the number of recyclables to be picked up separately basing on the specific features of the installation for processing segregated recyclables. Recyclables are collected by a fleet of garbage trucks (V ). The recyclables collection system is represented by a graph G = (N , E ), where N = {0, ..., n} is the set of nodes, and E is the set of edges. In set N nodes represent pickup points (BRB locations), except nodes 0 and n; node 0 is the depot where all the garbage trucks start and end their routes, while n represents a recycling center, which must be visited by every garbage truck to unload the collected recyclables.
The problem is deterministic, and for each recyclable r we know in advance the amount to be picked up from pickup point j during the whole planning horizon (d rj ). Each garbage truck v has its characteristic capacity c v . Note that parameters c v and d rj have the same unit based on the volume after compression. The RCRBP is less restrictive and allows the vehicle v to collect different types of recyclables during the planning horizon on condition that the vehicle is assigned to only one type each day. For each edge, (i, j) ∈ E , we know the travel time between nodes i and j (a ij ), so we can easily compute the distance between i and j. For each garbage truck v we know its maximum working time s v ; as working time is equal to driving time of garbage truck v, we can use s v to compute the maximum range of vehicle v.
Binary variable y vrijt yields value 1, if garbage truck v goes from node i to node j to collect recyclable r; otherwise, variable y vrijt yields value 0. The binary variable z vrt is equal to 1 if garbage truck v collects recyclable r in planning period t, otherwise z vrt is equal to 0. Non-negative variable x vjt represents the moment of the arrival of garbage truck v at pickup point j in planning period t.
Using the notation presented above, we can define the following MIP model for the Recyclables Collection Route Balancing Problem: ∑ v∈V ∑ i∈N ∑ t∈T y vrijt = 1, r ∈ R, j ∈ N : j > 0, j = n; ∑ i∈N ∑ j∈N :j>0,j =n,i =j ∑ j∈N :j>0 The objective function (2) minimizes the average working time (i.e., distance) in the system, so that the working time (i.e., kilometrage) assigned to each garbage truck (collector crew) can be balanced. Such an approach is typical for lot scheduling problems (such as the key performance indicator C sum ) but it is unusual for VRPs. Such an objective function suits the RCRBP as it is mandatory to use all vehicles in the fleet and exploit them to the same extent and intensity.
Constraints (3) assure that in the planning period, t, at most one garbage truck, v, collecting recyclables, r, traverses edge, (i, j), i.e., garbage truck v goes from pickup point i directly do pickup point j. Moreover, in the planning period, t vehicle v is allowed to collect only recyclable r from pickup point j, to which the garbage truck comes directly from pickup point i (4). Constraints (5) guarantee that recyclables of all types are collected from each node during the planning horizon, while constraints (6) assure that in each planning period vehicle v is allowed to collect recyclables of one type r top. In the planning period t, the amount of collected recyclables cannot exceed maximum capacity of vehicle v (c v ); it means that the capacity of a garbage truck A assigned to the route Z must be at least equal to the total capacity of BRBs to be collected along this rout. In other words, even if all BRBs belonging to the route Z are 100% filled, the capacity of garbage truck A is large enough to empty all bins. Constraint (7) guarantees that in planning period tdemand for collecting recyclables r (d tj ) cannot exceed capacity of vehicle v (c v ). This constraint requires the restrictive assumption that in the planning period t vehicle v can visit the recycling center and get unloaded only once. Such a simplification is allowed when the problem is adjusted to the specific requirements of a segregated solid waste collection system with BRBs in Cracow, Poland, where each garbage truck serves one type of recyclables in a given day and collects them from BRBs belonging to one route. The last location belonging to the route is the recycling center. For the rest of the day, the garbage truck is used for other tasks, e.g., collecting recyclables from detached houses in a residential district, but these tasks are out of the scope of the research to which this paper refers.
Constraints (8) assure that loops from node i to node i are avoided, while constraints (9) guarantee that loops i − j − i in the route of vehicle v are also forbidden. Constraints (10) guarantee that in the planning period t picking point i is served at most once no matter what type of recyclables is collected from this node; in other words, at most one type of recyclables can be picked up from node i in planning period t. Flow constraints (11) assure that if vehicle v enters picking point i, in the planning period, the t must also leave it.
Constraints (12)- (14) integrate the model by establishing the relationship between decision variables z vrt and y vr0jt . If in the planning period t, garbage truck v picks up recyclables r from pickup point i, then this garbage truck in this planning period can collect only recyclable r. Constraints (13) guarantee that each garbage truck has a mandatory stop at node n, where the vehicle is unloaded. To reinforce the mutual dependencies between decision variables, constraints (14) narrow the feasible region.
Constraints (15) and (16) enable the computation of the arrival time (i.e., distance) x vit of garbage truck v at pickup point i in planning period t. Arrival times are computed only if in the planning period t garbage truck v traverses edge (i, j), otherwise big constant M makes these contracts "inactive". Constraints (17) set arrival at node 0 to 0, which is interpreted as departure from depot and the beginning of the route of vehicle v in planning period t. Similarly, constraints (18) forbid the use of vehicle v after the end of working day s v .
The newly developed MIP model for RCRBP shows that it is possible to solve a real decision-making problem using an exact method after adjusting it to the requirements of the selected solid waste collection system. Here, an innovation is the adoption of an objective function representing balanced working time, i.e., balancing the routes assigned to garbage collectors teams. Therefore, the RCRBP can be used as a part of the decision support system in SWM for generating optimal weekly schedules of serving BRBs, despite the computation time for an MIP model potentially being relatively long. The weekly schedule of collecting segregated recyclables from BRBs is not frequently modified. Moreover, the locations of BRBs hardly change over a long time, so that the citizens may not change their habit of segregating recyclables with BRBs in fixed locations.
Due to the fact that this scheduling problem belongs to strategic decision making, we can invest some time in computations to obtain an optimal or suboptimal solution using an exact method. Therefore, despite RCRBP for large-size instances cannot be solved in polynomial time, it still can be a useful decision support tool, because after clustering we can obtain small-and medium-size instances, for which the optimal solution can be found in acceptable computational time. The MIP model can also be the point of reference for effective heuristics or meta-heuristics, which can find suboptimal solutions in a reasonable time.

RCRBP Adapted for Using Electric Garbage Trucks
Adjusting the RCRBP for the fleet partially composed electric garbage trucks some alterations are needed: (1) Removing the index t from decision variable z vrt , so it can be replaced with new decision variable z vr . Variable z vr yields value 1, if garbage truck v collects recyclables r during entire planning horizon, otherwise z vr is equal to 0. (2) Constraints need to be modified in the way presented in Formulas (19)- (22). ∑ r∈R z vrt ≤ 1, . . . ; −→ ∑ r∈R z vr ≤ 1, v ∈ V; (19) ∑ j∈N :j>0 Modified constraints (19) assure that during the entire planning horizon, vehicle v serves only one type of recyclables r. Constraints (20) were modified by introducing variables z vr and y vr0jt . If a garbage truck v departs from the depot to pickup point j to collect the recyclables r, then this vehicle is allowed to collect exclusively recyclables r during the entire planning horizon. Modified constraints (21) guarantee that every garbage truck travels to the recycling center (node n) where it can get unloaded. Constraints (22) still narrow the feasible region by strengthening the relationship between decision variables y vrijt and z vr .
When the minimum requested number of routes L (e.g., computed using Formula (1)), next two constraints- (23) and (24)-must be added to the MIP model. If each garbage truck is assigned to at most one type of garbage, it is sufficient to only add constraints (24).

Computational Experiments
Computational experiments were conducted with the aim to find the size of instance of the RCRBP for which the presented MIP model can be solved using a standard computer, because the model is supposed to be used in the company using typical available computational power. We used solver GUROBI 9.0.1 and a computer equipped with a two-core processor Intel Core i7-4710HQ CPU @2.50 GH and RAM 16 GB. There was no expectation to find a ready-to-implement route merely to test the RCRBP for an aggregated BRB network and its flexibly changed composition of the heterogeneous garbage truck fleet.
To conduct experiments with a clustering algorithm, we randomly generated a BRB network using an empirical probability distribution based on historical data. Randomly generated instances were used for validating and verifying the model, so that we could demonstrate its universal applicability for problems of the RCRBP family. In the recyclable collection system with BRBs, the distances between BRB locations were in the range 1-29 km, distances between with the average equal to 7.3 km (see Section 3). In Table 2, we present real distances, their standardized values, frequencies, and cumulative frequency; the data were used for developing the probability distribution function used for generating random instances. The obtained probability distribution was standardized and then the cumulative distribution function was computed (see Table 2). The empirical distribution of historical data is presented in Figure 1, while Figure 2 shows the probability distribution function fitted to the empirical distribution.   Similarly to the referential SWM system in the agglomeration of Cracow, Poland, the instance data are comprised of three types of recyclables to be collected during the five-day planning horizon. We adopt a strong simplification assuming that the demand for collection service is identical in every pickup point and it is equal to 1. Selective recyclables collection was to be performed by a fleet of 12 garbage trucks. We assume that 5 vehicles (i.e., ca. 40% of the fleet) are replaced by electric garbage trucks, whose range is up to 150 km, when the average route in the system under investigation is ca. 219 km. We assume the range of conventional-fuel-based vehicles has no upper limit. The total capacity of the fleet must (according to the problem description) be at least equal to the total demand. Therefore, there should not be a situation when the total demand is greater than the total vehicle capacity. Moreover, we assume that the BRB service frequency is adjusted to BRB fill level, so the amount of recyclables to be picked up never exceeds the BRB capacity. As it has already been mentioned, the filling level is monitored by another unit, so we can assume that the location and sizes of BRB answer properly the citizens' needs.
As the network of 639 BRB locations can be considered as a large-size instance, we used the clustering algorithm to aggregate BRB locations to get smaller instances of, respectively, 10, 15, and 20 pickup points (see Figure 3). Clustering is indispensable for getting instances solvable on a standard computer, as the VRPs are not solvable in polynomial time [34]. For clustering, we used the k-mean algorithm (the centroid algorithm). After clustering, the demand of clustered BRB points was aggregated. Certainly, any other clustering algorithm can be used for this problem with the same efficacy, because the goal is not to find routes but only to verify and validate the MIP model. A much more efficient approach is to obtain an optimal solution for a small-sized instance of clustered BRB location using an exact method (i.e., MIP) and then find a route within the cluster using an efficient heuristic algorithm. For further experiments with the MIP model, we selected three instances of the following size: 10, 15, and 20 clusters, i.e., aggregated BRB locations. Instances were generated using the probability distribution function described at the beginning of this subsection. We considered two types of the RCRBP: (A) when each garbage truck has a fixed assignment of one type of recyclables, and (B) when the minimum number of routes is defined. Note that the RCRBP(A) is far less restrictive than the way of generating pickup schedules adopted by the decision makers of the company (see Section 3), while the RCRBP(B) is an attempt to mirror the decision makers' approach with the predefined number of routes. The aim was to determine the recyclable collection routes, so that the total working time (based on distance) is minimized, i.e., the routes are balanced in terms of working time. Parameters used for the performed computations are presented in Table 3: we have instances of three sizes: 10, 15, and 20 aggregated pickup points to be served by maximum 12 garbage trucks including five ZEVs. In experiments with the RCRBP(B), we assume that the minimum number of routes is 20. Certainly, due to the clustered working time (distance) and aggregated pickup demand, we can expect that routes usually consist of one clustered pickup point. An interesting thing to observe will be the number of vehicles used to serve the system and their assignment to recyclables and routes. Here we can expect differences from the solutions generated manually by the decision makers. This is also the area where savings on vehicle usage and kilometrage can be made. Obtained results are presented in Table 4. The quality of solutions were assessed using GAP measurement, which is the gap between the obtained feasible MIP solution and the best estimation for the problem being solved for the given data instance. where: ObjBound-the best estimation for the instance, ObjVal-the best MIP solution for the instance. In Table 4, the obtained results are presented. The first and foremost observation is the exponential growth of the size of instances: for 10 clusters (instance I) the size of RCRBP(A) as the MIP problem is 37,525 × 18,780, for 15 clusters-82,915 × 41,520, and for 20 clusters-146,305 × 73,380. Moreover, the gap between MIP solutions obtained after 600, 1800, and 3600 s of computation and the optimal solution found for the linearized problem is great and range from 77.7% up to 93.7%. The gap increases with the growth of instance size. These observations confirm that the problem can be solved with exact methods only for small-sized instances and there is a great need for efficient algorithms for clustering and routing to be used to obtain applicable schedules based on the initial solutions found using MIP model.
We would like to emphasize that the assumption which helps the decision maker in generating recyclables pickup schedules using manual tools, i.e., computing the minimum number of routes (L ), may make the MIP model more almost imperceptibly less complicated. We can compare the size of instances for RCRBP(A) and RCRBP(B) for the same case: in case I (10 clusters), the number of columns of the RCRBP(A) is 18,780 while the the number of columns of the RCRBP(B) is 18,636; in case II (15 clusters) they are, respectively, 41,580 and 41,436; in case III (20 clusters) they are, respectively, 73,380 and 73,236. Nevertheless, introducing the required minimum number of routes to the problem can change significantly the obtained results, i.e., the number of routes and the assignment of collectors teams to routes and to recyclables. As a consequence, it may have a relevant impact on the total operational costs.

Conclusions
To summarize the obtained results, we must clearly state that research on the SWM systems with the fleet fully of partially consisting of ZEVs is scarce but urgently required. Thr presented example of the use of an exact method for decision system support is possible, despite the obvious limitations of them. However, on the level of strategic decision-making, these limitations can be overcome and the results may be used for such fleet routing, which results in decreasing energy usage and carbon footprint of segregated recyclables collection.The RCRBP presented in this paper has already been utilized for solving instances based on historical data, but future research should concentrate on testing the schedules obtained with the RCRBP in a real system.
The RCRBP and the developed MIP model can serve as a useful tool for assessing the fulfillment of the SWM's requirements and goals. Using the model, we cannot receive the ready-to-implement recyclables collection schedules, because of the complexity of the problem. However, optimal solutions for small-sized instances with clustered BRB locations can be used as the initial solutions for other efficient routing algorithms, so that recyclables collection routes and schedules can be easily modified if changes in pickup demand occurred. This paper justifies the need for new models for routing the fleet serving a segregated recyclable collection system. Due to the different requirements of SWMs, also atypical approaches to defining the problem can be successfully launched, such as the one used for the RCRBP.
Apart from developing tools based on exact methods, computer simulations are needed to enable experiments with real data on the fleet (number of vehicles, their range, and capacity), crew (different team sizes), operational costs (travel costs, staff costs, etc.) as well as stochastic recyclables' pickup demand.
The future research in this area should be focused on including more flexible garbage trucks in the problem, e.g., vehicles with compartments which can collect separately two types of recyclables in one go. Moreover, the BRBs in the SWM under consideration are not equipped with sensors monitoring the filling level, so for the time being, the collection system is not included in the Internet of Things. The fill level is visually inspected by collectors at every emptying service and reported to the management. These data are taken into consideration by decision makers during developing a routing strategy and changing BRBs locations or capacity. Data collected constantly could have brought interesting insights into the BRB filling dynamics and could have resulted in new routing strategies. Finally, the research should be concentrated on minimizing operational costs influenced by environmental decisions such as renewing the garbage trucks fleet with ZEVs.
Funding: This work was supported by the AGH University of Science and Technology with the grant for the maintenance of research potential. Research project supported by the program Excellence initiative-research university" at the AGH University of Science and Technology in Cracow, Poland.
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflict of interest. The funding institution had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript: