Next Article in Journal
On the Constrained Solution of RBF Surface Approximation
Next Article in Special Issue
Optimal Emergency Evacuation Route Planning Model Based on Fire Prediction Data
Previous Article in Journal
On the Double ARA-Sumudu Transform and Its Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Artificial Bee Colony Algorithm with Nelder–Mead Method to Solve Nurse Scheduling Problem

1
Department of Information Technology, Sri Manakula Vinayagar Engineering College, Pondicherry 605107, India
2
Department of Computer Science and Technology, Madanapalle Institute of Technology & Science, Madanapalle 517325, India
3
Department of Information Technology, College of Computer and Information Technology, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
4
Department of Computer Science & Engineering, Graphic Era Deemed to Be University, Dehradun 248002, India
5
Department of Computer Science and Engineering, Women Institute of Technology, Dehradun 248007, India
6
Uttaranchal Institute of Technology, Uttaranchal University, Dehradun 248007, India
7
Department of Project Management, Universidad Internacional Iberoamericana, Campeche 24560, Mexico
8
Department of Computer Engineering, Faculty of Science and Technology, Vishwakarma University, Pune 411048, India
*
Author to whom correspondence should be addressed.
Submission received: 17 June 2022 / Revised: 19 July 2022 / Accepted: 19 July 2022 / Published: 25 July 2022
(This article belongs to the Special Issue Combinatorial Optimization Problems in Planning and Decision Making)

Abstract

:
The nurse scheduling problem (NSP) is an NP-Hard combinatorial optimization scheduling problem that allocates a set of shifts to the group of nurses concerning the schedule period subject to the constraints. The objective of the NSP is to create a schedule that satisfies both hard and soft constraints suggested by the healthcare management. This work explores the meta-heuristic approach to an artificial bee colony algorithm with the Nelder–Mead method (NM-ABC) to perform efficient nurse scheduling. Nelder–Mead (NM) method is used as a local search in the onlooker bee phase of ABC to enhance the intensification process of ABC. Thus, the author proposed an improvised solution strategy at the onlooker bee phase with the benefits of the NM method. The proposed algorithm NM-ABC is evaluated using the standard dataset NSPLib, and the experiments are performed on various-sized NSP instances. The performance of the NM-ABC is measured using eight performance metrics: best time, standard deviation, least error rate, success percentage, cost reduction, gap, and feasibility analysis. The results of our experiment reveal that the proposed NM-ABC algorithm attains highly significant achievements compared to other existing algorithms. The cost of our algorithm is reduced by 0.66%, and the gap percentage to move towards the optimum value is 94.30%. Instances have been successfully solved to obtain the best deal with the known optimal value recorded in NSPLib.

1. Introduction

In a hospital, various operations are performed; nurse rostering is a resource allocation problem. Every day, work is divided into three periods: day shift, noon shift, and night shift. The process consists of allocating workload to the nurses periodically by considering hospital terms, namely constraints and requirements, for a scheduling period of one month. Constraints are classified as hard and soft constraints; hard constraints are needed to be satisfied when allocating the roster. The soft constraints are considered when the prerequisite is desirable but not obligatory. The soft constraint will estimate the quality of the roster; these constraints are considered as much as possible, but violating soft constraints leads to a penalty, while the roster is still considered feasible. Each nurse is allocated to a specific number of shifts with various constraints. Each shift is accomplished by the group of nurses based on their preferences and skills required for the particular operation. The NSP is a multi-objective NP-hard combinatorial problem [1,2,3,4].
Nurse scheduling is allocating nurses with different services and bonds in a subject satisfying the constraints to reduce the overall cost of the scheduling problem. The NSP is the placement of nurses with different skill patterns into shifts pattern in a subject satisfying the constraints. This problem focuses on generating a roster and representing slot allotment for the nurses to shift. The obtained roster solution should persuade hospital regulations, nurse preferences, and requests listed in terms of soft constraints. Two methods can enable the scheduling of nurses: cyclic rostering and non-cyclic rostering. In cyclic rostering, each nurse is assigned to a shift for the planning epoch. The same cyclic method is repeated without modifications in the shift pattern for consecutive days. In this approach, the workload among the nurses can easily be distributed among them within a stipulated time [5,6]. In non-cyclic rostering, the shift pattern allotted to the nurse keeps changing for various periods. The new shift pattern is scheduled based on the nurse’s preferences, requests, and constraints. The non-cyclic approach is more advantageous than cyclic rostering because updating the roster will cause the entire schedule to modify. In contrast, in the non-cyclic process, only the pattern part needs to be changed. The non-cyclic approach of rostering will effectively inscribe all constraints, preferences, and requirements required [7].
The objective is to increase the nurses’ priority and diminish the total penalty cost from desecrations of the soft constraints. Most researchers have considered only a single goal, and only a few seek multi-objectives to solve NSP [8,9,10]. NSP is the assignment of nurses to the shifts, and the main objective is to assign nurses by satisfying the hospital regulations. The NSP is solved using heuristic, meta-heuristic, and mathematical programming models [11,12,13].
The automatic planning of NSP [14] is to generate a roster for working nurses and rest days for a particular period concerning a set of constraints. The head nurse manually commences this roster; scheduling nurses for three shifts for one month hardly takes six to eight hours. This roster may be preliminary since the changes occur based on the nurse’s opinion. Some researchers developed approaches to solving NSP. Aickelin and Dowsland [15] used a genetic algorithm framework with Tabu search as a local search to exploit the solution. Burke et al. used simulated annealing multi-objective approach to solve NSP and generate a key in the Pareto front [16].
The mathematical programming model is the first approach for solving the NSP; although mathematical programming methods ensure that they will provide the optimum solution, they fail to do so within a rational amount of time when the solution search space is vast. This requires a constraint-based programming method to solve complicated and soft constraints [17]. A heuristic approach is ideal for solving the problem with benign conditions and functionality. The heuristic-based approaches, such as the local search approaches of simulated annealing [18], iterated local search [19], and Tabu search [8], are combined with other methods to solve the NSP effectively. The meta-heuristic approach is implemented over low-level heuristics to improve their performance. The hybrid system combines the heuristic method with any deterministic process, such as a mathematical programming model [20]. The population-based way is a type of heuristic such as a genetic algorithm [21], scatter search [22], particle search [23], and ant colony optimization [24]. These approaches generate a population of solutions from which a better solution can be achieved.
Berrada et al. proposed a lexico-dominance technique to solve NSP with several soft constraints. The soft constraints are prefixed with priority values and ordered to obtain the best solution [25]. Burke et al. applied weight values to the soft constraints, which are determined by the hospital regulations [26]. The conditions are in order based on the weight values, and the best solution is chosen. Nikola and Petrovic introduced a bee colony optimization algorithm to solve NSP, and unscheduled shifts are assigned to the nurse constructively, leading to higher penalties [27,28]. Meng Yang et al. introduced a multi-objective artificial bee colony algorithm to solve the multi-objective home healthcare routing and scheduling problem. Based on the enormous neighborhood search, a heuristic solution update strategy has been proposed to trade off the search balance and achieve the workload balance in the scheduling [29].
In this research work, the authors introduced the artificial bee colony algorithm with the Nelder–Mead method (NM-ABC) to solve nurse scheduling problems. The local search mechanism is incorporated in the onlooker bee phase of ABC to obtain an optimal solution for NSP. ABC algorithm produced good results in solving optimization problems and improved the efficiency of the global search, memory management, and search improvement process [30,31,32,33,34,35,36,37]. However, ABC follows the same neighborhood search strategy in the onlooker and employee bee phase. Improving neighborhood search on quality individuals enhances the probability of convergence towards global optimum compared to a random process. Thus, this paper introduces the Nelder–Mead Local search in the onlooker bee phase to improve the neighborhood search.
The main contributions of this work are illustrated as follows:
  • A hybrid meta-heuristic algorithm, namely artificial bee colony optimization with Nelder–Mead Method, is proposed.
  • The search capability of ABC is enriched with the aid of the Nelder–Mead method, which consists of search strategies such as midpoint, reflection, expansion, contraction, and shrinkage processes. These search strategies enhance the balance between exploration and exploitation.
  • NM-ABC is implemented and tested on the nurse scheduling problem (NSPLib).
  • The performance of NM-ABC is compared with that of some classical optimization algorithms.
The rest of the paper is structured as follows: Section 2 illustrates the scientific formulation of NSP. In Section 3, the proposed artificial bee colony algorithm with NM is discussed. Section 4 shows the applicability of NM-ABC to solve NSP with experimental analysis. Section 5 deliberates a detailed analysis of the performance outcome of NM-ABC to solve NSP. Finally, Section 6 concludes the research work.

2. Nurse Scheduling Problem

The NSP can be described as allocating a set of nurses to a group of shifts for a given time. The constraints are organized by hospital regulations, nurse preferences, nurse requests, and working practices. Generally, NSP is divided into two different types of constraints: they are hard and soft constraints. Hard constraints are the regulations that must be gratified to achieve a feasible solution. The roster pattern generated should satisfy all hard constraints; generally, the general hard constraint is allocating shifts to the restricted number of nurses. Soft constraints are desirable and not obligatory but must be satisfied as much as possible. The soft constraints will determine the quality of the roster formed. The violation of any soft constraints leads to the penalty of the solution. The general soft constraint in NSP is to balance the workload among all nurses and usage of nurse resources efficiently [38,39]. The visual representation of NSP is illustrated in Figure 1.
The problem of NSP consists of a set of nurses N = { 1 , 2 , , i } that is allocated a set of shifts S = { 1 , 2 , , j } for the scheduled period D = { 1 , 2 , , k } in days. The shift pattern r consists of the allocation of shifts for the particular nurse over the expected period for the 0/1 matrix. The main objective of the NSP is to identify a possible solution such that the total cost is minimized. The solution representation of the NSP is given as
X j , d , r = { 1 ,   p a t t e r n   c o v e r s   s h i f t   j   o n   d a y   d 0 ,                                                                             o t h e r w i s e
The nurse preference is the wish expressed by a particular nurse to work on a distinct shift for a specific day. The penalty cost is added to the solution if the nurse’s preference is not achieved. The preference cost C i , j , d is the wish of the particular nurse i to work on shift j on a particular day d . The total preference cost of the shift pattern r for the scheduled period can be calculated as
P i , r = j = 1 S d = 1 D X j , d , r C i , j , d  
The objective function of the NSP can be formulated as
M i n i m i z e i = 1 N r = 1 δ i X j , r P i , r  
where δ i is the set of a possible shift pattern for the nurse i to work on the scheduled period.
Subject to
r = 1 δ i X j , r = 1  
The feasible solution is for all nurses available in the schedule, and this constraint specifies that precisely one shift pattern is allocated to every nurse in the hospital. The smallest number M j , d of nurses essential for the shift j on the scheduled day d   can be restricted by using Equation (5).
i = 1 N r = 1 δ i X j , d , r M j , d  
X i , r { 0 , 1 }
The Equations (4)–(6) are the coverage constraints considered in solving the NSP.
All the hard constraints in the NSP are mandatory and should be considered when designing the schedule pattern. The soft constraints are included in the objective function of increasing the quality of the solution. Soft restrictions are not mandatory and enhance the quality of the schedule. Violation of soft constraints leads to penalty cost, and the hospital management suggests the penalty cost. Some of the soft constraints are listed in this work, and the soft constraints (SC) are as follows
SC1:
In this constraint, the nurse has restricted work on a specific day. The constraint can be evaluated as
i = 1 N X i , d 0
SC2:
The nurse can request to work on a specific day. The constraint can be calculated as
i = 1 N X i , d 1
SC3:
The nurse can request not to work on a specific shift. The constraint can be
i = 1 N X i , j 0
SC4:
The nurse can request to work on a specific shift. The constraint can be determined as
i = 1 N X i , j 1
SC5:
No single work shift between two days off. The constraint can be calculated as
X i , ( j 1 ) X i , j + X i , ( j + 1 ) 0
SC6:
The nurses are not allowed to work more than three consecutive days. The constraint can be
X i , 2 j 3 0

3. Proposed Algorithm

3.1. Artificial Bee Colony Algorithm

Karaboga and Bahriye developed an artificial bee colony algorithm (ABC), inspired by honey bees’ natural behavior. The intelligent behavior of honey bees helps to find the near-optimal solution for the optimization problem [32]. ABC is the population-based algorithm and consists of three groups of honey bees: employed bee, onlooker bee, and scout bee. The colony consists of an equal number of employed bees and onlooker bees. Each solution in the population is held by one employed bee. The employed search for the food source and share the direction of the food source with onlooker bees through waggle dance. Based on the probability calculation, the higher-quality food sources are selected by the onlooker bee phase, and the bees continue further searches. The low-quality food sources are rejected, and employed bees are converted to scout bees. The scout bee will search for a new food source or food position.

3.1.1. Initialization

The initialization of the population is created with food source FS for n-dimensional vectors. The population of the solution is represented as X i = { x i , 1 , x i , 2 , , x i , n } . The food source in the population is generated using Equation (13).
x i , j = x m i n , j + r a n d   ( 0 , 1 ) ( x m a x , j x m i n , j )
The food sources are randomly distributed to employed bees. The objective calculation for the solution is evaluated accordingly.

3.1.2. Employed Bee Phase

In the employed bee phase, the candidate solution is generated, and the position of the food sources is monitored. The candidate solution can be developed using
v i , j = x i , j + Ø i , j ( x i , j x k , j )
where j = 1, 2, …, S, and we choose the value of k as different from I, k = 1, 2, …, FS, and the value of Ø ranges (−1, 1). The greedy selection is made between v i and x i based on the fitness calculation. If the fitness value of v i is greater than x i , the solution x i is replaced by v i .

3.1.3. Probability Calculation

After fitness calculation, the employed bees share the direction of the food source with onlooker bees. The onlooker bees evaluate the fitness of the employed bee’s solution using the probability value p i . The solution’s probability and fitness calculation are shown in detail in Algorithm 1.
Algorithm 1: Probability calculation
1:Fori = 1, 2, …, FS, do
2: Calculate probability values P i j for the solution v i , j
3:   P i = f i t i j = 1 F P f i t j
4:   f i t i = { 1 1 + f i ,                               f i 0 1 + a b s ( f i ) ,                   f i < 0
5:End For

3.1.4. Onlooker Bee Phase

Based on the probability value p i , Each onlooker bee randomly chooses the food source x i with probability p i . Onlooker bee performs modification on x i using Equation (14). To select the best solution among x i and v i , the greedy method is used, which is similar to the employed bee phase.

3.1.5. Scout Bee Phase

After employed and onlooker bees search, the food source is abandoned if the solution cannot be improved and exhausted for the predefined number of iterations. The corresponding employed bee develops a scout bee and explores new food sources using Equation (13).

3.2. Nelder–Mead Method

The NM method is used to find the local minimum for the given function and is represented as a simple triangle with three vertices. The worst vertex is found among the triangle, rejected for the next iteration, and replaced with a new vertex. The search continues toward the best solution in the triangle sequence, reducing the triangle size. At last, the vertex with a minimum point is chosen as the best solution.
Let f ( x , y ) be the function to be minimized; start with three vertices of a triangle. The fitness function is evaluated for all three points of the triangle. The vertices are ordered based on the fitness value: best (I), good (J), and worst (K) vertices are ordered. NM method works on four operations: reflection, expansion, contraction, and shrinkage.

3.2.1. Midpoint (M)

The midpoint between the line joining vertices is calculated for the first two best and good vertices using
M = I + J 2

3.2.2. Reflection (R)

The function increases towards the side of the triangle when we move from the worst vertex to the best vertex and decreases when moving from the worst to the excellent vertex. The test point reflection point R is chosen along the side of I J ¯ . To find R, calculate the midpoint of the best and good vertex because the best solution is away from the worst vertex. Draw a line joining K and M of length d . Place R extending at a distance d . The formula to calculate R is
R = 2 M K

3.2.3. Expansion (E)

The fitness calculation at vertex R is calculated, and if it is less than K , then the search is moved towards a minimum value. Now, extend the line segment through the vertex R to E by the distance d . If the fitness value at vertex E is better than R , then it is towards the minimum value. The formula to calculate E is given by
E = 2 R M

3.2.4. Contraction (C)

When the fitness value of R and K are small, another vertex in the triangle is needed to continue the process. Then, contraction towards midpoint without replacement is performed. The contraction points C 1 and C 2 are drawn along the line joining K M ¯ and M R ¯ for the length of d 2 . The formula to calculate C is
C = R + M

3.2.5. Shrinkage (S)

The fitness value at C is calculated, and if it is not less than K , then the vertices J and K should be shrunk towards the best vertex I . The vertex J is replaced with M , and K is replaced with S . The point S is the midpoint of I   and K . The process is continued until the minimum value is found. The detailed description of the flow of the NM method is shown in Algorithm 2.
Algorithm 2: Nelder–Mead Method
1:Produce new food source v i using modified Nelder–Mead Method
2:Let v i denote list of vertices
3: ɽ , μ, λ, and ζ are the constants of likeness, extension, shrinkage, and contraction
4:f is the fitness method to be reduced
5:For i = 1, 2, …, n + 1 vertices, do
6: Order the vertices from deepest fitness function f(v_1) to maximum fitness function f(〖v〗_(n + 1))
7: f ( v 1 ) f ( v 2 ) f ( v n + 1 )
8: Calculate midpoint for best two vertices
9: v m = v i n , where i = 1, 2, …, n
10: Calculate reflection point v_r
11: v r = v m + ɽ ( v m v n + 1 )
12:if f ( v 1 ) f ( v r ) f ( v n ) then
13:   v n v r and go to stopping criteria
14:End if
15: Calculate expansion point v e
16:if f ( v r ) f ( v 1 ) then
17:   v n v e and go to stopping criteria
18:End if
19:if f ( v e ) < f ( v r ) then
20:   v n v e and go to stopping criteria
21:else
22:   v n v r and go to stopping criteria
23:End if
24: Calculate contraction point v c
25:if f ( v n ) f ( v r ) f ( v n + 1 ) then
26:  Compute outside contraction
27:   v c = λ   v n + 1 + ( 1 λ )   v m .
28:End if
29:if f ( v r ) f ( v n + 1 ) then
30:  Compute inside contraction
31:   v c = λ   v n + 1 + ( 1 λ )   v m .
32:End if
33:if f ( v r ) f ( v n ) then
34:  Shrinkage is done between v m and the best vertex among v r and v n + 1 .
35:End if
36:if f ( v c ) < f ( v r ) then
37:   v n v c and go to Stopping criteria
38:else go to Shrinkage phase
39:End if
40:if f ( v c ) f ( v n + 1 ) then
41:   v n + 1 v c and go to Stopping criteria
42:else go to the Shrinkage phase
43:End if
44: Calculate Shrinkage
45: Shrink close the best individual with new apices
46: v i = ζ v i + v 1 ( 1 ζ ) , where i = 2, …, n + 1
47:End for
48:Determine the new vertices of the simplex thus formed based on their fitness and continue with the process of the reflection phase

3.3. Nelder–Mead Method-Based ABC (NM-ABC)

ABC algorithm is lithe to improve and progress; the intricacy of the algorithm is reduced since it uses fewer parameters. The improved search ability aids in attaining optimal solutions with less computational time and increased convergence speed. ABC is good at examination but fails in exploitation [40]. An improved local search algorithm is incorporated in ABC to tradeoff the search [41].
The Nelder–Mead method is the famous local search algorithm, and it is simple and efficient. It is also proficient at embedding in other global search algorithms. However, it is entrapped in a local optimum solution. Thus, it is poor at exploration and has low convergence towards the initial position. Nelder–Mead (NM) method is good at the exploitation process but poor in the exploration process. This paper uses the NM method in ABC to improve the exploitation. The detailed description of the proposed algorithm’s pseudo code is shown in Table 1, Table 2 and Table 3. The workflow of the proposed NM-ABC is shown in Figure 2. The proposed NM-AMC algorithm is presented in Algorithm 3.
NM method is used in the onlooker bee phase of the ABC algorithm. The intention is to use the NM method in the onlooker bee phase instead of the employee bee phase since the individuals participating in the onlooker bee phase are selected based on probability. When the individual is chosen with probability, it is considered a quality individual in terms of fitness. Intensive search on quality individuals results in global optimum rather than searching on random individuals.
Algorithm 3: NM-ABC
1:Initialize the population
2:For i = 1, 2, …, FS, do
3:For j = 1, 2, …, S, do
4:  Generate x i , j solution
5:   x i , j = x m i n , j ± r a n d   ( 0 , 1 ) ( x m a x , j x m i n , j )
6:  Where x m i n , j and x m a x , j are the min and max limit of the dimension j
7:End for
8:  Compute the objective of the population
9:   i t e r = 1
10:  Repeat
11:  {
12:   Employed Bee Phase
13:   For each food source i do
14:     Generate candidate solution v i using Equation (14)
15:    Select between v i and x i
16:   End For
17:   Onlooker Bee Phase
18:   Set r = 0
19:   While (r <= FS)
20:    If rand(0,1) < P i using Algorithm 3
21:     Generate candidate solution v i by Algorithm 2
22:     Select between v i and x i
23:     r = r + 1
24:    End if
25:   End while
26:   Scout Bee Phase
27:   Abandon the food source x i , which cannot improve further using Equation (13)
28:   Remember the best individual obtained so far
29:   iter = iter + 1
30:  }
31:End for
32:Until iter = max FEs

4. Experimental Results

4.1. Experimental Setup

The NSP datasets are taken from the library NSP lib and consist of 25 to 100 instances; each contains “N” nurses, “D” days, and “S” shift patterns. The nurses’ N varies from 25 to 100 nurses, and the schedule is to be made for given D days with S shift patterns. NSP lib contains cases to describe the maximum and minimum utilization of resources for the health care unit. It consists of “N” nurses, “D” days, and “S” shift patterns and provides the coverage matrix for days with shift patterns. The nurses’ workload for the shift concerning a particular day relates to N and D for handling 25 to 100 instances.
The performance of NM-ABC for NSP is evaluated using the NSPLib dataset [42]. The dataset is accessed on 10 October 2021 from the specified link (https://www.projectmanagement.ugent.be/research/personnel_scheduling/nsp). The author summarizes the characteristic of the dataset in Table 1. The proposed technique to solve NSP is implemented with the help of MATLAB 2016a under Windows on an Intel i5 processor with 8 GB of RAM and 1 TB storage. Table 1 designates the cases utilized by NM-ABC to solve NSP.
The parameter settings of NM-ABC and compared algorithms are presented in Table 2. In this experimentation, algorithms M1 to M5 and NM-ABC consist of a population size of 100, and the maximum number of iterations is 1000. These algorithms will stop its execution when the maximum iterations or optimal solution are reached. The association of the results obtained by NM-ABC clearly specifies that the proposed technique is relatively better than the existing methods.

4.2. Performance Metrics

The performance of the proposed technique is evaluated by relating it to five different opponent methods. Here, eight performance metrics are utilized to evaluate the experimental results listed in this section.

4.2.1. Average Best Time (ABT)

The best time is to achieve the best value for a particular instance. The average best time is the type of the best times of all test cases taken from the dataset. In experimental analysis, the average best time (ABT) is computed using
ABT = i = 1 n Time   is   taken   to   achieve   the   best   value Total   number   of   instances  
where n is the number of test cases in the given data set.

4.2.2. Standard Deviation

Standard deviation (SD) is the portion of distribution among set values from its mean value. Average standard deviation determines the type of the standard deviation of all test cases taken from the dataset. The average standard deviation (ASD) can be measured using
ASD = i = 1 n ( value   obtained   in   each   instance i Mean   value   of   the   instance ) 2
where n is the number of instances in the given data set.

4.2.3. Least Error Rate

The least error rate (LER) is the variance between the actual optimal value and the obtained best value. The LER can be calculated using
LER = V a l u e N S P L i b B e s t   V a l u e i

4.2.4. Success Percentage

Success percentage is the number of instances that attain optimal value for the given instance. Average success percentage (ASP) is the average number of models that achieve optimal value to the total number of test cases taken from the dataset. The value of ASP can be computed using
ASP = i = 1 n Number   of   instances   succeed   to   attain   optimal   value   Total   number   of   instances   100
where n is the number of instances in the given data set.

4.2.5. Cost Reduction

Cost reduction is the variation between actual cost in NSPLib and the cost obtained from our approach. Average cost reduction (ACR) is the average cost reduction to the total number of test cases taken from the dataset. The value of ACR can be computed from
ACR = i = 1 n C o s t i C o s t N S P L i b Total   number   of   instances  
where n is the number of test cases in the given data set.

4.2.6. Gap

The value of gap is the distance to attain the best deal. The average gap (AGap) is the average distance to achieve the best value from all instances of the total number of cases taken from the dataset. The value of AGap is calculated using
ACR = i = 1 n ( C o s t i C o s t N S P L i b )   /   C o s t N S P L i b Total   number   of   instances   100
where n is the number of instances in the given data set.

4.2.7. #Both Feasible Solution

#Both feasible solution (#BFS) is the number of feasible solutions found to obtain optimal value concerning both NSPLib and our approach’s best value.

4.2.8. #Feasible Solution

#Feasible solution (#FS) is the number of feasible solutions found to obtain optimal value concerning known optimal values of NSPLib.

4.3. Experimental Result Analysis

The results obtained by NM-ABC with competitive methods are shown in Table 3. The performance is associated with exiting techniques; the value in the table defines the attained best value by proposed and other competitor’s algorithms. The objective of NSP is to reduce cost; the lowermost principles are the best solution obtained. In the evaluation of the performance of the algorithm, the authors have considered 15 different cases of diverse sizes with five other instances. It is proven proposed NM-ABC accomplished 44 best results out of 75 instances.
The experimental results of 15 cases are summarized, and the best results are achieved within the time. The best solutions for each case in all methods are highlighted in bold font. The best values are obtained by using our proposed Algorithm 1. It is noted from Table 4 that, experimentally, NM-ABC obtained 12 best results out of 25 instances in smaller-sized datasets from case1 to case5, 14 best results in medium-sized datasets from case6 to case10, and 17 best results out of 25 instances in larger-sized datasets case12 to case15. The experimental results of the proposed algorithm have a high potential in exploiting search space solutions towards better results in various ways.
Table 4 shows the best time, standard deviation, and the least error rate for each case recorded for ten runs. The mean value of the proposed algorithm is 1.75% reduced compared to that of other competitive methods, showing that our proposed algorithm attained a lesser worst value in addition to the best solution. The least error rate for the proposed algorithm can be calculated using Equation (21) with the known optimal value recorded in NSPLib. The standard deviation is increased by 10% compared to other competitive methods. The computational time to attain each best value is shown, and the time taken to reach the best solution for each case is tabulated under the best time in a table. Our proposed method yields 39.32% less computational time to attain the best results than other competitor methods.
Table 5 describes the number of feasible solutions obtained by NM-ABC and other methods, and the table shows our proposed algorithm produced a more feasible solution than the solutions recorded in NSPLib. The solution is viable if the hard and soft constraints listed in Table 4 and Table 5 are satisfied. NM-ABC satisfied several feasible solutions. #Both feasible is the number of possible solutions recorded in NSPLib or NM-ABC algorithm, and #feasible is the number of viable solutions obtained in NSPLib. Table 5 shows NM-ABC has attained 90% of a better value than known values reported in NSPLib. Compared with other methods, our algorithm outperforms with 87% of the best value, which is higher than other methods listed in Table 5. Thus, our algorithm contributes a new deterministic search and practical heuristic approach to solve NSP using NSPLib dataset.
Table 6 summarizes the comparison and assessment of average best time, average standard deviation, and average success percentage obtained by our proposed algorithm NM-ABC with another competitor method. In Table 4, the columns 5 and 6 describe the average best time, average standard deviation, and average success percentage of the proposed algorithm. Columns 7 to 21 describe the performance metrics of another competitor’s method.
Figure 3 compares the average best time of NM-ABC with other methods. The best time is the time taken to attain the best value, for instance, using Algorithm 2. The time is taken to allocate and schedule nurses for a particular time without overruling hard constraints and reducing violations of soft constraints. For smaller datasets, the computational time taken by NM-ABC is reduced to 56.72%. For medium datasets, the time taken is reduced by 36.40%. For larger datasets, the time taken is reduced by 34.31% compared to other competitor methods. The ABT is calculated using Equation (19) and shows the reduced computational time to schedule nurses for a particular shift on days for the scheduled period. The HSHH algorithm consumes more computational time to solve NSP, while MAPA and BCO achieve 50% of our proposed approach.
Figure 4 portrays the comparison of the average standard deviation of NM-ABC with another competitor method. Figure 4 shows that an increase in the standard deviation will increase the search space to obtain the best value.
Figure 5 compares our proposed algorithm NM-ABC algorithm’s average success percentage with other methods. Success percentage is the number of instances in attaining optimal value for the given instance. Average success percentage (ASP) is the average number of instances that obtains optimal value for the total number of cases taken from the dataset. ASP is calculated using Equation (22) and proves NM-ABC has increased the success percentage over other competing methods. In Table 7, it is observed that NM-ABC achieved a 100% success percentage except for case 3 and case 14. Our algorithm NM-ABC shows an overall 97.33% of success percentage. For smaller datasets, it is shown that the success percentage is 96% and 26% more than that of all other competitor methods. For the medium dataset, our algorithm achieved a 100% success percentage of 56.25% more than other competing methods. The larger dataset NM-ABC attained a 96% success percentage, which is 48.18% more than other methods. In the competitor method, the multi-objective ant colony optimization algorithm (M4) achieved the second-best success percentage with an overall 72%.
Table 7 summarizes the comparison and assessment of average cost reduction and average gap percentage obtained by NM-ABC with another competitor method. Table 7, the 4th and 5th columns describe the average cost reduction and average gap percentage of the proposed algorithm. Columns 6 to 15 illustrate the performance metrics of another competitor’s method.
Figure 6 portrays the analysis and comparison of the average cost reduction of our proposed NM-ABC with another competitor method. ACR is the difference between the best-known value observed in NSPLib and the cost obtained from our algorithm. Average cost reduction (ACR) is the average cost reduction from the dataset to the total number of instances and calculated using Equation (23). In Table 7, it is shown that NM-ABC minimized the cost of the NSP. The main objective of NSP is to reduce resource utilization, which is reflected by the cost reduction, as shown in Figure 6. Using our proposed algorithm, NM-ABC, the cost of NSP is reduced by 0.66%. For more minor instances, NM-ABC is reduced to 1.12% compared to other competing methods and is reduced to 0.11%. For medium instances, our proposed algorithm is reduced to 0.62% more than the original cost value recorded in NSPLib, and other methods reduced it to 0.70%. Compared to more significant instances, the proposed algorithm is reduced to 0.63%; compared with other competitor methods, it is reduced to 0.68%. In the proposed NM-ABC algorithm, the cost of resource utilization decreases with an increase in the dataset size.
Figure 7 compares the average gap percentage of NM-ABC with another competitor method. The gap percentage is the distance between the attained best value and the known optimal value recorded in NSPLib. The average gap (AGap) is the average distance to obtain the best and known value from all instances to the total number of instances. The value of AGap is calculated using Equation (23), and from Table 7, it is proven that NM-ABC attained a positive value, which shows the algorithm moved towards the best optimum value. Our algorithm achieved 94.30% of successfully solved instances to reach the best value concerning the known value recorded in NSPLib.

5. Discussions

The experiments to solve NP-hard combinatorial NSP are conducted by our proposed algorithm NM-ABC. Problem-based existing algorithms are chosen and compared with the proposed NM-ABC algorithm. The result of our proposed algorithm is compared with other competitor methods, and the best values are compared in Table 7. Various performance metrics are considered to evaluate the proposed algorithm’s performance. Table 3, Table 4, Table 5, Table 6 and Table 7 show the outcome of our proposed algorithm and other existing method’s performance. From the table and figure, it is evident that NM-ABC has more ability to attain the best value with less computation time compared to the known optimal value listed in NSPLib. The average number of function evaluations (NFEs) for proposed NM-ABC is observed concerning to number solutions updated using reflection, contraction, expansion, and shrinkage phase. We noticed that NFEs of proposed algorithm is 106 for all the test cases.
Compared with other existing methods, the mean value of NM-ABC is reduced by 1.75% compared to that of other competitive methods and attained a lesser worst value in addition to the best solution. The proposed method yields 39.32% less computational time to obtain best results compared to other competitor methods. The datasets are divided based on their size as smaller, medium, and large datasets; the computational time taken by NM-ABC is reduced to 56.72%, 36.40%, and 34.31%, respectively. The success percentage to attain the best value of our proposed approach is 97.33%. Compared with other methods with various-sized datasets, our algorithm achieves 26% for the smaller dataset, 56.25% for medium datasets, and 48.18% for larger datasets. The cost of our algorithm is reduced by 0.66%, and the gap percentage to move towards the optimum value is that 94.30% instances were successfully solved to obtain the best deal with the known optimal value recorded in NSPLib.
Our algorithm has proven significant performance in attaining the best solution with optimized resource utilization and nurse preferences by satisfying both hard and soft constraints. It is also shown that the existing approach solves NSP with higher utilization of resources and violation of soft constraints that lead to increased cost. The ability to distribute the workload among nurses with nurse performance and satisfaction are achieved in our algorithm. The proposed system is tested on larger datasets and works astoundingly well than the other techniques.

6. Conclusions

This paper solves NSP using a hybrid artificial bee colony algorithm with Nelder–Mead (NM-ABC) method. The proposed algorithm is evaluated using the NSPLib dataset, and the performance of the proposed algorithm is compared with the other five existing methods and evaluated in the NSPLib dataset. To assess the implementation of our proposed algorithm, 75 different cases of various-sized datasets were chosen for evaluation, and in that, 44 out of 75 instances achieved the best result. The evaluation of the proposed algorithm is compared with other existing techniques in terms of the best time, standard deviation, least error rate, success percentage, cost reduction, average gap, #both feasible solutions, and the number of #feasible solutions. When comparing the results of existing algorithms in metrics listed, the proposed NM-ABC outperforms in most instances of NSP. Future work of this research can be extended with more objectives in NSP for optimization.

Author Contributions

Conceptualization, R.M.; methodology, R.M. and R.R.; validation, S.S.A. and M.R.; formal analysis, A.G., R.S. and A.D.; writing—original draft preparation, R.M.; writing—review and editing, M.R. and R.S.; supervision, R.R. and D.G.; funding acquisition, S.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Deanship of Scientific Research, Taif University Researchers Supporting Project number (TURSP-2020/215), Taif University, Taif, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data in this research paper will be shared upon request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anwar, K.; Awadallah, M.A.; Khader, A.T.; Al-Betar, M.A. Hyper-heuristic approach for solving nurse rostering problem. In Proceedings of the 2014 IEEE Symposium on Computational Intelligence in Ensemble Learning (CIEL), Orlando, FL, USA, 9–12 December 2014. [Google Scholar]
  2. Awadallah, M.A.; Bolaji, A.L.; Al-Betar, M.A. A hybrid artificial bee colony for a nurse rostering problem. Appl. Soft Comput. 2015, 35, 726–739. [Google Scholar] [CrossRef]
  3. Constantino, A.A.; Landa-Silva, D.; de Melo, E.L.; de Mendonça, C.F.X.; Rizzato, D.B.; Romão, W. A heuristic algorithm based on multi-assignment procedures for nurse scheduling. Ann. Oper. Res. 2014, 218, 165–183. [Google Scholar] [CrossRef] [Green Version]
  4. Garey, M.R.; Johnson, D.S. Computers and Intractability: A Guide to the Theory of NP-Completeness; W.H. Freeman: San Francisco, CA, USA, 1979. [Google Scholar]
  5. Megeath, J.D. Successful hospital personnel scheduling. Interfaces 1978, 8, 55–60. [Google Scholar] [CrossRef]
  6. Musliu, N.; Gärtner, J.; Slany, W. Efficient generation of rotating workforce schedules. Discret. Appl. Math. 2002, 118, 85–98. [Google Scholar] [CrossRef] [Green Version]
  7. Millar, H.H.; Kiragu, M. Cyclic and non-cyclic scheduling of 12 h shift nurses by network programming. Eur. J. Oper. Res. 1998, 104, 582–592. [Google Scholar] [CrossRef]
  8. Burke, E.; de Causmaecker, P.; Berghe, G.V. A hybrid tabu search algorithm for the nurse rostering problem. In Simulated Evolution and Learning. SEAL 1998; Springer: Berlin/Heidelberg, Germany, 1998; pp. 187–194. [Google Scholar] [CrossRef]
  9. Legrain, A.; Omer, J.; Rosat, S. An online stochastic algorithm for a dynamic nurse scheduling problem. Eur. J. Oper. Res. 2020, 285, 196–210. [Google Scholar] [CrossRef] [Green Version]
  10. Özcan, E.; Bilgin, B.; Korkmaz, E.E. Hill climbers and mutational heuristics in hyperheuristics. In Parallel Problem Solving from Nature-PPSN IX; Springer: Berlin/Heidelberg, Germany, 2006; pp. 202–211. [Google Scholar]
  11. Cheang, B.; Li, H.; Lim, A.; Rodrigues, B. Nurse rostering problems—A bibliographic survey. Eur. J. Oper. Res. 2003, 151, 447–460. [Google Scholar] [CrossRef]
  12. Jaradat, G.M.; Al-Badareen, A.; Ayob, M.; Al-Smadi, M.; Al-Marashdeh, I.; Ash-Shuqran, M.; Al-Odat, E. Hybrid elitist-ant system for nurse-rostering problem. J. King Saud Univ. Comput. Inf. Sci. 2019, 31, 378–384. [Google Scholar] [CrossRef]
  13. Van den Bergh, J.; Beliën, J.; De Bruecker, P.; Demeulemeester, E.; De Boeck, L. Personnel scheduling: A literature review. Eur. J. Oper. Res. 2013, 226, 367–385. [Google Scholar] [CrossRef]
  14. Burke, E.K.; De Causmaecker, P.; Berghe, G.V.; Van Landeghem, H. The state of the art of nurse rostering. J. Sched. 2004, 7, 441–499. [Google Scholar] [CrossRef]
  15. Aickelin, U.; Dowsland, K.A. Exploiting problem structure in a genetic algorithm approach to a nurse rostering problem. J. Sched. 2000, 3, 139–153. [Google Scholar] [CrossRef] [Green Version]
  16. Burke, E.K.; Li, J.; Qu, R. Pareto-Based Optimization for Multi-objective Nurse Scheduling. Computer Science Technical Report No. NOTTCS-TR-2007-5. 2007. Available online: https://www.researchgate.net/publication/277291801_Pareto-Based_Optimization_for_Multi-objective_Nurse_Scheduling_Corresponding_author (accessed on 1 October 2021).
  17. Weil, G.; Heus, K.; Francois, P.; Poujade, M. Constraint programming for nurse scheduling. IEEE Eng. Med. Biol. Mag. 1995, 14, 417–422. [Google Scholar] [CrossRef]
  18. Brusco, M.J.; Jacobs, L.W. Cost analysis of alternative formulations for personnel scheduling in continuously operating organizations. Eur. J. Oper. Res. 1995, 86, 249–261. [Google Scholar] [CrossRef]
  19. Lourenço, H.R.; Martin, O.C.; Stützle, T. Iterated local search. In Handbook of Metaheuristics; Springer: New York, NY, USA, 2003; pp. 320–353. [Google Scholar]
  20. Warner, D.M.; Prawda, J. A mathematical programming model for scheduling nursing personnel in a hospital. Manag. Sci. 1972, 19 Pt 1, 411–422. [Google Scholar] [CrossRef]
  21. Kawanaka, H.; Yamamoto, K.; Yoshikawa, T.; Shinogi, T.; Tsuruoka, S. Genetic algorithm with the constraints for nurse scheduling problem. In Proceedings of the 2001 Congress on Evolutionary Computation (IEEE Cat. No.01TH8546), Seoul, Korea, 27–30 May 2001; Volume 2, pp. 1123–1130. [Google Scholar]
  22. Burke, E.K.; Li, J.; Qu, R. A hybrid model of integer programming and variable neighborhood search for highly-constrained nurse rostering problems. Eur. J. Oper. Res. 2010, 203, 484–493. [Google Scholar] [CrossRef]
  23. Kennedy, J. Particle swarm optimization. In Encyclopedia of Machine Learning; Springer: New York, NY, USA, 2011; pp. 760–766. [Google Scholar]
  24. Gutjahr, W.J.; Rauner, M.S. An ACO algorithm for a dynamic regional nurse-scheduling problem in Austria. Comput. Oper. Res. 2007, 34, 642–666. [Google Scholar] [CrossRef] [Green Version]
  25. Berrada, I.; Ferland, J.A.; Michelon, P. A multi-objective approach to nurse scheduling with both hard and soft constraints. Socio-Econ. Plan. Sci. 1996, 30, 183–193. [Google Scholar] [CrossRef]
  26. Burke, E.K.; Curtois, T.; Qu, R.; Berghe, G.V. A scatter search methodology for the nurse rostering problem. J. Oper. Res. Soc. 2010, 61, 1667–1679. [Google Scholar] [CrossRef] [Green Version]
  27. Todorovic, N.; Petrovic, S. Bee colony optimization algorithm for nurse rostering. IEEE Trans. Syst. Man Cybern. Syst. 2013, 43, 467–473. [Google Scholar] [CrossRef]
  28. Xu, Y.; Wang, X. An artificial bee colony algorithm for scheduling call centres with weekend-off fairness. Appl. Soft Comput. 2021, 109, 107542. [Google Scholar] [CrossRef]
  29. Yang, M.; Ni, Y.; Yang, L. A multi-objective consistent home healthcare routing and scheduling problem in an uncertain environment. Comput. Ind. Eng. 2021, 160, 107560. [Google Scholar] [CrossRef]
  30. Basturk, B.; Karaboga, D. An artificial bee colony (ABC) algorithm for numeric function optimization. IEEE Swarm Intell. Symp. 2006, 8, 687–697. [Google Scholar]
  31. Erkoc, M.E.; Karaboga, N. A novel sparse reconstruction method based on multi-objective Artificial Bee Colony algorithm. Signal Process. 2021, 189, 108283. [Google Scholar] [CrossRef]
  32. Karaboga, D.; Basturk, B. A powerful and efficient algorithm for numerical function optimization: Artificial bee colony (ABC) algorithm. J. Glob. Optim. 2007, 39, 459–471. [Google Scholar] [CrossRef]
  33. Kong, D.; Chang, T.; Dai, W.; Wang, Q.; Sun, H. An improved artificial bee colony algorithm based on elite group guidance and combined breadth-depth search strategy. Inf. Sci. 2018, 442, 54–71. [Google Scholar] [CrossRef]
  34. Naidu, K.; Mokhlis, H.; Terzija, V. Performance investigation of ABC algorithm in multi-area power system with multiple interconnected generators. Appl. Soft Comput. 2017, 57, 436–451. [Google Scholar] [CrossRef]
  35. Ozturk, C.; Hancer, E.; Karaboga, D. Dynamic clustering with improved binary artificial bee colony algorithm. Appl. Soft Comput. 2015, 28, 69–80. [Google Scholar] [CrossRef]
  36. Su, S.; Zhou, F.; Yu, H. An artificial bee colony algorithm with variable neighborhood search and tabu list for long-term carpooling problem with time window. Appl. Soft Comput. 2019, 85, 105814. [Google Scholar] [CrossRef]
  37. Zhao, H.; Pei, Z.; Jiang, J.; Guan, R.; Wang, C.; Shi, X. A hybrid swarm intelligent method based on genetic algorithm and artificial bee colony. In Advances in Swarm Intelligence: ICSI 2010; Springer: Berlin/Heidelberg, Germany, 2010; pp. 558–565. [Google Scholar] [CrossRef]
  38. Wickert, T.I.; Smet, P.; Berghe, G.V. The nurse rerostering problem: Strategies for reconstructing disrupted schedules. Comput. Oper. Res. 2019, 104, 319–337. [Google Scholar] [CrossRef]
  39. Wolbeck, L.; Kliewer, N.; Marques, I. Fair shift change penalization scheme for nurse rescheduling problems. Eur. J. Oper. Res. 2020, 284, 1121–1135. [Google Scholar] [CrossRef]
  40. Zhu, G.; Kwong, S. Gbest-guided artificial bee colony algorithm for numerical function optimization. Appl. Math. Comput. 2010, 217, 3166–3173. [Google Scholar] [CrossRef]
  41. Črepinšek, M.; Liu, S.-H.; Mernik, M. Exploration and exploitation in evolutionary algorithms: A survey. ACM Comput. Surv. (CSUR) 2013, 45, 35. [Google Scholar] [CrossRef]
  42. Vanhoucke, M.; Maenhout, B. NSPLib—A nurse scheduling problem library: A tool to evaluate (meta-)heuristic procedures. In Operational Research for Health Policy: Making Better Decisions, Proceedings of the 31st Annual Conference of the European Working Group on Operational Research Applied to Health Services; Peter Lang AG: Bern, Switzerland, 2007. [Google Scholar]
Figure 1. Visual Representation of Nurse Scheduling Problem.
Figure 1. Visual Representation of Nurse Scheduling Problem.
Mathematics 10 02576 g001
Figure 2. The workflow of NM-ABC.
Figure 2. The workflow of NM-ABC.
Mathematics 10 02576 g002
Figure 3. Analysis and Comparison of Average Best time.
Figure 3. Analysis and Comparison of Average Best time.
Mathematics 10 02576 g003
Figure 4. Analysis and Comparison of Average Standard Deviation.
Figure 4. Analysis and Comparison of Average Standard Deviation.
Mathematics 10 02576 g004
Figure 5. Analysis and Comparison of Average Success Percentage.
Figure 5. Analysis and Comparison of Average Success Percentage.
Mathematics 10 02576 g005
Figure 6. Analysis and Comparison of Average Cost Reduction.
Figure 6. Analysis and Comparison of Average Cost Reduction.
Mathematics 10 02576 g006
Figure 7. Analysis and Comparison of Average Gap.
Figure 7. Analysis and Comparison of Average Gap.
Mathematics 10 02576 g007
Table 1. Instances Taken from dataset with the values of Nurse, Day, and Shift.
Table 1. Instances Taken from dataset with the values of Nurse, Day, and Shift.
CaseTypeInstancesNurseDayShift
1N251, 7, 12, 19, 252574
2N252, 5, 9, 15, 272574
3N251, 3, 16, 27, 352574
4N255, 10, 25, 38, 412574
5N257, 11, 30, 42, 472574
6N501, 4, 12, 26, 295074
7N503, 6,12, 26, 365074
8N504, 9, 15, 40, 475074
9N505, 10, 23, 29, 405074
10N506,14, 20, 32, 415074
11N602, 8, 14, 20, 3260284
12N603, 12, 19, 23, 3460284
13N601, 4, 19, 29, 4060284
14N605, 9, 15, 30, 4360284
15N606, 15, 26, 35, 4460284
Table 2. Parameters setting of compared methods.
Table 2. Parameters setting of compared methods.
TypeMethodParameters and Values
M1Multi-Assignment Problem-based Algorithm (MAPA) [3]Number of Iterations—1000,
Penalty Violation Value—100
M2Hybrid Artificial Bee Colony Algorithm (HABC) [2]Limit—100, Spin Track—10,
Number of Population—100,
Number of Iterations—100
M3Bee Colony Optimization Algorithm (BCO) [27]Knowledge Base (b)—2
Number of Population—100,
Number of Iterations—100
M4Hybrid Elitist–Ant System (HEAS) [12]Population Size—100,
Number of Iterations—1000,
Pheromone Initial Values—0.01
Evaporation Rate—0.25
M5Harmony Search-Based Hyper-Heuristic Algorithm (HSHH) [1]Number of Population—100,
Number of Iterations—100
ProposedArtificial Bee Colony with Nelder–Mead (NM-ABC)Number of Bees—100,
Number of Iterations—1000,
Maximum Runs—20,
Reflection Coefficient—α > 0
Expansion Coefficient—γ > 1
Contraction Coefficient—0 > β > 1
Shrinkage Coefficient—0 < δ < 1
Table 3. Summary of the Best result obtained by proposed Artificial bee colony with Nelder Mead method (NM-ABC) and competitor methods.
Table 3. Summary of the Best result obtained by proposed Artificial bee colony with Nelder Mead method (NM-ABC) and competitor methods.
CaseTypeInstanceOptimal ValueNM-ABCM1M2M3M4M5
C-1N251307307307307307306307
7291287290292290292292
12296296297296296296296
19302302302302303302302
25308307307306307308308
C-2N252274269277279276273273
5303303303302302302301
9276276276276276276275
15296289299300299296300
27293293293293293293293
C-3N251333331333339337332332
3315315315315315317315
16323325322323326326323
27318316318318319318317
35333327330330331333333
C-4N255313285301319310312316
10284273275288277281279
25308294303300305306312
38294294294294294293299
41296296296296296299296
C-5N257293289299299293293296
11299299299299299298299
30311309310319315311310
42283283287283283283283
47310309309315309309309
C-6N501575569579577573575577
4641640641640641640649
12575571578580575575575
26566566569566566566565
29575575572575575573577
C-7N503590512590597595601599
6571571571571571571571
12606600609603608606612
26579574578578579580581
36630630630630630630630
C-8N504644640642645642648645
9571571571571571571571
15580573577589580583585
40562562565562562562561
47562562572562561565562
C-9N6053362329933703362336233683372
103114310731123119311731143123
233476345034793480347534753477
293061302530693061306930613060
402786278627862786278627852786
C-10N6062756275627562756275627562756
143394339033933399339933943400
203441344134413441344134403441
323398339833973400339833983401
413514350435103520351335143520
C-11N6023870387038703870387038703870
83598359835983598359835983598
143703370037053706370837113705
203646364636463646364936463646
323642363936393641364536523641
C-12N6032721272027302721272527292725
122988297629882988299029982992
192988298829882987298829982988
233432342734313432343234323431
343197319731963197319731973197
C-13N6013244324332443244324432493249
42988296929892985299629873001
193136312531413135313931363139
293103310031053103310731023103
402834283428342834283428342834
C-14N6053293327532913299329933003299
92959294529592959296929692959
153063297230633056305830693060
302935287329282934293429342934
432963296529622963296329632963
C-15N6063031299730313040303030353035
153383338333833383338333833383
263969387539783979396939693979
353496349234993496349635003500
443475332834813479347234753473
Table 4. Summary of standard deviation, least error, and best time obtained by proposed Artificial bee colony with Nelder Mead method (NM-ABC) and competitor methods.
Table 4. Summary of standard deviation, least error, and best time obtained by proposed Artificial bee colony with Nelder Mead method (NM-ABC) and competitor methods.
CaseTypeInstanceNM-ABCM1M2M3M4M5
SDLERBest
Time
SDLERBest
Time
SDLERBest
Time
SDLERBest
Time
SDLERBest
Time
SDLERBest
Time
C-1N2511.36013.141.28017.711.17023.141.17020.641.251321.18024.87
C-1N2572.40411.091.27113.821.45121.091.45126.091.85119.141.37119.91
C-1N25121.900.010.00110.050.00010.010.00020.510.0000.270.0002.11
C-1N25191.8000.51.47010.711.9600.51.961121.7901.51.8503.70
C-1N25251.9010.91.96111.012.44230.92.44111.051.4700.761.7302.83
C-2N2523.70519.362.36321.942.32522.52.32223.182.04124.091.9115.28
C-2N2551.7601.541.96012.332.15123.442.24134.212.16135.552.8626.50
C-2N2591.8500.051.99012.191.870301.87020.031.95015.011.81113.63
C-2N25154.43721.22.71321.712.49429.22.34312.82.10023.62.1444.84
C-2N25271.5500.041.9200.361.570251.57031.521.96032.262.1203.42
C-3N2512.5221.21.7001.671.916122.14422.62.14123.31.8914.60
C-3N2531.6000.341.6200.431.87010.341.87010.511.45210.62.1002.18
C-3N25161.79223.071.91130.322.060281.91321.541.62350.771.92049.14
C-3N25272.84215.731.64020.581.990202.10127.871.89033.932.50133.47
C-3N25354.3563.282.6934.322.95334.32.79225.941.7307.271.5608.40
C-4N2559.482832.134.861240.565.426344.32350.072.33159.032.94358.47
C-4N25105.611112.674.57919.344.274293.82735.342.24346.673.53542.74
C-4N25255.821422.53.01530.53.52833.63.35344.851.96256.032.69453.17
C-4N25381.74017.291.95022.051.95019.341.95027.992.07133.332.44533.45
C-4N25411.8901.541.6702.381.6603.671.6604.441.6635.891.7906.77
C-5N2573.12417.322.45610.642.50614.262.27017.921.83023.221.85322.37
C-5N25111.6400.072.1109.691.9002.891.9002.932.1914.352.1505.10
C-5N25302.30215.32.24123.642.2484.32.2446.951.7202.361.8319.41
C-5N25421.9703.291.9548.322.1804.262.1805.912.2407.211.8508.36
C-5N25472.34114.22.20118.282.625652.62167.12.33198.552.02181.70
C-6N5016.14629.46.46415.954.822294.65233.74.58045.855.23241.10
C-6N5044.49116.14.6809.313.88133.973.88027.025.17122.483.69821.41
C-6N50125.06418.454.85310.923.99525.33.56029.534.24017.295.29024.19
C-6N50265.0601.564.08385.624.26023.924.26032.783.72018.393.83117.78
C-6N50293.7701.434.7639.013.31031.783.31022.55.06253.034.2124.45
C-7N50330.487834.674.35044.553.10740.323.10557.663.741169.154.74967.32
C-7N5064.1003.334.0709.353.26034.23.26035.873.87037.133.5208.31
C-7N50125.76659.025.15372.394.58324.24.45218.713.52023.565.82623.15
C-7N50265.44518.454.47130.484.84138.054.84012.284.27114.194.38215.49
C-7N50364.8801.364.37011.844.450224.45042.684.15033.344.7304.68
C-8N5045.16433.94.75244.665.28112.95.28214.855.36445.334.9617.03
C-8N5093.96021.284.67032.313.88034.63.88015.244.22037.224.5807.76
C-8N50157.28732.44.87335.824.56923.94.04020.14.36323.953.71524.48
C-8N50405.07012.254.87323.924.47027.44.47028.534.47061.665.06111.61
C-8N50474.27020.654.051025.334.96037.235.08189.064.04341.764.64012.10
C-9N60528.316332.438.94843.297.79045.297.01061.515.75676.045.761072.16
C-9N60109.70712.346.14216.298.53526.368.63322.536.12027.637.25927.40
C-9N602316.402623.547.26331.135.80431.55.69143.277.34153.145.54151.17
C-9N602922.013643.678.98853.87.19039.986.04861.826.44070.896.53171.41
C-9N60405.4700.375.57010.685.15021.45.15011.596.15112.196.4103.47
C-10N6066.09020.115.75023.495.51032.45.51042.465.73023.636.3304.54
C-10N60147.71442.035.08143.136.50524.786.50535.86.46047.687.2868.35
C-10N60206.40024.656.62030.245.640425.64054.336.05139.166.89018.33
C-10N60327.49020.046.89131.316.96265.96.83025.926.33048.866.8738.70
C-10N60419.671055.098.22479.367.54637.437.17144.985.77059.925.55630.06
C-11N6026.96035.076.24046.66.85056.36.85058.847.07030.726.08011.70
C-11N6085.34027.36.10039.556.28069.266.28052.916.17055.726.60016.36
C-11N60148.05349.766.29251.945.70348.536.18553.417.17875.246.95216.72
C-11N60206.17026.26.56039.196.98032.946.98346.045.99090.967.15020.23
C-11N60327.03351.878.60365.447.83170.437.69376.376.001074.997.79132.05
C-12N6035.66120.236.26932.296.28032.96.28463.027.16844.416.9145.20
C-12N601212.671249.225.88052.386.20043.26.20237.816.551052.116.93422.05
C-12N60198.26021.286.87042.378.26154.898.26035.536.901057.665.5908.11
C-12N60238.72519.97.42134.686.53025.826.37051.277.02056.467.16126.03
C-12N60346.580397.54149.395.680725.50061.57.19072.757.34071.50
C-13N6017.73133.116.36044.336.13045.196.13056.756.08558.566.9659.38
C-13N60412.771947.939.10152.957.89350.457.61849.427.12135.167.111335.11
C-13N60199.931144.439.06549.087.30139.296.00356.516.00062.548.27331.96
C-13N60298.06335.97.35238.857.09061.456.78464.47.44140.656.58018.29
C-13N60406.54032.717.26044.965.95065.235.95057.097.49078.775.7309.73
C-14N60513.001818.437.87223.56.45643.96.80653.128.35735.467.08639.50
C-14N60911.521410.057.49020.67.43052.567.431042.597.151023.856.8504.70
C-14N601528.159149.027.63060.46.617537.44577.517.57679.637.21389.71
C-14N603022.136247.3410.30760.36.05152.65.08176.276.05190.746.38188.39
C-14N60436.02237.16.561346.900496.73028.555.76034.285.22034.15
C-15N60617.403420.4713.08025.4710.86959.97.59130.145.71434.977.24435.72
C-15N60155.85024.236.17035.865.09036.95.09039.025.81051.416.51012.00
C-15N602643.309479.3812.159100.59.611085.346.570125.035.180147.866.7310143.70
C-15N60357.13451.287.05364.026.14059.256.14079.894.77489.26.18413.34
C-15N604457.3414784.36.45696.957.474877.89387.156.81088.586.16297.75
Table 5. Comparison of the number of #both feasible and #feasible solutions obtained by the proposed Artificial bee colony with Nelder Mead method (NM-ABC) and competitor’s methods.
Table 5. Comparison of the number of #both feasible and #feasible solutions obtained by the proposed Artificial bee colony with Nelder Mead method (NM-ABC) and competitor’s methods.
CaseTypeInstancesNM-ABCM1M2M3M4M5
#BFS#FS#BFS#FS#BFS#FS#BFS#FS#BFS#FS#BFS#FS
1N251, 7, 12, 19, 25352434242444
2N252, 5, 9, 15, 27353323233514
3N251, 3, 16, 27, 35143534121335
4N255, 10, 25, 38, 41252523250413
5N257, 11, 30, 42, 47251322343524
6N501, 4, 12, 26, 29250223453512
7N503, 6,12, 26, 36253424223322
8N504, 9, 15, 40, 47351333352223
9N505, 10, 23, 29, 40151233232412
10N506,14, 20, 32, 41353522344522
11N602, 8, 14, 20, 32353434223334
12N603, 12, 19, 23, 34252445332223
13N601, 4, 19, 29, 40152235222422
14N605, 9, 15, 30, 43042524131224
15N606, 15, 26, 35, 44152243242322
Table 6. Comparison and assessment of ABT, ASD, and ASP obtained by proposed algorithm NM-ABC and competitor methods.
Table 6. Comparison and assessment of ABT, ASD, and ASP obtained by proposed algorithm NM-ABC and competitor methods.
CaseTypeInstancesNM-ABCM1M2M3M4M5
ABTASDASPABTASDASPABTASDASPABTASDASPABTASDASPABTASDASP
1N251, 7, 12, 19, 255.131.8710012.661.28017.131.48018.061.48010.731.278018.681.2380
2N252, 5, 9, 15, 278.442.6610013.712.196026.032.086024.352.076026.12.0410024.732.1780
3N251, 3, 16, 27, 358.722.628011.461.9110020.932.168021.692.164025.171.776031.561.99100
4N255, 10, 25, 38, 4117.234.9110022.973.2110023.923.366032.543.0210040.192.058038.922.6860
5N257, 11, 30, 42, 4710.042.2710014.112.196018.142.294020.162.248027.142.0610031.391.9480
6N501, 4, 12, 26, 2913.394.910026.164.974028.794.056029.113.9310031.414.5510033.194.4540
7N503, 6,12, 26, 3623.3710.1310033.724.488031.754.058033.444.024035.473.916037.794.6440
8N504, 9, 15, 40, 4724.15.1510032.414.646027.214.636033.564.5510041.984.494044.64.5960
9N505, 10, 23, 29, 4022.4716.3810031.047.384032.916.896040.146.56047.986.368045.126.340
10N506,14, 20, 32, 4132.387.4710041.516.5110040.56.434040.76.338043.856.0710045.26.5840
11N602, 8, 14, 20, 3238.046.7110048.546.768055.496.738057.516.84065.536.486059.416.9180
12N603, 12, 19, 23, 3429.938.3810042.226.798045.766.5910049.836.526056.686.964060.586.7960
13N601, 4, 19, 29, 4038.829.0110046.037.834052.326.8710056.836.494055.146.838058.896.9340
14N605, 9, 15, 30, 4332.3916.168039.767.9710050.216.698055.616.76052.796.984061.296.5580
15N606, 15, 26, 35, 4451.9326.210064.568.984065.687.836072.256.668082.45.666099.36.5640
Table 7. Comparison and assessment of ACR and AGap obtained by proposed algorithm NM-ABC and competitor methods.
Table 7. Comparison and assessment of ACR and AGap obtained by proposed algorithm NM-ABC and competitor methods.
CaseTypeInstancesNM-ABCM1M2M3M4M5
ACRAGapACRAGapACRAGapACRAGapACRAGapACRAGap
1N251, 7, 12, 19, 251.000.330.200.070.200.070.200.070.000.00−0.20−0.07
2N252, 5, 9, 15, 272.400.83−1.20−0.42−1.60−0.55−0.80−0.280.400.140.000.00
3N251, 3, 16, 27, 351.600.490.800.25−0.60−0.18−1.20−0.37−0.80−0.250.400.12
4N255, 10, 25, 38, 4110.603.555.201.74−0.40−0.132.600.870.800.27−1.40−0.47
5N257, 11, 30, 42, 471.400.47−1.60−0.53−3.80−1.27−0.60−0.200.400.13−0.20−0.07
6N501, 4, 12, 26, 292.200.38−1.40−0.24−1.20−0.200.400.070.600.10−2.20−0.38
7N503, 6,12, 26, 3617.802.99−0.40−0.07−0.60−0.10−1.40−0.24−2.40−0.40−3.40−0.57
8N504, 9, 15, 40, 472.200.38−1.60−0.27−2.00−0.340.600.10−2.00−0.34−1.00−0.17
9N505, 10, 23, 29, 4026.400.84−3.40−0.11−1.80−0.06−2.00−0.06−0.80−0.03−3.80−0.12
10N506,14, 20, 32, 412.800.081.200.04−2.60−0.08−0.80−0.020.200.01−3.00−0.09
11N602, 8, 14, 20, 321.200.030.200.01−0.40−0.01−2.20−0.06−3.60−0.10−0.20−0.01
12N603, 12, 19, 23, 343.600.12−1.40−0.050.200.01−1.20−0.04−5.60−0.18−1.40−0.05
13N601, 4, 19, 29, 406.800.22−1.60−0.050.800.03−3.00−0.10−0.60−0.02−4.20−0.14
14N605, 9, 15, 30, 4336.601.202.000.070.400.01−2.00−0.07−4.40−0.14−0.40−0.01
15N606, 15, 26, 35, 4455.801.61−3.60−0.10−4.60−0.130.800.02−1.60−0.05−3.20−0.09
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Muniyan, R.; Ramalingam, R.; Alshamrani, S.S.; Gangodkar, D.; Dumka, A.; Singh, R.; Gehlot, A.; Rashid, M. Artificial Bee Colony Algorithm with Nelder–Mead Method to Solve Nurse Scheduling Problem. Mathematics 2022, 10, 2576. https://0-doi-org.brum.beds.ac.uk/10.3390/math10152576

AMA Style

Muniyan R, Ramalingam R, Alshamrani SS, Gangodkar D, Dumka A, Singh R, Gehlot A, Rashid M. Artificial Bee Colony Algorithm with Nelder–Mead Method to Solve Nurse Scheduling Problem. Mathematics. 2022; 10(15):2576. https://0-doi-org.brum.beds.ac.uk/10.3390/math10152576

Chicago/Turabian Style

Muniyan, Rajeswari, Rajakumar Ramalingam, Sultan S. Alshamrani, Durgaprasad Gangodkar, Ankur Dumka, Rajesh Singh, Anita Gehlot, and Mamoon Rashid. 2022. "Artificial Bee Colony Algorithm with Nelder–Mead Method to Solve Nurse Scheduling Problem" Mathematics 10, no. 15: 2576. https://0-doi-org.brum.beds.ac.uk/10.3390/math10152576

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