Next Article in Journal
Oil Price Shocks to Foreign Assets and Liabilities in Saudi Arabia under Pegged Exchange Rate
Next Article in Special Issue
The Vehicle Routing Problem with Simultaneous Pickup and Delivery Considering the Total Number of Collected Goods
Previous Article in Journal
Wearable Inertial Sensor-Based Hand-Guiding Gestures Recognition Method Robust to Significant Changes in the Body-Alignment of Subject
Previous Article in Special Issue
Robust Appointment Scheduling in Healthcare
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metaheuristics in Business Model Development for Local Tourism Sustainability Enhancement

by
Pawnrat Thumrongvut
1,
Kanchana Sethanan
1,*,
Thitipong Jamrus
1,
Chuleeporn Wongloucha
2,
Rapeepan Pitakaso
3 and
Paulina Golinska-Dawson
4
1
Research Unit on System Modelling for Industry, Department of Industrial Engineering, Faculty of Engineering, Khon Kaen University, Khon Kaen 40002, Thailand
2
Faculty of Economics, Khon Kaen University, Khon Kaen 40002, Thailand
3
Metaheuristics for Logistic Optimization Laboratory, Department of Industrial Engineering, Faculty of Engineering, Ubon Ratchathani University, Ubon Ratchathani 34190, Thailand
4
Faculty of Engineering Management, Poznan University of Technology, 60-965 Poznań, Poland
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(24), 4750; https://doi.org/10.3390/math10244750
Submission received: 22 November 2022 / Revised: 8 December 2022 / Accepted: 12 December 2022 / Published: 14 December 2022
(This article belongs to the Special Issue Planning and Scheduling in City Logistics Optimization)

Abstract

:
This study focused on analyzing planning and scheduling services in the tourism industry. Because dealing with these issues necessitates consideration of several important factors and stakeholders in the tourism business, it is challenging to operate resources efficiently. The purpose of this research is to propose a novel approach that allows maximizing the profits of tourism-related service sectors while considering many real-life constraints, such as sequence-dependent travel time, tourist time windows, points of interest, and specific destination constraints. We test our mathematical model for solving first small-scale problems and then metaheuristics proposed for finding a solution for real-life size problems. Moreover, sensitivity analysis was used to analyze the case study’s worthiness when the total cost and the revenue factor were changed. A real case study from Thailand’s Khon Kaen and Kanchanaburi provinces were used to verify the proposed models. The results indicate that the proposed models can be applied to investment decisions and strategy development. Furthermore, the outputs of the proposed models (i.e., the mathematical and metaheuristics models) can be employed to enhance the sustainability of other supply chains.

1. Introduction

Tourism is an essential and growing hospitality sector. It is a vital industry that drives economic and social growth worldwide [1]. Consequently, connection and coordination are critical for the most effective utilization of the available resources within the tourism supply chain. In addition, the various stakeholders in the tourism industry’s supply chain must be considered. This approach promotes tourism sustainability, employment opportunities, and the equal distribution of direct and indirect benefits from businesses to local areas. Tourism is a major driver of economic growth in Thailand. It is a considerable contributor to the nation’s gross domestic product (GDP) and a major employer. In 2019, the tourism industry recorded a contribution of 21.9 percent to the country’s GDP and accounted for more than 16.2 percent of all employees before the COVID-19 pandemic [2]. As a result of the virus, there have recently been fewer international tourists visiting the country. When the pandemic subsides, a community’s tourism business can bounce back by addressing new and exciting visitor requirements.
An important tool for maintaining the tourism industry’s sustainable growth is sharing tourism revenues with local businesses. Sustainable development is development that satisfies present needs without causing future problems. It is a type of development in which a given society, economy, and environment must be interconnected and interdependent. Managing tourism profits within a given area can help communities achieve financial sustainability, allowing them to expand their existing businesses or invest in new projects such as food and beverage, agriculture, and market initiatives to support increasing numbers of tourists. Effective community tourism business plans or commercial models require coordination across the network supply chains of an area’s various tourism businesses. All tourism supply chain stakeholders, including lodgings, tourist attractions, souvenir stores, and restaurants, must collaborate to maximize the potential of a tourism benefit-sharing scheme.
Such a scheme covers a wide variety of activities, including private, group, event, and tourist destinations or points of interest (POIs). Tourists often choose points of interest based on personal preferences, such as goodwill, customer reviews, and their price ranges [3]. Due to the different things that any one location can offer, tourists require there to be a particular activity that can only be carried out in the one location. Depending on the location, tourists’ closing time and sightseeing time are limited for each activity. An itinerary’s completion time must fall within the traveler’s time window. Each tourism company offers a series of diverse activities to service every kind of need of each tourist, such as shopping, rest, eating, visiting, sleeping, etc.
This paper carried out a case study focusing on Thailand’s tourism industry. The study dealt with the service scheduling problem for tourists who visit multiple points of interest, interact with multiple companies, engage in multiple activities, and have multiple lengths of stay. This problem was therefore formulated as a multi-period and multi-visit service scheduling problem with a time window (MMSSPTW). The study aimed to maximize the profits of tourism-related service sectors while considering many real-life constraints, such as sequence-dependent travel time, tourist time windows, points of interest, and specific destination constraints. These aspects make a MMSSPTW a difficult issue to solve manually for optimality. To address this issue, a mixed-integer linear programming model (MILP) was devised to solve a small-sized problem, while metaheuristics approaches were developed to solve various practical problems. Applying the proposed models to various tourism businesses would assist in the time-consuming planning process, raw material preparation, supplier coordination, and resource optimization among different supply chain participants. The solutions obtained from these approaches (i.e., the profits of each tourism business) were used to determine the NPV, IRR, B/C ratio, and payback period using the sensitivity analysis method, which attempts to complement the growing demands of businesses in any given tourism supply chain with the capability of the service sectors, while maintaining maximum tourist satisfaction.
The remainder of the paper is organized as follows. Section 2 discusses the existing literature. Section 3 describes the problem framework and the mathematical model proposed to solve the problem. In Section 4, metaheuristics are developed. A case study and computational results are presented in Section 5. In Section 6, a sensitivity analysis is performed to validate the model. Finally, Section 7 provides the conclusion and future research directions.

2. Literature Review

Tourism supply chain management (TSCM) is an important topic of focus in the tourism industry, which defines the tourism supply chain as a network of tourism businesses engaged in a variety of operations, with the provision of tourism services/diverse product aspects also accounted for. These networks comprise many different people from the public and private sectors [4]. Components of TSCM include the suppliers of products and services for visitors, such as transportation firms, attraction providers, accommodation suppliers, travel agents, memento stores, and so on. In [5], the authors stated that for supply chain efficiency, it is essential for the many tourism service providers to coordinate and collaborate with one another. The tourism industry has two main characteristics: the supply side and the demand side. Firstly, tourism goods are multifaceted and blend commodities and services. These goods are varied and complex, comprising lodging, commuting, touring, eating, and purchasing [6]. Secondly, tourism demand is characterized by extreme unpredictability and fragility, and is sometimes more volatile than demand in other sectors.
To provide tourism services to visitors, tourism business owners must plan and schedule services. Efficient tourist sequencing to minimize lost opportunity time is thus one of the essential variables in boosting the profitability of tourism businesses. In this study, this problem was defined as a multi-period and multi-visit service scheduling problem with a time window (MMSSPTW) to maximize the total profits of tourism-related service sectors while considering many real-life constraints, such as sequence-dependent travel time, tourist time windows, points of interest, and specific destination constraints. To the best of our knowledge, there are currently no existing publications on MMSSPTWs. These problems are classified as service scheduling problems (SSP), which involve scheduling arriving customers who demand services across a finite time horizon [7].
Numerous studies on SSPs have been undertaken, with these focusing on various attributes, including customer preference, travel time, planning horizon, time windows, time duration, servicing capacity, and location constraints. These attributes may result in distinct SSP characteristics. The research conducted in [8] is one example of a problem-related investigation. In their investigation of a multi-period service scheduling problem, their study examined two factors: the total number of ahead-of-time periods and operators. The objective was to determine the duration of service for each customer. In 2020, the authors of [7] focused on a logistics service scheduling problem to enhance the overall job delivery time. Recently, [9] addressed the patient transportation problem using a stochastic mixed-integer program, considering the constraints of customers, vehicles, and real-time handling to dynamically update planned routes and schedules. A sensitivity analysis is an analysis of a project’s worthiness when factors change, regardless of whether they increase or decrease; it assesses whether the project is worth further investment [10] if any changes in factors affect the Benefit/Cost Ratio (B/C Ratio), Internal Rate of Return (IRR) or Net Present Value (NPV). Investment value analysis has been applied in many industrial sectors, including consumer products, the agro-food industry, technology, property and construction, services, finance, and resources [11]. The authors of [12] designed a model for machine tool selection through utilizing various methods. A sensitivity analysis was used in their research to demonstrate how the model results varied in response to alterations in the relative importance of cost-related, technical, operational, and quality-related factors. Therefore, our study applied sensitivity analysis to analyze the opportunity cost incurred by each tourism business operator in the study area.
Based on the literature just reviewed, mathematical programming, heuristics, and metaheuristics were determined to be the most prevalent methods used to solve SSPs. However, the mathematical model has limitations in terms of its computational time, and heuristics may need help finding the optimal solution. Therefore, the use of heuristics is appropriate for simple issue features and a few linked components. In contrast, complicated issues or real-world scenarios need the employment of a metaheuristic, since the primary benefit of metaheuristic approaches is their ability to achieve near-optimal answers in an acceptable amount of time.
Therefore, this paper used the hybridization of a differential evolution algorithm and variable neighborhood search algorithm (HDEVNS), and also the hybridization of a traditional genetic algorithm and k-variable moves algorithm (HGAKV), which were firstly developed for solving the SSP with large-scale instances in order to maximize the profits of service operations in the tourism industry. Differential evolution (DE) is a population-based search algorithm that Storn and Price [13] initially proposed. The DE mechanism’s merits include its simplicity and ease of implementation, making it extremely popular among academics and practitioners, yet it frequently becomes trapped in the local optimum [14]. Consequently, the original structure was required to improve its performance. Various studies have investigated the DE approach, and it has been used extensively in multiple problems, such as a workforce scheduling and routing problem in a sugarcane mill [15], a multi-trip vehicle routing problem with backhauls and a heterogeneous fleet in the beverage logistics industry [16], a large-scale global black-box optimization problem [17], a cyclical multiple parallel machine scheduling problem in sugarcane unloading systems [18], and an employee transportation problem [19,20]. The variable neighborhood search algorithm is a metaheuristic that uses the idea of neighborhood change, which has more than one type of neighborhood structure, to systematically explore the solution space [18,21]. A genetic algorithm (GA) is a model that replicates natural evolution to identify approaches that combine stochastic and directed search to achieve a balance between exploitation and exploration of the search area [22]. Holland [23] introduced the GA concept in 1975. The benefits of GAs are their rapid calculation and simplicity. As a result, in this study, integer encoding was used to solve this problem since it is a simple procedure for genetic algorithms [24]. Still, its disadvantage is that local search is insufficient to identify the optimal solution. This paper used the K-Variable moves algorithm (KV) to enhance the performance of the traditional GA. The KV is a local search strategy [18]. Despite the availability of newer algorithms, GA are still widely used, particularly when combined with another algorithm [25].
The main contributions of this paper are threefold. Firstly, to the best of our knowledge, the literature has yet to discuss MMSSPTWs, which involve several visits and periods. Furthermore, our work considers the time windows of tourists, points of interest, location restrictions, and sequence-dependent travel time. Secondly, even though the mathematical model has limitations in terms of its computational time, it is often used to solve other kinds of problems, such as scheduling [8], allocation [26], manufacturing [27,28,29], planning [30], and routing problems [31,32,33]. MMSSPTWs have never been solved using such a paradigm. In this study, the mathematical model derived from [5] was presented as a solution to the MMSSPTW. Lastly, DE and GAs have never been used to solve MMSSPTWs, despite their extensive use in solving a variety of issues, as evidenced by a literature review of studies conducted on tourism-related difficulties. This research employed DE, a GA, hybrid differential evolution (HDEVNS), and a hybrid genetic algorithm (HGAKV) to solve a realistic-sized MMSSPTW.

3. Problem Description and Mathematical Model

We analyzed a multi-period and multi-visit service scheduling problem with a time window (MMSSPTW) in the local tourism industry. The problem framework is illustrated in Figure 1. Each tourist group n has a specific activity, and each activity o can be selected or determined based on which activity in their itinerary a tourist chooses. Each business k has a limited capacity Qk to serve tourists. Each group’s activities o are carried out in the time frame set by the tourist, such as the starting time, finishing time, earliest service time of activity j at the business k, and latest service time of activity j at the business k.
The following assumptions were considered to develop a mathematical model:
  • At time zero, all tourist groups and locations are available.
  • Each location can only accommodate one tourist group at a time.
  • Each activity has only one location constraint.
  • Each tourist group visits only one location at a time.
  • The tourists’ sightseeing time correlates.
  • Activities and locations are assigned ahead of time, and no preemption of activities is permitted.
  • Activities and locations are announced in advance, and there is a strict prohibition on any form of activity preemption; attendees are not permitted to disrupt ongoing activities until they have been completed.

Mathematical Model

This section describes the proposed MILP for the MMSSPTW. The indices, sets, parameters, and decision variables used in the MILP model are defined as follows.
 Indices and sets:
i,pTourist group index;
NSet of tourist groups;
j,lActivity index;
OiSet of activities for tourist group i;
k,hBusiness index;
KSet of businesses.
 Parameters:
A k i j Sightseeing time of business k for activity Oij (unit: hour);
T h k Travel time between the business h to the business k (unit: hour);
V Positive large number;
M R k i j Tourist group i activity j performs at the business k;
Q k Maximum capacity of the business k (unit: person);
D i Number of tourists in group i (unit: person);
B i Starting time of tourist group i (unit: hour);
U i Finishing time of tourist group i (unit: hour);
E k i j Earliest service time of business k for activities Oij (unit: hour);
F k i j Latest service time of business k for at activities Oij (unit: hour);
R V j k The revenue (collection fee per visit) of business k for activities Oij (unit: THB/person);
C C j k The cost (staff cost, service cost and administration cost) of business k for activities Oij (unit: THB/person).
 Decision variables:
S T k i j Starting time of tourist group i activity j at the business k (unit: hour);
C T k i j Completion time of tourist group i activity j at the business k (unit: hour);
X k i j = 1, if tourist group i activity j is assigned at the business k;
= 0, otherwise;
Y k i j p l = 1, if tourist group i activity j precedes activity l at the business k;
= 0, otherwise;
W k h i j = 1, if tourist group i activity j is processed at business k and activity j-1 at the business h;
= 0, otherwise.
 Objective function:
M A X       Z = i N j O i k K D E M i ( R V j k C C j k ) X k i j         (1)
The objective function is used to maximize businesses’ profits with tourist satisfaction.
 Subject to:
k K X k i j 1 i , j (2)
Equation (2) ensures that the activity is assigned to at least one business.
S T k i j   +   C T k i j     X k i j V i , j , k (3)
C T k i j       S T k i j + A k i j ( 1 X k i j )   V i , j , k (4)
Equations (3) and (4) ensure that the amount of time between the starting and finishing times must at least equal the amount of time spent sightseeing at business k.
S T k i j     C T k p l   Y k i j p l V i , j , p i , l , k (5)
S T k p l     C T k i j   ( 1 Y k i j p l )   V i , j , p i , l , k (6)
Y k i j p l + Y k p l i j       1 i , j , p i , l , k (7)
Equations (5)–(7) ensure that tourist i and activity j and tourist p and activity l cannot be carried out simultaneously at any business.
M R k i j X k i j 0 i , j , k (8)
i N j O i X k i j 0 k (9)
Equations (8) and (9) guarantee the feasibility of the businesses that are demonstrated for any activity.
k K S T k i 1     k K k > 1 W k 1 i 1 T 1 k i (10)
k K S T k i j     k K k > 1 C T k i ( j 1 ) + k K k > 1 h K T h k W k h i j i , j > 1 (11)
Equations (10) and (11) are sequencing constraints that take travel time into account.
k K h K h > 1 W k h i j 1 i , j > 1 (12)
k K W k 1 i 1 1 i (13)
Equations (12) and (13) guarantee that no action can include more than one of the businesses.
h K W k h i j = X k i j i , j , k (14)
k K W k h i j   k K W k h i ( j 1 ) i , j > 1 , h (15)
Equations (14) and (15) ensure that each activity’s business is chosen from the available alternatives.
D i X k i j   Q k i , j , k (16)
Equation (16) ensures that the number of tourists is fewer than the limit of business k.
B i k K S T k i 1 i (17)
U i k K C T k i j i , j (18)
E k i j X k i j   S T k i j i , j , k (19)
F k i j X k i j   C T k i j i , j , k (20)
Equations (17)–(20) ensure that the completion and starting times fall within the time windows of the business’s and the tourist group’s goals.
S T k i j , C T k i j 0 i , j , k (21)
X k i j , W k h i j , Y k i j p l { 0 , 1 } i , j , p , l , k (22)
Equation (21) is a continuous decision variable constraint, while Equation (22) is a binary decision variable constraint.

4. Development of Metaheuristics

When the size and complexity of problems grow too great for exact solution methods, metaheuristics are often used. This study used a metaheuristic employing a Genetic Algorithm (GA) and Differential Evolution (DE). To increase the exploitation and exploration capabilities of the search and control strategy, we designed a hybrid approach based on a GA (called HGAKV) and a hybrid based on DE (called HDEVNS). Exploration and exploitation are global and local searches from a latent perspective. The traditional DE and GA have strong global optimizing exploring capabilities. However, preserving motions to a local optimum still lacks exploitation potential. Therefore, the K-Variable moves algorithm (KV) and Variable Neighborhood Search algorithm (VNS) were introduced to increase our ability to search for the best solutions using the local search approach.

4.1. The Traditional Differential Evolution (DE) and Traditional Genetic Algorithm (GA)

4.1.1. The Traditional DE

Four processes comprise the DE algorithm: initial solution, mutation, recombination, and selection. Xji,G represents the vector containing dimension (D) variables. The population size or the number of vectors utilized in a differential evolution iteration is determined by the number of samples in the population. The first iteration will generate a random population size, and each vector will be identical. Then, the vectors will be changed into mutant and trial vectors using mutation and recombination processes.

4.1.2. The Traditional GA

The overall structure of a GA includes operators comprising an initial solution, crossover, mutation, and selection. In this study, GA operators were utilized to deal with this issue, including weight mapping crossover (WMX), swapping mutation, and the elitism method for selection.

4.2. Hybrid Differential Evolution Algorithm with Variable Neighborhood Search Algorithm (HDEVNS)

In this study, a Variable Neighborhood Search algorithm (VNS) was employed to develop a hybrid with the DE algorithm. The VNS is a metaheuristic that uses the idea of neighborhood change, which has more than one type of neighborhood structure, to systematically explore the solution space. The conceptual properties of the neighborhood descend to local minima and escape from valleys that contain them. The insert and swap neighborhoods are two different neighborhood structures that were applied in this study. Swap or insert operators are random at any position in the sequence. In Algorithm 1, the VNS pseudocode is illustrated. Nk is the neighborhood structures (k = 1, 2, …, kmax). The collection of solutions inside the kth neighborhood of s is represented by Nk(s), which is a solution set. A random solution, s′, is generated for each iteration using the Nk(s). A second solution, designated s″, is then produced via a local search for solution s′. If a solution is not produced, the procedures are repeated using the following neighborhood, k = k + 1. The old solution is updated, and the value of k is set to 1 if the new solution, s″, is superior to the previous one. The process continues until the termination condition is met.
Algorithm 1. Variable Neighborhood Search algorithm (VNS)
1Initialize solution s
2For k = 1 to kmax
3Shaking: generate at random a starting solution s’ ∈ Nk(x)
4Local search: apply a local search from the starting solution x′ using the base random neighborhood structure s
 (optional) insert algorithm
 (optional) swap algorithm
5Improve or not:
6 If s″ is better than s
7  s = s
8 End If
9Neighborhood change: k = k + 1
10End for
11Output: s
The structure of the hybrid DE algorithm with the VNS algorithm consists of an initial solution, mutation process (Equation (23)), recombination process (Equation (24)), and selection process (Equation (25)). The VNS algorithm was used after the crossover operator because we wanted to increase our chances of obtaining better results in the search area according to [5]. Algorithm 2 illustrates the pseudocode for the HDEVNS algorithm.
Algorithm 2. Hybrid Differential Evolution algorithm with Variable Neighborhood Search algorithm (HDEVNS)
1Input: MMSSPTWdata, DE parameters (NP, CR, F)
2Output: the best solution
3Begin Generate the initial population of NP
4while termination condition is not satisfied do
5for i = 1 to NP //NP is the predefined number of population
6  Perform mutation process using Equation (23)
7  Perform recombination process using Equation (24)
8  Perform VNS algorithm
9  Evaluate objective value of vectors
10  Perform selection process using Equation (25)
11  Update best solution
12end for
13end while
14end

4.2.1. Initial Solution

DE, an approach for solving some tourism problems, was used by Thumrongvut et al. [5]. They used two elements in each vector to denote the activity sequence and the place selection. The population number is the population size or the number of vectors, where each vector consists of the activity of a tourist group vector and a place vector. Firstly, each vector is a randomly generated real number equal to the number of activity counts. The vector’s dimension corresponds to the group’s serial number, and its sequence corresponds to the group’s activity number. Then, the decoding of each dimension sorts the rank order value (ROV) of each vector in ascending order so that it obtains the sequence of groups and activities to be operated. In addition, the place selection, or k selection (KS), is an alternative set in individual groups and activities. The internal KS number represents the serial number of the optional places, corresponding to the place selected by the matching activity sequence operation.

4.2.2. Mutation Process

The mutation process transforms the target vector into the mutant vector. Equation (23) is used in the generating process.
V j i , G + 1 = X r 1 , G + F ( X r 2 , G X r 3 , G )
where Vji,G+1 is called the mutant vector; i represents the vector number (i = 1, 2, …, NP); j is the position of a vector (j = 1, 2, …, D); F is a scaling factor, which is set as 0.8 [34]; r1, r2, and r3 are the indices of randomly selected vectors.

4.2.3. Recombination Process

The recombination process transforms the mutant vector into the trial vector. Equation (24) is used in the recombination or crossover process. Each position value in a vector may correlate to a target or trial vector depending on the random number (Rand) generated for that position and compared to the crossover rate (CR), a preset parameter set to 0.8 in this study [34].
U j i , G + 1 = { V j i , G + 1     i f   ( R a n d ( j ) C R ) X j i , G       i f   ( R a n d ( j ) > C R )
where Xji,G is the target vector; Uji,G+1 is the trial vector; Vji,G+1 is the mutant vector.

4.2.4. Selection Process

The selection process is used to select the new target vector. We compared the current trial vector fitness function and the current target vector fitness function with the highest objective value selected for the next iteration (Xji,G+1) using Equation (25).
X j i , G + 1 = { U j i , G           i f       f ( U j i , G + 1 ) < f ( X j i , G ) X j i , G         o t h e r w i s e

4.3. Hybrid Genetic Algorithm with K-Variable Moves Algorithm (HGAKV)

The KV algorithm is a local search-based heuristic that extends the swap method, in which k travels from one position to the next k positions until the final k in the chromosome arrives at the first k. We set the light and strong motion percentages to 20% and 80%, respectively. We used k = 20% in this study, which indicates that k should be the random number for transferring up to roughly 20% of the tourist groups. Figure 2 shows a KV example where k equals 3.
The hybrid GA with KV algorithm structure consists of an initial solution, crossover operator, mutation operator, and selection using elitism. Because we wanted to eliminate duplicate algorithm solutions, we employed the KV algorithm instead of the mutation operator. This method considerably reduces computational time while producing superior results. The detailed HGAKA algorithm is shown in Algorithm 3.
Algorithm 3. Hybrid Genetic Algorithm with K-Variable moves algorithm (HGAKV)
1Input: MMSSPTWdata, GA parameters (popSize, maxGen, pm, pc)
2Output: the best solution
3Begin t ← 0
 Generate the initial populations
 Evaluate objective value of populations
4while tmaxGen
5 Generate F(t) from P(t) by weight mapping crossover (WMX)
6 Generate F(t) from P(t) by swap mutation operation or KV algorithm
7 Evaluate objective value of F(t)
8 Select P(t+1) from P(t) and F(t) by elitist selection operation
9tt + 1
10end while
11end

4.3.1. Initial Solution

The initial population is crucial for the GA since it directly affects the convergence rate of fitness values and the quality of solutions. The initial chromosomal population is generated randomly as an integer number based on the activity number in each chromosome, since the genetic algorithm’s conventional operators easily convey it. There are two elements to consider: the Activity sequence (AS) and the k selection (KS). Figure 3 displays an example of the encoding procedure. In the decoding procedure, chromosomes are decoded to evaluate fitness under the constraints of the multi-period and multi-visit service scheduling problem with a time window (MMSSPTW) and taking into account conditions such as the time windows of tourists and locations, sequence-dependent travel time, and destination restrictions. The operation can be described as follows.
Step 1:
Construct a random initial solution based on the activity number in the chromosome, which is composed of the group numbers.
Step 2:
Consider the following tourist and location constraints:
  • Examine the time window restrictions of existing tourists and locations.
  • Evaluate the present capacity restrictions of the location.
  • Update the sequence of group I in the set number population.
  • Repeat until the sequence is completed.
Step 3:
Check the cumulative tourists in the location. If the operation locations are not in service, the defined sequence must be repaired, and then step 2 must be returned to.
Step 4:
Update capacity in locations.
Step 5:
Update the time window of tourists.
Step 6:
Evaluate fitness using the objective function.
Step 7:
Select the best solution.

4.3.2. Crossover Operator

The crossover operator facilitates the transfer of genetic information across two chromosomes from one generation to the next. This study employed a precedence weight-mapping crossover (WMX) for the operation sequence. According to [24,25], the WMX has significant effects when used for crossover operations. The WMX determines a one cut-point randomly and then exchanges substrings between the parents for both the tourist group priority and activity tourism type. To avoid the repairing process after a crossover, the weight of the ascending number segment is mapped, and an offspring is formed through a mapping relationship, as shown in Figure 4.

4.3.3. Mutation Operator

The mutation operator can expand the population’s variety by providing additional variability to avoid the local optima phenomenon. A swap mutation was applied in this study. A cut-point is selected randomly and assigned as a mutation point for two allele values in a string, similar to the crossover technique. If the random value from a uniform number falls within the mutation rate range, two allele values inside the chromosome will swap. Figure 5 shows an example of a swap mutation.

4.3.4. Selection Operator

The selection operation aims to develop a new population with greater fitness values than the current population. In this study, an elitism method was used, due its incredible effect on selection operation, according to [16]. A limited number of individuals with the highest fitness values are selected to continue to the next iteration in the elitist selection strategy.

4.4. Current Practice (CP)

Thailand’s tourism industry uses a constructive heuristic based on the first-come, first-served (FCFS) principle for scheduling travelers’ routes. FCFS is an inefficient method that leads to poor management practices in the tourism industry, negatively impacting income and tourist satisfaction. Due to significant constraints (such as time windows, location eligibility, and travel time between areas), the issue has become too complex for the authorities to resolve.

5. Case Study

In this paper, the issue under consideration was a real-life case study of a network supply chain of various tourism businesses. The Khon Kaen and Kanchanaburi provinces of northeastern and central Thailand were selected for the case study. Khon Kaen tourism is characterized by many tourism service providers, many cultural and historical sights, high-quality silk, leisure shopping, beautiful nature, and local cuisine, and various unique characteristics of the province. Likewise, Kanchanburi is a popular vacation spot for nature enthusiasts. As a place to unwind, the area allows visitors to view the water from the river’s banks, either inside a raft house or from a riverside restaurant. Along with witnessing stunning waterfalls, lush woods, and three of Thailand’s largest reservoirs, white water rafting, golf, and elephant trekking are activities popular with travelers. Furthermore, many people enjoy relaxing beside the river, which provides pleasant relief from the heat and congestion of the surrounding capital town.
Table 1 shows the details of the Khon Kaen tourists in each group and the value in parentheses is the business number. For example, as seen in Table 1, tourist group 1 had a sightseeing time (Aijk) of 2.00 h for activity 1 at business no. 1, 1.00 h for activity 2 at business no. 7, 2.50 h for activity 3 at business no. 12, and 0.50 h for activity 4 at business no. 15. The group’s time window ([Bi, Ui]) was between 2.00 and 8.00 h. Each business had a limited capacity (Qk) to serve tourists. The revenue (collection fee or income per visit per person) and various business costs (staff costs, service costs, and administration costs) (THB/person) of the businesses are shown in Table 2. As shown in Table 2, the case study featured four restuarants, two hotels, and nine tourist attractions or POIs. The travel time (Thk) between any two points was not a symmetric function; the time between points k and h differed from the time between points h and k. In addition, the number of tourists in each group was 15, 30, 45, 10, 4, and 20 persons, respectively. Figure 6a,b show various businesses in the Khon Kaen and Kanchanaburi provinces, respectively.
According to the data collected, restaurant income ranged between 70 and 250 THB per person, accommodation income ranged between 200 and 700 THB per person, and attraction income ranged between 10 and 200 THB per person. Therefore, if its service was fully utilized, a given business generated a certain amount of income per person, as shown in Table 3. Similarly, Table 4 presents the total costs, which were the proportions of the revenue of the restaurant, accommodation, and tourist attraction or POI businesses, which were 30–60%, 30–40%, and 30–60%, respectively.
Making a profit is the primary goal of any business in any industry, as it indicates that limited resources, such as humans, materials, and facilities, are being used to their full potential. The total profits of the businesses in the case study in Khon Kaen province were computed by multiplying their revenue by the number of customers minus the total costs by the number of customers. The profits of each business are shown in Table 5. The objective value aimed to maximize the total profits, which, for all the activities of all the tourist groups and businesses, was calculated to be 24,031.8 THB. When business owners increase their income through effective planning and management, increased profits and returns are obtained, contributing to better long-term stability.

6. Computational Results

The proposed methods were coded in MATLAB software (R2021a), and the mathematical model was executed using LINGO (v.18) for MILP on an Intel® Core™ i7-1165G7 processor 2.80 GHz, with 16.0 GB of RAM, running on Windows 10. The input and algorithm parameters were set according to the preliminary test, as summarized in Table 6. The experiments were tested with 20 randomly generated data sets, which comprised small sized problems (instance no. 1–5), medium sized problems (instance no. 6–12), and large sized problems (instance no. 13–20). All details of the test instances are shown in Table 7. The computational time of the mathematical model can be divided into two categories. For the first category, we ran LINGO until it discovered the optimal solution, and then reported the computational time. The second computational time was employed when LINGO could not identify the optimal solution within 28,800 s. At this point, we halted the process and reported the best solution obtained by LINGO within 28,800 s. Due to LINGO’s inability to find the optimal solution within the acceptable computing time for of the medium and large sized problems, the best and upper-bound solutions were used for comparison with the proposed methods.
All algorithms’ termination conditions were determined according to the problem size. The number of iterations was used as the stopping criteria for the small sized problems (set to be 300 iterations). In contrast, for the medium and large sized problems, the computational time was used, which was fixed to be 60 and 120 s, respectively. The experiment was repeated five times, and the average solution was reported. Computational results are shown in Table 8 and Figure 7. It can be seen that the mathematical model was able to find optimal solutions for instances no. 1 to 5, while for instances no. 6 to 12 and 13 to 20, the best objective was found and the upper bound solution was found during the limited computation time, respectively. The numerical experiment results show that our proposed HGAKV method outperformed the current practice, the traditional DE, the traditional GA, and the HDEVNS method. The results were statistically tested and are displayed in Table 9 using the Wilcoxon sign rank test with a 95% confidence interval, obtained using IBM SPSS Software (V.28). The statistical test indicates that the solutions obtained from the HGAKV method were significantly different at a p-value ≤ 0.05 compared to those obtained from all the proposed methods.
Based on the results shown in Table 8, the heuristic performance (HP) was measured by comparing the solution obtained from all the compared algorithms with the solution obtained from the mathematical model. Another factor examined in this study was the relative improvement (RI), which was compared to the total profit of the solution from all the algorithms.
H P ( % ) = S o l u t i o n a lg S o l u t i o n M I L P × 100
R I ( % ) = S o l o t i o n a lg - S o l o t i o n o l d S o l o t i o n a lg × 100
where
  • SolutionMILP = the solution of the mixed integer linear programming;
  • Solutionold = the solution of the current practice;
  • Solutionalg = the solution of the current practice, DE, HDEVNS, GA, and HGAKV.
The HP of all the proposed methods was calculated using Equation (26), and the RI using Equation (27). As shown in Table 10, the HP revealed that the HGAKV obtained a near-optimal result, recording an average of 89.24%. The RI was compared to the total profit of the solution obtained using the current practice, DE, HDEVNS, GA, and HGAKV methods. As shown in Table 10, these algorithms had average RIs of 29.04%, 30.80%, 29.48%, and 35.20%, respectively. Additionally, the HGAKV improved the solution quality by 9.75%, 6.82%, and 8.62% with regard to the DE, GA, and HDEVNS, respectively.
The hybrid algorithm (HGAKV)’s advantages are the combined advantages of extensive exploration in the search space, deriving from the GA algorithm, and high exploitation, deriving from the KV algorithm. Consequently, in most cases, the HGAKV produced the best solution, outperforming the other approaches and significantly improving a given company’s total profit. In contrast, the current practice approach recorded the lowest quality result among all the approaches.

Sensitivity Analyses

The profits of a tourism supply chain will increase if more demand is achieved, which is accomplished when available points of interest (restaurants, accommodation, shops, tourist attractions, etc.) can accommodate an increase in tourist demand. One way to increase the demand response is by decreasing the total costs of venues or by reducing the opportunity costs of businesses by planning and scheduling servicing effectiveness. Governments can boost a nation’s GDP and foreign exchange by investing in the tourism industry, which is economically significant for developing nations. Further, governments and the private sector must collaborate to develop profitable and sustainable growth plans. To encourage the participation of the private sector in financing, constructing, and managing infrastructure and other related development projects, a combination of policy reforms, institutional support, incentives, and financing modalities are used to foster public–private partnerships in the tourism industry. According to the references discussed above, an alternative to budgetary limitations enables the tourism industry to build assets and maximizes the use of private sector skills. The potential benefits of a public–private partnership for the private sector include, by the partnership’s very nature, risk sharing; a public–private partnership involves the sharing of project risks and cost reduction, which encourages private investment in tourism.
In sensitivity analysis approaches, the specific parameters of the objective function are given the most weight overall. The majority of these have a high to extremely high sensitivity. Sensitivity analysis determines how much the decision-making evaluation index changes when one or more uncertainties are changed. Through this, the degree to which the adjustments impact the achievement of the targeted outcomes is determined. Then, judgments about the plan’s capacity can be formed when the external environment changes negatively. In this study, sensitivity analysis was employed to analyze the case study’s worthiness when the total cost and the revenue factor were changed in terms of the NPV, IRR, B/C Ratio, and payback period. In the method, two parameters (revenue and cost) were varied simultaneously, but the range of variation was defined differently. The similar outcomes obtained show that the parameter sensitivity was able to be measured without the chosen method’s impacting the outcomes.
Table 11 presents the effects on a business’s profit when their revenue and costs were changed by 5% and 10%. The three examples shown are business no. 3, a restaurant, business no. 6, a hotel, and business no.10, a tourist attraction, with each business example having eight scenarios. For scenario 1 of business no. 3, increasing its revenue and costs by 5% was found to increase its profit by 5% (3543.75THB). Thus, profit was found to be significantly sensitive to changes in the revenue and cost values. Another way to formulate this is that the business’s profit was found to grow by 1% for every 1% increase in revenue and costs, or fall by 1% for every 1% decrease in revenue and cost (see scenario 4 of the restaurant). Similarly, in scenario 2, if the business’s revenue was increased by 5% and its costs were decreased by 5%, its profit was increased by 15% (3881.25THB). On the other hand, if its revenue was decreased by 5% and its costs were increased by 5%, its profit was decreased by 15% (scenario 3). Scenarios 5 to 8 for each business show the impact of changing the business’s revenue and costs by 10%. For example, scenario 7 of business no. 3 shows that decreasing its revenue by 10% and increasing its cost by 10% resulted in a 30% (3712.50THB) decrease in its profit. Consequently, profit was found to be likewise quite sensitive to shifts in these parameters. To state this another way, every 1% decrease in these parameters was found to result in a 3% decrease in profit, or every 1% increase in these parameters was found to result in a 3% gain in profit. Comparing the changed factors in the four scenarios found that the profits of each of the example businesses were more sensitive to changes that involved revenue increases, which lead to a positive objective value or increases in profits than to changes that involved revenue decreases. However, all changes were found to have a significant impact on profit.
Based on the sensitivity analysis, company owners or stakeholders may need to devise strategies to boost income and eliminate expenditures to ensure a business’s sustainability. In addition, the strength of sensitivity analysis is that it determines how changes in variables will affect a business’s profits, and this approach can also determine economic indicators. For example, the economic indicators of business no. 3 were as follows: its Net Present Value (NPV) was 13,542,512.92 THB, its Internal Rate of Return (IRR) was 59.44%, its Benefit/Cost Ratio (B/C Ratio) was 2.86 %, and its payback period was nine months, indicating that the case study was worth the investment. Since the NPV and the B/C Ratio values were greater than zero, the IRR was greater than the Minimum Attractive Rate of Return (MARR) (set to be 20%), and the payback period was faster than the target payback period.

7. Conclusions

This paper addresses a multi-period and multi-visit service scheduling problem with a time window (MMSSPTW) in local tourism businesses. The study intended to maximize the total profits of service operations in the tourism industry within service hours by considering the time windows of tourists and locations, sequence-dependent travel time, and destination restrictions. The proposed methods were developed to solve a real-world case study of a local tourism business network supply chain. This case study featured Khon Kaen and Kanchanaburi, Thailand. A mixed-integer linear programming model (MILP) was devised to solve the small-sized problems, while metaheuristics approaches were developed to obtain solutions for the practical problems. In the study, the DE, GA, HDEVNS, and HGAKV algorithms were employed.
Numerical results from these metaheuristics were investigated and compared to a set of constructed instances to illustrate the algorithms’ benefits in terms of their solution quality. For a significant number of instances and real-world issues, the HGAKV algorithm outperformed the DE, GA, and HDEVNS algorithms in terms of its solution quality under identical experimental conditions. The results obtained in the present work show the service schedule of each tourism business with the sequence of tourists, including the number of tourists, time duration, activity details, the time windows of tourists, and the tardiness time of a given location. The solution obtained from these approaches (i.e., the profit of each tourism business) was used to determine the NPV, IRR, B/C ratio, and payback period using the sensitivity analysis method, which attempts to complement the growing demands of businesses in a tourism supply chain with the capability of the service sectors while maintaining maximum tourist satisfaction.
In future research, we plan to consider multiple objective functions and vehicle conditions in transport modes. We believe that the proposed model can solve complicated combinatorial optimization problems in other industries (e.g., hospitality, hotel, transportation, and food and beverage industries). Although the proposed method is beneficial, additional work should be conducted to discover solutions by implementing other hybrid methods to compare the merits of other approaches in tackling issues of this nature.

Author Contributions

Conceptualization, P.T. and K.S.; methodology, P.T.; software, T.J., P.G.-D. and C.W.; validation, P.T. and P.G-D.; formal analysis, K.S., C.W. and R.P.; writing—original draft, P.T., K.S. and T.J.; writing—review and editing, K.S. and R.P.; supervision, K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by the Research and Graduate Studies Khon Kaen University, Thailand (grant number RP66-1-004) and the Research Unit on System Modeling for Industry, Department of Industrial Engineering, Faculty of Engineering, Khon Kaen University, Thailand.

Data Availability Statement

Not applicable.

Acknowledgments

We also thank Somnuk Theerakulpisut for his critical review of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, K.; Wu, D.; Zou, Y. Tourist Arrival Forecasting Using Multiscale Mode Learning Model. Mathematics 2022, 10, 2999. [Google Scholar] [CrossRef]
  2. Department of Tourism the Tourism Development Strategic Plan 2018–2021. Available online: https://checkinthailand.net/wp-content/uploads/2021/02/Plan-department-of-tourism_compressed-1.pdf (accessed on 10 July 2022).
  3. Díaz, M.R.; Espino-Rodríguez, T.F. Determining the Sustainability Factors and Performance of a Tourism Destination from the Stakeholders’ Perspective. Sustainability 2016, 8, 951. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, X.; Song, H.; Huang, G.Q. Tourism Supply Chain Management: A New Research Agenda. Tour. Manag. 2009, 30, 345–358. [Google Scholar] [CrossRef] [Green Version]
  5. Thumrongvut, P.; Sethanan, K.; Pitakaso, R.; Jamrus, T.; Golinska-Dawson, P. Application of Industry 3.5 Approach for Planning of More Sustainable Supply Chain Operations for Tourism Service Providers. Int. J. Logist. Res. Appl. 2022. [Google Scholar] [CrossRef]
  6. Ahmadimanesh, F.; Paydar, M.M.; Asadi-Gangraj, E. Designing a Mathematical Model for Dental Tourism Supply Chain. Tour. Manag. 2019, 75, 404–417. [Google Scholar] [CrossRef]
  7. Zhou, L.; Zhang, L.; Fang, Y. Logistics Service Scheduling with Manufacturing Provider Selection in Cloud Manufacturing. Robot Comput. Integr. Manuf. 2020, 65, 101914. [Google Scholar] [CrossRef]
  8. Fernández, E.; Kalcsics, J.; Núñez-del-Toro, C. A Branch-and-Price Algorithm for the Aperiodic Multi-Period Service Scheduling Problem. Eur. J. Oper. Res. 2017, 263, 805–814. [Google Scholar] [CrossRef] [Green Version]
  9. Yu, X.; Shen, S.; Wang, H. Integrated Vehicle Routing and Service Scheduling Under Time and Cancellation Uncertainties with Application in Nonemergency Medical Transportation. Serv. Sci. 2021, 13, 172–191. [Google Scholar] [CrossRef]
  10. Frey, H.C.; Patil, S.R. Identification and Review of Sensitivity Analysis Methods. Risk Anal. 2002, 22, 553–578. [Google Scholar] [CrossRef]
  11. Dhami, I.; Deng, J.; Strager, M.; Conley, J. Suitability-Sensitivity Analysis of Nature-Based Tourism Using Geographic Information Systems and Analytic Hierarchy Process. J. Ecotourism 2017, 16, 41–68. [Google Scholar] [CrossRef]
  12. Özgen, A.; Tuzkaya, G.; Tuzkaya, U.R.; Özgen, D. A Multi-Criteria Decision Making Approach for Machine Tool Selection Problem in a Fuzzy Environment. Int. J. Comput. Intell. Syst. 2011, 4, 431–445. [Google Scholar] [CrossRef]
  13. Storn, R.; Price, K. Differential Evolution—A Simple and Efficient Heuristic for Global Optimization over Continuous Spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  14. Zheng, L.; Luo, S. Adaptive Differential Evolution Algorithm Based on Fitness Landscape Characteristic. Mathematics 2022, 10, 1511. [Google Scholar] [CrossRef]
  15. Punyakum, V.; Sethanan, K.; Nitisiri, K.; Pitakaso, R.; Gen, M. Hybrid Differential Evolution and Particle Swarm Optimization for Multi-Visit and Multi-Period Workforce Scheduling and Routing Problems. Comput. Electron. Agric. 2022, 197, 106929. [Google Scholar] [CrossRef]
  16. Sethanan, K.; Jamrus, T. Hybrid Differential Evolution Algorithm and Genetic Operator for Multi-Trip Vehicle Routing Problem with Backhauls and Heterogeneous Fleet in the Beverage Logistics Industry. Comput. Ind. Eng. 2020, 146, 106571. [Google Scholar] [CrossRef]
  17. Vakhnin, A.; Sopov, E.; Semenkin, E. On Improving Adaptive Problem Decomposition Using Differential Evolution for Large-Scale Optimization Problems. Mathematics 2022, 10, 4297. [Google Scholar] [CrossRef]
  18. Kusoncum, C.; Sethanan, K.; Pitakaso, R.; Hartl, R.F. Heuristics with Novel Approaches for Cyclical Multiple Parallel Machine Scheduling in Sugarcane Unloading Systems. Int. J. Prod. Res. 2021, 59, 2479–2497. [Google Scholar] [CrossRef]
  19. Pitakaso, R.; Sethanan, K.; Srijaroon, N. Modified Differential Evolution Algorithms for Multi-Vehicle Allocation and Route Optimization for Employee Transportation. Eng. Optim. 2020, 52, 1225–1243. [Google Scholar] [CrossRef]
  20. Supattananon, N.; Akararungruangkul, R. Modified Differential Evolution Algorithm for a Transportation Software Application. J. Open Innov. Technol. Mark. Complex. 2019, 5, 84. [Google Scholar] [CrossRef] [Green Version]
  21. Fleszar, K.; Osman, I.H.; Hindi, K.S. A Variable Neighbourhood Search Algorithm for the Open Vehicle Routing Problem. Eur. J. Oper. Res. 2009, 195, 803–809. [Google Scholar] [CrossRef]
  22. Cho, D.H.; Moon, S.H.; Kim, Y.H. Genetic Feature Selection Applied to Kospi and Cryptocurrency Price Prediction. Mathematics 2021, 9, 2574. [Google Scholar] [CrossRef]
  23. Holland, J.H. Adaptation in Natural and Artificial Systems; University of Michigan Press: Ann Arbor, MI, USA, 1975. [Google Scholar]
  24. Jamrus, T.; Chien, C.F. Extended Priority-Based Hybrid Genetic Algorithm for the Less-than-Container Loading Problem. Comput. Ind. Eng. 2016, 96, 227–236. [Google Scholar] [CrossRef]
  25. Kamsopa, K.; Sethanan, K.; Jamrus, T.; Czwajda, L. Hybrid Genetic Algorithm for Multi-Period Vehicle Routing Problem with Mixed Pickup and Delivery with Time Window, Heterogeneous Fleet, Duration Time and Rest Area. Eng. J. 2021, 25, 71–86. [Google Scholar] [CrossRef]
  26. Ma, J.; Webb, T.; Schwartz, Z. A Blended Model of Restaurant Deliveries, Dine-in Demand and Capacity Constraints. Int. J. Hosp. Manag. 2021, 96, 102981. [Google Scholar] [CrossRef]
  27. Rahman, M.S.; Das, S.; Manna, A.K.; Shaikh, A.A.; Bhunia, A.K.; Cárdenas-Barrón, L.E.; Treviño-Garza, G.; Céspedes-Mota, A. A Mathematical Model of the Production Inventory Problem for Mixing Liquid Considering Preservation Facility. Mathematics 2021, 9, 3166. [Google Scholar] [CrossRef]
  28. Qiu, S.; Ming, X.; Sallak, M.; Lu, J. Joint Optimization of Production and Condition-Based Maintenance Scheduling for Make-to-Order Manufacturing Systems. Comput. Ind. Eng. 2021, 162, 107753. [Google Scholar] [CrossRef]
  29. Srichok, T.; Pitakaso, R.; Sethanan, K.; Sirirak, W.; Kwangmuang, P. Combined Response Surface Method and Modified Differential Evolution for Parameter Optimization of Friction Stir Welding. Processes 2020, 8, 1080. [Google Scholar] [CrossRef]
  30. Zhang, H.; Qiu, J.; Wang, Y. Planning Strategy of Fast-Charging Stations in Coupled Transportation and Distribution Systems Considering Human Health Impact. Int. J. Electr. Power Energy Syst. 2021, 133, 107316. [Google Scholar] [CrossRef]
  31. Delorme, M.; Iori, M.; Mendes, N.F.M. Solution Methods for Scheduling Problems with Sequence-Dependent Deterioration and Maintenance Events. Eur. J. Oper. Res. 2021, 295, 823–837. [Google Scholar] [CrossRef]
  32. Ferreira, K.M.; de Queiroz, T.A.; Toledo, F.M.B. An Exact Approach for the Green Vehicle Routing Problem with Two-Dimensional Loading Constraints and Split Delivery. Comput. Oper. Res. 2021, 136, 105452. [Google Scholar] [CrossRef]
  33. Chokanat, P.; Pitakaso, R.; Sethanan, K. Methodology to Solve a Special Case of the Vehicle Routing Problem: A Case Study in the Raw Milk Transportation System. AgriEngineering 2019, 1, 75–93. [Google Scholar] [CrossRef] [Green Version]
  34. Qin, A.K.; Huang, V.L.; Suganthan, P.N. Differential Evolution Algorithm with Strategy Adaptation for Global Numerical Optimization. IEEE Trans. Evol. Comput. 2009, 13, 398–417. [Google Scholar] [CrossRef]
Figure 1. Illustration of the MMSSPTW.
Figure 1. Illustration of the MMSSPTW.
Mathematics 10 04750 g001
Figure 2. Example of KV algorithm.
Figure 2. Example of KV algorithm.
Mathematics 10 04750 g002
Figure 3. Example of encoding procedure.
Figure 3. Example of encoding procedure.
Mathematics 10 04750 g003
Figure 4. Example of weight mapping crossover (WMX).
Figure 4. Example of weight mapping crossover (WMX).
Mathematics 10 04750 g004
Figure 5. Example of swap mutation.
Figure 5. Example of swap mutation.
Mathematics 10 04750 g005
Figure 6. Locations of various businesses in Thailand’s central and northeastern regions: (a) Kanchanaburi province; (b) Khon Kaen province.
Figure 6. Locations of various businesses in Thailand’s central and northeastern regions: (a) Kanchanaburi province; (b) Khon Kaen province.
Mathematics 10 04750 g006
Figure 7. Scatter plot of the total profit obtained from Table 8.
Figure 7. Scatter plot of the total profit obtained from Table 8.
Mathematics 10 04750 g007
Table 1. The sightseeing time and time windows of tourist groups (hour).
Table 1. The sightseeing time and time windows of tourist groups (hour).
GroupAijkBiUi
j = 1j = 2j = 3j = 4
12.00(1)1.00(7)2.50(12)0.50(15)2.008.00
27.00(5)1.00(4)2.00(11)1.50(9)0.0014.00
31.50(10)1.00(14)1.50(3)7.00(6)0.0018.00
42.50(8)1.15(2)2.50(12)2.00(13)2.0012.00
51.50(7)1.00(4)2.50(8)1.50(9)1.508.00
61.50(2)6.50(6)2.50(8)1.50(10)0.0015.00
Table 2. Characteristics of the businesses in Khon Kaen.
Table 2. Characteristics of the businesses in Khon Kaen.
BusinessDetailQk (Person)
1Restaurant50
2Restaurant70
3Restaurant30
4Restaurant50
5Accommodation500
6Accommodation300
7Tourist attraction or POI200
8Tourist attraction or POI50
9Tourist attraction or POI100
10Tourist attraction or POI60
11Tourist attraction or POI80
12Tourist attraction or POI90
13Tourist attraction or POI60
14Tourist attraction or POI150
15Tourist attraction or POI120
Table 3. The revenue type of activity j for business k (THB/person).
Table 3. The revenue type of activity j for business k (THB/person).
RVjk123456789101112131415
1128.00600.00181.00200.001249.001482.0027.00497.0098.00327.00359.00409.00106.00201.0042.00
2117.00195.00425.00616.00913.002399.0055.00136.0053.00207.00138.00179.00455.00381.00376.00
3288.00381.00550.00698.003203.005455.00139.00253.00267.0098.00116.00482.0083.00441.00231.00
4389.00593.00144.00116.006284.00822.00329.0034.00167.00379.00300.00320.0023.00105.0038.00
Table 4. The total cost type of activity j for business k (THB/person).
Table 4. The total cost type of activity j for business k (THB/person).
CCjk123456789101112131415
176.80240.0063.3560.00374.70592.8010.80149.1034.30147.15179.50163.6047.70120.6023.10
270.2078.00148.75184.80273.90959.6022.0040.8018.5593.1569.0071.60204.75228.60206.80
3172.80152.40192.50209.40960.902182.0055.6075.9093.4544.1058.00192.8037.35264.60127.05
4233.40237.2050.4034.801885.20328.80131.6010.2058.45170.55150.00128.0010.3563.0020.90
Table 5. The profit each business obtained (THB).
Table 5. The profit each business obtained (THB).
Profit123456789101112131415
1960.001000.000.000.002988.000.0040.80174.000.001215.000.000.000.000.000.00
20.00475.000.001972.000.00 3192.00108.000.000.000.000.000.000.00972.000.00
30.000.003375.000.000.000.000.00374.400.000.00648.00510.000.000.000.00
40.000.000.000.000.003996.000.000.001161.60456.000.000.000.000.00414.00
SUM960.001475.003375.001972.002988.007188.00148.80548.401161.601671.00648.00510.000.00972.00414.00
TOTAL24,031.8
Table 6. The input data and parameter values used in all proposed methods.
Table 6. The input data and parameter values used in all proposed methods.
DetailParameterInterval Values
Input data:
Travel time between locations (hour)ThkU [0.05–8.00]
Sightseeing time (hour)AijkU [1.00–12.00]
Number of tourists in each group (person)DiU [2.00–50.00]
Servicing capacity (person)QkU [2.00–100.00]
Earliest arrival time of activity for any business (hour)EkijU [0.00–2.00]
Latest arrival time of activity for any business (hour)FkijU [0.5–72.00]
Starting time (hour)BiU [0.00–2.00]
Finishing time (hour)UiU [0.05–72.00]
Costs of business (THB/person):
 Restaurant;
 Accommodation;
 Attraction
CCjk
U [70–250];
U [200–3000];
U [10–300]
Revenue of business (THB/person):
 Restaurant;
 Accommodation;
 Attraction
RVjk
U [100–700];
U [300–7000];
U [20–500]
Parameter values used in all proposed methods:
Mutation ratepm0.2
Crossover ratepc0.8
Number of population sizeNP25
Termination condition:
 Small-size instances;
 Medium-size instances;
 Large-size instances
maxGen
300 iterations;
60 s;
120 s
Scaling factorF2
Crossover rateCR0.8
Table 7. Details of the problem instances.
Table 7. Details of the problem instances.
Instance No.Number of Locations
(K)
Number of Tourists
(i)
Number of Activities
(j)
Total Number of Tourists
(Person)
133360
2443153
3544125
4554135
5564197
6684169
7884231
88105277
99104136
109126395
1110123364
1210146368
1312153389
1415154372
1515203556
1620205635
1730305792
1830306803
1940406941
20504061085
Table 8. Computational results of the problem instances.
Table 8. Computational results of the problem instances.
Instance No.Lingo V.18Proposed MethodComp. Time (s)
Total Profit (THB)Comp. Time (s)Total Profit (THB)
CPDEHDEVNSGAHGAKVCPDEHDEVNSGAHGAKV
123,400.002.0023,400.0023,400.0023,400.0023,400.0023,400.0016.1322.0820.3223.4025.10
259,670.002.0036,780.0059,670.0059,670.0059,670.0059,670.0031.0822.2228.1024.2031.50
365,000.003.0055,400.0065,000.0065,000.0065,000.0065,000.0039.6542.6344.2043.5050.10
470,200.005.0054,890.0070,200.0070,200.0070,200.0070,200.0029.3134.9836.2037.1039.30
5102,440.008.0078,290.00102,440.00102,440.00102,440.00102,440.0054.9035.2344.2037.4049.20
697,880.00 BOF28,800.0057,030.0078,520.0081,640.0078,510.0081,710.0060.0060.0060.0060.0060.00
7120,120.00 BOF28,800.0093,290.00105,960.00109,120.00106,200.00115,180.0060.0060.0060.0060.0060.00
8130,050.00 BOF28,800.0065,030.0087,750.0097,500.0087,900.00111,530.0060.0060.0060.0060.0060.00
990,720.00 BOF28,800.0060,010.0063,440.0064,480.0063,440.0064,480.0060.0060.0060.0060.0060.00
10168,100.00 BOF28,800.0074,930.00108,600.00115,440.00114,600.00131,480.0060.0060.0060.0060.0060.00
11191,960.00 BOF28,800.0091,120.00140,010.00141,960.00139,120.00152,690.0060.0060.0060.0060.0060.00
12107,120.00 BOF28,800.0049,320.0069,640.0070,200.0070,550.0095,660.0060.0060.0060.0060.0060.00
13161,710.00 UB28,800.0083,380.00149,370.00150,540.00148,910.00155,490.00120.00120.00120.00120.00120.00
14193,440.00 UB28,800.0090,120.00141,440.00151,840.00144,310.00165,720.00120.00120.00120.00120.00120.00
15216,840.00 UB28,800.00113,900.00185,640.00197,730.00189,930.00200,100.00120.00120.00120.00120.00120.00
16302,750.00 UB28,800.00108,240.00190,450.00198,250.00197,120.00223,530.00120.00120.00120.00120.00120.00
17284,800.00 UB28,800.00122,370.00206,700.00231,400.00226,500.00242,370.00120.00120.00120.00120.00120.00
18226,340.00 UB28,800.0099,810.00153,660.00161,460.00156,000.00199,290.00120.00120.00120.00120.00120.00
19233,980.00 UB28,800.00132,080.00204,360.00189,540.00200,900.00221,880.00120.00120.00120.00120.00120.00
20366,300.00 UB28,800.00159,020.00247,260.00280,800.00256,600.00313,670.00120.00120.00120.00120.00120.00
Note: BOF = Best Objective Found; UB = Upper Bound.
Table 9. The p-value of the statistical test of differences of the solutions of the proposed method.
Table 9. The p-value of the statistical test of differences of the solutions of the proposed method.
DEHDEVNSGAHGAKV
CP0.0010.0010.0010.001
DE-0.0080.0350.001
HDEVNS--0.0110.001
GA---0.001
Table 10. The heuristic performance and relative improvement of all proposed methods obtained from Table 8.
Table 10. The heuristic performance and relative improvement of all proposed methods obtained from Table 8.
Instance No.HP (%)RI (%)RI (%)
The Proposed Methods Compared to the MILPThe Proposed Methods Compared to
the CP Algorithm
The Proposed Methods Compared to
the HGAKV Algorithm
CPDEHDEVNSGAHGAKVDEHDEVNSGAHGAKVDEHDEVNSGA
1100.00100.00100.00100.00100.000.000.000.000.000.000.000.00
261.64100.00100.00100.00100.0038.3638.3638.3638.360.000.000.00
385.23100.00100.00100.00100.0014.7714.7714.7714.770.000.000.00
478.19100.00100.00100.00100.0021.8121.8121.8121.810.000.000.00
576.43100.00100.00100.00100.0023.5723.5723.5723.570.000.000.00
658.2780.2283.4180.2183.4827.3730.1427.3630.203.900.093.92
777.6688.2190.8488.4195.8911.9614.5112.1619.018.005.267.80
850.0067.4774.9767.5985.7633.3033.3026.0241.6921.3212.5821.19
966.1569.9371.0869.9371.085.416.935.416.931.610.001.61
1044.5764.6068.6768.1778.2231.0035.0934.6243.0117.4012.2012.84
1147.4772.9473.9572.4779.5434.9235.8134.5040.328.307.038.89
1246.0465.0165.5365.8689.3029.1829.7430.0948.4427.2026.6226.25
1351.5692.3793.0992.0896.1544.1844.6144.0146.383.943.184.23
1446.5973.1278.4974.6085.6736.2840.6537.5545.6214.658.3812.92
1552.5385.6191.1987.5992.2838.6442.4040.0343.087.231.185.08
1635.7562.9165.4865.1173.8343.1745.4045.0951.5814.8011.3111.81
1742.9772.5881.2579.5385.1040.8047.1245.9749.5114.724.536.55
1844.1067.8971.3468.9288.0535.0438.1836.0249.9222.9018.9821.72
1956.4587.3481.0185.8694.8335.3730.3234.2640.477.9014.589.46
2043.4167.5076.6670.0585.6335.6943.3738.0349.3021.1710.4818.19
MAX56.16100.00100.00100.00100.0044.1847.1245.9751.5827.2026.6226.25
MIN55.8962.9165.4865.1171.080.000.000.000.000.000.000.00
AVG58.2580.8983.3581.8289.2429.0430.8029.4835.209.756.828.62
Table 11. The effects of changing the parameters of the case study.
Table 11. The effects of changing the parameters of the case study.
Scenario 1Scenario 2Scenario 3Scenario 4Scenario 5Scenario 6Scenario 7Scenario 8
Revenue + 5%;
Cost + 5%
Revenue + 5%;
Cost − 5%
Revenue − 5%;
Cost + 5%
Revenue − 5%;
Cost − 5%
Revenue + 10%;
Cost + 10%
Revenue + 10%;
Cost − 10%
Revenue − 10%;
Cost +10%
Revenue − 10%;
Cost − 10%
Business no. 3: Restaurant
Revenue157.50157.50142.5142.5165.00165.00135.00135.00
Cost78.7571.2578.7571.2582.5067.5082.5067.5
Profit3543.753881.252868.753206.253712.504387.502362.503037.50
% Change in profit5.00%15.00%(15.00%)(5.00%)10.00%30.00%(30.00%)(10.00%)
Business no. 6: Hotel
Revenue652.05652.05589.95589.95683.10683.10558.90558.90
Cost391.23353.97391.23353.97409.86335.34409.86335.34
Profit7547.408625.605750.406828.607906.8010,063.204312.806469.20
% Change in profit5.00%20.00%(20.00%)(5.00%)10.00%40.00%(40.00%)(10.00%)
Business no. 10: Attraction
Revenue87.1587.1578.8578.8591.3091.3074.7074.70
Cost34.8631.5434.8631.5436.5229.8836.5229.88
Profit1754.551865.951476.051587.451838.102060.901281.101503.90
% Change in profit5.00%11.67%(11.67%)(5.00%)10.00%23.33%(23.33%)(10.00%)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Thumrongvut, P.; Sethanan, K.; Jamrus, T.; Wongloucha, C.; Pitakaso, R.; Golinska-Dawson, P. Metaheuristics in Business Model Development for Local Tourism Sustainability Enhancement. Mathematics 2022, 10, 4750. https://doi.org/10.3390/math10244750

AMA Style

Thumrongvut P, Sethanan K, Jamrus T, Wongloucha C, Pitakaso R, Golinska-Dawson P. Metaheuristics in Business Model Development for Local Tourism Sustainability Enhancement. Mathematics. 2022; 10(24):4750. https://doi.org/10.3390/math10244750

Chicago/Turabian Style

Thumrongvut, Pawnrat, Kanchana Sethanan, Thitipong Jamrus, Chuleeporn Wongloucha, Rapeepan Pitakaso, and Paulina Golinska-Dawson. 2022. "Metaheuristics in Business Model Development for Local Tourism Sustainability Enhancement" Mathematics 10, no. 24: 4750. https://doi.org/10.3390/math10244750

APA Style

Thumrongvut, P., Sethanan, K., Jamrus, T., Wongloucha, C., Pitakaso, R., & Golinska-Dawson, P. (2022). Metaheuristics in Business Model Development for Local Tourism Sustainability Enhancement. Mathematics, 10(24), 4750. https://doi.org/10.3390/math10244750

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

Article Metrics

Back to TopTop