A new human-based metaheuristic algorithm for solving optimization problems based on preschool education

In this paper, with motivation from the No Free Lunch theorem, a new human-based metaheuristic algorithm named Preschool Education Optimization Algorithm (PEOA) is introduced for solving optimization problems. Human activities in the preschool education process are the fundamental inspiration in the design of PEOA. Hence, PEOA is mathematically modeled in three phases: (i) the gradual growth of the preschool teacher's educational influence, (ii) individual knowledge development guided by the teacher, and (iii) individual increase of knowledge and self-awareness. The PEOA's performance in optimization is evaluated using fifty-two standard benchmark functions encompassing unimodal, high-dimensional multimodal, and fixed-dimensional multimodal types, as well as the CEC 2017 test suite. The optimization results show that PEOA has a high ability in exploration–exploitation and can balance them during the search process. To provide a comprehensive analysis, the performance of PEOA is compared against ten well-known metaheuristic algorithms. The simulation results show that the proposed PEOA approach performs better than competing algorithms by providing effective solutions for the benchmark functions and overall ranking as the first-best optimizer. Presenting a statistical analysis of the Wilcoxon signed-rank test shows that PEOA has significant statistical superiority in competition with compared algorithms. Furthermore, the implementation of PEOA in solving twenty-two optimization problems from the CEC 2011 test suite and four engineering design problems illustrates its efficacy in real-world optimization applications.

Alongside exploration and exploitation, achieving the desired performance of metaheuristic algorithms relies on striking a balance between these two aspects during the search process 6 .
Due to the random nature of their search process, metaheuristic algorithms do not guarantee finding the global optimum for optimization problems.Consequently, the solutions provided by metaheuristic algorithms for optimization problems are termed quasi-optimal.The pursuit of improved quasi-optimal solutions for optimization problems has spurred the development of numerous metaheuristic algorithms 7 .These algorithms are employed to handle optimization tasks in various sciences such as combined heat and power economic dispatch 8 , solving general systems of nonlinear equations 9 , numerical optimization problems 10 , semi-submersible platform boom 11 , dynamic positioning system (DPS) 12 , auto drum fashioned brake design 13 , search the optimal parameters for a bucket wheel reclaimer (BWR) 14,15 , medica applications 16 , and Feature Subset Selection (FSS) 17,18 .
The main research question is, according to the countless metaheuristic algorithms introduced so far, what is the primary motivation for introducing newer algorithms based on the need?The No Free Lunch (NFL) theorem 19 provides a definitive explanation to address this question.The NFL theorem states that achieving acceptable performance with a metaheuristic algorithm for a specific set of optimization problems does not guarantee similar performance for other optimization problems.An algorithm that has shown success in solving particular optimization problems may fail when applied to others.The NFL theorem highlights that no single metaheuristic algorithm can claim to be the best optimizer for all optimization problems.The NFL theorem serves as a catalyst for ongoing research in the field of metaheuristic algorithms, inspiring researchers to continually innovate and devise more efficient solutions for optimization problems by developing novel algorithms.According to this, the author of this paper, by motivation from the NFL theorem and based on the simulation of human activity in the process of preschool education, has designed a new metaheuristic algorithm to deal with optimization tasks in science.
Our extensive literature review shows no metaheuristic algorithm inspired by the concept of preschool education has yet been developed.This is even though creating an educational and fostering environment in preschool has typical characteristics of intelligent decision-making and optimizing the process.In this study, we aim to bridge this research gap by introducing a novel metaheuristic algorithm named the Preschool Education Optimization Algorithm (PEOA), which draws inspiration from the concept of preschool education.
As mentioned, to provide an effective search process in the problem-solving space, a metaheuristic algorithm must have a high ability in exploration, exploitation, and balancing between them during the search process.In the design of PEOA, by taking separate phases of updating the position of population members to manage exploitation, exploration, and balancing exploitation and exploration, an effort has been made to achieve a powerful and effective search process in the problem-solving space to achieve suitable solutions for optimization problems.
In the design of PEOA, the exploitation ability to manage local search is modeled based on the simulation of the gradual growth of the preschool teacher's educational influence.In this process, the child gradually develops under the teacher's influence.Modeling this gradual learning process by making small changes in the position of PEOA members in the problem-solving space leads to an increase in the exploitation ability of PEOA to manage the local search in the accurate scanning of the problem-solving space near the discovered solutions and promising areas to find better solutions.According to this, PEOA is expected to be effective in exploitation for local search in the problem-solving space.
In the design of PEOA, the exploration ability to manage the global search is modeled based on the simulation of individual knowledge development guided by the teacher.In this process, based on imitating the teacher, the child tries to learn the lesson taught by the teacher.Modeling this learning process by making extensive changes in the position of PEOA members in the problem-solving space leads to an increase in the exploration ability of PEOA to manage the global search in the comprehensive scan of the problem space to prevent the algorithm from getting stuck in local optima and identifying the region containing the global optimum.According to this, PEOA is expected to effectively explore global search in the problem-solving space.
In PEOA design, the simulation of Individual increase of knowledge and self-awareness positively affects the ability to exploit the algorithm for local search.Modeling this process by making small changes in the position of PEOA members leads to improving the algorithm's exploitability to manage the local search.
On the other hand, in the design of PEOA, to manage exploration and exploitation and establish a balance between them during the search process, priority has been given to exploration in the initial iterations so that by making extensive changes in the position of population members, the problem-solving space can be scanned well and the promising areas be identified.Then, by increasing the iterations of the algorithm, priority has been given to exploitation so that by shrinking the range of changes in the position of the population members in the problem-solving space, the algorithm can achieve more effective solutions for the given problem by accurately scanning promising areas.Therefore, the proposed PEOA approach is expected to perform well in exploration, exploitation, and balancing during the search process in the problem-solving space to achieve suitable solutions for optimization problems by managing an effective search process.
The aspects of innovation and novelty of this paper are in the introduction and design of a new human-based metaheuristic algorithm named the Preschool Education Optimization Algorithm (PEOA), which draws its inspiration from the preschool education process.The key contributions of this research are outlined as follows: • The development of PEOA is grounded in the concept of preschool education.
• PEOA is mathematically modeled through three distinct phases: (i) the gradual growth of the preschool teacher's educational influence, (ii) individual knowledge development guided by the teacher, and (iii) individual increase of knowledge and self-awareness.
• The efficacy of PEOA in solving optimization problems is assessed using fifty-two standard benchmark func- tions encompassing unimodal, high-dimensional multimodal, and fixed-dimensional multimodal types, as well as the CEC 2017 test suite.• A comprehensive comparative analysis is carried out to assess the performance of PEOA with ten widely recognized algorithms.• The practical applicability of PEOA is demonstrated by applying it to twenty-two optimization problems from the CEC 2011 test suite and four engineering design problems, showcasing its effectiveness in real-world scenarios.
The subsequent sections of the paper are thoughtfully structured to present the literature review in the "Literature Review" section, followed by the theoretical framework and mathematical model of the proposed optimizer in the dedicated "Preschool Education Optimization Algorithm" section.The "Simulation Studies and Results" section provides a concise summary of the simulation studies conducted and the corresponding outcomes.The implementation of PEOA in solving real-world applications is presented the "PEOA for real-world applications" section.Conclusions and several proposals for further research are provided the "Conclusion and future works" section.

Literature review
Metaheuristic algorithms have been developed by taking inspiration from a variety of sources, such as natural phenomena, animal behaviors, biological sciences, physical laws, human interactions, and game rules.These algorithms can be categorized into five main groups based on their fundamental design principles: swarm-based, evolutionary-based, physics-based, human-based, and game-based approaches.Each category represents a distinct approach to problem-solving, leveraging different concepts and techniques.
Swarm-based algorithms draw inspiration from the collective behavior of various organisms in nature, including birds, animals, aquatic creatures, insects, and more.Prominent examples of swarm-based approaches extensively employed for solving optimization problems include Ant Colony Optimization (ACO) 20 , Particle Swarm Optimization (PSO) 21 , Artificial Bee Colony (ABC) 22 , and Firefly Algorithm (FA) 23 .PSO is designed based on the swarming movement observed in flocks of fish and birds as they search for food sources in their environment.ACO leverages the ability of ants to find the optimal route between a food source and their nest.The foraging activities of honey bee colonies inspired the design of ABC.The flashing light behavior exhibited by fireflies, which serves to attract mates and prey through bioluminescence, forms the basis of FA's strategy.
Optimization (ASO) 31 , Reptile Search Algorithm (RSA) 32 , Conscious neighborhood-based Crow Search Algorithm (CCSA) 33 , Quantum-based Avian Navigation optimizer Algorithm (QANA) 34 , and Starling Murmuration Optimizer (SMO) 35 .These algorithms emulate the swarm intelligence exhibited by living organisms in their search for and acquisition of food resources.
Evolutionary-based algorithms are inspired by biological sciences, genetics, concepts of natural selection, and stochastic operators.Genetic Algorithm (GA) 36 and Differential Evolution (DE) 37 are among the most wellknown evolutionary-based approaches that are inspired by the reproduction process, the concepts of Darwin's theory of evolution, and the evolutionary operators of selection, crossover, and mutation.
Human-based algorithms are inspired by the mutual communication and interactions of humans in social and individual life.Tabu Search (TS) 54 creates a tabu list to keep track of recently explored solutions and prevent revisiting them, promoting diverse exploration.The algorithm iteratively generates neighboring solutions, evaluates their fitness, and updates the tabu list accordingly.By incorporating aspiration criteria, Tabu Search can escape local optima.
Teaching Learning Based Optimization (TLBO) 55 is a popular human-based approach inspired by the teaching and learning dynamics between teachers and students in a classroom setting.
In Queueing Search (QS) 56 algorithm, typical occurrences involve customers actively choosing fast-service queues, where individual customer service is primarily impacted by staff or the customer themselves.Additionally, others may influence customers during service when the queue order lacks strict adherence.Some other human-based metaheuristic algorithms are: Poor and Rich Optimization (PRO) 57 , Human Mental Search (HMS) 58 , Multi-Leader Optimizer (MLO) 59 , Following Optimization Algorithm (FOA) 60 , Teamwork

Inspiration and main idea of PEOA
Preschool education plays a crucial role in a child's early development and lays the foundation for their future learning journey.Attending nursery school provides young children with numerous benefits that contribute to their overall growth and well-being 66 .
One of the key advantages of preschool education is the enlargement of the opportunity for social interaction.Children at this age are naturally curious and eager to explore their surroundings.Preschool offers a nurturing environment where they can engage with peers and develop essential social skills.Children are exposed to different subjects, ideas, and challenges through various activities and play-based learning.Conversations with teachers and peers help them articulate their thoughts and emotions effectively, boosting their self-confidence and self-expression ability 67,68 .
In the realm of preschool education, the role of a teacher extends beyond mere instruction, encompassing a dynamic interplay of intelligent processes that shape young minds.Preschool teachers must solve intricate interactions and involve adaptive strategies to enable the full complexity of education in fostering holistic development.
A preschool teacher acts as a guiding force, steering children's curiosity and exploration toward constructive paths.Through structured activities and open-ended play, the teacher creates an environment where children can interact intelligently with their peers, stimulating cognitive growth and social adeptness.
The teacher, attuned to the unique needs of each child, facilitates this journey by encouraging self-expression, supporting decision-making, and promoting autonomy.By fostering such intelligent processes, the preschool teacher empowers children to embrace their individuality and develop a strong sense of identity.
Furthermore, the preschool teacher guides children through challenges, triumphs, and fails and mainly nurtures resilience and problem-solving skills.This echoes the notion that intelligent processes are at play, enabling children to overcome obstacles and emerge stronger.
In sum, the role of a preschool teacher extends far beyond conventional instruction.The teacher cultivates an environment where children learn, explore, and grow by orchestrating an intricate symphony of intelligent processes.As the guiding force behind young learners' development, the preschool teacher empowers children to embark on their own lifelong journey of intellectual curiosity, self-discovery, and personal achievement.
Mathematical modeling of these intelligent interactions in preschool education is the fundamental inspiration in PEOA design.

Mathematical model of PEOA
The proposed PEOA approach is a population-based technique that can provide suitable solutions for optimization problems in a repetition-based process based on the search power of its members.The PEOA population is formed by the members of a community so that the position of each of these members in the search space suggests values for the decision variables of the problem.Each population member is a candidate solution for the problem, which can be represented using a vector from a mathematical point of view.The PEOA population consisting of these vectors can be represented using a matrix according to Eq. (1).
where X(t) is the PEOA population matrix, − → X i (t) is the i th PEOA's member, x i,j (t) is the value of the j th vari- able determined by the i th PEOA's member, N is the number of PEOA population members, m is the number of problem variables, t ∈ {1, 2, . . ., T} is the iteration counter (i.e., the number of the actual population) and T is the total number of iterations.At the beginning of the algorithm, the initial position of the PEOA population in the search space is generated randomly using Eq. ( 2).
where r is a random number from the uniform distribution in the interval [0, 1] , lb j and ub j are the lower and upper bound of the j th problem variable respectively.
Since each PEOA member is a candidate solution for the problem variables, the objective function of the problem can be calculated based on the proposed values of each PEOA member.Therefore, the calculated values for the objective function of the problem can be represented using a vector according to Eq. ( 3). (1)

where
− → F (t) is the objective function vector and F i (t) is the objective function value based on the i th PEOA's member.
Based on the comparison of the calculated values for the objective function, the member that provides the best value for the objective function is known as the best population member − → X best (t) .Considering that in each iteration of PEOA, the position of the population members in the search space is updated, new values for the objective function are calculated.Based on the new values evaluated for the objective function, the best member should also be updated in each iteration.
The process of updating the PEOA population in the search space is perform in three phases (i) the gradual growth of the preschool teacher's educational influence, (ii) individual knowledge development guided by the teacher, and (iii) individual increase of knowledge and self-awareness.

Phase 1: The gradual growth of the preschool teacher's educational influence (exploitation phase)
It is evident that the role of the teacher changes significantly with the child's age and thus depending on the school education level.At the beginning of the educational process, in the nursery, i.e., at the child's age from 0 to 2 years, the teacher has primarily an upbringing role, and the educational one is insignificant.In the age of child from 2 to 5 years, i.e., in kindergarten, the educational role gradually increases over the caregiving role, and in kindergarten, i.e., in the age of children from 5 to 6 years, the level of teaching influence of the preschool teacher is almost at the same level as during the following primary education.
In the design of PEOA, the best member is considered as the preschool teacher.Because in preschool education, the teacher's influence increases with the passage of time 69 .To simulate this phase of PEOA, first, based on the teacher's impact, a new position is calculated for each PEOA member using Eq. ( 4).Then, if the value of the objective function is improved in the new position, this new position replaces the previous position of the corresponding member according to Eq. ( 5).where − → X P1 i (t + 1) is the new calculated position for the i th PEOA member based on first phase of PEOA, x P1 i,j (t + 1) is its j th dimension, � K(t) := � X best (t) is the preschool teacher (i.e., the kindergarten teacher), t is the iteration counter, T is the total number of iterations.

Phase 2: Individual knowledge development guided by the teacher (exploration phase)
In this phase of PEOA, population members are updated based on the modeling children's activities, as children try to imitate the work and take on the teacher's experience to be more successful than their classmates.To simulate this phase of the PEOA, first a new position is calculated for each member of the population based on following the preschool teacher using Eq. ( 6).This process leads to large shifts in the position of population members, which has a positive effect on exploration and global search in different areas of the problem-solving space.According to Eq. ( 7), the new position calculated for each member of the population is acceptable if it improves the value of the objective function.Equation ( 7) is a criterion for performing or not performing the process of updating the position of the PEOA member.Hence, Eq. (7) states that the new position is acceptable for a population member if the value of the objective function is improved in the new position, as the movement of population members in the problem-solving space aims to achieve better solutions and prevent the algorithm from moving toward inappropriate solutions.where − → X P2 i (t + 1) is the new calculated position for the i th PEOA member based on second phase of PEOA, x P2 i,j (t + 1) is its j th dimension, − − → rand is a random vector of the dimension m drawn from the uniform distribution . . .

Phase 3: Individual increase of knowledge and self-awareness (exploitation phase)
In addition to the influence of the kindergarten teacher, each child tries to increase their self-awareness through different ways such as playing games, analyzing the possibilities, expectations, etc. Increasing self-awareness leads to achieving an ideal state of themselves.
In the third phase of PEOA, the population members are updated based on modeling children's efforts to raise self-awareness.To simulate this phase of PEOA, first, a new position is randomly generated near each member of the population using Eq. ( 8).This process leads to small changes in the position of population members, which plays an influential role in increasing the PEOA local search and exploitation ability in finding possible better solutions around the discovered solutions.According to Eq. ( 9), the proposed calculated position for each member of the population is acceptable if it improves the value of the objective function.
where − → X P3 i (t + 1) is the new calculated position for the i th PEOA member based on third phase of PEOA, x P3 i,j (t + 1) is its j th dimension and − − → rand is a random vector of the dimension m drawn from the uniform distri- bution in the interval [0, 1].

Repetition process, pseudo-code, and flowchart of PEOA
After updating all PEOA members based on the first to third phases, the first iteration of PEOA is completed.After completing each iteration, the best candidate solution for the problem is updated.Then, based on the new values calculated for the position of the population members and their corresponding objective function, the algorithm enters the next iteration.The process of updating PEOA members using Eqs.(4) to (9) continues until the full implementation of PEOA.At the end, the best candidate solution found during the iterations of the algorithm is presented as the solution to the problem.The pseudo-code of different steps of PEOA implementation is presented in Algorithm 1.

Computational complexity of PEOA
In this subsection, the PEOA computational complexity analysis is discussed.PEOA initialization for an optimization problem has a complexity equal to O(N • m), where N is the number of population members and m is the number of decision variables of the problem.In each iteration, PEOA population members are updated in three phases.The PEOA update process has a complexity equal to O(3N • m • T) , where T is the total number For = 1 to N 7.
Phase 1: The gradual growth of the preschool teacher's educational influence 8.
Calculate the new position for the th PEOA member using Eq.(4).9.
Update position of the th PEOA member using Eq. ( 5).11.
Calculate new position of the th PEOA member using Eq. ( 6). 13.
Update position of the th PEOA member using Eq. ( 7).14.
Phase 3: Individual increase of knowledge and self-awareness (exploitation) 15.
Calculate the new position for the th PEOA member using Eq. ( 8).16.
Update position of the th PEOA member using Eq. ( 9).

Simulation studies and results
In this section, the performance of PEOA in solving optimization problems is evaluated.For this purpose, a set of twenty-three standard benchmark functions of unimodal, high-dimensional multimodal, and fixed-dimensional multimodal types are employed.Full details and explanations of these functions are provided in 70 .In addition, the performance of PEOA in handling the CEC 2017 test suite is also evaluated.Complete information and a detailed description of the CEC 2017 test suite information are available at 71 .The results of PEOA have been compared with the performance of ten famous algorithms GA, PSO, GSA, TLBO, MVO, GWO, WOA, MPA, TSA, and RSA.The rationale behind selecting these ten metaheuristic algorithms from the plethora of options in the literature can be summarized as follows.The first group, encompassing GA and PSO, are well-known and widely used algorithms.The second group consists of GSA, TLBO, GWO, and MVO, the most cited methods.The third group comprises recently published and widely used methods: WOA, MPA, TSA, and RSA.The values of the control parameters of competing algorithms are provided in Table 1.It should also be mentioned that due to the main advantage of the proposed PEOA approach, which lacks control parameters in its mathematical model, it does not need any parameter tuning process.In addressing the twenty-three standard benchmark functions F1 to F23, the PEOA and competing algorithms are each employed in twenty independent runs where each iteration contains 1000 iterations to optimize each benchmark function.In addressing the CEC 2017 test suite, the proposed PEOA approach and competitor algorithms are employed in fifty-one independent implementations where each execution contains function evaluations (FEs).Simulation results are reported using six indicators: mean, best, worst, standard deviation (std), median, and rank.

Qualitative analysis of PEOA
The qualitative analysis results of the proposed PEOA approach in solving some unimodal and multimodal benchmark functions are shown in Fig. 1.In this analysis, four metrics are considered: search history, the Table 1.Assigned values to the control parameters of competitor algorithms.The search history metric shows that PEOA searches the problem-solving space well at both global and local levels to discover the original optimal area and converge to the optimal solution.The trajectory metric shows that PEOA creates large changes in the position of the population members in the initial iterations with high exploration ability.Then, with increasing iterations, based on the exploitation ability with minor changes in the position of the population members, it converges towards solutions close to the global optimum.The average www.nature.com/scientificreports/fitness metric shows that the population of the algorithm moves towards better solutions during the iterations of the algorithm.The convergence curve metric shows that PEOA has a high convergence speed in solving the problem with a descending trend during the iterations of the algorithm, which indicates the high ability of the proposed algorithm to balance exploration and exploitation.

Evaluation of unimodal benchmark functions
In order to analyze the exploitation ability of PEOA and competitor algorithms in local search, seven unimodal functions of F1 to F7 are selected.The optimization results of unimodal functions of F1 to F7, using PEOA and competitor algorithms, are reported in Table 2. Based on the obtained results, PEOA with high exploitation ability has converged to the global optimum in solving functions F1, F2, F3, F4, and F6.In solving functions F5 and

Evaluation of high dimensional multimodal benchmark functions
In order to investigate the exploration capability of PEOA and competitor algorithms in solving problems that have a large number of local optima, six high-dimensional multimodal functions of F8 to F13 have been selected.The implementation results of PEOA and competitor algorithms on functions F8 to F13 are presented in Table 3.
The simulation results show that PEOA, with its high exploration capability, has provided the global optimum in solving F9 and F11 functions by discovering the main optimal area in search space.In solving functions F8, F10, F12, and F13, PEOA is the first best optimizer by providing optimal global search.The analysis of the high-dimensional multimodal simulation results shows that PEOA has an acceptable ability in exploration and global search, and compared to competitor algorithms, it has provided superior efficiency in optimizing F8 to F13 functions.

Evaluation of fixed dimensional multimodal benchmark functions
In order to evaluate the ability of PEOA and competitor algorithms, in creating a balance between exploration and exploitation during the search process, ten fixed-dimensional multimodal functions of F14 to F23 have been selected.The results of using PEOA and competitor algorithms in optimizing functions of F14 to F23 are released in Table 4.The simulation results show that PEOA is the first best optimizer in solving functions F14, F15, F18, F22, and F23.In solving functions F16, F17, F19, F20, and F21, PEOA and some competitor algorithms have provided similar results for the "mean" index.However, PEOA has provided better performance in solving these functions by providing better values in the "std" index.The analysis of the results of fixed-dimensional multimodal functions, shows that PEOA has a high ability to balance exploration and exploitation, and by providing better results for these functions, it has superior performance compared to competitor algorithms.Boxplot diagrams resulting from the performance of PEOA and competitor algorithms in optimizing functions F1 to F23 are presented in Fig. 2.

Evaluation of the CEC 2017 test suite
This subsection evaluates the performance of PEOA and competing algorithms in handling the CEC 2017 test suite.CEC 2017 test suite has thirty standard benchmark objective functions consisting of (i) three unimodal functions of C17-F1 to C17-F3, (ii) seven multimodal functions of C17-F4 to C17-F10, (iii) ten hybrid functions of C17-F11 to C17-F20, and (iv) ten composition functions of C17-F21 to C17-F30.From this test suite, function C17-F2 has been excluded from simulation studies due to its unstable behavior.The results of implementing PEOA and competing algorithms on the CEC 2017 test suite for problem dimensions equal to 10 are reported in Table 5.The boxplot diagrams obtained from the metaheuristic algorithms are drawn in Fig. 3. Based on the optimization results, PEOA is the first best function optimizer: C17-F1, C17-F3 to C17-F21, C17-F23, C17-F24, and C17-F26 to C17-F30.The simulation results show that PEOA has superior performance in handling the CEC 2017 test suite by achieving better results for most of the benchmark functions than competing algorithms.

Statistical analysis
Reporting optimization results using mean, best, worst, standard deviation (std), median, and rank indices provides valuable information about the performance of metaheuristic algorithms.However, even with a very low probability, the superiority of one algorithm over several others may be coincidental.Therefore, in this subsection, a statistical analysis is presented on the performance of PEOA and competing algorithms to show whether the superiority of PEOA from a statistical point of view has a significant difference compared to competing algorithms.For this purpose, the Wilcoxon signed-rank test 72 , a non-parametric statistical test, is used.This test uses an index called " p-value" to determine whether there is a significant difference between the average of two data samples.
The results of employing the Wilcoxon signed-rank test on the performance of PEOA and competitor algorithms in optimizing the mentioned benchmark functions are reported in Table 6.Based on the statistical analysis results, in cases where " p-value" is calculated to be less than 0.05, PEOA has a statistically significant superiority compared to the corresponding competitor algorithm.The statistical analysis shows that PEOA has a significant superiority in solving unimodal benchmark functions, high-dimensional multimodal, fixed-dimensional multimodal, and the CEC 2017 test suite, from a statistical point of view, compared to competitor algorithms.

Discussion
Metaheuristic algorithms are from the group of stochastic approaches to solve optimization problems that can provide suitable solutions for optimization problems based on random search in the problem-solving space in an iterative process.To have an effective search process, metaheuristic algorithms must have the appropriate power in exploitation, exploration, and balancing during the search process.
Unimodal functions do not have any local optima except the global optimum.For this reason, they are suitable options for measuring the ability to exploit metaheuristic algorithms to manage the local search in the problemsolving space to achieve solutions close to (even matching) the global optimum.Benchmark functions F1 to F7, C17-F1, and C17-F3 are selected from the unimodal type.Based on the optimization results, PEOA has provided the global optimum for the functions F1 to F4, F6, C17-F1, and C17-F3, with high ability in exploitation and powerful local search.Also, PEOA is the first best optimizer for F5 and F7 functions.These results confirm and Table 4. Evaluation results of fixed-dimensional multimodal functions.
guarantee the high exploitation power of PEOA to manage local search in problem-solving space.What is evident from the simulation results is that PEOA, by providing better results for unimodal functions and obtaining the rank of the first best generator, has provided a superior performance in competition with competing algorithms, which indicates the exceptional ability of PEOA in exploitation compared to competing algorithms.In addition, the statistical analysis results obtained from the Wilcoxon signed-rank test confirm that PEOA has a significant statistical superiority in competition with the compared algorithms in dealing with unimodal functions and exploitation to manage local search.High-dimensional multimodal functions of F8 to F13 and C17-F4 to C17-F10 have several local optima in addition to the global optimum.For this reason, these functions are suitable options for measuring the quality  Based on the optimization results, PEOA, with high ability in exploration for global search management, has provided the global optimum for functions F9, F11, C17-F4, C17-F6, and C17-F9.In addition, PEOA is the first best optimizer for functions F8, F10, F12, F13, and C17-F4 to C17-F10.What is evident from the optimization results is that PEOA has a high ability in exploration to manage the global search to identify the region containing the global optimum in the problem-solving space.Based on the simulation results, PEOA has provided superior performance compared to competing algorithms by achieving better results for highdimensional multimodal functions and C17-F4 to C17-F10 and getting the rank of the first best optimizer, which shows the superior exploration power of the proposed approach in global search management.In addition, the results of the statistical analysis obtained from the Wilcoxon signed-rank test confirm that this superiority of PEOA is significant from a statistical point of view.Fixed-dimensional multimodal functions F14-F23 have a number of local optima lesser than high-dimensional multimodal functions.For this reason, they are suitable options for evaluating the power of metaheuristic algorithms in balancing exploration and exploitation.Also, benchmark functions C17-F11 to C17-F30 from the CEC 2017 test suite are complex functions that challenge the power of metaheuristic algorithms to balance exploration and exploitation.What is evident from the optimization results is that PEOA, with a high ability to balance between exploration and exploitation, has identified the region containing the global optimum based on the global search and converged to suitable solutions close to the global optimum based on the local search.
Based on the obtained results, PEOA is the first best optimizer in handling F14-F23, C17-F11 to C17-F21, C17-F23, C17-F24, and C17-F27 to C17-F30 compared to competing algorithms.Based on the simulation results, PEOA has provided superior performance compared to competing algorithms to balance exploration and exploitation by achieving better results for most of the benchmark functions and ranking as the first-best optimizer in total.In addition, the statistical analysis results obtained by employing the Wilcoxon signed-rank test confirm that PEOA has a significant statistical superiority compared to competing algorithms to balance exploration and exploitation to deal with high-dimensional multimodal functions of F14 to F23, and functions of C17-F11 to C17-F30 from CEC 2017 test suite.
The main findings from the analysis of the simulation results are that PEOA has high ability in exploitation to manage the local search based on the optimization results of F1 to F7, C17-F1, and C17-F3 functions, has high ability in exploration to manage the global search based on the optimization results of the functions F8 to F13, and C17-F4 to C17-F10, and has a high ability to balance exploration and exploitation based on the optimization results of F14 to F23, and C17-F11 to C17-F30 functions.

PEOA for real-world applications
In this section, the effectiveness of PEOA in solving optimization problems in real world applications has been evaluated.For this purpose, the ability of PEOA and competitor algorithms in optimizing twenty-two real-world optimization problems from CEC 2011 test suite and four engineering design problems has been challenged.In addressing the CEC 2011 test suite, the proposed PEOA approach and each of the competitor algorithms is implemented in twenty-five independent implementations where each implementation contains 150,000 FEs.In addressing the four engineering design problems, the proposed PEOA approach and each of the competitor algorithms is implemented in twenty independent implementations where each implementation contains 1000 iterations.
In this section, the effectiveness of PEOA in solving optimization problems in real-world applications has been evaluated.For this purpose, the ability of PEOA and competitor algorithms to optimize twenty-two realworld optimization problems from the CEC 2011 test suite and four engineering design problems has been www.nature.com/scientificreports/challenged.In addressing the CEC 2011 test suite, the proposed PEOA approach and each competitor algorithm are implemented in twenty-five independent implementations, each containing 150,000 FEs.In addressing the four engineering design problems, the proposed PEOA approach and each of the competitor algorithms are implemented in twenty independent implementations, where each implementation contains 1000 iterations.

Evaluation of the CEC 2011 test suite
This subsection evaluates the efficiency of PEOA and competing algorithms in handling the CEC 2011 test suite.
The CEC 2011 test suite consists of twenty-two constrained optimization problems from real-world applications.Full description, details, and detailed information of the CEC 2011 test suite is available at 73 .The results of employing PEOA and competing algorithms to optimize the CEC 2011 test suite are reported in Table 7.The boxplot diagrams obtained from the performance of the metaheuristic algorithms are plotted in Fig. 4. The optimization results show that PEOA is the first best optimizer for optimization problems C11-F1 to C11-F22.
The simulation results show that PEOA had better results than all competitors for all twenty-two optimization problems and was the first-best optimizer to handle the CEC 2011 test suite.In addition, the results obtained from the statistical analysis of the Wilcoxon signed-rank test confirm that PEOA has significant statistical superiority over the competing algorithms to optimize the CEC 2011 test suite.

Pressure vessel design problem
Pressure vessel design is a real-world application where the main goal in this design is to minimize the construction cost.The schematic of this design is presented in Fig. 5.The mathematical model of pressure vessel design is as follows 74 :

With
The results of pressure vessel design optimization are presented in Table 8.Based on the simulation results, PEOA has presented the optimal design of this problem with the values of the design variables equal to (0.778027, 0.384579, 40.31228, 200) and the objective function equal to 5882.9013.The convergence curve of PEOA in the optimization of pressure vessel design is presented in Fig. 6.Analysis of the simulation results shows that PEOA has provided superior performance in pressure vessel design optimization compared to competitor algorithms.

With
The results of using PEOA and competitor algorithms in optimizing speed reducer design are presented in Table 9.Based on the simulation results, PEOA has presented the optimal design of this problem with the values of the design variables equal to (3.5, 0.7, 17, 7.3, 7.8, 3.3502147, 5.2866832) and the objective function equal to Subjectto :

With
The implementation results of PEOA and competitor algorithms in the optimization of welded beam design are presented in Table 10.Based on the simulation results, PEOA has presented the optimal design of this problem with the values of the design variables equal to (0.20573, 3.470482, 9.036637, 0.20573) and the objective function equal to 1.724856 .The convergence curve of PEOA in the optimization of welded beam design is presented in Fig. 10.What is clear from the analysis of simulation results is that compared to competitor algorithms, PEOA has performed better in optimizing welded beam design.Minimize : f (x) = (x 3 + 2)x 2 x 2 1 .

With
The results of tension/compression spring optimization using PEOA and competitor algorithms are presented in Table 11.Based on the simulation results, PEOA has presented the optimal design of this problem with the values of the design variables equal to (0.051606, 0.354725, 11.40679) and the objective function equal to 0.012665 .The convergence curve of PEOA in tension/compression spring optimization is presented in Fig. 12.Based on the analysis of simulation results, it is concluded that PEOA has provided more effective performance in tension/ compression spring optimization compared to competitor algorithms.

Conclusion and future works
In this paper, a new human-based metaheuristic algorithm called Preschool Education Optimization Algorithm (PEOA) was introduced.The design of PEOA draws its primary inspiration from the dynamic processes within preschool education.The proposed PEOA was explained in three phases (i) the gradual growth of the preschool teacher's educational influence, (ii) individual knowledge development guided by the teacher, and (iii) individual increase of knowledge and self-awareness, and then its mathematical model was presented.A set of fifty-two standard benchmark functions representing unimodal, high-dimensional multimodal, fixed-dimensional multimodal, and CEC 2017 test suite were utilized to assess the optimization capabilities of PEOA.The obtained optimization results highlighted PEOA's proficiency in effectively balancing global exploration and local exploitation during the search.Additionally, PEOA's performance was benchmarked against ten established metaheuristic algorithms.Analyzing the simulation outcomes revealed that PEOA consistently outperforms competing algorithms by delivering enhanced solutions across a majority of benchmark functions.Furthermore, when applied to address twenty-two real-world optimization problems from CEC 2011 test suite and four distinct engineering design challenges, PEOA demonstrated remarkable efficacy in resolving real-world problems, underscoring its efficiency in practical applications.
The proposed PEOA approach has several advantages for global optimization problems.Against these advantages, PEOA also has several disadvantages.The main advantage of PEOA is that its mathematical model has no control parameters, which must be adjusted by the user (of course, exclude the population size N and maximal number of iterations T ).For this reason, the proposed approach does not need the parameter tuning process.The second advantage of PEOA is its high efficiency in dealing with various optimization problems in many sciences and complex high-dimensional problems.The proposed method's third advantage is its excellent ability to balance exploration and exploitation in the search process, which allows high-speed convergence to provide suitable values for decision variables in optimization tasks, especially in complex problems.The fourth advantage of the proposed PEOA is its robust performance in handling real-world optimization applications.The proposed PEOA approach is a stochastic-based solving method.So, the main disadvantage of PEOA, similar to all stochasticbased optimizers, is there is no guarantee that PEOA will achieve the global optimal solution.In addition, PEOA may fail to address some optimization applications because, according to the NFL theorem, there is no a priory presumption that any metaheuristic algorithm will be successful or not.Another disadvantage of PEOA is that it is always possible to develop newer algorithms that perform better than existing algorithms and PEOA.
The introduction of the PEOA approach paves the way for many potential avenues of future investigation and development.One notable prospect is the development of binary and multi-objective adaptations of the PEOA, which holds significant promise for further research.Additionally, exploring the application of PEOA in tackling optimization challenges across diverse disciplines and real-world scenarios offers an enticing direction for future studies.

Figure 2 .
Figure 2. Boxplot of performance of PEOA and competitor algorithms in solving F1 to F23.

Figure 3 .
Figure 3. Boxplot of performance of PEOA and competitor algorithms in solving the CEC 2017 test suite.

Figure 5 .
Figure 5. Schematics of the pressure vessel design.

Figure 6 .
Figure 6.Convergence curves of PEOA on the pressure vessel design.

Figure 7 .
Figure 7. Schematics of the speed reducer design.

Figure 8 .
Figure 8. Convergence curves of PEOA on the speed reducer design.

Figure 9 .
Figure 9. Schematics of the welded beam design.
www.nature.com/scientificreports/Tension/compression spring design problem Tension/compression spring design is a real-world application aimed at minimizing tension/compression spring weight.The schematic of this design is presented in Fig. 11.The mathematical model of tension/compression spring is as follows 23 : Consider : X = [x 1 , x 2 , x 3 ] = [d, D, P].

Figure 10 .
Figure 10.Convergence curves of PEOA on the welded beam design.

Figure 11 .
Figure 11.Schematics of the tension/compression spring design.

Figure 12 .
Figure 12.Convergence curves of PEOA on the tension/compression spring design.

Table 2 .
Evaluation results of unimodal functions.F7, PEOA is the first best optimizer.Analysis of the simulation results shows that PEOA has a high exploitation ability in local search and compared to competitor algorithms, it has provided superior performance.

Table 3 .
Evaluation results of high-dimensional multimodal functions.
Vol.:(0123456789) Scientific Reports | (2023) 13:21472 | https://doi.org/10.1038/s41598-023-48462-1www.nature.com/scientificreports/ of the exploration ability of metaheuristic algorithms to manage the global search in the problem-solving space to prevent the algorithm from getting stuck in local optima and discovering the region containing the global optima.

Table 5 .
Evaluation results of CEC 2017 test suite.

Table 7 .
Evaluation results of CEC 2011 test suite.

Table 8 .
Evaluation results of the pressure vessel design ( SIs are statistical indicators and DVs are design variables).

Table 9 .
Evaluation results of the speed reducer design ( SIs are statistical indicators and DVs are design variables).

Table 10 .
Evaluation results of the welded beam design ( SIs are statistical indicators and DVs are design variables).

Table 11 .
Evaluation results of the tension/compression spring design ( SIs are statistical indicators and DVs are design variables).