A Fuzzy Crew Rostering Model Based on Crew Preferences and Seniorities considering Training Courses: A Robust Optimization Approach

Crew scheduling problem is divided into crew pairing problem (CPP) and crew rostering problem (CRP). In this paper, a rostering model is presented to assign crew to pairings in such a way that total weighted preference is maximized. Crew members declare which parings they wish to be assigned and which ones are undesirable for them. A score is calculated in the objective function if a crew member is assigned to his/her preferred pairing, and a penalty is considered if he/she is assigned to an undesirable pairing. Moreover, crew seniorities are considered in calculating total preference. In addition, the model considers standard rules and regulations as well as crew attendance at the required training courses. The model is formulated in such a way that inconsistent crew members are not assigned to a flight. Due to the uncertainty in determining of the seniority weight, this parameter is considered as fuzzy. At the end, the robust counterpart of the nominal model is developed due to the uncertainty of time away from the base (TAFB). In this research, the issue of inconsistent crew in rostering problem is considered for the first time. Moreover, a new scoring mechanism is introduced to calculate desirable and undesirable assignments in the objective function. The proposed CRP is solved using the genetic algorithm (GA), and its performance is verified in comparison with GAMS in some test problems. On average, the optimality gap in GA is only 0.5 percent. Finally, the proposed model is examined with real-world data from Air India Airline. In comparison with the previous research studies, the suggested model (scoring mechanism) reduced the number of undesirable rosters by 61.59%.


Introduction
Using the mathematical models in airline operations is an interest field for researchers. Because of the complexity of optimization problems in the airlines, they are usually divided into several smaller problems: flight scheduling (FS), fleet assignment (FA), aircraft maintenance routing (AMR), and crew scheduling (CS) [1,2]. Flight scheduling is the starting point of airline operations [3]. Flight scheduling output is utilized as fleet assignment input. e purpose of fleet assignment is to assign fleet types to flights scheduled in the previous step. Following FA problem, individual aircraft routing is planned so that maintenance requirements are satisfied [4].
One of the major topics in airlines is crew scheduling problem (CSP). Due to its complexity, CSP is separated into two phases: crew pairing problem (CPP) and crew rostering problem (CRP) [5]. Pairings are generated in CPP based on carrier rules and regulations [6]. A pairing includes a set of sequential duty days separated by layovers. e purpose of CPP is to cover all flights in such a way that crew cost is minimized. e pairings created in the CPP are assigned to crew in CRP, and each (co)pilot's roster is determined [7]. Although the CPP focuses on cost minimization, the aim of CRP is to maximize crew satisfaction by assigning balanced workload and meeting crew preferences.
ere are three different approaches in CRP: bidline, personalized with strict seniority, and personalized with global objective [8]. In bidline approach, which was common in North American airlines, anonymous monthly schedules, called rosters, are generated. en, crew members bid for the rosters according to their preference based on seniority. In a personalized approach with strict seniority, which is an emerging approach in North America, crew preferences are met sequentially in decreasing order of seniority. In personalized approach with global objective, crew assignment is done in such a way that total satisfaction is maximized without considering seniority [9].
Crew rostering problem is one of key the phases of airline scheduling. If crew members are assigned to their favorite rosters, a high degree of satisfactions is resulted. Using mathematical models in rostering problem leads to crew satisfaction. In this paper, crew members are assigned to rosters in such a way that maximum satisfaction is reached. To calculate the satisfaction, a new scoring system is used. On the other hand, the issue of travel safety is one of the most important issues for airlines. Assignment of inconsistent crew to a single flight threatens flight safety considerably. As a main motivation, we develop the proposed model to prevent assignment of inconsonant crew to a single flight.
A new model for CRP is formulated in this paper. e objective function maximizes total weighted preference. Instead of using simple bidline or personalized approach, a new mechanism for calculating preference is used in this paper. At first, each crew member declares if he/she prefers to be assigned to a given pairing or not. If they are assigned to one of their preferred pairings, one positive score is calculated in the objective function. Similarly, a negative score (penalty) is considered if a crew member is assigned to an undesirable pairing. e objective function maximizes crew preferences considering seniority. On the other hand, the model considers a crucial qualitative factor in crew assignment to flights. As some of the crew members may be inconstant and have bad work relations to each other, the model prevents to assign them to a single flight. Assignment of the inconsistent crew to a flight may lead to irrecoverable damage. Moreover, crew attendance at the required training courses and some of the airline regulations such as maximum time away from the base (TAFB) and maximum flying hours are considered in assigning crew to rosters. As TAFB is an uncertain parameter in real-world cases, robust counterpart of the proposed model is formulated using the approach of Bertsimas and Sim [10]. A genetic algorithm is used to solve the model for an Indian airline (Air India). GA efficiency in solving the proposed model is examined by an exact solution obtained by CPLEX in GAMS for a variety of small and medium-size examples. is paper benefits from two main contributions. e first one is using a scoring mechanism instead of bidding for preferred pairings. Against previous research studies, both willingness and unwillingness of crew in assigning to each pairing are considered. e second one is considering the issue of inconsistent crew. is important qualitative factor should be considered to ensure a safe travel for passengers. e organization of the paper is as follows. e literature on the crew rostering problem is discussed in Section 2. In Section 3, problem statement and mathematical model are described in detail. Section 4 discusses the solution approach, application of the CRP model in a real-world case study, and the obtained results. Finally, Section 5 concludes with findings and ideas for future studies.

Literature Review
In this section, a comprehensive review on the CRP is presented. Different research studies have been classified based on if their contributions are in solution approach (heuristic, metaheuristic, and exact) or development of the mathematical model. Although many studies have concentrated on developing solution approaches, only a few ones have contribution in mathematical formulation.

Development of the Solution Approach.
Boubaker et al. [11] proposed two heuristic algorithms to solve the CRP. e first approach was the standard branch-and-price algorithm, and the second one was a combination of Elhallaoui et al. [12] with the branch-and-price. Computational findings indicated that the combined method for large-scale problems performed better than standard ones. Maenhout and Vanhoucke [13] considered the use of the additional crew and developed a scatter search algorithm to solve the CSP. e result of comparing the proposed approach with the branch-and-price algorithm and variable neighborhood search demonstrated that this technique could result in highquality solutions. Saddoune et al. [14] discussed the weaknesses of the traditional approaches, solved the CRP using the dynamic column generation method, and implemented it in an American airline. An integrated approach for aircraft and crew schedule recovery problems was proposed by Zhang et al. [15]. To solve the real problem, they used a twolevel metaheuristic algorithm and showed that it could lead to high-quality solutions. Zhang et al. [16] used the Tabu search algorithm to solve the CRP and used the concept of search memory, which accelerated the convergence. By testing the samples, the findings indicated that their method achieves near to optimal solutions in a reasonable time. Doi et al. [17] addressed obtaining a fair roster for the crew by applying a two-level decomposition approach. Zhang et al. [18] developed three multiobjective evolutionary algorithms to search for Pareto solutions. e experimental results on a Chinese airline proved that the proposed technique was capable to improve the rosters' quality by balancing the crew workload. Zhang et al. [19] presented an implicit mathematical model using bidline approach and developed a new variable neighborhood search to solve it. e results indicated that the proposed approach could provide high-quality solutions for two different scenarios of the crew workload. Chutima and Arayikanon [20] formulated a CSP in a low-cost airline and solved it using a multiobjective evolutionary optimization algorithm (MOEA). Comparing the proposed algorithm results with the honeybee mating optimization (HBMO) indicated that MOEA was superior in terms of convergence criteria. e improved discrete particle swarm optimization (IDPSO) method was proposed by Zheng [21] to solve CSP. e findings showed that IDPSO's performance was very efficient. Zhou et al. [22] improved an ant colony algorithm to solve the CRP in a biobjective model optimizing fairness and satisfaction. ey proposed a hybrid complementary heuristic strategy and a local search strategy to solve the model adequately. ey showed that the novel method led to higher-quality solutions than the greedy algorithm and other prominent multiobjective optimization techniques. Chutima and Krisanaphan [23] considered crew cost minimization as well as balancing the workload of pilots. ey solved the proposed model with the adaptive nondominated sorting differential evolution algorithm III (ANSDE-III) and compared the results with the NSGA-III and the multiobjective differential evolutionary approach (MODE).
e results showed that ANSDE-III performed better than other algorithms. [24] formulated pilot roster recovery problem with the aim of minimizing the maximum flight time of each pilot after a disruption and reduced exchange costs. ey solved the model with the nondominated sorting genetic algorithm II (NSGA-II) to find Pareto solutions and applied the constraint-loosening mechanism (CLM). e results indicated that the proposed method led to high-quality solutions to the CR problem. Crew changeability in case of delay and disruption was offered by Ilagan and Sy [25]. According to the results, although the proposed model had increased costs, the delay had been significantly reduced. Kasirzadeh et al. [26] formulated the personalized CR problem with undetermined pairings and fixed pilots. e model aimed to optimize crews' preferences and their cost. ey used the approach of Saddoune et al. [14] and obtained an acceptable level of crew satisfaction. Armas et al. [27] proposed the CRP with the aim of fair distribution of workload among crew members, and they applied regulations that had not been considered in previous studies. To solve this model, they used the multistart randomized heuristic algorithm and showed that this algorithm has the ability to solve the problem with high quality. An integrated model for CSP was presented by Zeighami and Soumis [8].

Development of CR Models. Chen and Chou
ey considered pilots' vacation requests in the CP to generate better pairings.
e authors formulated the model by minimizing pairing costs and maximizing the desired number of vacations.
e experiment results showed that combining Benders' decomposition algorithm and column generation could significantly save CPU time. Quesnel et al. [28] considered rewards for the pairings that matched the crews' preferences in the objective function. ey solved the model via column generation and implemented it to a North American airline. e findings demonstrated that the suggested model could improve CRP solutions compared to previous studies. Mirjafari et al. [29] presented an integrated model for aircraft routing and CSP with the aim of minimizing crew, aircraft replacement, maintenance, and deadhead flight costs. ey suggest a new frame for maintenance operations based on flight hours. Because the Lagrangian relaxation method results have a lower gap to an optimum solution, this method is better than the particle swarm optimization (PSO) algorithm. A new integrated model for fleet assignment and crew scheduling problems was proposed by Rashidi Komijan et al. [1]. ey considered closed routes for crews and fleets simultaneously. To compare the effectiveness of the two methods, PSO and vibration damping optimization (VDO) algorithm solved ten large-scale examples. e findings showed that VDO leads to optimal results in a reasonable time. Saemi et al. [30] suggested an integrated CSP to minimize crew costs. e approach presented a better answer than computing CPP and CRP consecutively. It was also possible to add a day(s) off in a pairing to allow crew members to attend training courses and complete health checks. A multiobjective personalized model for the airline multiskilled was proposed by Baradaran and Hosseinian [31]. e multiskilled crew used in this model was assigned to flights and aircraft due to the crew's experience. ey used multiobjective differential evolution (MODE) and nondominated sorting genetic algorithm II (NSGA-II) to analyze the model. According to the results, the suggested model and algorithms could generate adequate schedules for CS problems, according to comparisons between algorithms. Ben Ahmed et al. [32] assigned the aircraft and crew to each flight simultaneously and included maintenance constraints in their model. ey solved the robust mixed-integer programming model using the decomposition approach and achieved good results. Table 1 briefly compares the proposed CRP model with the previous research studies.
To our best knowledge and Table 1, previous research studies mainly focused on cost minimization, fair workload distribution, and use of bidline and personalized approaches for crew rostering. is paper includes two main contributions. e first one is related to the issue of inconsistent (co)-pilots. is critical qualitative component should be evaluated to ensure passenger safety. e model prevents to use inconsistent crew in a single flight. e second contribution is that instead of suggesting preferred pairings, a scoring method is used. In contrast to previous studies, this strategy takes into account both the crew's willingness and unwillingness to assign to each pairing. In other words, one positive score is calculated for each favorable allocation and one negative score is considered for each undesirable allocation in the objective function. Total weighted preference is maximized by considering crew seniority.

Problem Statement
In this paper, a new mathematical model for the CRP is proposed. e objective function calculates the total weighted preference of the crew. We assume that each crew member declares his/her preferred and undesirable pairings.

Computational Intelligence and Neuroscience
By assigning a crew member to his/her preferred pairing, a positive score is calculated in the objective function. Similarly, if a crew member is assigned to his/her undesirable pairing, a negative score (penalty) is calculated. e objective function maximizes total weighted preference by considering crew seniority.
In each airline, some of the crew members may have bad relations with each other. Assigning inconsistent crew to a single flight may threaten flight safety. To cover this qualitative factor, the model should prevent to assign inconsistent crew to a flight. Also, several rules and regulations which are common in airlines are considered in the model. Some of these rules are as follows: (i) Time away from the base should not exceed a predefined value (ii) Maximum and minimum flying hours of crew are restricted (iii) An inexperienced co-pilot must be accompanied by an experienced pilot (iv) Some of the crew members should attend at training courses Due to the uncertainty of time away from the base (TAFB), the approach of Bertsimas and Sim [10] is used to formulate the robust counterpart. e advantage of using the approach of Bertsimas and Sim [10] is that the robust model will be linear. Also, the most well-known Indian airline, Air score cp , penalty cp ≥ 0.
e objective function maximizes total weighted preference of the crew. e way of calculating score and penalty is described in the constraints. Equations (2) and (3) ensure that each pairing is handled by one pilot and one co-pilot. Constraint (4) prevents assignment of inconsistent crew to a pairing. According to (5), the TAFB for each crew should not exceed the maximum allowed value. Relation (6) indicates that the total flying hours for a crew member should fall in the defined interval. Some pairings cannot be handled sequentially due to some reasons such as the destination of the first one is not the same as the origin of the next one. Constraint (7) prevents the assignment of a crew member to such pairings. Constraint (8) states that in each pairing, an inexperienced copilot must be accompanied with an experienced pilot. Relation (9) is about attendance of crew at the training course or meetings. Constraint (10) prevents concurrence of attendance at the training course and handling a pairing. Constraints (11) and (12) show the way of calculating score and penalty, which have been used in objective function. According to (11), one score is calculated if a crew member is assigned to his/her preferred pairing. Similarly, if he/she is assigned to an undesirable pairing, a penalty is considered. Constraint (13) ensures that a crew member can only be assigned to the pairings that origin from his/her home base. Finally, Constraints (14) and (15) indicate decision variable types. e seniority weight of pilots is in the range [0, 1]. Since the nature of this parameter is fuzzy, the Jimenez et al. [33] method has used to convert this uncertain coefficient in the objective function to definite. For more information, refer the study by Jimenez et al. [33]. For this matter, the triangular fuzzy distribution has chosen because of its simplicity and efficiency. After performing the calculations, the following equation replaces the indefinite objective function (1).
As TAFB is an uncertain parameter, the robust counterpart of the nominal model is formulated using the approach of Bertsimas and Sim [10]. TAFB includes total flying time, sit times, and layovers. TAFB p and TAFB p have assumed the nominal value and tolerance of time away from the base for pairing p. According to Bertsimas and Sim [10], protection level (Γ) is defined to make the model robustness. In other words, the solution's robustness is regulated by defining a protection level, which specifies a maximum deviation from the nominal TAFB. e main purpose of this parameter is to make the solutions generated workable for all uncertain scenarios. If at most Γ i of technical coefficients of the constraint i change in the defined interval, the model will definitely be robust. Also, a technical coefficient can be changed up to(Γ i − Γ i without violating the relation. Moreover, σ cp and θ cp are variables defined to formulate robust counterpart. To formulate robust counterpart, constraint (5) is replaced by the following constraints.

Solution Approach
According to NP-hard nature of the crew rostering problem, an appropriate optimization strategy is required [34]. Many rules and policies may deploy for crew rostering with Computational Intelligence and Neuroscience thousands of variables in a real-world problem. Furthermore, the presence of binary variables makes the solution approach more complicated.
Various papers in the literature use genetic algorithm to solve the airline crew scheduling problem. e GA is the second commonly applied solution technique and first metaheuristic approach for the CSP after column generation [35]. One of the benefits of the GA is the ability to manage numerous solution search areas [36]. As a consequence, the genetic algorithm is deployed as solution strategy for this study. In this section, problems are solved by a genetic algorithm in different sizes to evaluate the efficiency of the proposed approach. e GA outcomes are assessed and compared to the CPLEX solver in GAMS. Also, the solution method presented in this work and the components of the algorithm utilized are highlighted in Section 4.1.

Presentation of Solution.
e genetic algorithm uses the concept of chromosome to set parameters that offers a proposed problem solution. In practice, the chromosome is a string of solutions to solve the problem. e proposed chromosome can be designed as a string of discrete variables, binary numbers, and a continuous vector based on type of the problem. As a result, choosing an appropriate display form of the chromosome is an essential aspect of algorithm design. e designed chromosome for the proposed CRP is made up of the two-row vector. e number of columns in this vector is equal to the number of pairings |P|. e first and second rows specify the pilot number and the co-pilot number, respectively. Pilot Co-pilot Figure 1: Chromosome with 5 co-pilots and 6 pairings.

Signal-to-noise: Smaller is better
Maxlt npop pc pm

Generating the Initial Solution.
In this paper, the copilots are randomly assigned to each pairing. erefore, each compatible co-pilot gets selected to operate in pairing P. As an example, 6 pairings, 5 pilots, and 5 co-pilots are available (pilots are numbered from 1-5, and co-pilots are numbered from 6-10). An initial solution for the co-pilots assigned to the pairing is shown in Figure 1. To avoid assigning of a crew member to nonconsecutive pairings, Constraint (7) is applied to prevent the generation of infeasible solutions.
Pilot assignment process begins after all co-pilots get placed in their pairings. Each eligible pilot may get assigned to pairing P randomly, if the following conditions met: (1) e pilot should not be in conflict with the co-pilot (2) An inexperienced co-pilot must be accompanied by an experienced pilot in each pairing (3) Set of flights which may be assigned to a specific pilot in a sequence e pilot row is determined according to Figure 2 for the chromosome.

Generating a New Solution.
A suitable operator for generating a random neighborhood solution depends on to the nature of the problem. e CRP is an allocation problem; therefore, we considered the appropriate operators to generate new solutions. ese operators are defined in the following section.
(1) Mutation Operator. By using the swap operator, two pilots or two co-pilots are selected from two pairings and get swapped if the following conditions are met: (i) e co-pilot and the pilot are inconsistent (ii) If an inexperienced co-pilot is not paired with an experienced pilot (iii) Pairings that may not be done sequentially by one pilot or co-pilot For the mutation operator, a random number between 0 and 1 is generated. If the number is less than 0.5, a swap operation is performed on the co-pilots, and if it is greater than 0.5, it is performed on the pilots. Swapping the co-pilots of pairings P 5 and P 3 is presented in Figure 3.
(2) Crossover Operator. Figure 4 shows the single crossover operator for the CRP. Selecting the pilot/co-pilot row is similar to Section 4.1.3.1 for crossover operator as well.
e proposed chromosome in this paper complies with all constraints of the model except for the maximum/ minimum flight hours for crews and maximum time away from the base for crews. To tackle this problem, a penalty function strategy is applied. erefore, the average violation 1400 1200 100 -C 10 1400 1200 100 - Table 5: Crew assignment to pairings.

Pairings Crews Pilot
Co-pilot Computational Intelligence and Neuroscience 7 of these constraints is multiplied by a big number and then subtracted from the objective function.

Genetic Parameter Tuning.
In the genetic algorithm, maximum iterations, population size, crossover rate, and mutation rate parameters affect the performance of the algorithm. In this study, the Taguchi [37] is applied to get the best value for mentioned parameters. Each of the GA parameters is considered at three levels, as presented in Table 2.
Not all but only part of the factors may reach to their best level in Taguchi's [37] method. Each obtained result from Taguchi [37] method experiment is converting to a "signalto-noise" ratio. At this rate, the optimal value (average) is called the signal, and the undesirable value (standard deviation) is called the noise. According to the objective function which is to maximize the crew's satisfaction, the higher S/N means more desirability of results. erefore, the maximum point of the graph is selected for each parameter. Also, the corresponding level is considered as the optimal  e estimated values for GA parameters using the Taguchi [37] method are illustrated in Table 3 and Figure 5.

Computational Results for Test Problems.
e model is solved for a real case: Air India airline. Before solving the case study, 30 test problems in different scales are generated and solved using GAMS and GA to justify GA's efficiency in solving the proposed model. is operation is performed on a laptop with 16 GB RAM and Windows 10 64 bit. To illustrate the model results, a problem with 10 crew members is explained in detail. In this example, C 1 ′ , C 2 ′ , and C 3 ′ are assumed experienced pilots, while C 6 ′ and C 7 ′ are inexperienced co-pilots. Table 4 shows the data for this test problem. By solving the model, crew assignment is done. According to Table 5, one pilot and one co-pilot who are not inconsistent are assigned to each pairing. Crew's rosters can be easily concluded from Table 5. To illustrate, Figure 6 displays the first pilot's monthly roster. As the planning horizon is one    Table 6. In solving numerical examples by GAMS, the CPU time has been set to 10,800 seconds. According to Table 5, the performance of the GA to solve the proposed model is justified due to two reasons: (1) the maximum optimality gap between the genetic algorithm and exact solution is only 1.482%, and (2) in medium and large-sized problems, the CPU times in the GA are considerably less than GAMS.
A comparison between the CPU times in GA and GAMS is shown in Figure 7. As maximum solution time has been set to 10,800 seconds in GAMS, computations in problems 14-30 are terminated after three hours.
A comparison of the results between the results of GAMS and genetic algorithm is shown in Figure 8. e results of the objective function are shown in Figure 8. e green and orange points represent the objective function values resulted from GA and GAMS, respectively. It is worth noting that GAMS solution time has been limited to three hours. As a result, GAMS did not reach the optimum solution in problems 14 to 30.

Case Study.
Air India is a well-known airline in India, headquartered in New Delhi. As one of the biggest international airlines, Air India has 11% of the domestic market share [38]. Air India's main hub is Indira Gandhi International Airport in New Delhi, with a secondary hub at Chhatrapati Shivaji Maharaj International Airport in Mumbai. According to reports in 2019, Air India fleets include Airbus and Boeing that serve 170 destinations all over the world [32].
To solve the proposed model for the case study, a genetic algorithm was applied and the Taguchi [37] approach was used to adjust its parameters.
is case included 6190 pairings and 1340 (co)-pilots. e model was solved in 761.188 seconds, and the optimum value of the objective function was 5034.002. e convergence of the GA results is shown in Figure 9. Converging after 400 iterations indicates good solution quality.
As discussed in Section 3, the objective function includes two terms: score and penalty for assigning crew to desirable/ undesirable pairings. According to the Air India results, crew members were assigned to their preferred pairings in 72% of assignments. Table 7 shows the monthly rosters of 20 experienced pilots, and 20 co-pilots as a small sample. Moreover, the attendance days of 90 co-pilots in the required training courses are presented in Table 8.
To compare the performance of the proposed model with the previous research studies in which undesirable assignment penalty had not been considered, penalty term was removed from the objective function and constraint (12) was also disabled. By removing penalty concept, 2065 undesirable pairings are assigned to crews. is means that each crew will have an average of at least two undesirable pairings per month. However, our proposed model has reduced the number of undesirable rosters to 739, which is about 61.59% less than the previous models.

Robust Optimization.
In this section, Bertsimas and Sim [10] approach is used to formulate the robust counterpart of the proposed model. e approach of Bertsimas and Sim leads to a linear robust model, while it is independent of knowing the distribution of data and control of the level of protection through parameter Γ [39].
TAFB is subject to several factors such as flight time, flight duration, and connection time, which makes it an uncertain parameter. By considering different values for Γ (0 ≤ Γ ≤ 1), the robust model was solved for numerical

Computational Intelligence and Neuroscience
′ × example 3, and the results of GAMS and GA were compared. After justifying the performance of the GA, the robust model was applied to the case study and the results were analyzed ( Table 9).
As represented in Table 9, the optimality gap of GA for numerical example 3 is negligible. erefore, the performance of GA to solve the robust model is acceptable. According to the robust model applied to the case study, the model acts more conservatively by increasing of the Γ. In the case of Γ � 1, it guarantees that feasible solution to be at its maximum value and the objective function to be at its minimum value. Figure 10 shows the total weighted desirability obtained from the genetic algorithm for 400 iterations.

Conclusions
Crew scheduling problem is often divided into CPP and CRP and solved consecutively to reduce complexity. In this article, a novel MIP model was introduced for CRP. e difference between the proposed model and previous research is considering a scoring mechanism instead of bidding for preferred pairings, crew seniorities, and inconsistent crew. In addition, the model considers common rules and regulations as well as crew attendance at the required training courses. With this advantage, obviously, better solutions are obtained. Moreover, the robust counterpart of the nominal CRP model was presented to the uncertainty of mission time of the crews (TAFB). In addition, the suggested model's efficiency was assessed using real-world data from Air India Airline. Because the model was NP-hard, the CPLEX solver in GAMS was used for 13 first test problems, and the GA was utilized for 30 test problems with different scales and the case study. e results show that GA has a 1.482% average gap in optimum solution compared to exact solution. e proposed CRP is computed utilizing Air India data, which is one of the largest international carriers. Air India provides services to 170 destinations in 31 countries. e proposed model was tested for 6190 pairings and 1340 crews. For this real-world data, the objective function value was 5034.002 with 761.188 seconds of computational time. is research has two limitations that can be improved in future research studies. e first one is that (co)-pilots' medical checks have not been considered. e second one is that reserve crew have not been included in the model. Moreover, the following suggestions can be considered for future research [40][41][42]: (i) Considering a fair distribution of working time among crew members (ii) Providing an integrated model for CRP with other steps of the airline scheduling problem, and considering disruption for the crew rostering (iii) Developing the cooperative game between airlines to exchange crews among them

P:
Set of all pairings p: Index of pairing (p � 1, 2, . . .) C: Set of all crew members (pilots and co-pilots) c, c ′ : Index of pilots and co-pilots (c, c ′ �1, 2, . . .) C 1 : Set of all pilots c 1 : Index of pilots (c 1 � 1, 2, . . ., C 1 ) C 2 : Set of all co-pilots c 2 : Index of pilots (c 2 � 1, 2, . . ., C 2 ) C 1 ′ : Set of all experienced pilots c 1 ′ : Index of experienced pilots (c 1 ′ �1, 2, . . ., C 1 ′ ) C 2 ′ : Set of all inexperienced co-pilots c 2 ′ : Index of inexperienced co-pilots (c 2 ′ � 1, 2, . . ., C 2 ′ ) C inc : Set of inconsistent crew C inst : Set of crew who must participate at the training course P d : Set of all pairings that include day d NA: Set of all pairings that cannot be performed sequentially by a (co)-pilot D: Set of all days in which the training course is held d: Index of training course (d �1, 2, . . .) Binary parameter that equals 1 if pairing p is desirable for (co)-pilot c, and 0 otherwise Base cp :

Parameters
Binary parameter that equals 1 if pairing p begins and ends at the home base of (co)-pilot c, and 0 otherwise w c : Seniority weight of (co)-pilot c Decision variables X cp : Binary variable that equals 1 if pairing p is assigned to (co)-pilot c, and 0 otherwise Y cd : Binary variable that equals 1 if (co)-pilot c attends at training course held on day d, and 0 otherwise score cp : Score of assigning (co)-pilot c to his/her preferred pairing p penalty cp : Penalty for assigning (co)-pilot c to the undesirable pairing p. Computational Intelligence and Neuroscience 13 Data Availability e data used to support this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.