SciELO - Scientific Electronic Library Online

 
vol.13 número4Control con rechazo activo a perturbaciones para procesos con retardo de transporteFactores limitantes en la implantación de MoProSoft. Revisión sistemática índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Artigo

Indicadores

  • Não possue artigos citadosCitado por SciELO

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Revista Cubana de Ciencias Informáticas

versão On-line ISSN 2227-1899

Rev cuba cienc informat vol.13 no.4 La Habana oct.-dez. 2019

 

Artículo original

Heurísticas de Construccción Voraces para el Problema de Máxima Cobertura Dinámico.

Heuristic Constructive Greedy for Dynamic Maximal Covering Location Problem.

Yadira Rosabal Alfonso1  *  , Cynthia Porras Nodarse1  , Jenny Fajardo Calderín2 

1 Universidad Tecnológica de La Habana José Antonio Echeverría, La Habana, Cuba. Correo electrónico: yadirarosabalalfonso@gmail.com, cporrasn@ceis.cujae.edu.cu

2 Universidad de Deusto, Bilbao, España. Correo electrónico: fajardo.jenny@deusto.es

Resumen

Los problemas de localización de instalaciones se aplican a una gama amplia de situaciones reales ganando un gran interés en el ámbito de la investigación. El problema de localización de máxima cobertura (MCLP) es uno de los modelos clásicos de la literatura y su objetivo es maximizar la cobertura sobre la demanda de una población con recursos limitados. De este problema se conoce la variante dinámica (DMCLP) que a diferencia del clásico maximiza la cobertura en múltiples períodos. Los métodos aproximados han sido muy utilizados en la resolución del MCLP y sus variantes. Las heurísticas constructivas, clasificadas dentro del grupo de las heurísticas, son muy conocidas por la rapidez en encontrar soluciones de buena calidad de los problemas de localización, utilizadas en ocasiones como solución inicial en los algoritmos metaheurísticos y obteniendo buenos resultados. Para dar solución al DMCLP se proponen en este trabajo cuatro heurísticas constructivas voraces. Se realiza un primer experimento para conocer la heurística más adecuada para el DMCLP. Finalmente, se comparan las mejores heurísticas con el algoritmo metaheurístico Escalador de Colinas empleado en la literatura para resolver el DMCLP.

Palabras-clave: Palabras claves: Algoritmos Metaheurísticos; Heurísticas Constructivas; Algoritmos Voraces; Problema Máxima Cobertura Dinámico

Abstract

The problems of localization of facilities are applied to a wide range of real situations winning a great interest in the environment of the investigation. The maximal covering location problem (MCLP) it is one of the classic models of the literature and their objective is to maximize the population's covering with limited resources. Of this problem the dynamic variant is known (DMCLP) that maximizes the covering in multiple periods contrary to the classic. The approximate methods have been very used in the resolution of the localization problems and their variants. The heuristic constructive classified inside the group of the heuristic ones they are very well-known for the speed and quality in finding solutions to the localization problems, used in occasions as initial solution in the algorithms metaheuristics, obtaining good results this combination. To give solution to the dynamic pattern they intend in this work four heuristic constructive greedy. they are carried out a first experiment to know which of the heuristic ones outlined it is the most appropriate in the proposed pattern. Finally, it is compared the best heuristic with the algorithm metaheuristics Hill Climbing used in the literature to solve the DMCLP.

Key words: Algorithms Metaheuristics; Heuristics Constructive; Algorithms Greedy; Dynamic Maximal Covering Location Problem

INTRODUCCIÓN

Gran cantidad de problemas de optimización pueden ser encontrados en muchas situaciones reales y modeladas como un Problema de Cobertura (CLP). De este problema se derivan dos variantes: Problema de Localización de Máxima Cobertura (MCLP por sus siglas en inglés) y el Problema de Localización de Conjuntos (SLCP por sus siglas en inglés). En el SCLP se tiene como objetivo cubrir todos los nodos de demanda con un número mínimo de instalaciones, mientras que en el MCLP el objetivo es maximizar la cobertura sobre los nodos de demanda con un número conocido de instalaciones (Zarandi et al., 2013). El MCLP fue presentado por primera vez por Church y ReVelle en 1974 (Church and ReVelle, 1974) y desde entonces se han hecho numerosas aplicaciones y teorías al modelo clásico del MCLP (Murray, 2013). En cada problema la decisión se basa en buscar la mejor combinación de las p instalaciones que tenga un máximo de cobertura sobre la demanda de los nodos. Durante muchos años el MCLP ha sido objeto de varias investigaciones aplicándose en muchas situaciones reales tales como: localización de las estaciones de la policía y patrullaje policial (Curtin et al., 2010; Guarín et al., 2015; Chawathe, 2007), distribución de redes de celulares (Tutschku, 1998; Erdemir et al., 2008), localización de antenas (Farahani et al., 2012; Wei and ReVelle, 2015), localización de estaciones de taxis (Rezazadeh et al., 2018) entre otras.

Algunas clases de problemas de optimización son relativamente fáciles de resolver, sin embargo, muchos otros tipos de problemas de optimización son muy difíciles de resolver. La idea intuitiva de problema “difícil de resolver” queda reflejada en el término científico “NP-Duro”. Un problema de optimización difícil es aquel para el que no se puede garantizar la mejor solución posible en un tiempo polinomial. La existencia de una gran cantidad y variedad de problemas difíciles, que aparecen en la práctica y que necesitan ser resueltos de forma eficiente, impulsó el desarrollo de procedimientos eficientes para encontrar buenas soluciones, aunque no fueran óptimas (Martí, 2006).

El modelo MCLP es un problema NP-Duro (Tutschku, 1998) y varios autores (Church and ReVelle, 1974; Murray and Church, 1996; Adenso-Díaz and Rodríguez, 1995), han utilizado técnicas de solución aproximadas, tanto heurísticas como metaheurísticas para resolver el problema, lográndose muy buenos resultados y alcanzado así un alto prestigio. De este problema se conoce la variante de dinamismo o multiperíodo (DMCLP por sus siglas en inglés), la cual es un modelo más real del problema, ya que considera varios períodos de tiempo y ayuda a los decisores a elaborar planes en un rango de tiempo prolongado teniendo en cuenta las variaciones de la demanda (Fajardo, 2015).

Los métodos heurísticos de construcción son procedimientos iterativos que, en cada paso, añaden un elemento hasta completar una solución. Usualmente son métodos deterministas y están basados en seleccionar en cada iteración, el elemento con mejor evaluación (Martí, 2006). Los algoritmos heurísticos de construcción son normalmente muy rápidos y devuelven soluciones razonablemente buenas. La base de los algoritmos constructivos es formada por los llamados algoritmos voraces (Greedy) (Stutzle and Ruiz, 2018). Algunos algoritmos metaheurísticos inician con una fase constructiva en la generación de solución inicial para resolver los problemas de optimización. Los métodos usados en la fase constructiva en una metaheurísticas juegan un rol importante en la eficiencia del algoritmo. Los enfoques aleatorios y voraces son dos estrategias principales usadas en la generación de soluciones iniciales (Talbi, 2009).

Los métodos heurísticos constructivos también han sido empleados para la resolver el problema MCLP (Church and ReVelle, 1974; Kuehn and Hamburger, 1963; Feldman et al., 1966). Algunos autores (Downs and Camm, 1996; Xia et al., 2019) han utilizado heurísticas constructivas basadas en algoritmos voraces en la generación de solución inicial de los algoritmos metaheurísticos para resolver el MCLP. Sin embargo, al revisar la literatura se puede apreciar que es escasa la variedad de heurísticas constructivas aplicadas al modelo dinámico del MCLP, fundamentalmente en la fase constructiva de los algoritmos metaheurísticos. La variante DMCLP es un modelo más real del MCLP y el uso de heurísticas en combinación con las metaheurísticas brindaría una mejor calidad en la solución. En este trabajo se propone una modificación de cuatro heurísticas de construcción dinámicas, basadas algoritmos voraces estáticos, para revolver el DMCLP, con el objetivo que exploren el espacio de soluciones como estrategia de inicialización en las técnicas metaheurísticas.

El artículo presente está estructurado de la siguiente manera: Sección 1 Materiales y métodos o Metodología computacional se realiza revisión de la literatura relacionada con la investigación, definición del problema DMCLP y descripción de las heurísticas voraces propuestas. En la Sección 2 Resultados y discusión realizan el marco de experimentación para conocer qué solución es la más adecuada para el modelo propuesto. Finalmente, en la Sección 3 Conclusiones se presentan las conclusiones finales que guían la línea a seguir para la continuidad de la investigación.

MATERIALES Y MÉTODOS O METODOLOGÍA COMPUTACIONAL

Revisión de la Literatura

Desde su surgimiento, diversos son los métodos aplicados al MCLP como son las heurísticas de construcción basadas en algoritmos Voraces (Greedy). Los algoritmos voraces son algoritmos simples de diseñar, sin embargo, no es fácil de entender su capacidad y limitaciones. Un algoritmo voraz está diseñado para dar soluciones de buena calidad en poco tiempo y resuelve un problema de optimización global haciendo una sucesión de decisiones localmente óptima (Dzatora and Dzatorb, 2015). En esta sección se realiza una revisión del estado del arte de las heurísticas Voraces aplicadas a las variantes del modelo MCLP.

En el año 1963 Kuehn y Hamburger (Kuehn and Hamburger, 1963) propusieron uno de las heurísticas más conocidas para los problemas de Localización de Instalaciones (FLP por sus siglas en inglés) conocido como el Procedimiento Adición (ADD). Según plantea esta heurística, la ubicación que produce la mejora máxima del valor objetivo será seleccionada para colocar una instalación, y este procedimiento se repite hasta que todas instalaciones hayan sido ubicadas. Luego en 1966, Feldman y colaboradores (Feldman et al., 1966) propusieron una heurística similar, a la anterior, pero en sentido inverso llamada Eliminar (GD) para el FLP. Este método inicia con instalaciones localizadas en todos los sitios de instalaciones potenciales y elimina aquellas instalaciones que tienen el menor impacto en la función objetivo. En 1974, Church y ReVelle (Church and ReVelle, 1974) extendieron el procedimiento ADD para el MCLP y lo renombraron como Adición Voraz (GA). La idea básica de esta heurística es seleccionar un sitio potencial que maximice la mejora de la demanda cubierta total en cada ronda recalculando la demanda no cubierta en cada iteración. El procedimiento termina cuando se abren las p instalaciones. Otra heurística desarrollada por Church y ReVelle (Church and ReVelle, 1974) fue la heurística de Adición Voraz con Sustitución (GAS), desarrollado para obtener una solución más satisfactoria que la solución obtenida por el procedimiento GA.

En 1996, Haghani (Haghani, 1996) utilizó también la técnica GA para resolver el MCLP con capacidad (CMCLP por sus siglas en inglés). Se basa en que cuando una facilidad es añadida la demanda asignada es recalculada, similar a la propuesta realizada por Church y ReVelle (Church and ReVelle, 1974). En 1998 además de los métodos básicos GA y GAS, Hochbaum y Pathria (Hochbaum and Pathria, 1998) diseñaron una heurística con enfoque voraz para el MCLP, que selecciona el subconjunto que brinda la mejora máxima en cuanto a término de la cobertura del conjunto. Los resultados obtenidos fueron muy positivos, lográndose una buena solución. Ese mismo año, Saldanha y Captivo (Saldanha and Captivo, 1998) utilizaron heurística GD y búsqueda local, implementado todo ello en CPLEX, para resolver un problema de localización discreta dinámica con establecimientos y eliminación de costos. Luego en el 2008, Indriasari (Indriasari, 2008) propuso una variante del MCLP, llamado Problema de Área de Servicio Máximo (MSAP por sus siglas en inglés). El objetivo del modelo es maximizar el área de servicio total, dado un número fijo de instalaciones. Dos algoritmos voraces: GA y Adición Voraz con Evaluación Tiempo de Viaje (GAT), fueron propuestos para solucionar el problema de MSAP. Los algoritmos GA y el GAT produjeron una mejor cobertura que la cobertura existente, teniendo el mismo número de estaciones de bomberos.

En el 2009, Ahmad Rodzi Mahmud y Vini Indriasari (Rodzi and Indriasari, 2009) propusieron dos algoritmos heurísticos voraces, la heurística GA y GAS para resolver el MSAP, desarrollado para maximizar el servicio total de un número fijo de facilidades, resultando ambas muy eficientes en dar solución al problema, sin tener mucha diferencia en sus resultados. En el 2015, Pratap (Pratap, 2015) propuso una heurística voraz llamada Gran Paso (GBS) para el modelo MCLP. Se realizó una comparación del algoritmo propuesto con un algoritmo voraz clásico y un algoritmo aleatorio, mostrando resultados satisfactorios en comparación con los dos últimos. Por último, en 2016, Bulut y Szymanski (Bulut and Szymanski, 2016) diseñaron una variante del MCLP con restricciones adicionales específicas. Adoptan la variante de despliegue de AP WIFI con la demanda máxima satisfecha usando un número definidos de AP, donde pueden ser considerados algunas restricciones adicionales. Para resolver este modelo hacen uso del GAS, el cual mostró muy buenos resultados.

Definición y Representación del DMCLP

Los problemas de localización de instalaciones se aplican a una gama muy amplia de situaciones reales y por ello cada día ganan un gran interés en el ámbito de investigación de operaciones (Murray, 2016). Entre los problemas de localización se encuentra el problema de máxima cobertura (MCLP), conocido como uno de los modelos clásicos de la literatura (Church and ReVelle, 1974).

El objetivo de su modelo era maximizar la cobertura de la población con recursos. Desde su surgimiento hasta la actualidad se han desarrollado distintas variantes del modelo clásico de MCLP. En el año 1980 Schilling (Schilling, 1980) propone una variante dinámica del modelo multiobjetivo de localización (MODL por sus siglas en inglés). Gunawardane (Gunawardane, 1982) sugiere una extensión dinámica del modelo MCLP (DMCLP) a diferencia del clásico maximiza la cobertura en múltiples períodos. Zarandi en (Zarandi et al., 2013) propone una variante dinámica donde se conoce el número de instalaciones a ubicar sin preferencia sobre el número de instalaciones que se ubicarán en cada período. A continuación, se describen los parámetros y variables del modelo del problema DMCLP (Zarandi et al., 2013).

Parámetros:

: posición y conjunto de nodos demanda.

j, J: posición y conjunto nodo de instalación.

: Índice y cantidad de periodos de tiempo.

: demanda del nodo en el periodo .

: Distancia de un nodo de demanda en el período t.

cantidad de instalaciones a ubicar teniendo en cuenta todos los periodos de tiempo T.

: Radio de cobertura de las instalaciones.

: { j| }es el conjunto de instalaciones cuya distancia al nodo es menor que .

Variables:

: variable binaria que para el valor 1 significa que la instalación j en el periodo t está abierto.

: variable binaria que para el valor 1 significa que el nodo de demanda en el periodo t está siendo cubierto, 0 en caso contrario.

La función objetivo del problema es:

Maximizar la suma de la demanda de los nodos para aquellos = 1 en todos los periodos de tiempo.

Sujeto a:

Donde la restricción (2) establece que la cantidad de instalaciones abiertas en T periodo de tiempo debe ser igual al valor p y la (3) establece que un nodo de demanda en un periodo t solo pude ser cubierto si se encuentra dentro del radio de cobertura S de al menos una instalación abierta en el mismo periodo t.

Aunque muchos enfoques diferentes pueden ser empleados para representar una solución de DMCLP, en este empleamos una representación de cadena binaria en el cual cada bit representa el estado del nodo (Zarandi et al., 2013), donde el valor 1 significa que esa instalación se encuentra abierta o disponible y 0 en caso contrario (Fajardo, 2015). La Figura 1 se muestra un esquema de representación de una solución para T períodos (Zarandi et al., 2013).

Fig.1 Representación de una solución para el problema DMCLP. El valor 1 significa cuales son las instancias abiertas en cada periodo de tiempo. En este ejemplo se tiene J = 3, T = 3 y p = 5 (Fajardo, 2015) 

Algoritmos Heurísticos y Metaheurísticos

La necesidad de encontrar soluciones buenas de problemas difíciles en un tiempo razonable es otra de las razones para usar métodos heurísticos. El método heurístico es más flexible que el método exacto, permitiendo, por ejemplo, la incorporación de condiciones que son difíciles plantear. Las heurísticas se usan como parte de un procedimiento global que garantiza la búsqueda de la solución óptima de un problema. Hay muchos métodos heurísticos que son muy diferente en la naturaleza. Por consiguiente, es difícil de proporcionar una clasificación exacta. Además, muchos de ellos se han diseñado para resolver un problema específico sin la posibilidad de generalización o aplicación a otros problemas similares (Martí and Reinel, 2011). En (Martí and Reinel, 2011) plantean los métodos constructivos como uno de las clasificaciones de métodos de resolución mediante heurísticos. Consisten en construir literalmente paso a paso una solución del problema. Usualmente son métodos deterministas y suelen estar basados en la mejor elección en cada iteración.

El término metaheurística fue introducido por Fred Glover en 1986 (Glover, 1986) quien lo define de la siguiente manera:

«Una metaheurística es referida a la estrategia maestra que guía y modifica otras heurísticas para producir soluciones más allá de aquellas que son normalmente generadas en la búsqueda de un óptimo local. Las heurísticas guiadas por una estrategia meta pueden ser un procedimiento de alto nivel o pueden contener una descripción de movimientos permitidos para trasformar una solución en otra, conjuntamente con una regla asociada de evaluación».

Heurísticas Voraces para resolver el DMCLP

A continuación, se describirán cuatro heurísticas voraces para resolver el DMCLP.

Adición Dinámica Voraz (GAD)

El procedimiento GA aplicado al problema de localización de máxima Cobertura (MCLP) fue presentado por Church y ReVelle (Church and ReVelle, 1974). La idea básica de la heurística es seleccionar un sitio potencial que maximice la mejor demanda cubierta total en cada iteración. Este proceso es continuo hasta que se hayan seleccionado p instalaciones o toda la población esté cubierta (Church and ReVelle, 1974).

El funcionamiento del algoritmo Adición Dinámica Voraz (GAD) es muy parecido a la variante estática, la única diferencia es que se deben abrir las instalaciones a lo largo del tiempo. El algoritmo comienza seleccionando el sitio potencial que cubre la mayor cantidad demanda y escoge el período donde tiene su mayor cubrimiento en la demanda, ubicando así la primera instalación en ese período. Los criterios de la segunda instalación es determinar el sitio que cubre la mayor cantidad de demanda, no cubierta por la primera instalación. Para ello reevalúa nuevamente la cobertura total de cada sitio potencial por período, no teniendo en cuenta las demandas cubiertas del período seleccionado de la primera instalación. El procedimiento se repite hasta que el conjunto de elementos seleccionados constituya una solución, es decir se haya seleccionado p instalaciones.

Adición con Sustitución Voraz Dinámico (GASD)

El procedimiento GAS aplicado al MCLP fue presentado también por Church y ReVelle (Church and ReVelle, 1974). El algoritmo GAS para la variante dinámica (GASD) comienza ubicando sus instalaciones en cada iteración del mismo modo que lo realiza el algoritmo GAD. La diferencia se establece en que las instalaciones fijas obtenidas como las mejores no son añadidas aún al conjunto de solución final, antes se verifica si las propuestas de instalaciones fijas a ubicar pueden ser reemplazadas por otros sitios potenciales. El GASD realiza la posible sustitución mediante un análisis en cuanto a la ganancia en la cobertura general, y el sitio potencial que cumpla la mayor ganancia es la que sustituye a la instalación fija, ubicando la nueva instalación en el mejor periodo de dicho sitio potencial, en caso contrario, de existir pérdida en la cobertura general, la instalación fija analizada se mantiene. El algoritmo es continuo hasta que se haya analizado todas las instalaciones fijas, ubicando finalmente las instalaciones resultantes al conjunto de solución. En otras palabras, el algoritmo GASD determina nuevas ubicaciones en cada iteración como el algoritmo GAD y además trata de mejorar la solución en cada iteración mediante la sustitución, o intercambio de instalaciones.

Eliminación Voraz Dinámica (GDD)

La heurística GD fue presentado en 1966 por Feldman, Lehrer y Ray (Feldman et al., 1966). El algoritmo GD comienza con todos los sitios potenciales instalados, y en cada iteración elimina la instalación que tenga el menor decremento en la función objetivo. La variante Eliminación Voraz Dinámica (GDD) funciona del mismo modo. A diferencia del algoritmo clásico, el GDD verifica en cada iteración que los nodos de demanda que cubría la instalación eliminada, estén cubiertos en otras instalaciones, y aquella instalación que contenga la mayor demanda cubierta es la que se le asigna los nodos de demanda de la instalación eliminada. Posteriormente se recalcula la cobertura total de cada instalación. Este procedimiento es continuo hasta que solo queden abiertas las p instalaciones. Finalmente, el algoritmo GDD ubica las p instalaciones en los mejores t periodos de tiempo de cada instalación seleccionadas como las mejores.

Gran Paso Voraz Dinámico (GBSD)

El algoritmo GBS para el problema de localización de máxima cobertura fue presentado por Pratap (Pratap, 2015). El algoritmo Gran Paso Voraz para la variante dinámica (GBSD) se le incorpora un nuevo elemento denominado período de tiempo. En cada iteración el algoritmo evalúa todas las posibles K combinaciones de sitios de instalaciones potenciales formadas del conjunto y busca la combinación de instalaciones que en su unión contengan el mayor número de elementos cubiertos. Se añade después al conjunto las K instalaciones en el período donde su sumatoria es la de mayor incremento en la demanda cubierta. Este proceso es continuo hasta que se haya seleccionado p instalaciones. Este algoritmo tiene como restricción que las K combinaciones de sitios potenciales que se formen tienen que ser menor o igual que las p a ubicar, así mimo la p tienen que ser mayor que 1. Cuando las K combinaciones de sitios potenciales a formar sea menor que p, ejemplo K=2 y p=3, el algoritmo realiza combinaciones de dos en la primera corrida y en la próxima corrida K tomaría valor 1, ubicando así las 3 instalaciones.

RESULTADOS Y DISCUSIÓN

En esta sección se realiza una descripción experimental para evaluar la efectividad de los algoritmos voraces propuestos aplicados al DMCLP. Con el fin de evaluar el rendimiento de los algoritmos y sus resultados se aplican pruebas estadísticas no paramétricas. Para verificar si existen diferencias significativas en los algoritmos propuestos y para obtener el algoritmo de mejor rendimiento se aplica test Friedman (Friedman, 1937). Para evaluar si el mejor algoritmo en cada caso tiene un rendimiento significativamente diferente a los demás, se aplicaron los test post-hoc: Holm (Holm, 1979), Finner (Finner, 1993) y Li (Li, 2008), con un nivel de significación de α = 0.05. La herramienta Keel (García et al., 2009) fue la seleccionada para realizar los test estadísticos.

Inicialmente se describen las características de las instancias utilizadas en los experimentos y la configuración de los algoritmos. Con el objetivo de evaluar el comportamiento de los algoritmos voraces, se ejecuta de manera independiente los algoritmos propuestos y se realiza un análisis comparativo para seleccionar la variante que mejor resultado obtuvo. Por último, las mejores heurísticas se comparan con el algoritmo metaheurístico Escalador de Colinas con solución inicial por ruleta y con solución inicial Voraz.

Descripción de la experimentación

Se seleccionaron tres instancias de la librería TSP-Lib (Reinelt, 1991): ar6723, pm4951 y fi10639. De ellas se seleccionó aleatoriamente un conjunto de puntos para los nodos de demanda con una 𝐼 de 1800 a 2500 puntos y para las instalaciones una 𝐽 de 100 a 200. Se realizaron diferentes combinaciones de radios de coberturas y cantidad de instalaciones a ubicar. En la Figura 2 se muestra una visualización de las instancias originales de TSP-Lib, de las cuales se crearon las instancias para la experimentación.

Fig.2 Representación de las instancias originales 

La p representa el 30, 40, 50, 60 y 70 % de J. En la Tabla 1 se muestra la descripción de los diferentes elementos de las instancias. Para la experimentación con el algoritmo Escalador de Colinas se realizaron 30 ejecuciones y 10000 iteraciones, para todas las instancias.

Tabla 1 Descripción de las instancias para la experimentación con la variante DMCLP. 

Análisis y comparación de las Heurísticas Voraces de manera independiente

A continuación, se explican los resultados de comparar las heurísticas propuestas entre ellas. En los resultados del test de Friedman, ilustrado en la Figura 3 el p-valor en todas las instancias y de forma global fue igual a 0, por lo que es menor que 0.05 y se rechaza la hipótesis nula. Entre las heurísticas sí existen diferencias significativas. Según el ranking que ofrece el test de Friedman, el algoritmo GASD fue el mejor algoritmo en las instancias fi1800, ar2000 y pm2500 siendo así el mejor algoritmo de manera global.

Fig.3 Gráfica de Friedman. Comportamiento de las Heurísticas Voraces para las distintas instancias 

En la Tabla 2. Resultado de la Prueba no paramétrica de Friedman post-hoc Holm, Finner y Li. para las distintas Instancias se muestran los resultados de los post-hoc Holm, Finner y Li. El signo “>” significa que el algoritmo A se comportó mejor que el B, y “-” quiere decir que no hay diferencias significativas, o sea que no hay evidencias suficientes para rechazar la hipótesis nula. Se compara el algoritmo que mejor ranking obtuvo de forma global en la prueba de Friedman, el GASD, con los restantes algoritmos.

Tabla 2 Resultado de la Prueba no paramétrica de Friedman post-hoc Holm, Finner y Li. para las distintas Instancias 

El resultado obtenido en la Tabla 2 muestra que el GASD se comportó mejor que el GDD y GBSD de forma global ya que el p-valor es menor que 0.05, por tanto, se rechaza la hipótesis nula y sí existen diferencias significativas entre los comportamientos de los algoritmos. Hay que destacar que el algoritmo GASD en esta prueba tuvo un comportamiento similar al algoritmo GAD en todas las instancias y de forma global donde no hubo diferencias significativas.

Análisis y comparación del Escalador de Colinas con Solución Inicial por Ruleta y Heurísticas Voraces

En esta sección se realiza un análisis con el Escalador de Colinas (Bagherinejad et al., 2017) utilizando en su solución inicial algoritmos voraces y ruleta propuesto en (Zarandi et al., 2013). Los algoritmos voraces seleccionados fueron lo que obtuvieron mejores resultados en el análisis anterior, en este caso el GASD y el GAD. Se pretende con este experimento demostrar la calidad de la solución de las heurísticas voraces sobre el problema propuesto.

En la Figura 4 se muestra una gráfica de los resultados obtenidos al aplicar el Test de Friedman. El p-valor en todas las instancias y de forma global fue igual a 0, por lo que es menor que 0.05 y se rechaza la hipótesis nula. Ello implica que sí existen diferencias significativas entre los comportamientos de los algoritmos. Según el ranking que ofrece el test de Friedman, la combinación EC_GASD fue el mejor algoritmo en las instancias fi1800, ar2000 y EC_GAD fue mejor en la instancia pm2500. De forma global el algoritmo EC_GASD fue el de mejor rendimiento.

Fig.4 Gráfica de Friedman. Comportamiento de las Heurísticas Voraces para las distintas instancias 

En la Tabla 3 se muestran los resultados de los post-hoc Holm, Finner y Li para todas las instancias, donde se compara el algoritmo que mejor ranking obtuvo de forma global en la prueba de Friedman, el EC_GASD, con los restantes algoritmos. Se puede concluir que no hay diferencias significativas entre el algoritmo EC_GASD y EC_GAD ya que tuvieron un comportamiento similar por lo que no existen evidencia para rechazar la hipótesis nula. El algoritmo EC_GASD se comportó mejor que el EC_Ruleta.

Tabla 3 Resultado de la Prueba no paramétrica de Friedman post-hoc Holm, Finner y Li. para las distintas Instancias 

CONCLUSIONES

Los algoritmos voraces han sido aplicados a una gama amplia de problemas de optimización entre los que se encuentra el MCLP, dado a la rapidez y calidad en la solución que brindan. La variante dinámica del MCLP es un modelo más real del clásico, ya que al considerar varios períodos de tiempo se pueden tomar decisiones en un rango de tiempo prolongado teniendo en cuenta fluctuaciones de la demanda. Los algoritmos voraces propuestos en esta investigación para la variante dinámica obtuvieron resultados satisfactorios. De los algoritmos voraces propuestos, el que mejor comportamiento obtuvo de manera global fue el algoritmo GASD, aunque no hubo diferencias significativas con el algoritmo GAD. En el último análisis comparativo, los algoritmos GAD y GASD como solución inicial del Escalador de Colinas brindaron soluciones de mejor calidad en comparación con la solución inicial utilizando ruleta dentro del algoritmo metaheurístico Escalador de Colinas. Las heurísticas voraces utilizadas como estrategia de inicialización en los algoritmos metaheurísticos obtienen buenos resultados para resolver el modelo DMCLP.

REFERENCIAS

Zarandi, M. H. F., et al. "The large-scale dynamic maximal covering location problem". Mathematical Computer Modelling 2013, 57 (3-4): p. 710-719. [ Links ]

Church, R. and C. ReVelle. "The maximal covering location problem". Papers of the Regional Science Association, 1974, 32(21): p. 101-118. [ Links ]

Murray, A. T. "Maximal Coverage Location Problem: Impacts, Significance, and Evolution". Internacional Regional Science Review, 2016, 39(1): p. 5-27. [ Links ]

Curtin, K., et al. " Determining Optimal Police Patrol Areas with Maximal Covering and Backup Covering Location Models". Netw Spat Econ, 2010, 10(1): p. 125-145. [ Links ]

Guarín, A., et al. "Fast reaction police units in Medellín: A budget-constrained maximal homicide covering location approach". Colombia, Banco de la República, 2015. [ Links ]

Chawathe, S. S. "Organizing Hot-Spot Police Patrol Routes". Intelligence and Security Informatics, 2007: p. 79-86. [ Links ]

Tutschku, K. "Demand-based Radio Network Planing of Cellular Mobile Communication System". Seventeenth Annual Joint Conference of the IEEE Computer and Communications Societies Proceedings, 1998: p. 1054-1061. [ Links ]

Erdemir, E. T., et al. " Location coverage models with demand originating from nodes and paths: application to cellular network design". European Journal of Operational Research, 2008, 190(3): p. 610-632. [ Links ]

Farahani, R. Z., et al. "Covering problems in facility location: A review". Computers & Industrial Engineering, 2012, 62: p. 368-407. [ Links ]

Wei, R. and A. T. Murray "Continuous space maximal coverage: insights, advances and challenges". Computers & Operations Research, 2015, 62: p. 325-336. [ Links ]

Rezazadeh, H., et al. " Robust cooperative maximal covering location problem: a case study of the locating Tele-Taxi stations in Tabriz, Iran". International Journal of Services and Operations Management, 2018, 29(2): p. 163-183. [ Links ]

Martí, R. "Procedimientos Metaheurísticos en Optimización Combinatoria". Matemátiques, 2006, 60: p. 3-62. [ Links ]

Murray, A. T. and R. L. Church. " Applying Simulated Annealing to Location-planning Models". Journal of Heuristics, 1996, 2: p. 31-53. [ Links ]

Adenso-Díaz, B. and F. Rodríguez "A Simple Search Heuristic for the MCLP: Application to the Location of Ambulance Bases in a Rural Region". Omega, 1995, 25(2): p. 187-187. [ Links ]

Fajardo, J. "Soft Computing en Problemas de Optimización Dinámicos". Tesis de Doctorado, Universidad de Granada, España, 2015. [ Links ]

Stutzle, T. and R. Ruiz. "Iterated Greedy", IRIDIA-Technical Report Series ISSN 1781-3794, Belgium, 2018. 36p. [ Links ]

Talbi, E. G. "Metahueristics: from design to implementation". John Wiley and Sons, Inc, 2009, 74: p. 593. [ Links ]

Kuehn, A. A., and Hamburger, M. J "A heuristic program for locating warehouses". Management Science, 1963, 9(4): p. 643-666. [ Links ]

Feldman, E., et al. "Warehouse location under continuous economics of scale." Management Sciences, 1966, 12(9): p. 670-684. [ Links ]

Downs, B. T. and J. D. Camm "An exact algorithm for the maximal covering problem." Naval Research Logistics Quarterly, 1996, 43(3): p. 435-461. [ Links ]

Xia, L., et al. "An empirical comparison of five efficient heuristics for maximal covering location problems". IEEE/INFORMS International Conference on Service Operations, Logistics and Informatics, 2019: p. 747-753. [ Links ]

Dzatora, M. and J. Dzatorb. "An Efficient Modified Greedy Algorithm for the P-Median Problem." Proceedings of the 21st International Congress on Modelling and Simulation. Australia: Gold Coast, 2015: p. 1855-1861. [ Links ]

Haghani, A. "Capacitated maximum covering location models: Formulations and solution procedures". Journal of Advanced Transportation, 1996, 30(3): p. 101-136. [ Links ]

Hochbaum, D. S. and A. Pathria. "Analysis of the greedy approach in problems of maximum k-coverage". Naval Research Logistics Quarterly, 1998, 45(6): p. 615-627. [ Links ]

Saldanha, F. and M. E. Captivo "Heuristic approach for the discrete dynamic location problem". Location Science, 1998, 6(1-4): p. 1211-1223. [ Links ]

Indriasari, V. "Integration of Travel Time Zone for Optimal Siting of Emergency Facilities”. Thesis Masters, University Putra Malaysia, 2008. [ Links ]

Rodsi, M. and V. Indriasari "Facility Location Models Development to Maximize Total Service Area”. Theoretical and Empirical Researches in Urban Management, Abril 2009: p. 2087-2100. [ Links ]

Pratap, D. "Big Step Greedy Heuristic for Maximum Coverage Problem". International Journal of Computer Applications (0975 - 8887), 2015, 125(7): p. 19-24. [ Links ]

Bulut, E. and B. Szymanski "Rethinking Offloading WiFi Access Point Deployment from User Perspective". . Workshop on Smart Environments & Urban Networking, (SEUNet) in Proc. IEEE WiMob, 2016: p. 59-65. [ Links ]

D. Schilling "Dynamic Location Modeling for Public-Sector Facilities: A Multicriteria Approach", D. S., 1980,11(4): p. 714-724. [ Links ]

Gunawardane, G. "Dynamic versions of set covering type public facility location problems". European Journal of Operational Research, 1982, 10(2): p. 190-195. [ Links ]

Marti, R. and G. Reinel. "The Linear Ordering Problem Exact and Heuristics Method in Combinatorial Optimization". Springer, 2011, 12: p. 172. [ Links ]

Glover, F. " Future paths for integer programming and links to artificial intelligence". Computers & Operations Research, 1986, 13: p. 533-549. [ Links ]

Friedman, M. "The use of ranks to avoid the assumption of normality implicit in the analysis of variance". Journal of the American Statistical Association, 1937, 32(200): p. 675-701. [ Links ]

Holm, S. "A simple sequentially rejective multiple test procedure". Scandinavian Journal of Statistics, 1979, 6: p. 65-70. [ Links ]

Finner, H. "On a monotonicity problem in step-down multiple test procedures". Journal of the American Statistical Association, 1993, 88: p. 920-923. [ Links ]

Li, J. "A two-step rejection procedure for testing multiple hypotheses". Journal of Statistical Planning and Inference, 2008, 138: p. 1521-1527. [ Links ]

García, S., et al. "A study on the use of nonparametric tests for analyzing the evolutionary algorithms behaviour: a case study on the cec2005 special session on real parameter optimization". Journal of Heuristics, 2009, 16(6): p. 617. [ Links ]

Reinelt, G. "TSPLIB.A traveling salesman problem library". ORSA Journal on Computing, 1991, 3(4): p. 376-384. [ Links ]

Bagherinejad, J., et al. "Developing dynamic maximal covering location problem considering capacitated facilities and solving it using hill climbing and genetic algorithm". Engineering Review, 2017, 37 (2): p:178-193. [ Links ]

Recibido: 25 de Septiembre de 2019; Aprobado: 26 de Noviembre de 2019

* Autor para correspondencia: yadirarosabalalfonso@gmail.com

Creative Commons License