Extended Maximal Covering Location and Vehicle Routing Problems in Designing Smartphone Waste Collection Channels: A Case Study of Yogyakarta Province, Indonesia

: Most people will store smartphone waste or give it to others; this is due to inadequate waste collection facilities in all cities/regencies in Indonesia. In Yogyakarta Province, there is no electronic waste collection facility. Therefore, an e-waste collection network is needed to cover all potential e-waste in the province of Yogyakarta. This study aims to design a collection network to provide easy access to facilities for smartphone users, which includes the number and location of each collection center and the route of transporting smartphone waste to the ﬁnal disposal site. We proposed an extended maximal covering location problem to determine the number and location of collection centers. Nearest neighbor and tabu search are used in forming transportation routes. The nearest neighbor is used for initial solution search, and tabu search is used for ﬁnal solution search. The study results indicate that to facilitate all potential smartphone waste with a maximum distance of 11.2 km, the number of collection centers that must be established is 30 units with three pick-up routes. This research is the starting point of the smartphone waste management process, with further study needed for sorting, recycling, repairing, or remanufacturing after the waste has been collected.


Introduction
Developing countries such as Indonesia currently have the problem of handling large volumes of electronic waste (e-waste) [1]. It is associated with rapid technological and economic developments, leading to the production of a wider selection of electronic products at more affordable prices [2], thereby increasing public consumption and potential for electronic waste. The Global E-waste Monitor 2017 Quantities, Flows, and Resources ranked Indonesia ninth among the global producers of electronic waste, with smartphones being observed to contribute significantly. It is, however, important to note that the use of smartphones started increasing in 2020 due to the emergence of coronavirus, which prompted people to work and learn from home using online platforms. Records show smartphones are the technological devices with the highest consumption rate (70%), followed by laptops and personal computers [3], but there is no appropriate waste management process [4]. This is indicated by the absence of regulations for the collection and transportation of electronic wastes in Indonesia, with those implemented observed to be limited to informal initiatives. This, therefore, led to the low ranking of the country in the waste management system by the United Nations University. This is one of the major differences between Indonesia and developed countries [5]. A previous study also showed that improper handling of waste is dangerous for environmental sustainability [6].
About 80% of the materials composing smartphones can be recycled effectively [7]. Smartphones contain valuable materials, such as gold, silver, and palladium [8]. Metals in electronic waste, especially smartphones, are present in higher concentrations than in primary ore found in the ground. As an illustration, 300-350 g of secondary gold can be extracted from one ton of smartphones, while every ton of soil in ordinary gold mines only produces 5 g of primary gold [9]. Resource extraction from e-waste is more economical than extracting metal ores from the ground [10]. Thus, smartphone recycling is done because the economic benefits outweigh the costs [11]. Proper management of e-waste is necessary to reduce the problem of metal scarcity [8]. The potential for smartphone waste in Indonesia is quite significant. The total population of Indonesia in 2020 was 270,203,917 people [12]. If 63.53% are smartphone users [13], then the total number of smartphone users is 171,660,549 people. With an average smartphone lifetime of 4.7 years [14], these users produce 36,523,521 units of smartphone waste per year. When this waste is appropriately managed, in addition to minimizing the environmental impact, it can also provide significant economic benefits by producing 5.48-6.39 tons of secondary gold and saving natural resources.
However, so far, the amount of secondary metal recovered through e-waste recycling has been limited [15]; this is due to the limited supply of e-waste. A preliminary study conducted on smartphone users in Indonesia showed that 59% save non-functioning smartphones; 21% dispose of them; and the rest give them out to other people, sell them, and use them in other ways. This is because the public does not know what to do with these items. Meanwhile, Yogyakarta is one of the barometer provinces in Indonesia with an improper electronic waste management system through the formal channel. According to previous studies, government drivers are the factor with the most influence on consumers' intentions to participate in smartphone waste collection programs, followed by facility accessibility [16]. This means that the government needs to develop and implement a formal e-waste management system, starting with the e-waste collection process. One of the alternative electronic waste collection programs applicable to Indonesia is the use of a dropbox [17], but Yogyakarta Province does not currently have any collection points for smartphone waste. Therefore, there is the need to provide a convenient collection channel for the consumers, which is expected to be a major starting point for a formal channel to waste management in the area.
This study aimed to design a collection channel by determining the number and location of the collection center facilities followed by a transportation route from the collection center to the final disposal site. Facility location is related to the finding of a solution that covers customers using a number of facilities. It is, however, important to note that covering problems are fundamental facility location problems [18], which are often categorized as location set covering problems (LSCPs) and maximal covering location problems (MCLPs). The classic MCLP involves looking for the location of several facilities on the network in such a way that the population covered can be maximized [19]. Church and ReVelle first introduced this model in 1973 at a North American Regional Science Council [20]. The purpose was to maximize the demand covered by a particular service distance by placing a certain number of distribution facilities [21]. Therefore, customers or clients are declared covered when they are within a certain coverage distance from at least one facility [22]. The model is also important in the decision-making of the supply chain process, making it relatively important for practical use [23].
Several previous studies have used MCLP to design models and approaches in determining locations. MCLP is used in both the public and private sectors. In the public sector, it has been applied to determine the spread of an ambulance in emergency services [24], the location of emergency warning sirens [25], the location of medical equipment supply centers [26], the location of treatment centers in the event of a disease outbreak [27], appropriate locations for shelters for those temporarily displaced by floods [28], and the location of a waste cooking oil collection center [29]. Its use in the private sector involves determining the location of bank branches [30]. Several researchers have developed MCLPs. For example, Davari et al. [31] developed a MCLP with fuzzy travel times; Arana-Jiménez et al. [32] developed a fuzzy MCLP; Vatsa and Jayaswal ( [33,34]) modeled a capacitated multiperiod MCLP with server uncertainty; and Cordeau et al. [9] introduced the MCLP algorithm to determine a subset of facilities, maximizing customer requests by considering budget constraints. A continuous MCLP was also developed by Yang et al. [35] to optimize a continuous location of the cellular network's communication centers for natural disaster rescue. ReVelle et al. [36] also solved the MCLP with heuristic concentration which is used to determine a prominent case solution to maximum coverage locations with a high coverage percentage. Ibarra-Rojas et al. [37] developed a MCLP with accessibility indicators for when facilities have limited service areas, while Alizadeh and Nishi [38] used the hybrid covering location problem for strategic and tactical decisions. Alizadeh and Nishi [39] also developed a multiperiod maximal coverage location problem with different facility configurations as an extension of the classic MCLP. Zhang et al. [40] addressed the issue of locating multimodal facilities in emergency medical rescue.
The classical MCLP is used to determine the minimum number of facilities to maximize the demand covered by a given service distance. The model does not consider costs; it assumes that the number of facilities is minimal and that the investment costs are also minimal. Because each alternative location is assumed to have the same investment costs, it is necessary to develop a model that considers the difference in investment costs between potential locations. In this study, the collection center to be built is an intermediary facility, so it is also necessary to consider transportation costs to the final facility. Therefore, in this study, we develop the MCLP method by considering the investment and transportation costs, hereinafter referred to as an extended maximal covering location problem (e-MCLP). With this development, in addition to minimizing the number of facilities, it will also minimize total costs, including investment costs and transportation costs. Thus, the developed model is expected to provide an affordable facility location to consumers with a minimum total investment and transportation costs from the collection center to the final disposal facility.
The selection of the number and location of the collection centers was followed by the transportation route scheduling plan from the collection center to the final disposal site to determine the optimal route for efficient product distribution. It is defined as the route with the shortest distance and is considered important due to its ability to reduce transportation costs [41]. The vehicle route optimization problem, however, is known as the vehicle routing problem (VRP), which was introduced by Dantzig and Ramser in 1959 to solve the problem of gasoline distribution [42]. VRP is a common discrete optimization problem in transportation and logistics [43]. It is generally an integral part of the vehicle route with the exact delivery location visited once while all the routes start and end at the warehouse [44]. VRP focuses on the distribution of goods from the company's depot to customers and aims to minimize global transportation costs related to distance, fixed costs associated with vehicles and balance routes, and the number of vehicles required to serve consumers [45]. There are three methods of solving VRPs: the exact, heuristic, and metaheuristic methods [46]. However, the exact method is not applicable to a problem with a large input size and a limited time.
The methods used in this study include the heuristic and metaheuristic methods. The heuristic method involved the application of the nearest neighbor (NN) method, which has been widely used to solve VRP. Solomon introduced it in 1987 based on the idea of visiting the closest location from every other location visited [47], and it has been observed to be significantly better and to have more realistic performance in route formation than other methods [48]. This led to its wide application in solving the traveling salesman problem [49], determining routes from one city to another [50], designing waste transportation routes [51], and minimizing travel time and fuel consumption for transportation of agricultural products [52]. The nearest neighbor method is quite effective in its application due to its ability to look for consumers based on the closest distance from the vehicle's last location. It is, important to note that the nearest neighbor method produces the route with the shortest distance compared to other heuristic methods [41]. It is also easy to implement and execute the algorithm, but it does not guarantee the best resulting solution [53], so in this study, nearest neighbor was used to determine the initial solution.
This research applied the tabu search (TS) method, an algorithm considered to have the ability to produce an optimal solution. It was first introduced by Glover [54] based on the idea that allowing uphill motion helps to prevent the solution from becoming stuck in local optimal conditions [55]. The strength of this method lies in its flexible memory structure [54]. This makes its solutions very similar every time it is applied and makes it better than the other methods, such as simulated annealing and genetic algorithm [54]. Several studies have used tabu searches to solve VRP [56], classical VRP, periodic VRP, multidepot VRP, site-dependent VRP [57], heterogeneous fleet VRP [58], VRP with discrete split deliveries and pickups [59], multicompartment VRP [60], heterogeneous multitype fleet VRP with time windows and an incompatible loading constraint [61], multidepot open VRP [62], VRP with cross docks and split deliveries [63], VRP with private fleet and common carrier [57], time-dependent VRP with time windows on a road network [64], consistent VRP [65], and heterogeneous VRP on a multigraph [66]. Shi et al. [67] also used the heuristic solution method for the problem of multidepot vehicle routing-based waste collection and compared the results with the tabu search. Khan et al. [68] presented a sustainable closed-loop supply chain framework that uses a metaheuristic approach, tabu search, and simulated annealing. Tebaldi et al. [69] determined the best route to visit a set of customers, considering vehicle capacity and time constraints. This result underlies the use of the nearest neighbor approach to obtain an initial solution and the use of the metaheuristic tabu search approach to determine the final solution.

Materials and Methods
This research was conducted in two main stages: determining the number and location of collection centers and determining the smartphone waste transportation route. The location, number, and capacity of collection centers were determined by developing a maximal covering location problem hereinafter referred to as the extended maximal covering location problem (e-MCLP). The focus of the MCLP is to minimize the number of facilities while ensuring all consumers are covered, but the e-MCLP was developed to consider the costs involved. The model's objective was, therefore, to minimize the total costs, including those associated with investment and transportation from the collection facility to the final disposal site. The costs associated with collecting smartphones are not as high as those for other large volumes of e-waste, but the developed model can be used for other types of waste. The reason for choosing this type of waste is because it has a higher economic value (containing precious metals such as gold, silver, and palladium) than others, with components that allow up to 80% recycling and a large potential for smartphone waste. Meanwhile, for now, informal actors dominate the practice of recycling smartphone waste, which harms the environment. The low collection cost and high economic and environmental benefits are expected to motivate the government to implement the proposed scenario.
The development scenario involves two levels of collection center (CC) facilities, namely the primary collection center (PCC) and the secondary collection center (SCC). Consumers collect their waste at PCC. Instead, local governments carry out transportation from PCC to SCC. Transportation routes are needed in this study because smartphones are products with small volumes, so the capacity of the collection center is not as large as vehicle capacity. If one trip only picks up from one PCC, it becomes inefficient because the vehicle's utility is low, and transportation costs will be higher due to many trips being needed. For this reason, it is necessary to consider the route determination in this study. Routing is expected to increase vehicle utilization and save transportation costs. The output of determining the transportation route is expected to be an input for local governments to schedule waste collection.
Yogyakarta, one of the provinces in Indonesia, is located on Java Island and has an area of 3178.79 km 2 . It has a municipality and four regencies: Yogyakarta city and Gunung Kidul, Bantul, Sleman, and Kulon Progo regencies, with respective areas of 32.5, 1485.36, 506.85, 574.82, and 579.26 km 2 . These areas contain 14, 18, 17, 17, and 12 districts, respectively [70], as indicated in Appendix A, for a total of 78 districts. These districts were used as candidates for primary collection centers (PCCs) in this study. The parameters used as input in the mathematical model include the distance between the PCCs, the distance expected by consumers, and the distance from PCC to SCC.
Yogyakarta Province currently has 3 locations serving as final disposal sites (TPAs). The first is the Regional TPA, commonly called Piyungan TPA in Ngablak, Sitimulyo Village, Piyungan District, Bantul Regency. It is an integrated waste disposal site created to serve Yogyakarta City, Bantul Regency, and Sleman Regency [71]. The second location is Wonosari TPA in Wukirsari, Baleharjo, Wonosari, Gunung Kidul Regency, and the third is the Banyuroto TPA in Dlingo, Banyuroto, Nanggulan, Kulon Progo Regency. The Piyungan TPA has the largest capacity and most strategic location among the three, and this makes it suitable to be used as the secondary collection center (SCC). The candidates for the PCCs are district offices, which means the distances between PCCs are the same as those between district offices, and the distance from the PCC to SCC is the distance from the district office to the TPA Piyungan.
The PCC is provided by the government for consumers in the form of a dropbox, while SCC is a waste collection point for all the PCCs in a province. For this research, one SCC was located at the final disposal site in one province while the PCCs were built at the minimum number required to minimize investment costs incurred but with the ability to reach all consumers. Further, a survey conducted on smartphone users, with a total of 325 valid questionnaires, showed the consumers are willing to bring their smartphone waste to a collection facility with a maximum distance of 11.2 km. This means the PCC to be established is based on the number of districts to accommodate the interests of the consumers. Meanwhile, the PCC with the closest distance to the SCC was selected for this research to accommodate government interests by minimizing transportation costs. The PCC is located in the district office, a government-owned facility, and this means it does not require large investment costs since there is no need to procure land and a building, as only the dropbox needs to be prepared. This collection center has the capacity to accommodate all the smartphone waste supplies in the area due to the small product volume. It is important to determine the transportation route to optimize vehicle utility due to the relatively small volume of waste.
The location and capacity of the PCC were used to determine the transportation routes by joining the nearest neighbor approach and the tabu search model (NN-TB). The application of the NN was initiated from the starting point, which is the depot/SCC, and directed towards the PCC with the closest distance, which has not been visited due to several restrictions. The solution obtained at this stage is limited to determining the best route and the consumers to be served next based on the nearest point to the vehicle's last location [72]. It has been previously stated that the nearest neighbor algorithm is easy to implement and execute but does not guarantee the maximum resulting solution [53], and this was the reason it was used in this study to determine only the initial solution. Afterward, the tabu search method was used to search for the optimal route. The metaheuristic method is usually applied to solve combinatorial optimization problems, where the combinations are usually used to calculate the number of exchanges to be made in each iteration [73]. The tabu search algorithm is also a mathematical optimization method that guides the iterative search for solutions by providing tabu status for solutions found [74].

Collection Center Determination Steps
The parameters used as input in the mathematical model are the distance between PCC candidates, the distance expected by consumers, and the PCC candidate's distance to the final disposal site or secondary collection center (SCC). The distance between PCC candidates and distances between each PCC and SCC were based on Google Maps. The distance matrices between PCC candidates and from the PCC candidates to the SCC are shown in Appendix B. The distance value is essential to determine the number and location of PCCs to be built in the area.
The notation used in the mathematical model of e-MCLP is as follows: The number of the district (m = 1, 2, . . . , |m|) k The number of SCC (k = 1) Distance requirements (fulfilled or not) The capacity of PCC at point j Y i Coverage of smartphone waste supply at point i (covered or not) Supply of smartphone waste at point i The basic model was developed from the MCLP [75] in the form of e-MCLP, and its functional objective was to minimize the total cost of the number of facilities to be established within the range wanted by the consumer, as shown in Equation (1). The costs considered include those associated with the investment and transportation from PCC to SCC. Furthermore, the PCC was established in a district office, a government facility, which means there was no need to invest money in land acquisition. Therefore, the only investment needed was the procurement of dropbox, and the value is the same for all candidate locations. It is important to note that the PCC locations selected were those with the lowest investment costs and closer to the SCC. The decision variable X j has a value of 1 or 0, where a value of 1 indicates the point j is selected as a PCC and a value of 0 indicates the point j is not selected as a PC. Dropbox procurement costs are USD 350.37 (USD 1 is equivalent to IDR 14,270.75), and the dropbox service life is 5 years; using the straight-line depreciation method, the annual depreciation cost is USD 70.07 per dropbox. Thus, the investment cost per year is USD 70.07 per dropbox. The vehicle's fuel consumption is 10 km/L at USD 0.67 per liter; therefore the transportation cost is USD 0.067 per kilometer.
Equations (2)-(6) are constraint functions. Equation (2) is a limiting function that requires a ij X j to be 1, and this means a minimum of one PCC needs to be established within the range of the consumers' point. Meanwhile, Equations (3)-(5) state that X j , a ij , and Y i are binary, while Equation (6) states that the PCC capacity at point j is the accumulation of the waste supply multiplication at point i by 1 or 0, where 1 means the waste supply at point i is covered and 0 means it is not covered. Smartphones are, however, usually in small volume and not too large a supply due to the estimation of lifespan at two years. Therefore, the PCC capacity value used in this research is 1, which indicates that the entire waste supply was accommodated.

Steps to Determine the Transportation Route
The method used in this research was the nearest neighbor and tabu search (NN-TS) method, where the results obtained from the nearest neighbor were used as input in the tabu search. It is important to note that the tabu search was initiated by approaching a local minimum and noting recent movements in a tabu list that forms an adaptive memory to explore better solutions, with its size indicating the degree of diversification and intensification [76]. The mathematical model was, however, first determined before the calculations, and this was based on several assumptions and limitations, which include the following: (1) the vehicle has enough capacity to accommodate smartphone waste; (2) the distance from location j(a) to j(b) is the same as the distance from location j(b) to j(a) due to symmetry; (3) collection activities to PCCs start from 08:00-16:00 WIB with a rest time of 1 h, and this means the planning time horizon for a day is 7 h; (4) one vehicle visits more than one PCC but each PCC is only visited by one vehicle; (5) the average vehicle speed is 45 km/h; (6) the loading time at a PCC is 10 min; (7) the unloading and administration time at the SCC is 30 min. The notation used in the mathematical model of VRP is as follows: V The set of all vertices with 0 is a SCC {0, 1, 2, . . . , v} P The There is a trip from PCC at point j(a) to j(b) on trip t or not The objective function of the VRP mathematical model is to minimize the total distance traveled from the route as shown in Equation (7). The decision variable X t j(a)j(b)c has a value of 1 or 0; 1 indicates the selected route when vehicle c travels from PCC at point j(a) to j(b) on the trip t, and 0 indicates when the situation is otherwise. Equations (8)- (14) are constraint functions, with Equations (8) and (9) used to show that the route starts from and returns to SCC. Equations (10) and (11) state that each PCC is served exactly once on one route. Hereinafter, the vehicle's load capacity on a trip is the accumulation of the PCC capacity served, and its maximum capacity is not exceeded, as shown in Equation (12). This is because the supply is not large and the product volume is small, which allows the vehicle to carry the entire supply of smartphone consumers at once. Meanwhile, Equation (13) shows the vehicles going to the SCC to unload. However, the route completion time was calculated from the vehicle's total time plus the service time, which is loading-unloading time, and observed not to have exceeded the planning time horizon in a day, which is 7 h, as shown in Equations (14) and (15).
The steps to determine the initial solution using the nearest neighbor method [72] are as follows: a.
Select the center point as the starting point of transport, which is the SCC in this study. b.
Determine the point with the smallest distance from SCC and move to the PCC point. c.
The last point visited is the starting point; therefore, determine the point with the closest distance from the point. d.
Repeat the process until the vehicle does not have sufficient capacity for transportation; but because there is always enough capacity of the vehicle used in this research, the repetition is conducted until it meets the planning time horizon for a day but does not exceed it. e.
Drag this point to a line which is called a route with the working hours used as a constraint to form a freight route.
The Tabu search algorithm used in this study is based on [77,78] and includes the following steps: a.
Determine solution representation. This is a sequence of nodes where each is only visible once in the sequence. These nodes represent PCC and SCC. b.
Formulate initial solution formation, S. c.
Determine the neighborhood solution. This is an alternative solution obtained by moving the nodes such that each move produces a neighborhood solution and the number of solutions is calculated using the following Equation (16): where n is the number of PCCs visited in a route. d.
Create a tabu list. This list contains the moved attribute previously found, and its length increases with the size of the issue and also corresponds to the number of PCCs to be visited. e.
Find the best solution, S*. f.
Fix the tabu list. g.
Determine aspiration criteria. This is a method of overturning the tabu status. h.
Determine termination criteria. These are used after all predetermined iterations have been fulfilled. The number of iterations selected is the same as the number of points visited because the maximum number of iterations is the same as the length of the tabu list [79].

Number and Location of Collection Centers
The number and location of the PCCs were determined using the e-MCLP method. The solver software was used to determine the optimal solution. The calculations showed that 30 PCCs are to be built as shown in Figures 1 and 2 with a distribution of 1 unit in Yogyakarta city (Y6), 13 units in Gunung Kidul Regency (G1, G2, G4, G5, G6, G7, G8, G9, G13, G14, G15, G16, and G17), 6 units in Bantul Regency (B4, B10, B12, B13, B15, and B16), 6 units in Sleman Regency (S2, S6, S11, S12, S13, and S17), and 4 units in Kulon Progo Regency (K5, K8, K10, and K12). The selected PCC numbers and locations are shown in Appendix C.  The nearest neighbor method's search for initial solutions started with the 7 h obtained for planning horizon time, a loading time of 10 min for each PCC, and 30 min of unloading and administration time at SCC. This was followed by the determination of the depot as the starting location, which is the SCC. The vehicle has the capacity to accommodate the entire PCC because the supply is not large and the product volume is small; therefore, the planning time horizon was considered. The next step was the determination of the PCC with the closest distance, and this was discovered to be Pleret PCC, which has a distance of 4.3 km from the SCC. The distance matrices between the selected PCCs and from the selected PCC to the SCC are shown in Appendix D. It is important to note that the retrieval process was continued to the next PCC when the completion time (CT) was less than or equal to the planning time horizon but canceled when the completion time was greater than the planning time horizon. Furthermore, the next PCC was determined based on the closest distance with the initial steps implemented when it was discovered not to have been served. It is also important to point out that just one type of vehicle was used.
The number of trips or tours required to make the collection was calculated to be 3 with a total distance of 659.1 km, travel time of 14.65 h, and a completion time of 21.16 h, as shown in the sequence presented in Table 1 It was discovered that Route 3 has a longer travel time than the planning time horizon, and it was used as an initial solution in the tabu search method with the expectation that it will improve and provide shorter distances and times for the optimal solution. The tabu search method was applied using the initial solution calculated from the nearest neighbor method. Route 1 was found to be SCC-Pleret PCC-Kotagede PCC-Sewon PCC-Pandak PCC-Bambanglipuro PCC-Sedayu PCC-Dlingo PCC-Playen PCC-Patuk PCC-Ngglipar PCC-Ngawen PCC-SCC. This was followed by the input of the number of elements to be searched, which was found to be in accordance with the points to be visited, i.e., 11 PCCs. The number of neighborhood solutions was later determined using Equation (14), and 55 lines were recorded. Furthermore, the tabu list length was discovered to be in line with the number of PCCs to be visited, which was 11 customer locations. This was followed by the maximum number of iterations, which was recorded to be 11 iterations in line with the number of PCCs. These steps were repeated for the other routes, and the determination of the best route produced three routes with a total distance of 602.2 km, a travel time of 13.4 h, and a completion time of 19.89 h. The time for a shipment was found to be 3 days. Furthermore, the best sequences for Routes 1, 2, and 3 had total distances of 178.3, 198.3, and 224.5 km; travel times of 3.98, 4.41, and 5.01 h; and completion times of 6.63, 6.58, and 6.68 h, respectively, as shown in Figure 3 and Table 2.

Discussion
The results showed that the city/regency with the fewest PCCs is Yogyakarta city due to the short distance between its districts, with the one PCC established in Kota Gede district being found to have the ability to reach 13 other districts. The farthest is the Tegalrejo district, which is 9.5 km away, and this is also considered to be within the distance desired by the consumers. Meanwhile, most of the PCCs were built in Gunung Kidul Regency due to its large area relative to the other cities and regencies, and this caused quite a long distance between the districts. The area is 47% of the total area of Yogyakarta Province, as shown in Figure 2. Therefore, there is a need to build 13 PCCs in the existing 18 districts to cover all consumers, and the remaining 5 will be accessible because they are less than 10 km from the built locations. For example, Playen PCC covers Paliyan District while Semanu PCC covers Ponjong and Karangmojo districts. It is also possible for the waste from Wonosari District to be transported to Playen or Semanu PCC, while Ngawen PCC covers the Semin district.
This problem, if solved using MCLP as done by Church and Davis [70], Murray [17], Boonmee et al. [18], and Hartini et al. [26] to minimize the number of collection centers that must be built, results in the same number of PCCs that must be built as with e-MCLP, namely 30 PCCs; the number of PCCs established in each city/regency is the same, but there are several different locations. The comparison of selected CCP locations from each method is shown in Figure 4. The different locations are Kotagede, Semanu, Berbah, and Kalibawang when using the e-MCLP method. When using the MCLP method, the selected locations are Gondomanan, Ponjong, Minggir, and Samigaluh, as shown in Figure 4. The difference between the four locations will have implications for saving transportation costs from PCC to SCC because of the shorter distance, while the investment costs, in this case, are the same for each selected PCC. Comparison of the distance from PCC to SCC between the two methods is shown in Table 3 Table 4 shows that the distance between the selected PCCs and SCC is shorter in e-MCLP than MCLP. This indicates the numbers and locations calculated using the two approaches were able to accommodate the range expected by the consumers, but e-MCLP considered the investment costs and the distance between the PCC and SCC, unlike the MCLP. Therefore, MCLP provided a greater total PCC to SCC distance, which is indicated by 989.6 km with transportation costs of USD 131.76 when the collection is at only one PCC and the total cost required per year is USD 8426.27. Meanwhile, e-MCLP provided a shorter total distance of 40.9 km with 4.13% savings in transportation costs at USD 261.6 per year. e-MCLP is very suitable for PCCs with large waste volumes because vehicle capacity is filled faster when the volume of waste is large so that there are fewer pick-up points on one route. When there are fewer pick-up points in one route, the more routes there will be, and the development of this method is suitable for implementation. This model's savings in transportation costs will be felt when the number of routes increases because vehicles will depart and return to SCC more often. That is, the closer CCP distance to the SCC is very beneficial for the vehicle. In this study, the selected location does not affect the investment cost because each candidate location requires a procurement cost of the same amount. However, the developed model can accommodate each candidate location requiring a different investment cost. Later, the selected location will provide a minimum total cost, including investment and transportation costs.
The best route was determined using the tabu search method to improve the results of the nearest neighbor. This is in line with the opinion found in [77,80,81] that metaheuristics are popular optimization problem-solving techniques to overcome the weaknesses of the heuristic method due to their ability to avoid being trapped in a local optimum solution [82]. The first route was found to be better than the original solution due to its ability to reduce the distance traveled by 15.6 km and the travel time by 0.33 h, thereby reducing the distance and travel time by 8%. Meanwhile, the optimal solution in the second route is the same as the initial solution, but the route completion sequence is reversed such that the first PCC visited using the nearest neighbor was the last in the tabu search method. This shows the nearest neighbor method also has the ability to provide the best solution, and this is in accordance with the findings of [41] that the nearest neighbor method produces the shortest route compared to other heuristic methods. The NN algorithm was able to minimize distribution costs [83] and could easily and quickly resolve problems for several small cities [84]. Furthermore, the initial solution was observed to be infeasible for the third route because the completion time, which was recorded to be 7.6 h, exceeds the planning time horizon, which is 7 h. The continuation of the iteration using the tabu search method changed the initially infeasible solution to feasible as indicated by the shortening of the completion time (CT) to 6.68 h with a total distance (D) of 225.5 km and a travel time (TT) of 5.01 h, saving 8.53% of travel time. This means the tabu search was able to reduce the distance and travel time by 15% as indicated by the 41.3 km and 0.92 h results when compared to the nearest neighbor method, as shown in Table 4.
The tabu search method was generally able to provide better performance than the nearest neighbor method. The metaheuristic approach gives better performance results than the heuristic approach [70]. The results showed the possibility of collecting all the smartphone waste in Yogyakarta Province using three routes. This can be completed in a day through the use of three vehicles or in three days through the use of one vehicle. The total distance required to be covered is 602.2 km with a travel time of 13.4 h and a total completion time of 19.89 h. This means the tabu search method generally saved 56.9 km (8.6%) distance and 1.25 h (8.5%) travel time.
Determination of smartphone waste collection routes in the province of Yogyakarta with one route picking up at several PCC points managed to save a mileage of 346.5 km compared to one route only picking up at one PCC point and a total of 30 pick-up points. If the smartphone waste collection is done once a week, this shorter distance can provide transportation cost savings of USD 2214.39 per year. The area of Yogyakarta Province is only 0.16% of the territory of Indonesia; if this model is implemented nationally, the estimated transportation cost savings will be more than USD 1 million.
This research is expected to be the initial framework in formulating e-waste management policies for the national formal channel. If this proposal is successfully implemented in Yogyakarta, it is likely to be implemented in other provinces in Indonesia. The developed model can also be used for other solid waste collection scenarios. The proposed e-MCLP model is very suitable for large e-waste because there is no need to proceed to route determination, considering that the supply from PCC may already meet vehicle capacity. One trip only picks up from a PCC and then returns to the SCC again. However, to use the proposed model, it is necessary to consider whether the community is willing to bring their large size/volume e-waste to the provided PCC. This research is also the first step in electronic waste management, which will then be followed by the next stage of management, which includes separation, repair, recycling, remanufacturing, or disposal.

Conclusions
Smartphone waste has a high economic value and has great potential. The tendency of people to store and dispose of smartphone waste is due to the absence of waste collection facilities and government regulations that specifically regulate electronic waste management mechanisms. With the public's willingness to bring smartphone waste to a collection point with a maximum reach of 11.2 km and the benefits that will be obtained, this is a challenge and an opportunity for the government to design an optimal collection channel. The design of the collection channel involves consumers as suppliers of electronic waste, primary collection centers (PCCs), and Secondary Collection Centers (SCC). Due to the small area of Yogyakarta Province, 1 SCC is sufficient to accommodate the supply of smartphone waste from all selected PCCs. Based on the results of calculations using e-MCLP, as many as 30 PCCs should be built, with a distribution of 1 PCC in Yogyakarta City, 13 PCCs in Gunung Kidul Regency, 6 PCCs in Bantul Regency, 6 PCCs in Sleman Regency, and 4 PCCs in Kulon Progo Regency. e-MCLP can produce the minimum number of primary collection facilities required to cover all consumers with the shortest distance from secondary collection facilities to minimize total costs, including investment and transportation costs, with a total cost of USD 3617.92 per year.
The best transportation route from PCC to SCC was determined using the nearest neighbor and tabu search method (NN-TB). The pick-up route starts and ends at SCC, and the result shows three routes to use in smartphone waste collection. These routes take three days to complete by using one vehicle or one day using three vehicles with a total time required of 19.89 h and a distance of 602.2 km.
Further research can expand the study of e-waste with a large volume because the large volume will affect the willingness of consumers to bring their e-waste and the need to calculate the capacity of the collection center. This research is expected to be the initial framework in formulating e-waste management policies for a formal national channel. Research can also be continued with the design of management following the collection of e-waste in a final disposal site, such as separation, repair, recycling, remanufacturing, or disposal.

Data Availability Statement:
The data used to support the findings of this study are available from the corresponding author on request.