An Interactive Differential Evolution Algorithm Based on Backtracking Strategy Applied in Interior Layout Design

: The Interior layout model is to optimize the arrangement position of each room to maximize the comfort and quality of life of residents. Due to the complexity of the Interior layout problem, the computation of ﬁtness function costs lots of time. To reduce the high computational cost while maintaining the solution performance. An interactive differential evolution algorithm based on Backtracking operator (IDE-BO) is proposed as the solver of the Interior layout model. The human-computer interaction mechanism of IDE beneﬁts the automatic adjustment of ﬁtness parameters that best meet the user’s subjective preferences to achieve the optimal solution. At the same time, the backtracking strategy can also help jump out when the algorithm falls into local optimization. The IDE is compared to other two conventional optimization methods based on two different layout scenarios. The experimental results show that in interior layout model IDE-BO is better than conventional interactive genetic algorithm (IGA) and IDE which do not use BO strategy, the super-performance of IDE-BO in complex situations in terms of execution time and convergence rate


Introduction
With the rapid development of the real estate industry, the demand for house interior design has experienced remarkable growth in the past several decades. Given certain house type and height of storey, indoor space arrangement is the main factor affecting house comfort. Different areas and spaces determine the quality of furniture functions, people's work and rest life is also closely related to a good spatial layout. Thus, it is important to determine the location of the room to maximize the quality of people's life, which is usually referred to as layout optimization problem.
In the traditional interior design, the designer understands the user's needs in advance, draws the floor plan based on a fixed structure, and finally provides the design plan for the user to choose. In the floor plan design process, the user cannot participate in it. Unfortunately, the vast majority of users are lack of professional domain knowledge, they cannot accurately describe their preferences or describe their favorite design styles. They can only use relatively colloquial, vague, and emotional language to describe them, such as "atmosphere", "warm", "western", These abstract terms represent their desired decorating style, but they are unable to give specific requirements or standards. In the end, there will often be situations where the designer spends a lot of energy in designing several sets of solutions, but they do not meet the requirements of the users, and the users cannot tell where they are not satisfied with. The design need to be adjusted in time. This has caused a great waste of social resources. Therefore, in the early stage of house design, especially during the interior floor plan design stage, it is necessary to consider users' individual needs. Encouraging users to timely participate in the interior floor plan design process is the key to solving this problem.
Interior graphic design can be regarded as a spatial layout problem, which is a typical optimization problem. Early literature results can be divided into two major research areas: space scheduling and site layout planning. Thabet et al. [1] developed the SCaRC (Space Constraint and Resource Constraint) scheduling system to quantify the workspace parameters, and the workspace availability during the scheduling process was defined and merged. Heng Li [2] et al. used genetic algorithms to transform the construction sitelevel facility layout problem into a problem of assigning scheduled facilities to scheduled locations. The latest methods for automated interior-space design, such as [3,4], relied on custom mathematical formulations of geometry rules, which are implemented and compiled by using general-purpose programming languages. Merrell mainly studies the baseline of automated interior layout based on design guidelines. The empirical weighted sum of the criteria for each mathematical expression is used [3]. When each rule is implemented in the software process, a large amount of user input is required to define object relations and properties. Other methods are used to generate internal layout, such as have centered around learning relations based on existing interior layouts [4].
In this paper, we propose a method for easily creating indoor space layout only through the user's selection operations. The proposed method uses the cost functions and interactive differential evolution algorithm (IDE), establish a man-machine collaborative work mode. The algorithm is responsible for generating layout map, and the user is responsible for online evaluation. The evaluation information is fed back to the algorithm, and a new generation of drawings are generated again. In this way, it can effectively solve the defects of the previous design model that users are not able to participate in the design process due to their lack of professionalism, reasonably introduce users' subjective preferences, the template that meets the needs is finally obtained, providing blueprints for users and designers. Section 1 introduces the background and significance of the research in the field of interior layout, analyzes the existing design defects in the field of interior layout, illustrates the urgent need to improve the existing design methods in the context of "high demand" and "low efficiency", and finally proposes Finally, a new design approach combining interior layout and interactive differential evolution algorithm is proposed. Section 2 presents a detailed introduction to the layout field, illustrates the current state of contemporary research in the layout field, and points out the shortcomings of modern design software. Section 3 proposes the interactive evolutionary algorithm, introduces the interactive genetic algorithm and the interactive differential evolutionary algorithm, focuses on the relevant contents of the interactive differential evolutionary algorithm, and proposes a backtracking strategy to help the algorithm get rid of the local optimum dilemma in view of the shortcoming that the traditional heuristic algorithm is easy to fall into the local optimum. Section 4 introduces the application model of the algorithm in indoor layout. According to the demand conditions of population in different household stages, the corresponding design framework is given, and the constraint function as well as the model of genetic coding are designed according to the relevant rules, the algorithm design process is introduced in detail, and the whole operation process from inputting parameters to generating the final target solution is completed through the user interaction interface. Section 5 tests the convergence effect of the algorithm and the satisfaction of users using this system through simulation experiments and subjectivity experiments, and it can be seen that it is indeed feasible to apply the interactive differential evolution algorithm in the field of indoor layout.

Layout Problem
The problem of layout design is the research on the reasonable arrangement and placement of specific objects. There are a large number of layout design problems in modern industry and engineering, such as container placement and shipment in marine terminal [5], equipment and material flow in warehouse [6] and workshop scheduling [7], architectural layout design in construction industry, etc. Layout design has strong demand and huge market space. The research on its layout design has extensive and profound practical significance. Most layout design problems are inseparable from space, such as the division of one-dimensional space, the indoor plan design of two-dimensional space, and the cockpit layout of three-dimensional space [8,9]. Therefore, the space layout problem accounts for the main part of the layout design.
Interior spatial layout is not only the main problem involved in interior design, but also an important research direction in layout design. Compared with other types of layout problems, interior spatial layout has more research value and prospect interior design [10], virtual reality [11], computer games [12] and other fields have been widely used in indoor space layout in the field of architecture, interior space layout includes layout graphic design (space division) and scene optimization layout (home decoration). Traditional interior design often lacks personalized needs due to the differences in communication between designers and users, resulting in the consumption of a lot of human resources.
In the early research, a large number of restrictive research methods were proposed for the interior graphic design of residential buildings [13,14]. Shami and Mirahmadi [15] designed an evolutionary system to generate multiple groups of configurations according to the given constraints, and rank their priorities according to the evaluation indicators to find the best configuration in the indoor layout. A semi-automatic modeling system was proposed by Rosser [16] to integrate building plan and real-time map data to generate building model. In addition, the expert system (Expert System) [17] implementation plan layout has also been the focus of scholars. Rio-Cidoncha [18] integrated the idea of artificial intelligence and expert system to design a kind of layout design model. Ocheol and Kwon [19] encapsulated the architectural design knowledge and relevant design specifications into the expert system, referenced and deduced the spatial layout based on the given BIM spatial information and existing knowledge, and compared the spatial layout schemes with the expert system to determine the best choice under different accessibility conditions.
Unfortunately, due to the complexity of the combinatorial problem, the traditional constraint methods and evolutionary algorithms can not fully match the constraints and evaluation indicators in the layout problem, such as room adjacent comfort, different indoor styles, natural factors, Feng Shui layout and other factors that can not be formulated [20][21][22]. In the user evaluation mode in IDE, the user's subjective evaluation replaces the complex data evaluation, so that ordinary people who lack professional knowledge can also replace the designer's perspective. Compared with others, users can understand their own needs, simplify complexity and improve efficiency.

Auxiliary Design Tools
Although the style of interior design will change with the fashion trend or social atmosphere, the work done by designers in the design layout has not changed much. They collect materials and on-site investigation according to the requirements of the owners, and design the corresponding sketch, and negotiate with party A to analyze and rectify the good sketch scheme, and finally make the final effect drawing and construction drawing. In recent years, with the development and popularization of computer field, more and more design software is gradually applied to people's daily life and work. The famous computer design software AutoCAD, 3dsmax and Archicad can be well applied to interior design and greatly reduce the workload. However, the above aided design systems provide great convenience to the design industry, but they are only one-way operation tools, they are lack of communication platform between users and designers, and the degree of user participation is limited. Moreover, the use of these tools requires professional design basis and software use knowledge, which is difficult for non-professionals to master, so it has certain limitations in the use process.

Interactive Evolutionary Algorithm
As a traditional multi-objective optimization problem, interior spatial layout can be solved by evolutionary algorithm. However, the user demand is subjective and difficult to be described by explicit mathematical functions, which is a typical implicit optimization problem. As a representative method to solve implicit optimization problems, interactive evolutionary algorithm does not solve many non-inferior solutions first, but works out the final solution step by step through the dialogue between analysts and decision-makers, which has attracted the attention of many scholars.

Interactive Genetic Algorithm
Early IEA applications were mostly implemented using interactive genetic algorithms [23][24][25], individual evaluation scores are given by users as a bridge to interact with the algorithm [26]. The iterative process of the algorithm is shown in Figure 1. Evolutionary individuals are provided by the interactive interface in explicit form for users to score, using a purely numerical or hierarchical system. The disadvantage of this evaluation method is that the scores given by users can only be relative results which need to be mapped to absolute fitness values, increasing the complexity of the algorithm. However, the above aided design systems provide great convenience to the design industry, but they are only one-way operation tools, they are lack of communication platform between users and designers, and the degree of user participation is limited. Moreover, the use of these tools requires professional design basis and software use knowledge, which is difficult for non-professionals to master, so it has certain limitations in the use process.

Interactive Evolutionary Algorithm
As a traditional multi-objective optimization problem, interior spatial layout can be solved by evolutionary algorithm. However, the user demand is subjective and difficult to be described by explicit mathematical functions, which is a typical implicit optimization problem. As a representative method to solve implicit optimization problems, interactive evolutionary algorithm does not solve many non-inferior solutions first, but works out the final solution step by step through the dialogue between analysts and decision-makers, which has attracted the attention of many scholars.

Interactive Genetic Algorithm
Early IEA applications were mostly implemented using interactive genetic algorithms [23][24][25], individual evaluation scores are given by users as a bridge to interact with the algorithm [26]. The iterative process of the algorithm is shown in Figure 1. Evolutionary individuals are provided by the interactive interface in explicit form for users to score, using a purely numerical or hierarchical system. The disadvantage of this evaluation method is that the scores given by users can only be relative results which need to be mapped to absolute fitness values, increasing the complexity of the algorithm.

Interactive Differential Evolution Algorithm
Although the introduction of "human-computer interaction mode" can produce better solutions, it also produces fatigue related problems, with the increase of the number

Interactive Differential Evolution Algorithm
Although the introduction of "human-computer interaction mode" can produce better solutions, it also produces fatigue related problems, with the increase of the number of iterations. The user's fatigue will gradually increase, and the given evaluation score will also be distorted. Therefore, it is difficult to find a balance between the number of iterations and fatigue. Compared with the genetic algorithm, the interactive differential evolution algorithm replaces the scoring mechanism with the selection operation. The user only needs to select a more satisfactory individual between two individuals, which simplifies a series of comparison scoring operations in the mind, which reduces the difficulty of the user's operation and helps to reduce the degree of fatigue in the operation process. The iterative process of the algorithm is shown in Figure 2. In recent years, interactive differential evolution algorithm has been applied in image retrieval [27], image enhancement [28], image filter [29] and other aspects.

Selection Operation
After the test vector is obtained, the selection operator is used to select a better one from the parent vector and the trial vector to make it survive to the next generation. The selection operator is expressed by the following formula:

Interactive Differential Evolution Algorithm Based on Backtracking Strategy
As a kind of heuristic algorithm, differential evolution algorithm is a greedy strategy in essence, which objectively determines that the optimal solution that does not conform to the greedy rule will miss, and the local optimal solution will inevitably fall into.
Generally, the greedy selection operator is excellent. However, when the individual is in a locally optimal state, if the greedy selection operator continues to be used, a better vector from the trial and parent vectors will be selected to survive to the next generation. Since the parent vector is a locally optimal value, the newly generated trial vector is likely to continue further in the locally optimal direction. When an individual falls into a local

Initialization
At the start of the algorithm, a specific method is used to generate the initial population, the method is expressed using the following formula:

Mutation Operation
During each iteration, a mutation vector is generated by the mutation operation. Mutation operation plays an significant role in the performance of differential evolution algorithm. Therefore, since the differential evolution algorithm was first proposed, a great deal of researchers have conducted research on mutation operators and proposed a large number of mutation operators. The most classic mutation operators are as follows.
(1) DE/rand/1: (2) DE/best/1: (3) DE/current-to-best/1: where X j,t refers to the jth individual in the population at the tth iteration, V j,t is the corresponding mutation vector, r1, r2, r3 are randomly selected from {1, 2, . . . , N} and are not equal to j. X best,t is the best individual at tth iteration, X pbest,t is chosen at random from the top p% individuals, F j is the scaling factor of the jth individual, meaning that the scale factors of each individual are independent of each other. At the same time, the adaptive strategy can also be applied to F j .
The three randomly selected individuals X r1,t , X r2,t and X pbest,t in the mutation operator are sorted from good to bad, corresponding to fitness f b , f w and f m . Where F l = 0.1 and F u = 0.9.

Crossover Operation
Then cross operation is carried out on the mutation vector and the parent vector to get the trial vector. The crossover operator is generally represented by the following formula: where U j,i,t is the ith component of the trial vector of the jth individual at tth iteration. CR is the crossover rate of jth individual, i rand is a random number from {1, 2, . . . , D}.

Selection Operation
After the test vector is obtained, the selection operator is used to select a better one from the parent vector and the trial vector to make it survive to the next generation. The selection operator is expressed by the following formula:

Interactive Differential Evolution Algorithm Based on Backtracking Strategy
As a kind of heuristic algorithm, differential evolution algorithm is a greedy strategy in essence, which objectively determines that the optimal solution that does not conform to the greedy rule will miss, and the local optimal solution will inevitably fall into.
Generally, the greedy selection operator is excellent. However, when the individual is in a locally optimal state, if the greedy selection operator continues to be used, a better vector from the trial and parent vectors will be selected to survive to the next generation. Since the parent vector is a locally optimal value, the newly generated trial vector is likely to continue further in the locally optimal direction. When an individual falls into a local optimal value, it is characterized by stagnation and is difficult to evolve to a better solution [30]. Literature [31] also proposed that greedy selection is difficult to help individuals to jump out of the local optimal value, so it is necessary to develop a new selection operator to help individuals jump out of the local optimal value.
To solve the above problems, an interactive differential evolution algorithm based on backtracking strategy is proposed, which helps the algorithm to return to a previous state by gradually returning to the upper layer, so as to jump out of the local optimal state. This method has been applied to differential evolution algorithm [32] and particle swarm optimization algorithm [33] respectively, and good experimental results have been obtained in high-dimensional problems.
As the greedy selection operator performs well when individuals do not fall into local optimum, greedy selection operator is used when individuals do not fall into local optimum, a new selection operator based on backtracking strategy is used when the individual falls into the local optimum. The first step is to determine whether the individual is trapped in a local optimum. An individual is considered to be in a stagnant state when its data is slowly updated within a certain range [34]. Although this method may have a subtle failure rate, it is an effective way to determine whether a target individual is stuck in a local optimum, and a strategy can be developed to reduce the impact of misclassification. The new selection operator borrows the idea of the backtracking algorithm, as shown in Figure 3. A spatial warehouse is set up at each iteration to store the individuals Rx (x denotes the number of layers of the current iteration) that failed to compete and were discarded in this round of comparison. When an individual is detected to be trapped in the local optimum at iteration y, vectors R y−1 , R y−2 , and R y−3 are successively removed from the storage space and one of them is randomly selected to survive to the next generation, and if the individual still does not jump out of the local optimum, vectors R y−4 , R y−5 , and R y−6 are successively removed from the storage space and removed from [R y−1 , R y−2 , R y−3 , R y−4 , R y−5 , R y−6 ] to choose a random vector to survive to the next generation, and so on. If the parent vector is replaced with a randomly generated vector, it is difficult to guarantee the adaptation value of the random vector and the adaptation value must be calculated, which increases the computational complexity. Therefore, replacing the parent vectors with discarded trial vectors is considered. Since these discarded trial vectors are not passed on to the next generation, individuals evolving based on these vectors are more likely to help individuals jump out of local optima, and Figure 4 shows the complete iterative process of the algorithm. where √ represents the better individual who wins in each selection.

Generation Method of Indoor Plan
The interior layout of houses can be solved as a multi-objective optimization problem [35]. Therefore, the relative importance of different objective functions should be considered. In the traditional multi-objective weighting method, each objective function is multiplied by a customized weight and integrated as a single objective, which is transformed into a single-objective optimization algorithm to obtain the optimal solution [36]. In interior design, multi-objective optimization is a big problem that perplexes decision makers. Some factors, such as orientation, daylighting, adjacency and feng shui, are generally difficult to set in advance. A set of Pareto optimal solutions without preference can be obtained by using genetic algorithm, and decision-makers can choose according to their own preferences. However, it will undoubtedly intensify decision-makers' sense of fatigue and greatly reduce the effect to select the appropriate results from the numerous feasible solutions, and it is difficult to determine the adaptive value with specific functions in design and art industries. The "human-computer interaction" mode in this paper can replace the

Generation Method of Indoor Plan
The interior layout of houses can be solved as a multi-objective optimization problem [35]. Therefore, the relative importance of different objective functions should be considered. In the traditional multi-objective weighting method, each objective function is multiplied by a customized weight and integrated as a single objective, which is transformed into a single-objective optimization algorithm to obtain the optimal solution [36]. In interior design, multi-objective optimization is a big problem that perplexes decision makers. Some factors, such as orientation, daylighting, adjacency and feng shui, are generally difficult to set in advance. A set of Pareto optimal solutions without preference can be obtained by using genetic algorithm, and decision-makers can choose according to their own preferences. However, it will undoubtedly intensify decision-makers' sense of fatigue and greatly reduce the effect to select the appropriate results from the numerous feasible solutions, and it is difficult to determine the adaptive value with specific functions in design and art industries. The "human-computer interaction" mode in this paper can replace the complex adaptive value function with the subjective choice of the user in a more complex situation, then evolve to the individual that the user is satisfied with, and gradually search to obtain the solution that meets the requirements.

Exterior Frame
Interior layout refers to the overall framework of the house and the internal space distribution, usually based on the external overall framework fixed in the case of interior design. At different stages of family development, the composition of family members and living functions are characterized by stages [37]. According to the research results, the number of family members and space demand of current home buyers are shown in Table 1. According to the structural design of the house, the house type map presents a variety of forms, which can be divided into one-bedroom, two-bedroom, three-bedroom, multibedroom, etc. Table 2 shows a schematic representation of the different types of frameworks.

Coding Mode
There are many encoding forms of chromosomes in differential evolution algorithm [38]. In this paper, binary coding is adopted, and the parameters of distribution space are expressed in the form of chromosome bit string based on character set {0,1}. In this paper, taking block diagram ④ as an example, the plan is divided into seven parts: room β1-β6 and information of distribution location of each room, the chromosomes of the population are composed of the coding arrangement of all variables (as shown in Figure 5), A-F correspond to 6 rooms respectively, and G represents the position information of their arrangement. The process of gene corresponding to individual is explained in detail ( Figure  6). When generating and arranging rooms, the proportion needs to meet certain constraints, but this relationship is by no means specific, and can only be determined according to the design style, room type and user requirements. The result of this algorithm is only the preliminary blueprint of the design, which provides the possibility for users and designers to discuss. Therefore, there may be small errors between the parameters and the actual situation.

Coding Mode
There are many encoding forms of chromosomes in differential evolution algorithm [38]. In this paper, binary coding is adopted, and the parameters of distribution space are expressed in the form of chromosome bit string based on character set {0,1}. In this paper, taking block diagram ④ as an example, the plan is divided into seven parts: room β1-β6 and information of distribution location of each room, the chromosomes of the population are composed of the coding arrangement of all variables (as shown in Figure 5), A-F correspond to 6 rooms respectively, and G represents the position information of their arrangement. The process of gene corresponding to individual is explained in detail ( Figure  6). When generating and arranging rooms, the proportion needs to meet certain constraints, but this relationship is by no means specific, and can only be determined according to the design style, room type and user requirements. The result of this algorithm is only the preliminary blueprint of the design, which provides the possibility for users and designers to discuss. Therefore, there may be small errors between the parameters and the actual situation.

Coding Mode
There are many encoding forms of chromosomes in differential evolution algorithm [38]. In this paper, binary coding is adopted, and the parameters of distribution space are expressed in the form of chromosome bit string based on character set {0,1}. In this paper, taking block diagram ④ as an example, the plan is divided into seven parts: room β1-β6 and information of distribution location of each room, the chromosomes of the population are composed of the coding arrangement of all variables (as shown in Figure 5), A-F correspond to 6 rooms respectively, and G represents the position information of their arrangement. The process of gene corresponding to individual is explained in detail ( Figure  6). When generating and arranging rooms, the proportion needs to meet certain constraints, but this relationship is by no means specific, and can only be determined according to the design style, room type and user requirements. The result of this algorithm is only the preliminary blueprint of the design, which provides the possibility for users and designers to discuss. Therefore, there may be small errors between the parameters and the actual situation.

Coding Mode
There are many encoding forms of chromosomes in differential evolution algorithm [38]. In this paper, binary coding is adopted, and the parameters of distribution space are expressed in the form of chromosome bit string based on character set {0,1}. In this paper, taking block diagram ④ as an example, the plan is divided into seven parts: room β1-β6 and information of distribution location of each room, the chromosomes of the population are composed of the coding arrangement of all variables (as shown in Figure 5), A-F correspond to 6 rooms respectively, and G represents the position information of their arrangement. The process of gene corresponding to individual is explained in detail ( Figure  6). When generating and arranging rooms, the proportion needs to meet certain constraints, but this relationship is by no means specific, and can only be determined according to the design style, room type and user requirements. The result of this algorithm is only the preliminary blueprint of the design, which provides the possibility for users and designers to discuss. Therefore, there may be small errors between the parameters and the actual situation. One-bedroom Two-bedroom Three-bedroom Multi-bedroom

Coding Mode
There are many encoding forms of chromosomes in differential evolution algorithm [38]. In this paper, binary coding is adopted, and the parameters of distribution space are expressed in the form of chromosome bit string based on character set {0,1}. In this paper, taking block diagram 4 as an example, the plan is divided into seven parts: room β1-β6 and information of distribution location of each room, the chromosomes of the population are composed of the coding arrangement of all variables (as shown in Figure 5), A-F correspond to 6 rooms respectively, and G represents the position information of their arrangement. The process of gene corresponding to individual is explained in detail ( Figure 6). When generating and arranging rooms, the proportion needs to meet certain constraints, but this relationship is by no means specific, and can only be determined according to the design style, room type and user requirements. The result of this algorithm is only the preliminary blueprint of the design, which provides the possibility for users and designers to discuss. Therefore, there may be small errors between the parameters and the actual situation. (a) Type Ⅰ ( b) Type Ⅱ ( c) Type Ⅲ Figure 6. Detailed explanation of coding.

Algorithm Flow
The specific idea of this interior layout design is as follows: according to the number (a) Type Ⅰ ( b) Type Ⅱ ( c) Type Ⅲ

Algorithm Flow
The specific idea of this interior layout design is as follows: according to the number of users' housing needs, choose the external frame of the house with suitable space size, according to the given external frame of the house, choose the regular space area to allocate the functional space required by the users, the shape of the functional space is chosen

Algorithm Flow
The specific idea of this interior layout design is as follows: according to the number of users' housing needs, choose the external frame of the house with suitable space size, according to the given external frame of the house, choose the regular space area to allocate the functional space required by the users, the shape of the functional space is chosen to represent the most common rectangle in the residential building in life, and the irregular polygon composed of multiple rectangular spaces constitutes the basic interior space layout room unit, the horizontal and vertical distances of functional spaces as well as the location and neighboring relationships are disassembled into the genetic algorithm population individual genetic display, and specific spatial rule constraints are organized according to the survey of a large number of layout maps in the market to prevent a large number of invalid solutions, which affect the efficiency of the algorithm.
The specific process of the algorithm is described as shown in Algorithm 1.

Algorithm 1 The algorithm flow of IDE with the proposed selection operator
(1) Initialization: At the start of the algorithm, a specific method is used to generate the initial population, the method is expressed using the following formula: During each iteration, a mutation vector is generated by the mutation operation. The variation criteria are as follows: where X j,t refers to the jth individual in the population at the tth iteration, V j,t is the corresponding mutation vector, r1, r2 are randomly selected from {1, 2, . . . , N} and are not equal to X pbest,t is chosen at random from the top p% individuals, F j is the scaling factor of the jth individual. (3) Crossover operation: Then cross operation is carried out on the mutation vector and the parent vector to get the trial vector. The crossover operator is generally represented by the following formula: where U j,i,t is the ith component of the trial vector of the jth individual at tth iteration. CR is the crossover rate of jth individual. i rand is a random number from {1, 2, . . . , D}.

(4) Selection operation:
After the trial vector is obtained, the selection operator is used to select a better one from the parent vector and the trial vector to make it survive to the next generation. The selection operator is expressed by the following formula: In this paper, the selection operation is handed over to the user through the interactive interface. The user can choose one of two operations on the interface, and put the unselected individuals into a specific storage space. If f U j,t > f X j,t , make X j,t equal to and place the unselected into the storage space. In addition, the competing failed vectors are put into storage for use in the backtracking strategy when trapped in a local optimum (5) After evolution, a new round of population K + 1 was obtained: 1 If an individual is found to be trapped in a local optimum after detection, the "backtracking" operation is adopted to take out the abandoned individuals of the previous three generations from the storage space, and randomly select one of them to survive to the next generation. Finally, return to step (5). 2 If the user has found the existence of an individual satisfying his preference, the algorithm operation will be terminated. Otherwise, go to step (2).
(6) Backtracking strategy to solve local optimum dilemma: When the algorithm is detected to be stuck in a local optimum for a long time, the predecessor variables are continuously removed from the storage space.
A vector is randomly extracted from the set of extracted vectors and assigned to the individual in the current generation that is caught in the local optimum.
When the individual is still in the local optimum after iteration, continue the above operation by taking out the vector from the previous 6 generations and assigning the value, and so on.

Interactive Interface
The user interaction interface is shown in Figure 7. According to the principle of simplicity and ease of use the interface can be roughly divided into three parts, the left part is the dominant display of individual genes of the population, from which the user can observe the individuals generated by each iteration of evolution and score them to give adaptation values, the upper half of the right part is the parameter setting area, where the user can adjust parameters such as algorithm crossover and variation probability by himself before the algorithm starts, and select the suitable furniture style, the lower half of the right half of the right side is the individual genetic details, you can select a layout drawing in the left population, and the detailed parameters of the drawing's room will be displayed on the interface. When the experiment starts, first enter the interface, the user selects the item "home style" in the parameter setting section, such as "Simplicity", "European", etc. Afterwards, the user can modify the values of the scale factor (default 0.65) and cross ratio (default 0.2) according to their wishes, and click the "Start" button. Wait until the initial population is generated on the left interface, then the user selects the individual, checks the detailed parameters of the individual and clicks the "Next" button. The algorithm evolves itself to generate the next generation of populations according to the user's choice, and so on, until we get a layout map that matches the user's subjective preference.
pean", etc. Afterwards, the user can modify the values of the scale factor (default 0.65) and cross ratio (default 0.2) according to their wishes, and click the "Start" button. Wait until the initial population is generated on the left interface, then the user selects the individual, checks the detailed parameters of the individual and clicks the "Next" button. The algorithm evolves itself to generate the next generation of populations according to the user's choice, and so on, until we get a layout map that matches the user's subjective preference.

Test in Simulated Environment
In order to quickly verify the feasibility of the algorithm, it is assumed that the user's final satisfactory solution is the plan view shown in Figure 8, that is O(o 1 , o 2 , · · · o D ). The distance between the evolutionary individual and the satisfactory solution is used as the fitness function instead of user evaluation. It is assumed that the user has obtained a satisfactory solution (as shown in Figure 8), the optimization objective has been determined, by the distance between the candidate solution and the target solution as the individual adaptation value. This represents to some extent the gap between the candidate solution and the user's expectation. The simulation experiments are mainly designed to verify the feasibility of user evaluation in the IDE algorithm framework and that the algorithm converges properly.

Test in Simulated Environment
In order to quickly verify the feasibility of the algorithm, it is assumed that the user's final satisfactory solution is the plan view shown in Figure 8, that is O( , , ⋯ ). The distance between the evolutionary individual and the satisfactory solution is used as the fitness function instead of user evaluation. It is assumed that the user has obtained a satisfactory solution (as shown in Figure 8), the optimization objective has been determined, by the distance between the candidate solution and the target solution as the individual adaptation value. This represents to some extent the gap between the candidate solution and the user's expectation. The simulation experiments are mainly designed to verify the feasibility of user evaluation in the IDE algorithm framework and that the algorithm converges properly. The euclidean distance between evolutionary individual X( , , ⋯ , ) and target feature vector can help users participate in simulation evaluation. The distance calculation formula is expressed as: The euclidean distance between evolutionary individual X(x 1 , x 2 , · · · , x n ) and target feature vector can help users participate in simulation evaluation. The distance calculation formula is expressed as:

Comparison of Mutation Operators
An excellent mutation operator has quite a significant impact on the performance of the algorithm. In this experiment, four classical mutation operators are compared in order to select the operator that best fits the interior layout model (Figure 9). Meanwhile, the traditional differential evolutionary algorithm uses the adaptation value function to calculate the population individual scores, and selects the top p% of individuals from the highest to the lowest according to the score ranking, from which individuals are randomly selected as the random vector in the variation operator DE/current-to-pbest/1, while the adaptation value function is replaced by the human selection operation in IDE, and the traditional adaptation value ranking is no longer applicable, in order to minimize the error, time of choice is used as the index of individual quality ranking, when the individuals of the population are compared two by two, the individuals that meet the user satisfaction and have higher scores are more likely to be selected by the users in the shortest time, although the time index may misjudge a few good target individuals in the high quality population, the error caused by such misjudgment is almost negligible in the case of having a large number of high quality target individuals. satisfactory solution as the evaluation criteria. The faster the distance decreases or the smaller the gap, the better the effect of the strategy in practical application. The experiment was divided into four control groups with different difference strategies. In order to make the experimental data more accurate, the experiment was conducted in 5 times, and the average value was finally selected as the data of the experiment. From the comparison results listed in Figure 9, we can observe that DE/current-to-pbest/1 and DE/current-tobest/1 attain the best performance since they show better results in convergence. At the same time, compared with DE/current-to-best/1, DE/current-to-pbest/1 is better in convergence speed and convergence, followed by DE/rand/1 and DE/best/1.

Experiment of Parameter Sensitivity
After comparative experiments, it can be found that DE/current-to-pbest/1 is more suitable for the model in this paper than the other three mutation operators. Therefore, in the following experiments, we compare the effects of different values of mutation parameter p on the experimental results, from the two aspects of convergence speed and convergence, the experiment was conducted in 5 times, and the average value was taken as the data of the experiment at the end.
As can be seen from Figure 10, in terms of convergence speed, when p ∈ 0,20 there will be a better result. However, if p is too small, the mutation operator will degenerate into DE/current-to-best/1 and lose convergence. From the two aspects, we can know that when p is 15, the algorithm not only has good convergence, but also has fast convergence speed. We set a final satisfactory solution (Figure 8), and set the corresponding fitness function, taking the Euclidean distance between the experimental individual and the final satisfactory solution as the evaluation criteria. The faster the distance decreases or the smaller the gap, the better the effect of the strategy in practical application. The experiment was divided into four control groups with different difference strategies. In order to make the experimental data more accurate, the experiment was conducted in 5 times, and the average value was finally selected as the data of the experiment. From the comparison results listed in Figure 9, we can observe that DE/current-to-pbest/1 and DE/current-to-best/1 attain the best performance since they show better results in convergence. At the same time, compared with DE/current-to-best/1, DE/current-to-pbest/1 is better in convergence speed and convergence, followed by DE/rand/1 and DE/best/1.

Experiment of Parameter Sensitivity
After comparative experiments, it can be found that DE/current-to-pbest/1 is more suitable for the model in this paper than the other three mutation operators. Therefore, in the following experiments, we compare the effects of different values of mutation parameter p on the experimental results, from the two aspects of convergence speed and convergence, the experiment was conducted in 5 times, and the average value was taken as the data of the experiment at the end.
As can be seen from Figure 10, in terms of convergence speed, when p∈ [0, 20] there will be a better result. However, if p is too small, the mutation operator will degenerate into DE/current-to-best/1 and lose convergence. From the two aspects, we can know that when p is 15, the algorithm not only has good convergence, but also has fast convergence speed.

Comparison of Algorithms
In IGA algorithm, in order to solve the requirement of non-negative probability in roulette selection, the calculation formula of adaptive value is revised as follows: where = √10. Figure 11 shows the average convergence curve variation of the three algorithms in finding the target value in ten simulation experiments. When t = 20, there is no significant difference between the three algorithms in the initial stage of iteration. In the later stage of iteration, IGA converges in advance, and IDE and IDE-BO have better effects. When t = 60, the convergence speed and convergence accuracy of IDE-BO are obviously better than the other two algorithms. In order to obtain accurate comparison results, the number of iterations in the simulation experiment is t = 150. In the real human-computer interaction environment, in order to avoid user fatigue, t ≤ 30 is generally appropriate.

Comparison of Algorithms
In IGA algorithm, in order to solve the requirement of non-negative probability in roulette selection, the calculation formula of adaptive value is revised as follows: (20) where c max = √ 10. Figure 11 shows the average convergence curve variation of the three algorithms in finding the target value in ten simulation experiments. When t = 20, there is no significant difference between the three algorithms in the initial stage of iteration. In the later stage of iteration, IGA converges in advance, and IDE and IDE-BO have better effects. When t = 60, the convergence speed and convergence accuracy of IDE-BO are obviously better than the other two algorithms. In order to obtain accurate comparison results, the number of iterations in the simulation experiment is t = 150. In the real human-computer interaction environment, in order to avoid user fatigue, t ≤ 30 is generally appropriate. Table 3 shows the number of iterations required to reach convergence for IDE-BO and the other two traditional algorithms in 10 simulation experiments, and it can be seen that IDE-BO has reached the convergence condition most of the time at around t = 125, while the number of iterations required to reach around 133 in IDE. Although IGA can converge significantly faster than IDE-BO, its convergence is not good as can be seen in Figure 11. To statistically compare the performance difference in convergence speed between IDE-BO and its competitors, a Wilcoxon signed-ranks test with a 0.05 significance level is conducted. Table 4 shows the p-values of the results when comparing IDE-BO with the other two algorithms. From the table, it can be seen that the p-value of IDE-BO is much lower than 0.05 for IDE which has a small difference in convergence and has a much better convergence than IGA which has a faster convergence. In general, the performance difference between IDE-BO and the other two algorithms is significant.
where = √10. Figure 11 shows the average convergence curve variation of the three algorithms in finding the target value in ten simulation experiments. When t = 20, there is no significant difference between the three algorithms in the initial stage of iteration. In the later stage of iteration, IGA converges in advance, and IDE and IDE-BO have better effects. When t = 60, the convergence speed and convergence accuracy of IDE-BO are obviously better than the other two algorithms. In order to obtain accurate comparison results, the number of iterations in the simulation experiment is t = 150. In the real human-computer interaction environment, in order to avoid user fatigue, t ≤ 30 is generally appropriate. Figure 11. Result of the simulation evaluation. selected the individuals in the population that best met their preferences and scored them on a scale from 0 to 100, and the average score given by the 25 testers in each iteration was calculated. The higher the score of the two indexes, the stronger the convergence ability of the algorithm. When the score of the final solution reaches 75 or more, the system can be considered to generate the final satisfactory solution for the users. Figure 12 shows how the satisfaction level varies with the user's subjective choice and the algorithm. judgment, and then affect the convergence of the algorithm. In this test, compared with "Simplicity", the word "European" has more complex meaning and thinking, which is difficult for users to judge. Therefore, the "Simplicity" style shows better results.
To statistically compare the performance difference between IDE-BO and traditional evolutionary algorithms, a Wilcoxon signed-ranks test with a 0.05 significance level is conducted. Table 5 shows the p-values of the results when comparing IDE-BO with the other two algorithms at different numbers of iterations. The p-values below 0.05 are shown in bold. According to the number of iterations, the experiments were divided into three stages: pre, intermediate, and post. From Table 4, it can be seen that there is little difference in the performance of the three algorithms at the pre stage, at the intermediate stage, comparing IGA, IDE-BO has a significant advantage, and at the post stage, comparing both algorithms, the p-value is less than 0.05. Overall, the difference between GODE and the other two algorithms is remarkable.
(a) Simplicity style (b) European style Figure 12. Subjective test. The p-values below 0.05 are shown in bold.

Summary and Future Work
In this paper, the interactive differential evolution algorithm is used to generate the interior layout plan, and the backtracking strategy is introduced to reduce the influence when the population falls into local optimum. The interior design mentioned is a local design adjustment (user's psychological expectation layout) under the satisfaction of certain constraints (domain expert's reasonableness assessment). The interior design discussed in this paper is not a question of the rationality of the interior layout, but whether It can be seen from the Figure 12 that the three algorithms have no obvious effect on the result when the iteration number t ∈ [0,12]. When t > 14, IDE will show stronger convergence than IGA, and the improved IDE is easier to jump out after falling into local optimization, so it will have higher satisfaction in the end. At the same time, we know that different subjective needs of users will have an important impact on the convergence of the algorithm, and different style choices will also have an impact on users' subjective judgment, and then affect the convergence of the algorithm. In this test, compared with "Simplicity", the word "European" has more complex meaning and thinking, which is difficult for users to judge. Therefore, the "Simplicity" style shows better results.
To statistically compare the performance difference between IDE-BO and traditional evolutionary algorithms, a Wilcoxon signed-ranks test with a 0.05 significance level is conducted. Table 5 Table 4, it can be seen that there is little difference in the performance of the three algorithms at the pre stage, at the intermediate stage, comparing IGA, IDE-BO has a significant advantage, and at the post stage, comparing both algorithms, the p-value is less than 0.05. Overall, the difference between IDE-BO and the other two algorithms is remarkable. The p-values below 0.05 are shown in bold.

Summary and Future Work
In this paper, the interactive differential evolution algorithm is used to generate the interior layout plan, and the backtracking strategy is introduced to reduce the influence when the population falls into local optimum. The interior design mentioned is a local design adjustment (user's psychological expectation layout) under the satisfaction of certain constraints (domain expert's reasonableness assessment). The interior design discussed in this paper is not a question of the rationality of the interior layout, but whether the interior layout design meets the user's expectations. This adjustment of house interior design varies from person to person, and the psychological expectations of different users vary. The interactive algorithm and its related strategies are introduced to realize the individual pursuit of different users, thus effectively bridging the gap between the "professionalism" of design software and the "unprofessionalism" of ordinary users. However, the number of iterations and population range of the algorithm should not be too large due to user interaction during algorithm execution, which limits the breadth of search space. Otherwise, it will aggravate user fatigue and aesthetic fatigue, thus affecting the final results. Therefore, if a screening mechanism is added in subsequent studies to select representative solutions from each generation and present them to users, user evaluation fatigue will be alleviated. How to maintain a balance between slowing down user fatigue and expanding search space can be further discussed in the follow-up research.