A hierarchical hyper-heuristic for the bin packing problem

This paper addresses the two-dimensional irregular bin packing problem, whose main aim is to allocate a given set of irregular pieces to larger rectangular containers (bins), while minimizing the number of bins required to contain all pieces. To solve the problem under study a dynamic hierarchical hyper-heuristic approach is proposed. The main idea of the hyper-heuristics is to search the space of low-level heuristics for solving computationally difficult problems. The proposed approach is “dynamic” since the low-level heuristic to be executed is chosen on the basis of the main characteristics of the instance to be solved. The term “hierarchical” is used to indicate the fact that the main hyper-heuristic can execute either simple heuristics or can run in a “recursive fashion” a hyper-heuristic. The developed solution strategy is evaluated empirically by performing extensive experiments on irregular packing benchmark instances. A comparison with the state-of-the-art approaches is also carried out. The computational results are very encouraging.


Introduction
This paper addresses the two-dimensional irregular bin packing problem (2DIBPP), an NP-hard combinatorial optimization problem (Pandey 2021). The aim is to pack a set of irregular pieces, into homogeneous rectangular stock sheets (bins) of fixed dimensions, in such a way that the number of used bins is minimized, that is the space utilization is maximized.
On the basis of the topology proposed in Wäscher and Schumann (2007), the problem under study falls into the general class of the two-dimensional single bin-size bin packing problem. However, additional features are considered, that is all pieces have irregular convex shape, and they cannot overlap with each other and cannot be flipped.
The 2DIBPP arises in several real-life applications, where it is necessary to cut irregular shape pieces from multiple homogeneous rectangular sheets. It is commonly addressed in the ceramic tile sector, shipbuilding industry, clothing, footwear and furniture production, logistics as well as glass, sheet metal processing and automobile industries (Okano 2002;Martinez-Sykora et al. 2017;Hu et al. 2020;Bennell et al. 2018;Luo et al. 2022).
Despite its practical importance, given the need of dealing with complex geometry, the 2DIBPP has attracted little attention in the scientific literature compared to other bin packing problems, and the majority of the scientific contributions have been focused on heuristic approaches.
In what follows, a concise overview of the literature related to the 2DIBPP is provided. The attention is focused on the scientific contributions most relevant to our study. State-of-the-art Hyper-heuristic methods, based on genetic algorithms, for solving two-dimensional regular (rectangular) and irregular (convex polygonal) have been presented in Terashima-Marín et al. (2010), Lopez-Camacho et al. (2013), López-Camacho et al. (2014).
In particular, in Terashima-Marín et al. (2010) the authors proposed several heuristics for selecting the pieces and the bins, and for placing the pieces into the bins. In addition, they built a benchmark set of 540 instances, used to evaluate the performance of the developed approaches. An extension of the Djang and Finch heuristics, originally designed for the one-dimensional bin packing problem, to the twodimensional case has been proposed in Lopez-Camacho et al. (2013). Very good results have been achieved on convex poly-gon instances. An evolutionary selection and constructive hyper-heuristic approach that combines single heuristics has been proposed by the same authors in López-Camacho et al. (2014).
In Okano (2002), a scanline-based algorithm is proposed to approximate the irregular shape of the pieces and a First-Fit Decreasing (FFD) algorithm is used to address the placement problem.
All the aforementioned papers do not allow pieces rotation. This possibility is instead exploited in Martinez-Sykora et al. (2017), Abeysooriya et al. (2018), Liu et al. (2020) and Zhang et al. (2022). In particular, in Martinez-Sykora et al. (2017) the authors considered free rotations of pieces and presented several variants of a construction algorithm. In order to address both the pieces assignment and the pieces allocation problems, some integer programming models are formulated and solved. In Abeysooriya et al. (2018), it is assumed that both free rotations and limited rotations (i.e., 90 • , 180 • , 270 • , 360 • ) are allowed and, unlike previous contributions, the allocation and placement problems are solved together, by means of a heuristic approach.
The variant of the 2DIBPP with limited rotations has been also studied in Liu et al. (2020). In particular, a constructive solution approach is proposed, in which the pieces are assigned to the bins by using the FFD strategy; the placement problem is solved by the bottom-left algorithm and the pieces exchange method. Furthermore, a greedy local search approach is executed to improve the solution quality.
An extension of this work is given in Zhang et al. (2022), where the authors introduce a waste least first decreasing strategy to assign the pieces to the bins, whereas an overlap minimization approach, based on the separation algorithm presented in Imamichi et al. (2009) andLeung Lin (2012), is used to address the placement problem. To improve the solution, a greedy local search approach that relies on pieces swapping between two bins is proposed. Contribution and organization of the paper This work focuses on the development of a hyper-heuristic approach to solve the 2DIBPP. The hyper-heuristic is used to define a high-level heuristic, which controls low-level heuristics. In particular, it should decide when and where to apply each single low-level heuristic, depending on the characteristic of the given problem state.
There are two main types of hyper-heuristics (Drake et al. 2020), categorized by whether they are used for selecting existing heuristic (base low-level heuristic) or generating new heuristics from components of above heuristic. The approach presented in this paper is of the first type. The aim is to choose base low-level heuristics, which allow to select and place the pieces into the bins in an efficient way.
The main contributions of this paper are the following: -development of a novel hyper-heuristic, able to provide better results than those of the state-of-the-art; -implementation of a benchmark system that exploits information related to the behavior of each low-level heuristic in terms of solution quality, to improve the performance of the hyper-heuristic in a dynamic fashion; -design and implementation of innovative and efficient feature-based selection and placement heuristics; -development and testing of a parallel version of the proposed hyper-heuristic.
The rest of the paper is organized as follows. Section 2 describes the main characteristics of the 2DIBPP addressed in this work. The proposed low-level heuristics are presented in Sect. 3, whereas the developed solution approach is detailed in Sect. 4. Section 5 is devoted to the discussion on the computational results. The paper ends with some concluding remarks given in Sect. 6.

Problem description
In this section, a formal description of the specific variant of the 2DIBPP addressed in this paper is provided. Let P = {p 1 , p 2 , . . . , p n } denote the set of n irregular pieces. Each piece p k , k = 1, . . . , n is represented as a polygon, whose area is denoted as a k . It is assumed that pieces rotations and flipping are not allowed.
Let B = {b 1 , b 2 , . . . , b N } denote the set of N identical rectangular bins, which can be used to pack the n pieces. Each bin b i , i = 1, . . . , N is a rectangular object, characterized by a fixed width W and fixed length L. It is assumed that the number of available bins N is large enough to contain all n pieces.
The objective of the 2DIBPP is to minimize the number of bins used to pack all the n pieces, by ensuring that: -the pieces are placed entirely within the bin; -the pieces do not overlap with each other.
It is worth observing that the number of bins used to pack the n items does not allow to discriminate among solutions characterized by the same number of bins. For this reason, the following two functions are considered to evaluate a solution.
The first one is the fitness function F, introduced in Lopez- Camacho et al. (2013), defined as follows: where U i represents the ratio utilization of each bin b i , i = 1, . . . , N , that is: where n i represents the number of pieces placed in the bin The other evaluation measure is the fractional number of bins K , introduced in Han et al. (2013) and Martinez-Sykora et al. (2017) to take into account the reusable part of a bin and relies on the use of a cut, either vertical or horizontal, which recovers the final non-utilized part of a bin (see Fig. 1). The parameter K is defined as follows: where R * represents the portion of the bin used once the reusable residual area has been removed by applying the vertical (horizontal) cut.
In this paper, the fitness function F is used to guide the proposed hyper-heuristic to explore the solution space, whereas the K is used to evaluate the performance of the proposed solution strategy in comparison with the state-ofthe-art approaches.

Low-level heuristics
The proposed low-level heuristics rely on the combination of "heuristic operators," which are essentially of the following three types: selection, placement and filling operators.
They offer basic functionalities for the exploration of the solution space. Thus, the developed hyper-heuristic is used to find an appropriate combination of these operators, with the aim of building good-quality solutions. The selection operators are used to extract the pieces to be allocated into the containers, by taking into account some specific features associated with the pieces (see Sect. 3.1). More specifically, the next piece to be handled is chosen by considering the presence or absence of these features, according to a certain cutoff value, experimentally determined.
The extracted piece is then placed in the best position by applying a placement operator. Four different operators have been implemented: Bottom Up (PBU), Top Down (PTD), Left Right (PLR) and Right Left (PRL).
In particular, in PBU the pieces are positioned from bottom to top in columns (Bottom Up), in PTD from top to bottom always in columns (Top Down), in PLR from left to right proceeding by rows (Left Right), and finally, in PRL from right to left (Right Left) always in a row. Recursive combinations of these four operators were also considered. Finally, the filling operators are used to fill the empty spaces according to a greedy strategy.
In what follows, the defined heuristic operators are described in details.

Selection operators and features
This section presents the main selection operators that are used to extract the pieces to be placed, by both the placement and the filling operators. It is worth observing that other selection strategies can be derived, by combining the main operators in an appropriate way.
For each piece, the following features have been defined (see Fig. 2

):
HorizontalityDown: it represents the ratio between the width of the base (the lowest horizontal segment of the piece with slope = 0) and the total width of the piece; in the example of Fig. 2, it is equal to 66%. HorizontalityTop: as above, but it takes into consideration the upper segment; VerticalityRight: it indicates the ratio between the vertical right side of the piece (the one with slope = infinite) and the overall height of the piece; in the example of Fig.  2, it is equal to 0%, VerticalityLeft: as above, but for the left side; HorizontalityDelta (VerticalityDelta): it represents the delta between the two extreme points of the figure, and this parameter is very useful for choosing pieces to be put together; Area: Area of the piece; Barycenter: The barycenter of the piece.
The main aim of the proposed operators is to select in a "smart way" the pieces to be placed, by using as criteria the features introduced above.
Since the outcome depends on the placement strategy used, the pieces already inserted in the bin and the remaining area, each selection operator is defined in such a way to work in symbiosis with the applied placement operator. The selection operators behave differently depending on whether If the solution is built by using the PBU placement strategy, the selection operator (referred to as bottom-up selection operator) handles the first pieces to be selected, in descending order of the HorizontalityDown, VerticalityLeft and Area features, while for the subsequent ones, the operator gives priority to the pieces that can be placed side by side, i.e., the pieces for which the difference between the Horizontality-Delta feature is less than a given threshold parameter.
In the case in which the PLR placement algorithm is used, the corresponding selection operator (i.e., left-right selection operator) works in similar way: in the case of the first piece, it gives priority to the same features Horizontality-Down, VerticalityLeft, Area as the bottom-up version. In the case of subsequent pieces, it selects the pieces in descending order of the VerticabilityLeft feature, since the aim is to select pieces that can be placed side by side on the current row.

Filling operators
This section describes the proposed filling operators, whose main aim is to place the pieces within the available space, using several heuristics, which are guided by local fitness functions.
The proposed operators can be viewed as grouped into three main classes (referred to as FA, FB and FC), which differ from each other for the specific fitness function used. Other operators can be derived, by modifying, in an appropriate way, the relevant parameters.
The FA operator inserts pieces into the bin, trying to optimize the position of their barycenter.
In particular, a fitness function is used to determine the coordinate (X, Y) of the barycenter of each piece. The highest priority is assigned to the pieces with the minimum value of the X coordinate, which are then selected first and, in the case of a tie, the pieces with the minimum value of the Y coordinate are processed.
More specifically, the first piece is placed at the bottom left. For the subsequent pieces, the positions tested are all and only those corresponding to the vertices of the area formed by the pieces already placed in the bin. If the vertices are too close to each other (for example when the polygon represents the approximation of a circle), these are discarded, and if they are too sparse, a virtual one is added. It is worth observing that considering only the vertices allows to improve the efficiency of the algorithm, but a deterioration in the solution quality could be observed. The proposed strategy is referred to as virtual vertex search (VVS), since the main idea is to either reduce or virtually increase the vertices to be considered, depending on whether they are too dense or too sparse.
The VVS differs from the no-fit polygon approach, generally used in the scientific literature (Dean et al. 2006;Domovic et al. 2014) since in the latter all possible placements of a polygon are evaluated. Figure 3 reports the configuration obtained by applying the FA operator. More specifically, for each vertex of the area made up of the pieces already placed, the operator tries to place the new one and then check if some constraints (i.e., the piece must be inside the bin and must not overlap with the already inserted pieces) are satisfied; in the affirmative case, it verifies whether the position is improved by using the fitness function described above. With reference to the example given in Fig. 3, the pieces are outside in 2, 3,4,5,6,8,9,10,11,12; the pieces overlap in configuration 7.
While the FA heuristics processes the pieces sorted in decreasing order of the area and checks all the vertices where it is possible to place the pieces, the FB operator instead tries to handle the pieces using a different fitness function. More specifically, it calculates the distance between the vertices of the piece to be placed and those of the polygon corresponding to the pieces already placed (i.e., the occupied area). If the determined distance is less than a predefined cutoff value (i.e., ), the fitness function is increased by 1 distance (i.e., Fitness := Fitness + 1 distance ). The piece for which the fitness function assumes the minimum value is selected.
The last filling operator (i.e., FC) is based on the same fitness function used in FB. Its purpose is not to find a solution for the current bin, but to create clusters of pieces as compact as possible and with a fairly rectangular shape, which will be handled as single larger pieces.

Placement operators
This section gives a detailed description of the considered placement operators (that is, PBU, PTD, PLR, PRL). As mentioned above, the PBU (PTD) places the extracted piece from bottom to up (from top to down) in columns, whereas in PLR (PRL) the piece is placed from left to the right (from right to left) proceeding by rows.
In what follows, the attention will be focused on the PBU and the PLR operators.
For the sake of comprehension, it is useful to consider the example reported in the upper box of Fig. 4, which gives the solution obtained by applying the PBU placement operator The PBU operator works by columns from bottom to up: in the first step it receives, from the selection operator, a piece that has the maximum height among all the pieces that have the greatest value of the HorizontabilityDown and Vertica-bilityLeft features.
This allows to extract the best piece to be placed in the lower left corner (indicated in Fig. 4 with the number 1). Successively, the selection operator is used to individuate a piece with the width as much as possible equal to the one just inserted and with HorizontabilityDown > 0, in order to be sure that the two pieces can be placed side by side (the piece numbered with 2 in the upper box of Fig. 4). The subsequent column is then considered and the approach tries to find the piece (numbered with 3 in Fig. 4), which occupies as much as possible the residual width. Since, after this case, it is not possible to find other pieces to be inserted above the one selected, the placement heuristic stops and as the columns are finished, the hyper-heuristic executes the filling operators, which fill the empty spaces, by selecting the pieces in decreasing order of the area.
The right-left placement operator works in a similar manner, but instead of operating by columns and from bottom to up, it operates by rows from right to left (see the lower box of Fig. 4).
Even though the placement operators are very simple, their combination allows to find different effective solutions and the obtained results are comparable with those of the stateof-the-art. This occurs especially when the area occupied by the pieces is commensurated with the area of the bin (as explained in Sect. 5).

The hyper-heuristic approach
In this section, the hyper-heuristic developed to address the 2DIBPP is described in details.
The proposed solution strategy relies on the execution of the low-level heuristics presented in the previous section. These operators have been defined by considering some basic elements derived from the solution strategies proposed in the scientific literature and are executed and combined in different ways, by the main hyper-heuristic approach. Thus, the main algorithm is not used to build a solution, but instead the aim is to select, in a hyper-heuristic way, these heuristic operators in order to obtain a good-quality solution.
It is worth observing that the basic operators behave differently, depending on the characteristics of the instance to be solved. Thus, in order to fully exploit the capability of the developed operators, since for some instances the mere combinations could not provide good-quality solutions, a hierarchical hyper-heuristic approach has been defined (detailed in Sect. 4.1).
The main idea is to divide the instance into blocks, to which the heuristic operators are applied recursively. Thus a hierarchical approach is obtained: the main hyper-heuristic relies on the execution of lower-level hyper-heuristics, applied on a single block.
In the proposed approach, the heuristic operators always work on only one bin at a time; indeed, it is not allowed the possibility of exchanging pieces between bins.
This choice was made to allow the algorithm to work even in parallel, i.e., the hyper-heuristic can explore in parallel the search space of the basic heuristics, in order to speed up the overall algorithm. This possibility is not considered in the algorithms proposed in the scientific literature, since generally the solution strategy either takes into consideration all the bins or, as a last step, performs a local search, by exchanging pieces between different bins. The flow diagram of the developed approach is depicted in Fig.5, from which it is evident that, until there are pieces to be processed, the following operations are executed: -A new bin is initialized and the low-level heuristics are loaded to be executed; -Several combinations of different low-level heuristics in the current bin (sequentially or in parallel) are evaluated; -The solution, which provides the best value, is selected and the pieces are placed into the current bin, according to the selected configuration.

The hierarchical heuristic scheme
In order to improve the performance of the placement operators, when small pieces need to be handled, a "hierarchical heuristic" has been devised. The main idea is to divide the bin area into blocks and to apply to each of them the heuristic operators, in a similar way to the main hyper-heuristic (Lopez-Camacho et al. 2013). Two different versions of the hierarchical heuristic have been defined, referred to as HHA and HHB.
HHA builds a solution by recursively applying placement operators. Figure 6 provides a graphical representation of the main operations executed by HHA.
The approach starts by considering an empty bin. In the first step (phase 1), the hyper-heuristic decides whether and which low-level heuristic operators to apply, by considering also the dimensions and the residual space of the bin (initially empty).
In the second step (phase 2), the previously selected lowlevel heuristics are executed and, thus, some initial solutions are determined. For each of them, the residual bound box is calculated and, therefore, if in the previous phase at least one piece has been inserted, the placement operators are applied recursively again (return to phase 1).
Finally, when the placement operators are not able to insert other pieces, due to the limited available free space, the filling operators are executed on the current solutions (phase 4) and the best one is selected.
The second variant of the proposed hierarchical heuristic (referred to as HHB) can be viewed as an enhanced version of HHA.
More specifically HHB, similarly to HHA, divides the space into blocks, but in addition, it executes a precomposition procedure, whose main aim is to combine the pieces in a block of fairly rectangular shape. These suboptimal solutions are determined by using the FC operator and considered in the hyper-heuristic as single pieces to be inserted in the bins.

Sequential and parallel versions
In the sequential version of the approach (Fig. 7), the order in which the low-level heuristics are executed plays an important role, since if an optimal solution is obtained for the current bin (i.e., the bin is fully filled), the other operators are not executed, with a consequent reduction in the total execution time.
Several approaches can be used to decide which low-level heuristic should be applied at each iteration. In the present work, the order has been determined experimentally. In par-ticular, to empirically determine the order in which these low-level heuristics should be executed, a benchmark system has been developed that reports, for each bin, the combination of the basic operators that provided the best results. Then, the best performing heuristics are selected to be executed first.
More specifically, a priority queue, initialized on the basis of the specific features (Lopez-Camacho et al. 2013) possessed by the pieces to be placed, was created. The priority value assigned to each low heuristic gives a measure of the quality of the solutions determined so far. This value is In particular, each time a heuristic is executed and a solution is determined, it receives a reward proportional to the fitness value obtained and the corresponding priority value is updated. This priority scheme is embedded into the hyperheuristic, and it is used to select the low-level heuristic to be executed at each iteration.
In the parallel version (Fig. 8), all basic heuristics are performed in parallel. In particular, all types of placement heuristics, all hierarchical hyper-heuristic and all filling operators. Each of these heuristics runs in parallel on a separate thread. If an optimal solution is found (e.g., all bins are filled up to 100%), this solution is returned immediately and the other threads are stopped. In the general case, however, all the low-level heuristics are performed and the most effective solution is selected upon completion.

Computational experiments
This section presents the computational results, collected by testing the proposed solution strategies. The developed algo-rithms have been implemented in the c# language, by using visual studio 2019 as development environment. The experiments have been carried out on a PC with a core i7 processor and 16 GB memory.
The computational experiments have been designed with the aim of assessing the performance of the single operators and their combination in the hyper-heuristic framework. A comparison with the state-of-the-art (i.e., results reported in Lopez-Camacho et al. (2013), Martinez-Sykora et al. (2017), Liu et al. (2020)) is also carried out.
The testings have been conducted on the benchmarks proposed in Terashima-Marín et al. (2010), which are grouped into 18 different sets (A . . . R), each of them containing 30 instances randomly generated. The characteristics of these 540 instances are given in Table 1, where, for each type, the object size, the number of pieces and the optimal solution are reported. It is worth observing that an optimal solution is determined only when all bins are filled up to 100%. From Table 1, it is evident that an optimal solution is available for all the instances, with the only exception of those of type G.
The pieces to be placed are made up of irregular convex polygons, with a number of sides between 3 and 8. The bins are square-shaped. These pieces have been studied in Lopez-  Camacho et al. (2013), and the values of some features (i.e., size, rectangularity, etc.), reported in Table 2, have been calculated . The performance of the proposed approaches has been evaluated by considering the following criteria: -the average number of used bins N ; -the percentage usage of each bin F (see Sect. 2); -the average fractional number of bins K (see Sect. 2); -the execution time T in seconds.
In order to evaluate the gain obtained by combining the developed operators in the proposed hyper-heuristic framework, in a first phase of the computational experiments, the heuristic procedures were tested separately. The related results are given in Tables 3-5. In particular, Table 3 contains the results obtained using the placement strategies (i.e., PBU, PTD, PLR, PRL) combined with the corresponding selection and filling operators; the results provided by the hierarchical heuristics (HHA and HHB) are reported in Table 4, whereas  Table 5 gives the computational results obtained by the filling operators (FA, FB and FC).    The numerical results given in Table 3 clearly underline that the placement strategies allow to obtain the optimal solutions only for the instances of type H . In addition, they showed a good behavior when solving instances of types B, O and Q, for which the value of the criterion F is greater than 0.96. Very poor performances (F ≤ 0.5) are observed on 5 out of 18 types of instances (i.e., D, E, F, L and N ). Thus, an average fitness value of 0.686 is obtained. The placement strategies are very efficient, requiring an average execution time of 11 seconds.
From Table 4, it is evident that the hierarchical heuristics outperform the placement strategies in terms of both solution quality and efficiency. Indeed, the average F increases to 0.720 and the average execution time reduces to about 9 seconds.
The results reported in Table 5, related to the filling operators, highlight that these operators are very efficient, but they are not able to bring satisfactory results in terms of solution quality. Indeed, an average value of F equal to 0.567 is obtained, whereas the average execution time is equal to 3.015 seconds. Thus, these operators allow to obtain reasonable results only if they are used in combination with the other strategies.
The computational results obtained by executing a combination of the operators in the order defined by the proposed hyper-heuristic are reported in Table 6. The results reported in Table 6 underline that the proposed hyper-heuristic outperforms the single operators in terms of solution quality; on the other hand, the combined approach is more time-consuming. Indeed, for instances of type B, H and Q, it is able to find the optimal solutions (100% filling of the bins); in addition, for other four types of instances the value of F is greater than 0.9 (i.e., C, K, O and R). The computational overhead increases to 21.421 seconds on average.
In order to assess the behavior of the implemented operators in terms of solution quality, Table 7 gives the number of times in which each operator has determined the best solution, on the considered instances. From Table 7, it is evident that the best performance is obtained by PBU, which has returned the best solution 1976 times, followed by HH, PRL, PTD and PLR.
In Table 8 the performance of the proposed hyper-heuristic is compared with the state-of-the-art approaches presented in Lopez-Camacho et al. (2013), Martinez-Sykora et al. (2017, Liu et al. (2020)). It is important to point out that the results of Table 8 are obtained with the sequential version of the proposed approach. In addition, in order to ensure a fair comparison in terms of computational overhead, the experiments Table 8 Comparison with the state-of-the-art approaches  have been carried out on a PC with an Intel i5-8250 processor and 8 GB memory.  Liu et al. (2020), where an average fitness of 0.732 has been determined, with an average execution time of 75 seconds.
On average, the proposed hyper-heuristic outperforms the state-of-the-art approaches both in terms of solution quality and efficiency. Indeed, an average fitness equal to 0.811 has been obtained, while keeping the average run time equal to about 38 seconds.
In order to investigate how the performance of the proposed approach is affected by the instance characteristics, it is useful to consider the features of the pieces, described in Lopez-Camacho et al. (2013) and reported in Table 2.
From the collected results, it is evident the implemented solution strategy behaves the best over than 60% of the types of instances.
In the remaining instances, comparable results have been obtained, with a slight deterioration of the performance for the instances for which the average area of the pieces is less than 0.05 and, therefore, the number of pieces is very high (i.e., 60). Table 9 reports the computational times of the sequential version and the parallel counterpart, using 21 threads.
The collected results underline that on average the parallel version is more efficient than the sequential counterpart. In particular, an average reduction of the 50% in the execution time has been observed. However, the advantage of using the parallel algorithm is more evident for those instances, which require relevant computational time. Indeed, in the easiest instances, the overhead due to the parallel algorithm setup hinders the advantages of the parallel approach.

Conclusions
In this paper, the two-dimensional irregular bin packing problem has been addressed. In particular, a hyper-heuristic approach has been developed and tested on benchmark instances. The proposed approach relies on the use of very simple low-level heuristics and exploits a hierarchical strategy, where the hyper-heuristic is executed recursively. The behavior of the proposed approach has been empirically tested, on benchmark instances. The obtained results are very encouraging: on average, the proposed solution strategy outperforms the state-of-the-art approaches in terms of both solution quality and efficiency. The present work has shown that simple low-level heuristics, which are not able to find good-quality solutions, if combined in a hyper-heuristic allow to find solution of good quality, in a limited amount of time.
The proposed approach allows to manage one bin at a time. This choice has reduced the possibility of implementing local search moves that consider the possibility of exchanging pieces between bins and to obtain better solutions. However, on the other hand, it was possible to define a parallel version of the approach, with a consequent reduction in the computational time.
The solution approach proposed in this work allows to address the bin packing problem when convex pieces are considered and the pieces cannot be rotated and flipped. The development of enhanced versions of the proposed hyperheuristic defined in such a way to handle the variants of the problem not considered in this work (i.e., non-convex, rotated and flipped pieces) will be the subject of future research.
Funding Open access funding provided by Universitá della Calabria within the CRUI-CARE Agreement.
Availability of data and material The data are available on request from the authors.

Declarations
Ethical approval No human participants or animals have been involved in this study.

Conflicts of interest/Competing interests
The authors have no conflicts of interest to declare that are relevant to the content of this article.

Informed Consent Not applicable.
Code availability The implemented codes are available on request from the authors.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.