在現實生活中，廠址選擇是一個在任何產業界皆非常重要而且常見的議題，這一類的問題在學術上稱為有限產能的分群廠址選擇問題（Capacitated Clustering Location Problem, CCP）。CCP主要是在一群已知的顧客點中，找尋已知所需數量之配銷中心 (Distribution center, DC)的最佳位置。廠址選擇需要考量的成本包含許多，但在過去研究中，最主要考量的目標為最小化從每個DC到每個所屬的顧客點的距離之加總，使得整體運送成本能降到最低。如此挑選DC的準則在於，每當顧客有需求時可以保證從DC到顧客點的運送距離為最小，實際的應用如：消防局、醫院等的設置，符合此種配送貨物或是提供服務的運送方式在本研究中我們稱之為「分批配送」，但並非所有產業或是公司之配送方式都如上述，某些行業必須在固定的時間一次出貨給所有的顧客才返回DC，實際的應用包含：超商生鮮商品的補給車、送報員、郵差等，此種首次被提出應用於廠址選擇問題的配送方式模型，在本研究中稱為「一次性配送」。在此概念下的場址選擇問題目標則為最小化從DC出發後，走完每個顧客點再回到DC的距離加總，一次性配送的運送方式將可以視為旅行銷售員問題做處理以求出最短迴路。
由於CCP是個NP-hard問題，本研究將提出一個混合式的分群演算法：簡化群體演算法結合調和平均數分群法 (SSOKHM)，在需求限制下來挑選特定數量DC的最佳位置，並在分群後，使用貪婪演算法來獲得每群之最短迴路以最小化整體之運送成本。這也是簡化群體演算法首次使用於廠址選擇問題上。
透過實驗結果與其他三種混合式分群法做驗證，顯示無論在分群的結果及整體運送成本上，SSOKHM演算法都具有相當穩定且優異的成效。且同時本研究在貪婪法後加入了一個檢查機制：交換局部搜尋(Exchange local search)，此機制可以透過移動顧客點所屬的分群，使廠址選擇更加精確藉此來節省公司之運送成本。
關鍵詞：分群場址選擇、簡化群體演算法、調和平均數分群、旅行銷售員問題

How to set up the distribution centers not only a very general and important issue but also a quite profound knowledge for many different industries in real life. This kind of problem called capacitated clustering location problem (CCP), then the previous literature, the target of the location problem is to minimize the sum of distances from each cluster centers to all customers in their cluster and the practical application of this type, for example: logistics centers, police stations, etc., this approach of shipping called “batch delivery”.
However, some types of business or the company shipping methods are different from the above. The distribution centers must ship goods to all customers in its cluster or cluster in a cycle time, such as postman, newsman, etc., this approach of shipping is called “one-time delivery”. One-time delivery can be modeled as a Travelling Salesman Problem (TSP). This type of problem is rarely discussed in the former literature.
In order to solve this practical problem, a new hybrid clustering algorithm, Simplified Swarm Optimization combining with K-harmonic means (SSOKHM) is proposed in this paper for the CCP with one-time delivery and using greedy algorithm for the TSP to find the shortest loop and obtain the target: minimize shipping costs. An exchange local search is an inspection mechanism which adds after greedy algorithm. Proposed method is applied to several facility location problems from OR library and compared with other hybrid clustering algorithm.
Numerical results show that the proposed approach (SSOKHM) performance is better than using other hybrid clustering algorithms in terms of shipping costs. Finally, we embed the exchange local search for each hybrid clustering algorithm. The results are presented that this inspection mechanism can be work in capacity constraint to enhance the performance of the solution and demonstrate the usefulness of this approach in CCP.
Keywords: Capacitated Clustering Location Problem, Simplified Swarm Optimization (SSO), K-harmonic Means (KHM), Travelling Salesman Problem

Abstract I
中文摘要 III
List of figures VI
List of table VII
Chapter 1. Introduction 1
Chapter 2. Literature review 6
2.1 Location Problem 6
2.2 Clustering 8
2.3 Traveling Salesman Problem (TSP) 9
Chapter 3. Model Formulation 11
3.1 Problem Statement 11
3.2 The Mathematical Model 12
3.3 Summary and Conclusion 13
Chapter 4. Methodology 14
4.1 K-harmonic Means (KHM) 14
4.2 Simplified Swarm Optimization (SSO) 16
4.3 The Proposed Hybrid Algorithm 17
4.4 Exchange Local Search 19
Chapter 5. Experimental Results 24
Chapter6. Conclusion 33
Reference 34

[1] Cooper, L., Location–allocation problems. Operations Research, 1963. 11: p. 331–343.
[2] Cooper L., The transportation–location problems. Operations Research, 1972. 20: p. 94–108.
[3] Cooper L. Heuristic methods for location–allocation problem. SIAM Review, 1964. 6: p. 37–53.
[4] Chen, C.H., Ting, C.J., Combining Lagrangian heuristic and Ant Colony System to solve the Single Source Capacitated Facility Location Problem, Transportation Research Part, 2008: p. 1099–1122.
[5] Osman, I.H., Christofides, N., Capacitated clustering problems by hybrid simulated annealing and tabu search, International Transactions in Operational Research. 1994. 1 (3): p. 317–336.
[6] Garey, M.R. & Johnson, D.S., A Guide to the Theory of NP - Comp1eteness. Computers and Intractability, (1979) San Francisco: Freeman.
[7] Mulvey, J.M., Beck, M.P., Solving capacitated clustering problems. North-Holland European Journal of Operational Research, 1984.18: p. 339-348.
[8] Scheuerer, S., Wendolsky, R., A scatter search heuristic for the capacitated clustering problem, European Journal of Operational Research, 2006. 169: p. 533–547.
[9] Garey, M.R., Johnson, D.S., Witsenhausen, H.S., The complexity of the generalized Lloyd–Max problem, IEEE Trans. Inform. Theory, 1982. 28 (2): p. 255–256.
[10] Güngör, Z., Ünler, A., K-harmonic means data clustering with simulated annealing heuristic, Applied Mathematics and Computation, 2007. 184: p. 199-209.
[11] Barcelo, J., Casanovas, J., A heuristic capacitated lagrangean algorithm for the plant location problem. North-Holland European Journal of Operational Research, 1984. 15: p. 212-226.
[12] Fisher, M., Jaikumar, R., A generalized assignment heuristic for the vehicle routing problem, Networks, 1981. 11 (2): p. 109–124.
[13] Koskosidis, Y.A., Powell, W.B., Clustering algorithms for consolidation of customer orders into vehicle shipments, Transportation Research, 1992. B-26(5): p. 365–379.
[14] Esnaf, S., Küçükdeniz, T., A fuzzy clustering-based hybrid method for a multi-facility location problem, J Intell Manuf, 2009. 20: p. 259–265.
[15] Xu, G., Xu, J., LP rounding algorithm for approximating uncapacitated facility location problem with penalties, Information Processing Letters, 2005. 94: p. 119–123.
[16] Sun, M., Solving the uncapacitated facility location problem using tabu search, Computers & Operations Research, 2006. 33: p. 2563 – 2589.
[17] Vinod, H.D., Integer programming and the theory of clustering. Journal of the American Statistical Association, 1969. 64: p. 506-519.
[18] Rao, M.R., Clustering analysis and mathematical programming. Journal of the American Statistical Association, 1971. 66: p. 622-626.
[19] Mulvey, J.M., Crowder, H.P., Cluster analysis: an application Lagrangian relaxation. Management Science, 1979. 25: p. 329-340.
[20] Jacobsen, S.K., Heuristics for the capacitated plant location model. North-Holland Publishing Company European Journal of Operational Research, 1983. 12: p. 253-261.
[21] Scheuerer, S., Wendolsky, R., A scatter search heuristic for the capacitated clustering problem, European Journal of Operational Research, 2006. 169: p. 533–547.
[22] Doring, C., Lesot, M., Kruse, R., Data analysis with fuzzy clustering methods. Computational Statistics & Data Analysis, 2006. 51: p. 192–214.
[23] Küçükdeniz, T., Baray, A., Ecerkale, K., Esnaf, S., Integrated use of fuzzy c-means and convex programming for capacitated multi-facility location problem, Expert Systems with Applications, 2012. 39: p. 4306–4314.
[24] Forgy, E.W., Cluster analysis of multivariate data: Efficiency versus interpretability of classifications, Biometrics, 1965. 21 (3): p. 768–769.
[25] Zhang, B., Hsu, M., Dayal, U., K-harmonic means — A data clustering algorithm, Technical Report HPL-1999-124, Hewlett-Packard Laboratories. 1999.
[26] Hammerly, G., Elkan, C., Alternatives to the k-means algorithm that find better clustering, Proc. of the 11th Int. Conf. on Information and Knowledge Management, 2002: p. 600–607.
[27] Güngör, Z., Ünler, A., K-harmonic means data clustering with simulated annealing heuristic, Applied Mathematics and Computation, 2007. 184: p. 199-209.
[28] Güngör, Z., Ünler, A., K-Harmonic means data clustering with tabu-search method, Applied Mathematical Modelling, 2008. 32: P. 1115–1125.
[29] Jiang, H., Yi, S., Li, J., Yang, F., Hu, X., Ant clustering algorithm with K-harmonic means clustering, Expert Systems with Applications, 2010. 37: p. 8679–8684.
[30] Yin, M., Hu, Y., Yang, F., Li, X., Gu, W., A novel hybrid K-harmonic means and gravitational search algorithm approach for clustering, Expert Systems with Applications, 2011. 38: p. 9319–9324.
[31] Alguwaizani, A., Hansen, P., Mladenović, N., Ngai, E., Variable neighborhood search for harmonic means clustering, Applied Mathematical Modelling, 2011. 35: p. 2688–2694.
[32] Hung, C.H., Chiou, H.M., Yang, W.N., Candidate clusters search for K-harmonic means data clustering, Applied Mathematical Modelling, 2013. 37: p. 10123–10128.
[33] Yang, F., Sun, T., Zhang, C., An efficient hybrid data clustering method based on K-harmonic means and Particle Swarm Optimization, Expert Systems with Applications, 2009. 36: p. 9847–9852.
[34] Güngör, Z., Ünler, A., Applying K-harmonic means clustering to the part-machine classification problem. Expert Systems with Applications, 2009. 36: p. 1179–1194.
[35] Flood, M. M., The traveling salesman problem, Operations Research, 1956. 4: p. 61-75.
[36] Althofer, I., Koschnick, K. U., On the convergence of threshold accepting, Applied Mathematics and Optimization, 1991. 24: p. 183-195.
[37] Larranaga, P., Genetic Algorithm for the Travelling Salesman Problem: A Review of Representations and Operators, Artificial Intelligence Review, 1999. 13(2): p. 129-170.
[38] Ray, S.S., Bandyopadhyay, S., Pal, S.K., New Operators of Genetic Algorithms for Traveling Salesman Problem, Proceedings of the 17th International Conference on Pattern Recognition, 2004. 2: p. 497-500.
[39] Yoichi, T., Nobuo, F., An Improved Genetic Algorithm Using the Convex Hull for Traveling Salesman Problem, IEEE Transactions on Systems, Man, and Cybernetics, 1998. 3: p. 2279~2284.
[40] Clerc M., Discrete Particle swarm optimization illustrated by the traveling salesman problem, New Optimization Techniques in Engineering Studies in Fuzziness and Soft Computing, 2004. 141: p. 219-239.
[41] Brecklinghaus, J., Hougardy, S., The approximation ratio of the greedy algorithm for the metric traveling salesman problem, Operations Research Letters, 2015. 43: p. 259–261.
[42] Karabulut, K., Tasgetiren, M.F., A variable iterated greedy algorithm for the traveling salesman problem with time windows, Information Sciences, 2014. 279: p. 383–395.
[43] Hassin, R., Keinan, A., Greedy heuristics with regret, with application to the heapest insertion algorithm for the TSP, Operations Research Letters, 2008. 36: p. 243 – 246.
[44] Yeh, W.-C., Study on Quickest Path Networks with Dependent Components and Apply to RAP, Report, NSC 97-2221-E-007-099-MY3, 2008/08/01-2011107/31.
[45] Yeh, W.-C., A Two-Stage Discrete Particle Swarm Optimization for the Problem of Multiple Multi-Level Redundancy Allocation in Series Systems, Expert Systems with Applications, 2009. 36(5): p.9192-9200.
[46] Yeh, W.-C., Optimization the End-of-Life Disassembly Sequencing Problem Based on the Use of a Novel Method based on Selt-adaptive Simplified Swarm Optimization and Precedence Preservative Operator, IEEE Transactions on Systems, Man, and Cybernetics (Part A), 2012. 42(1): p. 250-261.
[47] Yeh, W.-C., Simplified Swarm Optimization in Disassembly Sequencing Problems with Learning Effects, Computers & Operations Research, 2012. 42: p. 2168-2177.
[48] Yeh, W.-C., Wei, S.C., Simplified Swarm Optimization in Efficient Tool Assignment of Disassembly Sequencing Problem, IEEE Congress on Evolutionary Computation June 20-23, Cancún, México, 2013.
[49] Yeh, W.-C., New Parameter-Free Simplified Swarm Optimization for Artificial Neural Network Training and Its Application in the Prediction of Time Series, IEEE Transactions on Neural Networks and Learning Systems, 2013. 24(4).