Article Search
닫기

## Original Article

Split Viewer

Int. J. Fuzzy Log. Intell. Syst. 2018; 18(2): 135-145

Published online June 25, 2018

https://doi.org/10.5391/IJFIS.2018.18.2.135

© The Korean Institute of Intelligent Systems

## Genetic Algorithm-based Optimal Investment Scheduling for Public Rental Housing Projects in South Korea

Jae Ho Park, Jung Suk Yu, and Zong Woo Geem

1Graduate School of Urban Planning & Real Estate Studies, Dankook University, Yongin, Korea, 2School of Urban Planning & Real Estate Studies, Dankook University, Yongin, Korea, 3College of Information Technology, Gachon University, Seongnam, Korea

Correspondence to :
Zong Woo Geem (zwgeem@gmail.com)

Received: April 7, 2018; Revised: June 9, 2018; Accepted: June 18, 2018

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Declining birthrate is a serious problem that threatens the sustainability of Korean society. The main cause of this phenomenon is high living cost where housing cost accounts for the majority in household expenditure. South Korea has a very low supply rate in public rental housing when compared to other OECD countries. Because young people cannot afford to buy or lease a house for their new houses, some of them postpone or even give up marriage. As a countermeasure, Gyeonggi Province (surrounding area of Seoul) recently announced the supplying plan of 10,000 public rental houses by 2020. We expect this measure to alleviate the low birthrate problem and increase the demographic sustainability of the province. This study optimizes multi-annual investment scheduling for rental housing projects using genetic algorithm while satisfying the constraints such as budget, human resources, regional balance, etc. Through the optimal investment scheduling, we hope that public corporation will supply public rental houses more efficiently and more sustainably for the community.

Keywords: Public rental house, Optimal investment scheduling, Sustainable housing, Genetic algorithm

Low birthrate is a serious problem that threatens the sustainability of Korean society. The OECD data shows the total fertility rate of Korea is 1.21 (children per woman) and it is the lowest among developed countries (United States 1.86, United Kingdom 1.81, France 1.98, Japan 1.42, Germany 1.47) in 2014 as shown in Figure 1. Assuming no migration and unchanged mortality, a total fertility rate of 2.1 ensures a stable population.

Two main reasons for the declining birthrate in Korea are unstable employment and high housing cost (purchasing or rental). Recent studies by Oh and Choi [2] and Kim and Hwang [3] have shown that higher housing cost is an important factor for fertility decline and postponement of having children. According to Oh and Choi [1], young people are more influenced by rental cost rather than house-purchasing price. Thus, it is important to supply public rental houses for newlyweds. The survey in 2015 by Korean Ministry of Health and Welfare shows that the housing is the number one factor related to marriage, as seen in Figure 2.

Price to income ratio (PIR) in house has been widely used, both as an indicator of house price level and as that of affordability to acquire residential property, in different geographical areas or time periods. According to Lee et al. [4], PIR in Korea (4.4) is higher than that in the United States (3.5) and Canada (3.4), but lower than those in United Kingdom (5.2), Australia (6.1) and Hong Kong (11.4). Moreover, Korea’s long-term public rental housing ratio is considerably lower than that of major developed countries as observed in Table 1.

Since there are differences in conditions and social environments, it may not be adequate to simply compare the percentage of long-term public rental housing ratio. However, it is an undeniable fact that long-term public rental houses are very deficient in South Korea.

As a countermeasure to this problem, the government of Gyeonggi Province, metropolitan area containing Seoul, announced the plan to supply 10,000 public rental houses, called Dda Bok House (DBH; here, ‘Dda Bok’ means warm-hearted welfare), by 2020 on May 17, 2016.

DBH is a unique housing program that combines public rental housing and rental subsidy policies to significantly lower the burden of housing cost. In this program, the rental subsidy is proportional to the number of the babies. While the area of typical public rental house (two-room type) supplied to the newlyweds couple is 36 m2, the area of DBH is 44 m2, which is 22% wider than the current two-room type residence. Table 2 shows the comparison in details.

DBH plans to provide a safe nursing environment as well as autonomous community by creating various community-based programs. Thus, talent donators, who can give childcare, medical service, cooking, arts activity, physical activity, and local community movement, have a priority to live in the community.

The purpose of this study is to find the optimal investment scheduling for the multi-annual DBH project by maximizing the utility while meeting all the constraints.

So far, many researchers have performed useful research in the field of land-project optimization. Fei Xie et al. [6] introduced a multistage chance-constrained stochastic model for strategic planning of battery electric vehicle (BEV) inter-city fast charging infrastructure. The model was applied to a case study in California and solved by genetic algorithm. Liu et al. [7] performed land-use spatial optimization by coordinating the competitions between different land-use types in the example of Gaoqiao Town, Zhejiang Province, China. Li and Parrott [8], Haque and Asami [9], Taromi et al. [10] and Zhang and Huang [11] reached a consensus that evolutionary algorithm-based model is powerful tool in optimizing land use allocation in the cases of the Delaware, USA and Shenzhen, China, etc.

Other researchers have mainly been interested in private housing demand forecasting [12] and real estate appraisal forecasting with ridge regression [13]. Overall, these studies highlighted the suitability of evolutionary algorithm in the research regarding the land planning and housing demand forecasting. Unlike these studies, Park et al. [14] proposed an optimization model for deriving the optimal investment combination from the various real estate development projects. With these tools, public corporation could have reasonable investment portfolio that maximize the return of investment while meeting the multiple constraints.

In this study, an optimization model is further developed and applied to optimal investment scheduling for public rental housing projects with different groundbreaking years in South Korea.

Gyeonggi Urban Innovation Corporation (GICO), the subsidiary of Gyeonggi Province, is in charge of providing DBH effectively with the above-mentioned missions while considering the constraints such as human resources, annual budget, regional balance, etc., as follows:

• 1) Because DBH should be supplied in many regions simultaneously, GICO needs a lot of supervisors for construction and personnel for getting permission from local municipality. However, if GICO hires many employees at initial stage, it is a heavy burden to maintain those people after DBH project is finished. Therefore, it is desirable to maintain the proper number of personnel.

• 2) In order to maintain the financial soundness, GICO should make annually balanced investment rather than initial excessive investment as the DBH project is non-for-profit business.

• 3) In order to maintain regional balance, Gyeonggi Province is divided into 5 sub-regions such as Seoul-Busan (SB) area in orange color, West coastal (WC) area in sky-blue color, Seoul-Sinuiju (SS) area in blue color, Seoul-Wonsan (SW) area in grey color, and Eastern (ET) area in green color, as shown in Figure 3. Any sub-region cannot be excessively developed.

• 4) DBHs are built for newlywed, college students, junior employees and senior citizens. Thus, it is important to fairly

Genetic algorithm (GA) was first proposed by John Holland in the 1960s and was further developed by himself and his students and colleagues at the University of Michigan in the 1960s and the 1970s. Holland’s 1975 book “Adaptation in Natural and Artificial Systems” presented the GA as an abstraction of biological evolution and gave a theoretical framework for adaptation under the GA. GA is a method for moving from one population of “chromosomes” (e.g., strings of ones and zeros, or “bits”) to a new population by using a kind of “natural selection” together with the genetics-inspired operators of crossover, mutation, and inversion. Each chromosome consists of “genes” (e.g., bits), each gene being an instance of a particular “allele” (e.g., 0 or 1). The selection operator chooses those chromosomes in the population that will be allowed to reproduce, and on average the fitter chromosomes produce more offspring than the less fit ones. Crossover exchanges subparts of two chromosomes, roughly mimicking biological recombination between two single-chromosome (“haploid”) organisms; mutation randomly changes the allele values of some locations in the chromosome; and inversion reverses the order of a contiguous section of the chromosome, thus rearranging the order in which genes are arrayed.

Given a clearly defined problem to be solved and a bit string representation for candidate solutions, a simple GA works as follows:

• 1) Start with a randomly generated population of n l-bit chromosomes (candidate solutions to a problem).

• 2) Calculate the fitness f(x) of each chromosome x in the population.

• 3) Repeat the following steps until n offspring have been created:

• Select a pair of parent chromosomes from the current population, the probability of selection being an increasing function of fitness. Selection is done “with replacement,” meaning that the same chromosome can be selected more than once to become a parent.

• With probability Pc (the “crossover probability” or “crossover rate”), crossover the pair at a randomly chosen point (chosen with uniform probability) to form two offspring. If no crossover takes place, form two offspring that are exactly copies of their respective parents.

• Mutate the two offspring at each locus with probability Pm (the “mutation probability or “mutation rate”), and place the resulting chromosomes in the new population. If n is odd, one new population member can be discarded at random.

• 4) Replace the current population with the new population.

• 5) Go to step 2).

For this study, GA utilized parameter values: population size is 100, and mutation rate is 0.075.

For GA computation, we used Solver software from Frontline Systems, Inc. Average computational time ranges from 15 to 20 seconds per each run and termination condition is when convergence rate is equal to or less than 0.0001.

### 4. Optimization Formulation and Computational Results

To easily understand the optimal investment scheduling, 10 sample projects were selected as depicted in Table 3. Each project has different total project cost, resident type, land value and region. We want to find the optimal investment schedule by maximizing the effect of DBH supply while satisfying various constraints.

We think that the effect or utility of DBH can be maximized when more valuable houses are supplied earlier than less valuable ones. Although the housing value is not easily quantified, we instead use the appraised value of the project land. Therefore, the objective function of this scheduling problem is Max ∑ HiPiW(Yi), where Hi is the number of houses in project i; Pi is the appraised value of land in project i; and W(Yi) is weighting factor of groundbreaking year. To find the optimal investment schedule with the maximized objective function value, we set up weighting factor system. Because the earlier the better, if the groundbreaking year is 1, the weighting factor becomes 3 (highest value); if the groundbreaking year is 2, the weighting factor becomes 2 (medium value); and if the groundbreaking year is 3, the weighting factor becomes 1 (lowest value).

Figure 4 shows the structure of chromosome incorporating decision variables used in GA computation. The value of each decision variable is groundbreaking year (1st, 2nd or 3rd year).

### 4.1 Scenario 1

Here, we can consider several different scenarios. In Scenario 1, we set groundbreaking year (1st, 2nd, or 3rd year) of project i as decision variable. If Yi is decision variable, it has one of three candidate values, Yi ∈ {1, 2, 3}.

As previously mentioned, the objective function is Max ∑i HiPiW(Yi), where W(Yi) is weighting factor which becomes W(Yi) = 3 if Yi = 1; W(Yi) = 2 if Yi = 2; and W(Yi) = 1 if Yi = 3.

For the human resources limitation, we limit maximum four projects to launch in year j, N(Yi = j) ≤ 4, j = 1, 2, 3, where N(●) is number function which counts the total number of launched projects in year j.

For the project budget limitation, the maximum annual investment amount is 600 (unit= 100 million KRW), iCij600, j = 1, 2, 3, 4. Here, Cij is partial investment cost occurred in year j for project i, and 600 is about 36% of total 1,647 (100 million KRW).

For regional balance, any specific area cannot be excessively developed. In this example, among five areas (ET, SB, WC, SW, and SS), the maximum number of annually launched projects for SB area is limited up to two (the total number of SB projects in this example is four), N(Yi = j & Ai = SB) ≤ 2, j = 1, 2, 3. Here, Ai is sub-region where project i is included.

For resident-type balance, any specific class cannot be excessively developed. In this example, among three classes (junior employee, newlywed, and industrial worker), the maximum number of annually launched projects for junior employee is limited up to two (the total number of junior employee projects in this example is six), N(Yi = j & Ri =junior employee)≤ 2, j = 1, 2, 3. Here, Ri is resident type of project i.

The optimized results for this investment schedule example are shown in Table 4.

As observed in Table 4, projects B, C, I, and J will launch in 1st year; projects D, E, and G will launch in 2nd year; and projects A, F, and H will launch in 3rd year.

For the human resources limitation, this constraint is satisfied because the number of launched projects in 1st year is 4 (B, C, I, J), that in 2nd year is 3 (D, E, G), and that in 3rd year is 3 (A, F, H).

For the project budget limitation, this constraint is satisfied because the investment amount in 1st year is 286, that in 2nd year is 583, that in 3rd year is 531, and that in fourth year is 247. Every of them is less than 600.

For regional balance, this constraint is satisfied because the number of SB projects in 1st year is 2 (C, J), that in 2nd year is 2 (D, E), and that in 3rd year is 0. Every of them is less than two.

For resident-type balance, this constraint is satisfied because the number of junior employee projects in 1st year is 2 (C, J), that in 2nd year is 2 (D, E), and that in 3rd year is 2 (A, H). Every of them is less than two.

### 4.2 Scenario 2

In Scenario 2, two constraints are loosened from Scenario 1. While the maximum number of launched projects in each year was 4 in Scenario 1, that becomes 5 in Scenario 2. And, while the maximum number of junior employee projects in each year was 2 in Scenario 1, that becomes 4 in Scenario 2. Table 5 shows the detailed results from Scenario 2.

As observed in Table 5, projects B, C, H, and J will launch in 1st year; projects D, E, G, and I will launch in 2nd year; and projects A and F will launch in 3rd year. The differences between Scenarios 1 and 2 are projects H and I. The investment schedule of project H is moved forward from 3rd year to 1st year, and that of project I is moved back from 1st year to 2nd year for obtaining better effect of DBH project. The objective function value increased from 4,084 to 4,107.

For the human resources limitation, this constraint is satisfied because the number of launched projects in 1st year is 4 (B, C, H, J), that in 2nd year is 4 (D, E, G, I), and that in 3rd year is 2 (A, F). Every of them is less than 5.

For the project budget limitation, this constraint is satisfied because the investment amount in 1st year is 297, that in 2nd year is 583, that in 3rd year is 535, and that in 4th year is 232. Every of them is less than 600.

For regional balance, this constraint is satisfied because the number of SB projects in 1st year is 2 (C, J), that in 2nd year is 2 (D, E), and that in 3rd year is 0. Every of them is less than 2.

For resident-type balance, this constraint is satisfied because the number of junior employee projects in 1st year is 3 (C, H, J), that in 2nd year is 2 (D, E), and that in 3rd year is 1 (A). Every of them is less than 4.

### 4.3 Scenario 3

In Scenario 3, the maximum annual budget is tightened from Scenario 2. We apply maximum budget of 500 (100 million KRW) instead of 600. Table 6 shows the detailed results from Scenario 3.

As observed in Table 6, the schedules of two projects (G and I) have been changed from Scenario 2. The investment schedule of project G is moved back to 3rd year from 2nd year, and that of project I is moved forward to 1st year from 2nd year. Also, the objective function value decreased from 4,107 to 4,028.

For the human resources limitation, this constraint is satisfied because the number of launched projects in 1st year is 5 (B, C, H, I, J), that in 2nd year is 2 (D, E), and that in 3rd year is 3 (A, F, G). Every of them is less than 5.

For the project budget limitation, this constraint is satisfied because the investment amount in 1st year is 312, that in 2nd year is 488, that in 3rd year is 425, and that in 4th year is 362. Every of them is less than 500.

For regional balance, this constraint is satisfied because the number of SB projects in 1st year is 2 (C, J), that in 2nd year is 2 (D, E), and that in 3rd year is 0. Every of them is less than 2. For resident-type balance, this constraint is satisfied because the number of junior employee projects in 1st year is 3 (C, H, J), that in 2nd year is 2 (D, E), and that in 3rd year is 1 (A). Every of them is less than 4.

### 4.4 Scenario 4 (Real-World Example)

Now we consider a real-world investment schedule problem, which has 31 projects. Table 7 shows the details of the projects. Because the construction costs of some projects marked with asterisk are not determined yet, those costs are assumed by multiplying 130 million KRW and number of houses. Here, resident type consists of junior employee (JE), newlywed (NW), industrial worker (IW), senior citizen (SC), and to-be-determined (TBD).

Actually, 11 projects (A, B, D, H, J, K, L, M, N, O, and P) have already launch-ed in 2017 by GICO. Thus, the decision variables of those projects should have a value of one, which means 1st year, and those of other 20 projects should have a value of two or three, Yi ∈ {2, 3}. Here, the number of candidate investment schedule solutions becomes 220 (≥ 1, 000, 000) because 20 projects have two choices. As already utilized, the objective function is Max ∑i HiPiW(Yi).

For the human resources limitation, we limit maximum 11 projects to launch in year j, N(Yi = j) ≤ 11, j = 2, 3.

For the project budget limitation, the maximum annual investment amount is 2,000 (108 KRW), iCij2,000, j = 2, 3, 4. Here, 2,000 is around 1/3 of total project cost 6,102.

For regional balance, the maximum number of annually launched projects for SB area is limited up to 3, N(Yi = j & Ai = SB) ≤ 3, j = 2, 3.

For resident-type balance, the maximum number of annually launched projects for junior employee is limited up to 2, N(Yi = j & Ri = JE) ≤ 2, j = 2, 3.

The optimized results for this real-world investment schedule problem are shown in Table 8.

As observed in Table 8, 11 projects (C, F, Q, U, V, W, X, Z, AA, BB, CC) will launch in 2nd year; and 9 projects (E, G, I, R, S, T, Y, DD, EE) will launch in 3rd year.

For the human resources limitation, this constraint is satisfied because the number of launched projects in 1st year is 11, and that in 2nd year is 9. Every of them is less than 11.

For the project budget limitation, this constraint is satisfied because the investment amount in 2nd year is 1,886, that in 3rd year is 1,981, and that in 4th year is 1,149. Every of them is less than 2,000.

For regional balance, this constraint is satisfied because the number of SB projects in 2nd year is 3 (C, Q, W), and that in 3rd year is 3 (E, R, EE). Every of them is less than 3.

For resident-type balance, this constraint is satisfied because the number of JE projects in 2nd year is 2 (C, X), and that in 3rd year is 2 (E, R). Every of them is less than 2.

Public rental house is essential for low-income classes who cannot afford to purchase a house or pay a high security deposit (some property owners require 80% of house prices). South Korea has lowest fertility rate among OECD countries and the main reason is cost burden of house purchase or rent. Furthermore, public rental housing ratio in South Korea is considerably lower than that of major OECD countries. Thus, it is urgent to expand the supply of public rental houses for the housing sustainability of South Korea.

This study proposed a GA-based optimization model for multi-annual investment scheduling of public rental house project, obtaining good results. This cannot be easily obtained if we use only “engineering intuition”. The optimization model becomes very useful and practical tool for the massive DBH project while meeting various constraints such as limited human resources, limited annual budget, regional supply balance, and resident-type balance. We hope this optimization approach will be used for the investment scheduling of similar public projects in the future.

The proposed scheduling optimization model may be further applied to various cases such as financial investment, advisory service, and facility repair, as follows:

• - Financial investment: Financial investor may use the scheduling optimization model with various investment projects, which have different maturity, rate of return, and risk with limited financial resources to obtain multi-annual desirable cash flow.

• - Advisory service: Consulting firm, accounting firm, or law firm may use the scheduling optimization model with limited human resources to carry out more consulting projects simultaneously.

• - Facility repair: City officials may use the scheduling optimization model for lots of infrastructure facilities such as road, water supply and sewer pipe with limited human and financial resources to maintain facilities.

This research was supported by the Ministry of Science and ICT, Korea, under the Information Technology Research Center support program (No. IITP-2018-2017-0-01630) supervised by the Institute for Information & communications Technology Promotion (IITP).

Fig. 1.

Fertility rates in 2014. From OECD data [1].

Fig. 2.

Desirable measures to support marriage.

Fig. 3.

Major locations of DBH project in Gyeonggi Province. supply with respect to class balance.

Fig. 4.

Structure of chromosome incorporating decision variables.

Table. 1.

Table 1. Comparison of long-term public rental house ratio in various countries.

CountryKoreaUKFranceJapanOECD averageEU average
Long-term public rental house ratio5.5%17%19%6.1%11.5%9.4%
Base Year201420102007200820072007

From Ministry of Land, Infrastructure and Transport, 2012 [5]..

Table. 2.

Table 2. Comparison between current rental house and DBH.

Current two-room houseDBH
Area, exclusive (m2)3644
Floor plan

Table. 3.

Table 3. Samples of DBH projects.

ProjectProject districtTotal project cost (1st, 2nd, 3rd year) (KRW 100 million)Number of houses (H)Resident type (R)Appraised value of land (P) (KRW/m2)Region (A)
AGapyeong-eup190 (73+117+0)48Junior employee368,600ET
BNamyangju-Changhyeon61 (26+35+0)48Newlywed860,500ET
CSeongnam-Pangyo390 (200+190+0)300Junior employee2,897,000SB
DSuwon-Mangpo113 (38+75+0)100Junior employee827,000SB
ESuwon-Seodun130 (60+70+0)100Junior employee830,100SB
FSiheung-Sincheon85 (40+45+0)65Industrial worker409,200WC
GSiheung-Jeongwang380 (120+200+60)290Newlywed/Industrial worker399,300WC
HYangpyeong-Gongheung51 (26+25+0)49Junior employee610,800ET
IOsan-Gajang65 (15+40+10)50Industrial worker733,400WC
JYongin-Jukjeon182 (45+100+37)150Junior employee1,346,000SB

Table. 4.

Table 4. Results of optimal investment scheduling for Scenario 1.

Project nameResident typeAppraised value of land (KRW/m2)Start yearWeighted valueNumber of houses (H)Total project cost (KRW 100 million)
1st year2nd year3rd year4th year
AJunior employee368,6003148--73117
BNewlywed860,50013482635--
CJunior employee2,897,00013300200190--
DJunior employee827,00022100-3875-
EJunior employee830,10022100-6070-
FIndustrial worker409,2003165--4045
GNewlywed/Industrial worker399,30022290-12020060
HJunior employee610,8003149--2625
IIndustrial worker733,4001350154010-
JJunior employee1,346,000131504510037-
Total1,200286583531247

Table. 5.

Table 5. Results of optimal investment scheduling for Scenario 2.

ProjectResident typeAppraised value of land (KRW/m2)Start yearWeighted valueNumber of houses (H)Total project cost (KRW 100 million)
1st year2nd year3rd year4th year
AJunior employee368,6003148--73117
BNewlywed860,50013482635--
CJunior employee2,897,00013300200190--
DJunior employee827,00022100-3875-
EJunior employee830,10022100-6070-
FIndustrial worker409,2003165--4045
GNewlywed/Industrial worker399,30022290-12020060
HJunior employee610,80013492625
IIndustrial worker733,4002250154010
JJunior employee1,346,000131504510037-
Total1,200297583535232

Table. 6.

Table 6. Results of optimal investment scheduling for Scenario 3.

ProjectResident typeAppraised value of land (KRW/m2)Start yearWeighted valueNumber of houses (H)Total project cost (KRW 100 million)
1st year2nd year3rd year4th year
AJunior employee368,6003148--73117
BNewlywed860,50013482635--
CJunior employee2,897,00013300200190--
DJunior employee827,00022100-3875-
EJunior employee830,10022100-6070-
FIndustrial worker409,2003165--4045
GNewlywed/Industrial worker399,30031290-120200
HJunior employee610,80013492625
IIndustrial worker733,4001350154010
JJunior employee1,346,000131504510037-
Total1,200312488425362

Table. 7.

Table 7. Total DBH projects.

ProjectProject districtProject cost (1st, 2nd, 3rd year) (108 KRW)Number of housesResident typeLand value (KRW/m2)Region
AGapyeong-eup190 (73+117+0)48JE368,600ET
BNamyangju-Changhyeon61 (26+35+0)48NW860,500ET
C*Seongnam-Pangyo390 (200+190+0)300JE2,897,000SB
DSuwon-Mangpo113 (38+75+0)100JE827,000SB
ESuwon- Seodun130 (60+70+0)100JE830,100SB
F*Siheung- Sincheon85 (40+45+0)65IW409,200WC
G*Siheung-Jeongwang380 (120+200+60)290NW/IW399,300WC
HYangpyeong-Gongheung51 (26+25+0)49JE610,800ET
I*Osan-Gajang65 (15+40+10)50IW733,400WC
JYongin-Jukjeon182 (45+100+37)150JE1,346,000SB
KSuwon-Gwanggyo1393 (252+141+0)204NW2,211,000SB
LAnyang-Gwanyang171 (71+0+0)56NW2,428,000SB
MHwaseong-Jinan27 (27+0+0)31JE674,300WC
NSuwon-Woncheon196 (78+118+0)152SC2,145,000SB
OPaju-Byungwon134 (67+67+0)50SC966,000SS
PUiwang-Bugok56 (13+43+0)50JE1,151,000SB
Q*Anyang-Gwanyang2390 (200+190+0)300NW1,554,000SB
R*Yongin-Yeosung130 (70+60+0)100JE1,229,000SB
S*Icheon-Gwango130 (70+60+0)100SC819,600ET
T*Gwangju-Yeokdong650 (230+220+200)500NW318,500ET
U*Ansan-Wonsi299 (100+100+99)230IW240,000WC
V*Hanam-Deokpung2169 (80+89+0)130TBD2,448,000ET
W*Yongin-Yeongdeok130 (70+60+0)100TBD1,264,000SB
X*Uijeongbu-Howon130 (70+60+0)100JE777,300SW
Y*Yangpyeong-Changdae91 (45+46+0)70SC417,900ET
Z*Bucheon-Jungdong2260 (120+140+0)200TBD4,770,000WC
AA*Bucheon-Songnae130 (70+60+0)100TBD1,341,000WC
BB*Yangpyeong-Yanggeun444 (150+200+94)341NW1,631,000ET
CC*Gimpo-Sau65 (65+0+0)50TBD1,490,000SS
DD*Yangju-Nambang130 (70+60+0)100TBD275,500SW
EE*Suwon-Pajang430 (130+200+100)330TBD1,070,000SB

Table. 8.

Table 8. Results of optimal investment scheduling for Scenario 4.

 Project Resident type Land value (KRW/m2) Start year Number of houses Total Project Cost (108 KRW) 1st year 2nd year 3rd year 4th year A* JE 368,600 1 48 73 117 B* JE 860,500 1 48 26 35 C JE 2,897,000 2 300 200 190 D* JE 827,000 1 100 38 75 E JE 830,100 3 100 60 70 F IW 409,200 2 65 40 45 G NW/IW 399,300 3 290 120 200 H* JE 610,800 1 49 26 25 I IW 733,400 3 50 15 40 J* JE 1,346,000 1 150 45 100 37 K* NW 2,211,000 1 204 252 141 L* NW 2,428,000 1 56 71 M* JE 674,300 1 31 27 N* SC 2,145,000 1 152 78 118 O* SC 966,000 1 50 67 67 P* JE 1,151,000 1 50 13 43 Q NW 1,554,000 2 300 200 190 R JE 1,229,000 3 100 70 60 S SC 819,600 3 100 70 60 T NW 318,500 3 500 230 220 U IW 240,000 2 230 100 100 99 V TBD 2,448,000 2 130 80 89 W TBD 1,264,000 2 100 70 60 X JE 777,300 2 100 70 60 Y SC 417,900 3 70 45 46 Z TBD 4,770,000 2 200 120 140 AA TBD 1,341,000 2 100 70 60 BB NW 1,631,000 2 341 150 200 94 CC TBD 1,490,000 2 50 65 DD TBD 275,500 3 100 70 60 EE TBD 1,070,000 3 330 130 200 Total 4,494 716 1,886 1,981 1,149

1. OECD. Fertility rates. Available https://data.oecd.org/pop/fertility-rates.htm
2. Oh, CS, and Choi, SH (2012). The empirical study on the cause of low fertility-factors to impact on falling of nuptiality and rising of age at first marriage. Journal of Public Welfare Administration. 22, 91-125.
3. Kim, MY, and Hwang, JY (2016). Housing price and the level and timing of fertility in Korea: an empirical analysis of 16 cities and provinces. Health and Social Welfare Review. 36, 118-142.
4. Lee, CM, Kim, HA, and Cho, M (2012). A comparative analysis on house price to income ratios: focusing on measurement issues. Housing Studies Review. 20, 5-25.
5. Ministry of Land, Infrastructure and Transport (2012). Housing Business Handbook. Seoul: Ministry of Land, Infrastructure and Transport
6. Xie, F, Liu, C, Li, S, Lin, Z, and Huang, Y (2018). Long-term strategic planning of inter-city fast charging infrastructure for battery electric vehicles. Transportation Research Part E: Logistics and Transportation Review. 109, 261-276.
7. Liu, Y, Tang, W, He, J, Liu, Y, Ai, T, and Liu, D (2015). A land-use spatial optimization model based on genetic optimization and game theory. Computers, Environment and Urban Systems. 49, 1-14.
8. Li, X, and Parrott, L (2016). An improved genetic algorithm for spatial optimization of multi-objective and multi-site land use allocation. Computers, Environment and Urban Systems. 59, 184-194.
9. Haque, A, and Asami, Y (2014). Optimizing urban land use allocation for planners and real estate developers. Computers Environment and Urban Systems. 46, 57-69.
10. Taromi, R, DuRoss, M, Chen, B, Faghri, A, Li, M, and DeLiberty, T (2015). A multiobjective land development optimization model: the case of New Castle County, Delaware. Transportation Planning and Technology. 38, 277-304.
11. Zhang, W, and Huang, B (2015). Land use optimization for a rapidly urbanizing city with regard to local climate change: Shenzhen as a case study. Journal of Urban Planning and Development. 141.
12. Ng, ST, Skitmore, M, and Wong, KF (2008). Using genetic algorithms and linear regression analysis for private housing demand forecast. Building and Environment. 43, 1171-1184.
13. Ahn, JJ, Byun, HW, Oh, KJ, and Kim, TY (2012). Using ridge regression with genetic algorithm to enhance real estate appraisal forecasting. Expert Systems with Applications. 39, 8369-8379.
14. Park, JH, Geem, ZW, and Yu, JS (2016). The optimal investment portfolio for land development projects of public company. Korean Real Estate Review. 26, 23-38.

Jae Ho Park received his B.A. in economics from Yonsei University in 1997, M.A. in Graduate School of Real Estate & Construction from Dankook University in 2013, and currently pursuing Ph.D. in Graduate School of Urban Planning & Real Estate Studies at Dankook University. He passed the Uniform Certified Public Accountant Examination in California, U.S.A. in 2003. Now he works as a general manager at Gyeonggi Urban Innovation Corporation which develops new cities, housing lands, industrial complexes and constructs housing units to the vulnerable in Gyeonggi province, South Korea.

E-mail: jhpark1@gico.or.kr

Jung-Suk Yu is an associate professor in School of Urban Planning & Real Estate Studies at Dankook University. He teaches and researches in the field of real estate finance & investments. He holds a Ph.D. degree in financial economics from the University of New Orleans and served as a research fellow for Samsung Economic Research Institute.

E-mail: jsyu@dankook.ac.kr

ZongWoo Geem is an associate professor in College of Information Technology, Gachon University. He received his B.E. from Chung-Ang University, M.Sc. from Johns Hopkins University, and Ph.D. from Korea University. He has researched at Virginia Tech, University of Maryland -College Park, and Johns Hopkins University. He invented a music-inspired algorithm, Harmony Search, which has been applied to various optimization problems. His research interest includes phenomenon-mimicking algorithm and its applications to energy, environment and water fields.

E-Mail: geem@gachon.ac.kr

### Article

#### Original Article

Int. J. Fuzzy Log. Intell. Syst. 2018; 18(2): 135-145

Published online June 25, 2018 https://doi.org/10.5391/IJFIS.2018.18.2.135

## Genetic Algorithm-based Optimal Investment Scheduling for Public Rental Housing Projects in South Korea

Jae Ho Park, Jung Suk Yu, and Zong Woo Geem

1Graduate School of Urban Planning & Real Estate Studies, Dankook University, Yongin, Korea, 2School of Urban Planning & Real Estate Studies, Dankook University, Yongin, Korea, 3College of Information Technology, Gachon University, Seongnam, Korea

Correspondence to:Zong Woo Geem (zwgeem@gmail.com)

Received: April 7, 2018; Revised: June 9, 2018; Accepted: June 18, 2018

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

### Abstract

Declining birthrate is a serious problem that threatens the sustainability of Korean society. The main cause of this phenomenon is high living cost where housing cost accounts for the majority in household expenditure. South Korea has a very low supply rate in public rental housing when compared to other OECD countries. Because young people cannot afford to buy or lease a house for their new houses, some of them postpone or even give up marriage. As a countermeasure, Gyeonggi Province (surrounding area of Seoul) recently announced the supplying plan of 10,000 public rental houses by 2020. We expect this measure to alleviate the low birthrate problem and increase the demographic sustainability of the province. This study optimizes multi-annual investment scheduling for rental housing projects using genetic algorithm while satisfying the constraints such as budget, human resources, regional balance, etc. Through the optimal investment scheduling, we hope that public corporation will supply public rental houses more efficiently and more sustainably for the community.

Keywords: Public rental house, Optimal investment scheduling, Sustainable housing, Genetic algorithm

### 1. Introduction

Low birthrate is a serious problem that threatens the sustainability of Korean society. The OECD data shows the total fertility rate of Korea is 1.21 (children per woman) and it is the lowest among developed countries (United States 1.86, United Kingdom 1.81, France 1.98, Japan 1.42, Germany 1.47) in 2014 as shown in Figure 1. Assuming no migration and unchanged mortality, a total fertility rate of 2.1 ensures a stable population.

Two main reasons for the declining birthrate in Korea are unstable employment and high housing cost (purchasing or rental). Recent studies by Oh and Choi [2] and Kim and Hwang [3] have shown that higher housing cost is an important factor for fertility decline and postponement of having children. According to Oh and Choi [1], young people are more influenced by rental cost rather than house-purchasing price. Thus, it is important to supply public rental houses for newlyweds. The survey in 2015 by Korean Ministry of Health and Welfare shows that the housing is the number one factor related to marriage, as seen in Figure 2.

Price to income ratio (PIR) in house has been widely used, both as an indicator of house price level and as that of affordability to acquire residential property, in different geographical areas or time periods. According to Lee et al. [4], PIR in Korea (4.4) is higher than that in the United States (3.5) and Canada (3.4), but lower than those in United Kingdom (5.2), Australia (6.1) and Hong Kong (11.4). Moreover, Korea’s long-term public rental housing ratio is considerably lower than that of major developed countries as observed in Table 1.

Since there are differences in conditions and social environments, it may not be adequate to simply compare the percentage of long-term public rental housing ratio. However, it is an undeniable fact that long-term public rental houses are very deficient in South Korea.

As a countermeasure to this problem, the government of Gyeonggi Province, metropolitan area containing Seoul, announced the plan to supply 10,000 public rental houses, called Dda Bok House (DBH; here, ‘Dda Bok’ means warm-hearted welfare), by 2020 on May 17, 2016.

DBH is a unique housing program that combines public rental housing and rental subsidy policies to significantly lower the burden of housing cost. In this program, the rental subsidy is proportional to the number of the babies. While the area of typical public rental house (two-room type) supplied to the newlyweds couple is 36 m2, the area of DBH is 44 m2, which is 22% wider than the current two-room type residence. Table 2 shows the comparison in details.

DBH plans to provide a safe nursing environment as well as autonomous community by creating various community-based programs. Thus, talent donators, who can give childcare, medical service, cooking, arts activity, physical activity, and local community movement, have a priority to live in the community.

The purpose of this study is to find the optimal investment scheduling for the multi-annual DBH project by maximizing the utility while meeting all the constraints.

So far, many researchers have performed useful research in the field of land-project optimization. Fei Xie et al. [6] introduced a multistage chance-constrained stochastic model for strategic planning of battery electric vehicle (BEV) inter-city fast charging infrastructure. The model was applied to a case study in California and solved by genetic algorithm. Liu et al. [7] performed land-use spatial optimization by coordinating the competitions between different land-use types in the example of Gaoqiao Town, Zhejiang Province, China. Li and Parrott [8], Haque and Asami [9], Taromi et al. [10] and Zhang and Huang [11] reached a consensus that evolutionary algorithm-based model is powerful tool in optimizing land use allocation in the cases of the Delaware, USA and Shenzhen, China, etc.

Other researchers have mainly been interested in private housing demand forecasting [12] and real estate appraisal forecasting with ridge regression [13]. Overall, these studies highlighted the suitability of evolutionary algorithm in the research regarding the land planning and housing demand forecasting. Unlike these studies, Park et al. [14] proposed an optimization model for deriving the optimal investment combination from the various real estate development projects. With these tools, public corporation could have reasonable investment portfolio that maximize the return of investment while meeting the multiple constraints.

In this study, an optimization model is further developed and applied to optimal investment scheduling for public rental housing projects with different groundbreaking years in South Korea.

### 2. Problem Description

Gyeonggi Urban Innovation Corporation (GICO), the subsidiary of Gyeonggi Province, is in charge of providing DBH effectively with the above-mentioned missions while considering the constraints such as human resources, annual budget, regional balance, etc., as follows:

• 1) Because DBH should be supplied in many regions simultaneously, GICO needs a lot of supervisors for construction and personnel for getting permission from local municipality. However, if GICO hires many employees at initial stage, it is a heavy burden to maintain those people after DBH project is finished. Therefore, it is desirable to maintain the proper number of personnel.

• 2) In order to maintain the financial soundness, GICO should make annually balanced investment rather than initial excessive investment as the DBH project is non-for-profit business.

• 3) In order to maintain regional balance, Gyeonggi Province is divided into 5 sub-regions such as Seoul-Busan (SB) area in orange color, West coastal (WC) area in sky-blue color, Seoul-Sinuiju (SS) area in blue color, Seoul-Wonsan (SW) area in grey color, and Eastern (ET) area in green color, as shown in Figure 3. Any sub-region cannot be excessively developed.

• 4) DBHs are built for newlywed, college students, junior employees and senior citizens. Thus, it is important to fairly

### 3. Genetic Algorithm

Genetic algorithm (GA) was first proposed by John Holland in the 1960s and was further developed by himself and his students and colleagues at the University of Michigan in the 1960s and the 1970s. Holland’s 1975 book “Adaptation in Natural and Artificial Systems” presented the GA as an abstraction of biological evolution and gave a theoretical framework for adaptation under the GA. GA is a method for moving from one population of “chromosomes” (e.g., strings of ones and zeros, or “bits”) to a new population by using a kind of “natural selection” together with the genetics-inspired operators of crossover, mutation, and inversion. Each chromosome consists of “genes” (e.g., bits), each gene being an instance of a particular “allele” (e.g., 0 or 1). The selection operator chooses those chromosomes in the population that will be allowed to reproduce, and on average the fitter chromosomes produce more offspring than the less fit ones. Crossover exchanges subparts of two chromosomes, roughly mimicking biological recombination between two single-chromosome (“haploid”) organisms; mutation randomly changes the allele values of some locations in the chromosome; and inversion reverses the order of a contiguous section of the chromosome, thus rearranging the order in which genes are arrayed.

Given a clearly defined problem to be solved and a bit string representation for candidate solutions, a simple GA works as follows:

• 1) Start with a randomly generated population of n l-bit chromosomes (candidate solutions to a problem).

• 2) Calculate the fitness f(x) of each chromosome x in the population.

• 3) Repeat the following steps until n offspring have been created:

• Select a pair of parent chromosomes from the current population, the probability of selection being an increasing function of fitness. Selection is done “with replacement,” meaning that the same chromosome can be selected more than once to become a parent.

• With probability Pc (the “crossover probability” or “crossover rate”), crossover the pair at a randomly chosen point (chosen with uniform probability) to form two offspring. If no crossover takes place, form two offspring that are exactly copies of their respective parents.

• Mutate the two offspring at each locus with probability Pm (the “mutation probability or “mutation rate”), and place the resulting chromosomes in the new population. If n is odd, one new population member can be discarded at random.

• 4) Replace the current population with the new population.

• 5) Go to step 2).

For this study, GA utilized parameter values: population size is 100, and mutation rate is 0.075.

For GA computation, we used Solver software from Frontline Systems, Inc. Average computational time ranges from 15 to 20 seconds per each run and termination condition is when convergence rate is equal to or less than 0.0001.

### 4. Optimization Formulation and Computational Results

To easily understand the optimal investment scheduling, 10 sample projects were selected as depicted in Table 3. Each project has different total project cost, resident type, land value and region. We want to find the optimal investment schedule by maximizing the effect of DBH supply while satisfying various constraints.

We think that the effect or utility of DBH can be maximized when more valuable houses are supplied earlier than less valuable ones. Although the housing value is not easily quantified, we instead use the appraised value of the project land. Therefore, the objective function of this scheduling problem is Max ∑ HiPiW(Yi), where Hi is the number of houses in project i; Pi is the appraised value of land in project i; and W(Yi) is weighting factor of groundbreaking year. To find the optimal investment schedule with the maximized objective function value, we set up weighting factor system. Because the earlier the better, if the groundbreaking year is 1, the weighting factor becomes 3 (highest value); if the groundbreaking year is 2, the weighting factor becomes 2 (medium value); and if the groundbreaking year is 3, the weighting factor becomes 1 (lowest value).

Figure 4 shows the structure of chromosome incorporating decision variables used in GA computation. The value of each decision variable is groundbreaking year (1st, 2nd or 3rd year).

### 4.1 Scenario 1

Here, we can consider several different scenarios. In Scenario 1, we set groundbreaking year (1st, 2nd, or 3rd year) of project i as decision variable. If Yi is decision variable, it has one of three candidate values, Yi ∈ {1, 2, 3}.

As previously mentioned, the objective function is Max ∑i HiPiW(Yi), where W(Yi) is weighting factor which becomes W(Yi) = 3 if Yi = 1; W(Yi) = 2 if Yi = 2; and W(Yi) = 1 if Yi = 3.

For the human resources limitation, we limit maximum four projects to launch in year j, N(Yi = j) ≤ 4, j = 1, 2, 3, where N(●) is number function which counts the total number of launched projects in year j.

For the project budget limitation, the maximum annual investment amount is 600 (unit= 100 million KRW), $∑iCij≤600$, j = 1, 2, 3, 4. Here, $Cij$ is partial investment cost occurred in year j for project i, and 600 is about 36% of total 1,647 (100 million KRW).

For regional balance, any specific area cannot be excessively developed. In this example, among five areas (ET, SB, WC, SW, and SS), the maximum number of annually launched projects for SB area is limited up to two (the total number of SB projects in this example is four), N(Yi = j & Ai = SB) ≤ 2, j = 1, 2, 3. Here, Ai is sub-region where project i is included.

For resident-type balance, any specific class cannot be excessively developed. In this example, among three classes (junior employee, newlywed, and industrial worker), the maximum number of annually launched projects for junior employee is limited up to two (the total number of junior employee projects in this example is six), N(Yi = j & Ri =junior employee)≤ 2, j = 1, 2, 3. Here, Ri is resident type of project i.

The optimized results for this investment schedule example are shown in Table 4.

As observed in Table 4, projects B, C, I, and J will launch in 1st year; projects D, E, and G will launch in 2nd year; and projects A, F, and H will launch in 3rd year.

For the human resources limitation, this constraint is satisfied because the number of launched projects in 1st year is 4 (B, C, I, J), that in 2nd year is 3 (D, E, G), and that in 3rd year is 3 (A, F, H).

For the project budget limitation, this constraint is satisfied because the investment amount in 1st year is 286, that in 2nd year is 583, that in 3rd year is 531, and that in fourth year is 247. Every of them is less than 600.

For regional balance, this constraint is satisfied because the number of SB projects in 1st year is 2 (C, J), that in 2nd year is 2 (D, E), and that in 3rd year is 0. Every of them is less than two.

For resident-type balance, this constraint is satisfied because the number of junior employee projects in 1st year is 2 (C, J), that in 2nd year is 2 (D, E), and that in 3rd year is 2 (A, H). Every of them is less than two.

### 4.2 Scenario 2

In Scenario 2, two constraints are loosened from Scenario 1. While the maximum number of launched projects in each year was 4 in Scenario 1, that becomes 5 in Scenario 2. And, while the maximum number of junior employee projects in each year was 2 in Scenario 1, that becomes 4 in Scenario 2. Table 5 shows the detailed results from Scenario 2.

As observed in Table 5, projects B, C, H, and J will launch in 1st year; projects D, E, G, and I will launch in 2nd year; and projects A and F will launch in 3rd year. The differences between Scenarios 1 and 2 are projects H and I. The investment schedule of project H is moved forward from 3rd year to 1st year, and that of project I is moved back from 1st year to 2nd year for obtaining better effect of DBH project. The objective function value increased from 4,084 to 4,107.

For the human resources limitation, this constraint is satisfied because the number of launched projects in 1st year is 4 (B, C, H, J), that in 2nd year is 4 (D, E, G, I), and that in 3rd year is 2 (A, F). Every of them is less than 5.

For the project budget limitation, this constraint is satisfied because the investment amount in 1st year is 297, that in 2nd year is 583, that in 3rd year is 535, and that in 4th year is 232. Every of them is less than 600.

For regional balance, this constraint is satisfied because the number of SB projects in 1st year is 2 (C, J), that in 2nd year is 2 (D, E), and that in 3rd year is 0. Every of them is less than 2.

For resident-type balance, this constraint is satisfied because the number of junior employee projects in 1st year is 3 (C, H, J), that in 2nd year is 2 (D, E), and that in 3rd year is 1 (A). Every of them is less than 4.

### 4.3 Scenario 3

In Scenario 3, the maximum annual budget is tightened from Scenario 2. We apply maximum budget of 500 (100 million KRW) instead of 600. Table 6 shows the detailed results from Scenario 3.

As observed in Table 6, the schedules of two projects (G and I) have been changed from Scenario 2. The investment schedule of project G is moved back to 3rd year from 2nd year, and that of project I is moved forward to 1st year from 2nd year. Also, the objective function value decreased from 4,107 to 4,028.

For the human resources limitation, this constraint is satisfied because the number of launched projects in 1st year is 5 (B, C, H, I, J), that in 2nd year is 2 (D, E), and that in 3rd year is 3 (A, F, G). Every of them is less than 5.

For the project budget limitation, this constraint is satisfied because the investment amount in 1st year is 312, that in 2nd year is 488, that in 3rd year is 425, and that in 4th year is 362. Every of them is less than 500.

For regional balance, this constraint is satisfied because the number of SB projects in 1st year is 2 (C, J), that in 2nd year is 2 (D, E), and that in 3rd year is 0. Every of them is less than 2. For resident-type balance, this constraint is satisfied because the number of junior employee projects in 1st year is 3 (C, H, J), that in 2nd year is 2 (D, E), and that in 3rd year is 1 (A). Every of them is less than 4.

### 4.4 Scenario 4 (Real-World Example)

Now we consider a real-world investment schedule problem, which has 31 projects. Table 7 shows the details of the projects. Because the construction costs of some projects marked with asterisk are not determined yet, those costs are assumed by multiplying 130 million KRW and number of houses. Here, resident type consists of junior employee (JE), newlywed (NW), industrial worker (IW), senior citizen (SC), and to-be-determined (TBD).

Actually, 11 projects (A, B, D, H, J, K, L, M, N, O, and P) have already launch-ed in 2017 by GICO. Thus, the decision variables of those projects should have a value of one, which means 1st year, and those of other 20 projects should have a value of two or three, Yi ∈ {2, 3}. Here, the number of candidate investment schedule solutions becomes 220 (≥ 1, 000, 000) because 20 projects have two choices. As already utilized, the objective function is Max ∑i HiPiW(Yi).

For the human resources limitation, we limit maximum 11 projects to launch in year j, N(Yi = j) ≤ 11, j = 2, 3.

For the project budget limitation, the maximum annual investment amount is 2,000 (108 KRW), $∑iCij≤2,000$, j = 2, 3, 4. Here, 2,000 is around 1/3 of total project cost 6,102.

For regional balance, the maximum number of annually launched projects for SB area is limited up to 3, N(Yi = j & Ai = SB) ≤ 3, j = 2, 3.

For resident-type balance, the maximum number of annually launched projects for junior employee is limited up to 2, N(Yi = j & Ri = JE) ≤ 2, j = 2, 3.

The optimized results for this real-world investment schedule problem are shown in Table 8.

As observed in Table 8, 11 projects (C, F, Q, U, V, W, X, Z, AA, BB, CC) will launch in 2nd year; and 9 projects (E, G, I, R, S, T, Y, DD, EE) will launch in 3rd year.

For the human resources limitation, this constraint is satisfied because the number of launched projects in 1st year is 11, and that in 2nd year is 9. Every of them is less than 11.

For the project budget limitation, this constraint is satisfied because the investment amount in 2nd year is 1,886, that in 3rd year is 1,981, and that in 4th year is 1,149. Every of them is less than 2,000.

For regional balance, this constraint is satisfied because the number of SB projects in 2nd year is 3 (C, Q, W), and that in 3rd year is 3 (E, R, EE). Every of them is less than 3.

For resident-type balance, this constraint is satisfied because the number of JE projects in 2nd year is 2 (C, X), and that in 3rd year is 2 (E, R). Every of them is less than 2.

### 5. Conclusions

Public rental house is essential for low-income classes who cannot afford to purchase a house or pay a high security deposit (some property owners require 80% of house prices). South Korea has lowest fertility rate among OECD countries and the main reason is cost burden of house purchase or rent. Furthermore, public rental housing ratio in South Korea is considerably lower than that of major OECD countries. Thus, it is urgent to expand the supply of public rental houses for the housing sustainability of South Korea.

This study proposed a GA-based optimization model for multi-annual investment scheduling of public rental house project, obtaining good results. This cannot be easily obtained if we use only “engineering intuition”. The optimization model becomes very useful and practical tool for the massive DBH project while meeting various constraints such as limited human resources, limited annual budget, regional supply balance, and resident-type balance. We hope this optimization approach will be used for the investment scheduling of similar public projects in the future.

The proposed scheduling optimization model may be further applied to various cases such as financial investment, advisory service, and facility repair, as follows:

• - Financial investment: Financial investor may use the scheduling optimization model with various investment projects, which have different maturity, rate of return, and risk with limited financial resources to obtain multi-annual desirable cash flow.

• - Advisory service: Consulting firm, accounting firm, or law firm may use the scheduling optimization model with limited human resources to carry out more consulting projects simultaneously.

• - Facility repair: City officials may use the scheduling optimization model for lots of infrastructure facilities such as road, water supply and sewer pipe with limited human and financial resources to maintain facilities.

### Acknowledgements

This research was supported by the Ministry of Science and ICT, Korea, under the Information Technology Research Center support program (No. IITP-2018-2017-0-01630) supervised by the Institute for Information & communications Technology Promotion (IITP).

### Fig 1.

Figure 1.

Fertility rates in 2014. From OECD data [1].

The International Journal of Fuzzy Logic and Intelligent Systems 2018; 18: 135-145https://doi.org/10.5391/IJFIS.2018.18.2.135

### Fig 2.

Figure 2.

Desirable measures to support marriage.

The International Journal of Fuzzy Logic and Intelligent Systems 2018; 18: 135-145https://doi.org/10.5391/IJFIS.2018.18.2.135

### Fig 3.

Figure 3.

Major locations of DBH project in Gyeonggi Province. supply with respect to class balance.

The International Journal of Fuzzy Logic and Intelligent Systems 2018; 18: 135-145https://doi.org/10.5391/IJFIS.2018.18.2.135

### Fig 4.

Figure 4.

Structure of chromosome incorporating decision variables.

The International Journal of Fuzzy Logic and Intelligent Systems 2018; 18: 135-145https://doi.org/10.5391/IJFIS.2018.18.2.135

Comparison of long-term public rental house ratio in various countries.

CountryKoreaUKFranceJapanOECD averageEU average
Long-term public rental house ratio5.5%17%19%6.1%11.5%9.4%
Base Year201420102007200820072007

From Ministry of Land, Infrastructure and Transport, 2012 [5]..

Comparison between current rental house and DBH.

Current two-room houseDBH
Area, exclusive (m2)3644
Floor plan

Samples of DBH projects.

ProjectProject districtTotal project cost (1st, 2nd, 3rd year) (KRW 100 million)Number of houses (H)Resident type (R)Appraised value of land (P) (KRW/m2)Region (A)
AGapyeong-eup190 (73+117+0)48Junior employee368,600ET
BNamyangju-Changhyeon61 (26+35+0)48Newlywed860,500ET
CSeongnam-Pangyo390 (200+190+0)300Junior employee2,897,000SB
DSuwon-Mangpo113 (38+75+0)100Junior employee827,000SB
ESuwon-Seodun130 (60+70+0)100Junior employee830,100SB
FSiheung-Sincheon85 (40+45+0)65Industrial worker409,200WC
GSiheung-Jeongwang380 (120+200+60)290Newlywed/Industrial worker399,300WC
HYangpyeong-Gongheung51 (26+25+0)49Junior employee610,800ET
IOsan-Gajang65 (15+40+10)50Industrial worker733,400WC
JYongin-Jukjeon182 (45+100+37)150Junior employee1,346,000SB

Results of optimal investment scheduling for Scenario 1.

Project nameResident typeAppraised value of land (KRW/m2)Start yearWeighted valueNumber of houses (H)Total project cost (KRW 100 million)
1st year2nd year3rd year4th year
AJunior employee368,6003148--73117
BNewlywed860,50013482635--
CJunior employee2,897,00013300200190--
DJunior employee827,00022100-3875-
EJunior employee830,10022100-6070-
FIndustrial worker409,2003165--4045
GNewlywed/Industrial worker399,30022290-12020060
HJunior employee610,8003149--2625
IIndustrial worker733,4001350154010-
JJunior employee1,346,000131504510037-
Total1,200286583531247

Results of optimal investment scheduling for Scenario 2.

ProjectResident typeAppraised value of land (KRW/m2)Start yearWeighted valueNumber of houses (H)Total project cost (KRW 100 million)
1st year2nd year3rd year4th year
AJunior employee368,6003148--73117
BNewlywed860,50013482635--
CJunior employee2,897,00013300200190--
DJunior employee827,00022100-3875-
EJunior employee830,10022100-6070-
FIndustrial worker409,2003165--4045
GNewlywed/Industrial worker399,30022290-12020060
HJunior employee610,80013492625
IIndustrial worker733,4002250154010
JJunior employee1,346,000131504510037-
Total1,200297583535232

Results of optimal investment scheduling for Scenario 3.

ProjectResident typeAppraised value of land (KRW/m2)Start yearWeighted valueNumber of houses (H)Total project cost (KRW 100 million)
1st year2nd year3rd year4th year
AJunior employee368,6003148--73117
BNewlywed860,50013482635--
CJunior employee2,897,00013300200190--
DJunior employee827,00022100-3875-
EJunior employee830,10022100-6070-
FIndustrial worker409,2003165--4045
GNewlywed/Industrial worker399,30031290-120200
HJunior employee610,80013492625
IIndustrial worker733,4001350154010
JJunior employee1,346,000131504510037-
Total1,200312488425362

Total DBH projects.

ProjectProject districtProject cost (1st, 2nd, 3rd year) (108 KRW)Number of housesResident typeLand value (KRW/m2)Region
AGapyeong-eup190 (73+117+0)48JE368,600ET
BNamyangju-Changhyeon61 (26+35+0)48NW860,500ET
C*Seongnam-Pangyo390 (200+190+0)300JE2,897,000SB
DSuwon-Mangpo113 (38+75+0)100JE827,000SB
ESuwon- Seodun130 (60+70+0)100JE830,100SB
F*Siheung- Sincheon85 (40+45+0)65IW409,200WC
G*Siheung-Jeongwang380 (120+200+60)290NW/IW399,300WC
HYangpyeong-Gongheung51 (26+25+0)49JE610,800ET
I*Osan-Gajang65 (15+40+10)50IW733,400WC
JYongin-Jukjeon182 (45+100+37)150JE1,346,000SB
KSuwon-Gwanggyo1393 (252+141+0)204NW2,211,000SB
LAnyang-Gwanyang171 (71+0+0)56NW2,428,000SB
MHwaseong-Jinan27 (27+0+0)31JE674,300WC
NSuwon-Woncheon196 (78+118+0)152SC2,145,000SB
OPaju-Byungwon134 (67+67+0)50SC966,000SS
PUiwang-Bugok56 (13+43+0)50JE1,151,000SB
Q*Anyang-Gwanyang2390 (200+190+0)300NW1,554,000SB
R*Yongin-Yeosung130 (70+60+0)100JE1,229,000SB
S*Icheon-Gwango130 (70+60+0)100SC819,600ET
T*Gwangju-Yeokdong650 (230+220+200)500NW318,500ET
U*Ansan-Wonsi299 (100+100+99)230IW240,000WC
V*Hanam-Deokpung2169 (80+89+0)130TBD2,448,000ET
W*Yongin-Yeongdeok130 (70+60+0)100TBD1,264,000SB
X*Uijeongbu-Howon130 (70+60+0)100JE777,300SW
Y*Yangpyeong-Changdae91 (45+46+0)70SC417,900ET
Z*Bucheon-Jungdong2260 (120+140+0)200TBD4,770,000WC
AA*Bucheon-Songnae130 (70+60+0)100TBD1,341,000WC
BB*Yangpyeong-Yanggeun444 (150+200+94)341NW1,631,000ET
CC*Gimpo-Sau65 (65+0+0)50TBD1,490,000SS
DD*Yangju-Nambang130 (70+60+0)100TBD275,500SW
EE*Suwon-Pajang430 (130+200+100)330TBD1,070,000SB

Results of optimal investment scheduling for Scenario 4.

 Project Resident type Land value (KRW/m2) Start year Number of houses Total Project Cost (108 KRW) 1st year 2nd year 3rd year 4th year A* JE 368,600 1 48 73 117 B* JE 860,500 1 48 26 35 C JE 2,897,000 2 300 200 190 D* JE 827,000 1 100 38 75 E JE 830,100 3 100 60 70 F IW 409,200 2 65 40 45 G NW/IW 399,300 3 290 120 200 H* JE 610,800 1 49 26 25 I IW 733,400 3 50 15 40 J* JE 1,346,000 1 150 45 100 37 K* NW 2,211,000 1 204 252 141 L* NW 2,428,000 1 56 71 M* JE 674,300 1 31 27 N* SC 2,145,000 1 152 78 118 O* SC 966,000 1 50 67 67 P* JE 1,151,000 1 50 13 43 Q NW 1,554,000 2 300 200 190 R JE 1,229,000 3 100 70 60 S SC 819,600 3 100 70 60 T NW 318,500 3 500 230 220 U IW 240,000 2 230 100 100 99 V TBD 2,448,000 2 130 80 89 W TBD 1,264,000 2 100 70 60 X JE 777,300 2 100 70 60 Y SC 417,900 3 70 45 46 Z TBD 4,770,000 2 200 120 140 AA TBD 1,341,000 2 100 70 60 BB NW 1,631,000 2 341 150 200 94 CC TBD 1,490,000 2 50 65 DD TBD 275,500 3 100 70 60 EE TBD 1,070,000 3 330 130 200 Total 4,494 716 1,886 1,981 1,149

### References

1. OECD. Fertility rates. Available https://data.oecd.org/pop/fertility-rates.htm
2. Oh, CS, and Choi, SH (2012). The empirical study on the cause of low fertility-factors to impact on falling of nuptiality and rising of age at first marriage. Journal of Public Welfare Administration. 22, 91-125.
3. Kim, MY, and Hwang, JY (2016). Housing price and the level and timing of fertility in Korea: an empirical analysis of 16 cities and provinces. Health and Social Welfare Review. 36, 118-142.
4. Lee, CM, Kim, HA, and Cho, M (2012). A comparative analysis on house price to income ratios: focusing on measurement issues. Housing Studies Review. 20, 5-25.
5. Ministry of Land, Infrastructure and Transport (2012). Housing Business Handbook. Seoul: Ministry of Land, Infrastructure and Transport
6. Xie, F, Liu, C, Li, S, Lin, Z, and Huang, Y (2018). Long-term strategic planning of inter-city fast charging infrastructure for battery electric vehicles. Transportation Research Part E: Logistics and Transportation Review. 109, 261-276.
7. Liu, Y, Tang, W, He, J, Liu, Y, Ai, T, and Liu, D (2015). A land-use spatial optimization model based on genetic optimization and game theory. Computers, Environment and Urban Systems. 49, 1-14.
8. Li, X, and Parrott, L (2016). An improved genetic algorithm for spatial optimization of multi-objective and multi-site land use allocation. Computers, Environment and Urban Systems. 59, 184-194.
9. Haque, A, and Asami, Y (2014). Optimizing urban land use allocation for planners and real estate developers. Computers Environment and Urban Systems. 46, 57-69.
10. Taromi, R, DuRoss, M, Chen, B, Faghri, A, Li, M, and DeLiberty, T (2015). A multiobjective land development optimization model: the case of New Castle County, Delaware. Transportation Planning and Technology. 38, 277-304.
11. Zhang, W, and Huang, B (2015). Land use optimization for a rapidly urbanizing city with regard to local climate change: Shenzhen as a case study. Journal of Urban Planning and Development. 141.
12. Ng, ST, Skitmore, M, and Wong, KF (2008). Using genetic algorithms and linear regression analysis for private housing demand forecast. Building and Environment. 43, 1171-1184.
13. Ahn, JJ, Byun, HW, Oh, KJ, and Kim, TY (2012). Using ridge regression with genetic algorithm to enhance real estate appraisal forecasting. Expert Systems with Applications. 39, 8369-8379.
14. Park, JH, Geem, ZW, and Yu, JS (2016). The optimal investment portfolio for land development projects of public company. Korean Real Estate Review. 26, 23-38.