Guillermo González Vargas,¹ Felipe González Aristizábal²

¹ Ingeniero industrial, Universidad Nacional de Colombia. Aspirante a M.Sc. en Ingeniería Industrial, Universidad de los Andes, Colombia. gu-gonza@uniandes.edu.co
² Ingeniero industrial, Universidad Nacional de Colombia. Aspirante a M.Sc. en Ingeniería Industrial, Universidad de los Andes, Colombia. fe-gonza@uniandes.edu.co

RESUMEN

Palabras clave: distribución, ruteo de vehículos, formulación matemática.

ABSTRACT

This paper deals with VRP (vehicle routing problem) mathematical formulation and presents some methodologies used by different authors to solve VRP variation. This paper is presented as the springboard for introducing future papers about a manufacturing company’s location decisions based on the total distance traveled to distribute its product.

Keywords: distribution, vehicle routing problem, mathematical formulation.

Recibido: abril 18 de 2006
Aceptado: septiembre 3 de 2006

Introducción

Este artículo es el primero de una serie de tres, en los cuales se ilustrará la aplicación de técnicas metaheurísticas en la solución de un problema de ruteo de vehículos (VRP, por las siglas en inglés de Vehicle Routing Problem). Esta serie de artículos tiene su origen en un caso de aplicación desarrollado para una empresa manufacturera colombiana, la cual ha solicitado expresamente a los autores mantener en anonimato su identidad; por tanto, se ilustrará la metodología y los resultados obtenidos de la manera más abstracta posible, manteniendo los datos de entrada reales con los cuales se aplicaron las diferentes técnicas y la calidad y rigurosidad académica que merecen.

El artículo tiene como fin presentar al lector el problema que se abordará en los dos artículos siguientes, iniciando con una revisión bibliográfica en la que se ilustra la problemática del ruteo de vehículos y las técnicas utilizadas para resolverla, seguida por una breve descripción del problema; posteriormente se enuncia la formulación matemática del mismo y se finaliza con las conclusiones pertinentes.

Centrados en el problema de distribución, en el que se enmarca el presente artículo, es importante recurrir a la afirmación de Toth y Vigo (2000): “El problema de distribuir productos desde ciertos depósitos a sus usuarios finales juega un papel central en la gestión de algunos sistemas logísticos, y su adecuada planificación puede significar considerables ahorros. Esos potenciales ahorros justifican en gran medida la utilización de técnicas de investigación operativa como facilitadoras de la planificación, dado que se estima que los costos del transporte representan entre el 10% y el 20% del costo final de los bienes”. Dentro de este problema de transporte es necesario determinar el tipo de recurso a utilizar, la cantidad y las rutas a seguir, lo que se denomina problema de ruteo, y es tratado en la literatura como el problema del agente viajero (TSP, por las siglas en inglés de Traveling Salesman Problem), o en términos generales, para problemas con capacidad definida (Machado et al., 2002), es generalizado el VRP (Olivera, 2004).

El ruteo de vehículos (VRP) es un problema de optimización combinatoria complejo, considerado ya un paradigma en la literatura especializada (Hermosilla y Barán, s/f), que surgió, según Olivera (2004), desde 1959. Este tipo de situación, como se había mencionado anteriormente, es una generalización del problema del agente viajero, el mismo que puede ser explicado de la siguiente manera:

Existe un agente de ventas que debe visitar a sus clientes ubicados en diferentes ciudades y luego volver a su ciudad de partida, y dicha actividad debe ser llevada a cabo con el menor costo posible (Ahuja et al., 1993); según Hermosilla y Barán (s/f) el costo de la ruta puede estar dado por la duración total de la misma (en tiempo o distancia). El problema de ruteo de vehículos se representa en un grafo con nodos y arcos, los cuales representan la ubicación de los clientes y la red vial por la cual pueden circular los vehículos.

Una recopilación de técnicas exactas de solución existentes para los problemas de ruteo de vehículos puede encontrarse en Laporte (1992); no obstante los de gran dimensión resultan imposibles de solucionar en tiempo polinomial, por lo que el VRP se denomina NP-hard (Machado et al., 2000; Olivera, 2004), donde no es posible alcanzar una solución óptima, y, dependiendo de las características especiales de clientes, locaciones y producto/servicio, requiere la elaboración de una metodología de solución específica con la cual sea posible aproximarse lo mejor posible al óptimo. Las diferentes variaciones y restricciones del problema generan una “familia” de VRP (Medaglia, 2005) de la que vale la pena mencionar ocho casos típicos, los cuales al compartir características pueden dar lugar a todo un universo de problemas VRP. Los principales problemas de ruteo de vehículos se ilustran en la Figura 1 y pueden ser descritos así:

CVRP

(Capacited VRP), es el VRP más general y consiste en uno o varios vehículos con capacidad limitada y constante encargados de distribuir los productos según la demanda de los clientes (Olivera, 2004; Lee et al., 2002). Este problema ha sido resuelto mediante búsqueda Tabú (Olivera, 2004; Rego, s/f), algoritmos genéticos (Machado et al, 2002; Machado et al., 2003 (a); Olivera, 2004), algoritmos de colonias de hormigas (Olivera, 2004), Constraint programming (Shaw, 1998) y algoritmos híbridos de recocido simulado y algoritmos genéticos (Wendt y König, s/f).

]]>

MDVRP

(Multi-Depot VRP), o VRP con múltiples depósitos es un caso de ruteo de vehículos en el que existen varios depósitos (cada uno con una flota de vehículos independiente) que deben servir a todos los clientes, caso resuelto por Tansini et al., (s/f) mediante técnicas de cluster first – routen second, que serán descritas posteriormente.

PVRP

(Period VRP), contempla en su planteamiento un horizonte de operación de M días, periodo durante el cual cada cliente debe ser visitado una vez, problema propuesto por Francis et al., (2004) y resuelto por los mismos autores mediante relajación lagrangiana.

SDVRP

(Split Delivery VRP), o VRP de entrega dividida, donde se permite que un cliente pueda ser atendido por varios vehículos si el costo total se reduce, lo cual es importante si el tamaño de los pedidos excede la capacidad de un vehículo, (Lee et al., 2002; Archetti et al., 2001), resuelto en 2002 por Lee et al.

SVRP

(Stochastic VRP), se trata de un VRP en que uno o varios componentes son aleatorios; clientes, demandas y tiempos estocásticos son las principales inclusiones en este tipo de problemas. El SVRP ha sido resuelto por Bianchi et al., (s/f) a través de búsqueda Tabú, recocido simulado, algoritmos de colonias de hormigas, algoritmos genéticos y otros algoritmos evolutivos.

(VRP Pickup and Delivery), o VRP con entrega y recogida, es aquel en el que cabe la posibilidad de que los clientes pueden devolver determinados bienes, por tanto, se debe tener presente que estos quepan en el vehículo. Esta restricción hace más difícil el problema de planificación y puede causar una mala utilización de las capacidades de los vehículos, un aumento de las distancias recorridas o a un mayor número de vehículos (Volkan, 2005; Dethloff, 200; Halse, 1992; Gendreau et.al., 1994; Min, 1989). Una forma de solucionar el VRPPD mediante la utilización de algoritmos genéticos fue propuesta por Volkan en 2005, quien afirma que si este problema incluye la restricción de culminar todas las entregas antes de iniciar las recogidas se da lugar a un VRP con backhauls o VRPB, variación del VRP estudiada por Charlotte y Goetschalckx (1998).

MFVRP

(Mix Fleet VRP), es un VRP en el que se suponen vehículos con distintas capacidades o capacidad heterogénea, por lo que es necesario considerar estas capacidades en la ruta que seguirá cada recurso, ya que un camión más grande podrá realizar una ruta más larga o que tenga mayor concentración de demanda, lo cual fue estudiado inicialmente por Liu y Shen (1999) y posteriormente resuelto por Barchett y Campion mediante Búsqueda Tabú en 2002.

VRPTW

(VRP with Time Windows), es aquel en el que se incluye una restricción adicional en la que se asocia a cada cliente una ventana de tiempo, es decir, cada cliente sólo está dispuesto a recibir el bien o servicio durante un intervalo de tiempo predeterminado; este tipo de problema ha sido resuelto por diferentes autores, entre los que vale la pena mencionar a Olivera (2004), quien presenta una solución mediante búsqueda Tabú, Gendreau et al., (1998) proponen una heurística de inserción; Olivera (2004), Vacic (2002), Bräysy (2001), Zhu (2000) y Louis et al (1999) lo resuelven con algoritmos genéticos y Barán y Schaerer (2003) y Gambardella et al., (1999) presentan una propuesta a través de algoritmos de colonia de Hormigas.

Los diferentes problemas VRP, y básicamente los que utilizan múltiples vehículos (Restori, s/f; Olivera, 2004) y/o depósitos (Tansini et al., s/f), pueden reducir su complejidad acotando el universo de soluciones, disminuyendo el conjunto de clientes a ser visitados por cada vehículo o desde cada depósito, esto es, asignar a cada vehículo/depósito un conjunto de clientes para atender, lo que Medaglia (2005) llama set covering, o lo que otros autores conocen como clusterizar o asignar primero, rutear después cluster first – routen second (Olivera, 2004).

La clusterización de los clientes puede ser realizada a través de diferentes heurísticas, entre las que vale la pena mencionar:

Heurística de barrido o sweep, esta técnica propone establecer un punto de origen en el depósito y desde allí realizar un barrido para abarcar toda el área geográfica del problema, determinando así cada uno de los clusters (anónimo, s/f).³

Heurística de asignación generalizada de Fisher y Jaikumur, basa la generación de clusters en la solución de un problema de asignación generalizada (GAP) sobre los clientes, fue realizada por Fisher y Jaikumur (1981).

Una vez definido el conjunto de clientes a atender (cluster) se procede a realizar la asignación de la mejor ruta, este subproblema generalmente subyace en un caso de problema de agente viajero (TSP) y puede resultar tan pequeño como para ser solucionado mediante una técnica de optimización como la programación lineal (Ahuja et al., 1993).

Descripción del problema

La problemática tratada en este y los próximos artículos se basa en la necesidad que presenta una empresa manufacturera para decidir la localización de una bodega desde la cual sea posible distribuir su producto a 53 centros de consumo (en adelante se numerarán consecutivamente de 1 a 53), cada uno de los cuales tiene una demanda periódica constante,⁴ como se ilustra en la Tabla 1.

La empresa cuenta con máximo seis vehículos con capacidad constante y homogénea de 5.500 unidades, con los cuales debe entregar en cada periodo la totalidad de productos que se demandan (cada vehículo realiza un recorrido por periodo). El punto de partida de los vehículos es uno de tres centros de consumo que por sus características resultan opcionados para convertirse en bodegas de distribución (9, 28 y 49). La demanda del centro de consumo que se convierte en bodega se satisface in situ, por lo que no requiere desplazamiento ni consumo de la capacidad de los vehículos. La compañía realiza la distribución durante las 24 horas del día y se considera que el tiempo de entrega es despreciable.

Para el caso práctico del presente y los subsiguientes artículos se ilustrará la zona de influencia de la empresa con un grafo, en el cual los nodos representan los centros de consumo y los arcos las rutas directas existentes entre pares de centros (para el lector interesado se presenta la longitud de los arcos en el apéndice: “Matriz de distancias entre centros de consumo”, tabla en la cual se resaltan las distancias entre los arcos conectados de manera directa). El grafo de la Figura 2 ilustra la zona de influencia de la empresa y la ubicación geográfica de los centros de consumo, además se resaltan las tres posibles ubicaciones de la bodega de distribución.

Formulación del problema

Teniendo en cuenta la descripción del problema realizada anteriormente, es posible determinar que el problema planteado es un problema clásico de localización de una sola instalación (Schroeder, 1996) teniendo en cuenta un solo criterio de tipo tangible (Guerrero y Osorio, 2003) que puede ser denominado: “minimización de la distancia total a recorrer para la satisfacción de la demanda”⁵. Sabiendo que debido a criterios intangibles (Ballou, 1999) el universo de posibilidades para la localización de la bodega de la empresa queda reducida a tres posibilidades (centros de consumo 9, 28 y 49), el problema puede ser solucionado como tres subproblemas independientes de ruteo con múltiples vehículos, lo que Thomson y Orlin (1989) llaman multi-vehicle routing and scheduling, o lo que se ha descrito en este artículo como ruteo de vehículos con capacidad limitada CVRP (Olivera, 2004; García, 2000). La propuesta final de localización se ha de realizar teniendo en cuenta el mínimo costo de la ruta para cada una de las tres posibles localizaciones. Para este caso específico, y con base en lo dicho por Hermosilla y Barán (s/f), se asume que el costo de la ruta está dado por la longitud total de la misma, y se supone velocidad constante y unitaria para los vehículos.

Es importante aceptar que un problema similar al planteado se presentó en el concurso Whizzkids '96 (Applegate et al., 2001) y fue resuelto por Jan Karel Lenstra y Emile Aarts et al. (Hurkens, 1997) mediante la utilización de Recocido Simulado (Simulated Anealing) y que el mismo problema ha sido abordado por otros autores como Thomson y Orlin (1989) y Machado et al., (2003, a; 2003, b), entre otros.

Teniendo en cuenta que el problema se abordará como tres subproblemas de ruteo de vehículos, uno para cada posible localización de la bodega, en cada uno de los cuales se encontrará la menor distancia total a recorrer para atender la demanda de la empresa, el problema global puede ser definido como la minimización del mínimo de las distancias totales a recorrer para satisfacer la demanda de la zona de influencia de la empresa, partiendo desde cada uno de los tres centros de consumo opcionados (expresión 1).

Donde,

DT_i : Distancia total recorrida para satisfacer la demanda al localizar la bodega en el centro de consumo i (i= 9,28,49).

Formulación matemática de los subproblemas: VRP

Dado un número máximo de seis (6) vehículos con una capacidad homogénea determinada u=5500, que tienen como centro de operaciones único la bodega en el nodo 0 (centros de consumo 9, 28 y 49), y deben satisfacer a una cantidad de 52 clientes (centros de consumo), que se representan por j= 1,2,…,52, cada uno de los cuales tiene una demanda conocida d_j (véase Tabla 1), es posible realizar la siguiente formulación matemática:

Índices

Los índices del modelo son:

j = nodo de llegada j (1,2,…,52)

k = vehículo k (1,2,…, K)

Variables

Las variables que se definen en el modelo son:⁶

K = Número de recursos (vehículos) a utilizar.

Parámetros

Los parámetros del problema son:

d_i= Demanda en el nodo j

u= Capacidad del recurso k

n= Cantidad de clientes

Modelo matemático

El modelo matemático que representa cada uno de los tres subproblemas de ruteo se puede plantear según los aportes de Ahuja et al. (1993) y Olivera (2004) como sigue:⁷

Sujeto A:

]]>

]]>

El conjunto A se define como: A={(i,j) : y_ij=1}.

La variable y_ij presente en las restricciones (4) y (5) indica la activación del arco (i,j), lo que determina un recorrido entre los nodos i,j, además se asegura que todo cliente es un nodo intermedio de alguna ruta. Los grupos de restricciones (6) y (7) indican que k es la cantidad de vehículos utilizados en la solución y que todos los que parten del depósito deben regresar al mismo. La restricción (8) observa que cada vehículo no sobrepase su capacidad. La restricción (9) vigila que la solución no contenga ciclos usando los nodos 1,2,…n. De otra manera los arcos de A contendrían algún ciclo pasando a través de un conjunto de nodos Q y la solución violaría la restricción, porque el lado izquierdo de la restricción sería al menos |Q|. La restricción (10) limita el número máximo de vehículos a utilizar hasta una cantidad máxima. Las restricciones (10) y (11) indican que tanto la variable x^k_ij como la variable y_ij son binarias.

Consideraciones finales

Una vez definida la formulación matemática con la que se describe el problema al que se enfrenta el presente artículo y teniendo en cuenta que los parámetros d_j, u y n son conocidos, sólo resta determinar un input fundamental del CVRP, el costo de cada ruta (c_ij), que para este caso será la distancia de cada ruta.

Según lo expresado por Olivera (2004), puede suponerse que el grafo es completo, pues entre todo par de lugares de una red de transporte razonable, debería existir algún camino, concepto que necesariamente debe ser usado en el desarrollo del presente artículo; así, aunque existan nodos que sólo tienen una vía de acceso (véase por ejemplo nodo 50 en la Figura 2), es necesario suponer que entre cada par de nodos existe un arco que los une, por lo anterior los autores decidieron construir arcos ficticios que representen la conexión existente entre todo par de nodos del grafo, a través de la ruta más corta entre ellos.

La construcción de los arcos ficticios anteriormente mencionados supone la solución de problemas de ruta más corta (Shortest Path Problems), específicamente lo que Ahuja et al. (1993) llaman All-pair shortest path problem y que pretende encontrarla entre todo par de nodos de un grafo, problema que puede ser solucionado según los mismos autores mediante la aplicación de repeated shortest path algorithm, que consiste en la aplicación de un algoritmo para encontrar rutas más cortas para un solo origen n veces (donde n es el número de nodos), o bien, all-pairs lebel-correcting algorithm, que se basa en la aplicación del algoritmo de Floyd-Warshall. Los dos métodos mencionados pretenden alcanzar la condición de optimalidad representada por la expresión 13, donde d[i,j] representa la distancia más corta entre los nodos i, j.

Para crear la matriz de distancias entre todo par de centros de distribución del apéndice Matriz de distancias entre centros de consumo se decidió calcular la ruta más corta para cada municipio respecto de los demás (repeated shortest path algorithm) mediante la aplicación del algoritmo de Dijkstra (Taha, 1997).

Una vez hecho el recorrido bibliográfico pertinente, descrito y formulado el problema y recopilada la información necesaria para enfrentar el problema, los autores utilizaron diferentes técnicas aplicadas a encontrar la asignación que reduzca al máximo la distancia total recorrida por los vehículos encargados de satisfacer la demanda de los 53 centros de consumo, y con base en los resultados decidir la mejor localización para la empresa. Los resultados numéricos de la aplicación de cada una de las técnicas, al igual que la descripción de las mismas, se encontrarán en los dos artículos siguientes.

El de ruteo de vehículos es la generalización del problema del agente viajero y encierra una familia de que debe ser resuelta según las características específicas de cada caso.

La formulación matemática del ruteo de vehículos debe contener familias de restricciones que imposibiliten la construcción de ciclos o subtoures.

En los siguientes artículos los autores expondrán diferentes metodologías de solución para el caso de estudio ilustrado.

Bibliografía

AHUJA, R.; MAGNANTI, T. and ORLIN, J. Network flows: theory, algorithms, and applications, Englewood Cliffs, New Jersey: Prentice Hall, 1993.        [ Links ]

ANÓNIMO, Network. s/f., Disponible en: http://www.cs.tcd.ie/courses/baict/bass/4ict5/Networks2004.pdf. Consultado en: Marzo de 2004.        [ Links ]

APPLEGATE, D.; COOK, W.; DASH, S. and ROHE, A. Solution of a Min-Max Vehicle Routing Problem. August 15, 2001. Disponible en: http://www.research.ibm.com/ people/s/sanjeebd/reports/newsboys.pdf. Consultado en: Febrero de 2005.        [ Links ]

ARCHETTI, C.; MANSINI, R. and SPERANZA, M.G. The Split Delivery Vehicle Routing Problem with Small Capacity, November 7, 2001, Disponible en: http://fausto.eco.unibs.it/~www_mequ/ricerca/quaderni/201.pdf. Consultado en: Febrero de 2005.        [ Links ]

BALLOU, R. H. Business Logistics Management, Prentice-Hall, New Jersey, 4° edición, 1998, pp. 3-28, 481-592.         [ Links ]

BARÁN, B. and SCHAERER, M. A. Multiobjective Ant Colony System for Vehicle Routing Problem with Time Windows. Proceedings of the 21st IASTED International Conference APPLIED INFORMATICS, Innsbruck, Austria, 2003.        [ Links ]

BARCHETT D. and Campion, E., Mix fleet vehicle routing problem – An application of Tabu search in the grocery delivery industry; 2002, Disponible en: http://www.mang.canterbury.ac.nz/courseinfo/msci/msci480/WebPage.htm, Consultado en: Febrero de 2005.        [ Links ]

BIANCHI, L.; METAHEURISTICS for the Vehicle Routing Problem with Stochastic Demands; S/F., Disponible en: http://www.idsia.ch/idsiareport/IDSIA-06-04.pdf., Consultado en: Febrero de 2005.        [ Links ]

BRAMEL, J. and SIMCHI-LEVI, D. A location based heuristic for general routing problems, Operations Research, Vol. 43, 1995, pp. 649–660.        [ Links ]

BRÄYSY, O. Genetic Algorithms for the Vehicle Routing Problem with Time Windows. 2001, Disponible en: http://osiris.tuwien.ac.at/~wgarn/VehicleRouting/Braysy.pdf, Consultado en: Febrero de 2005.        [ Links ]

CHARLOTTE, J. and Goetschalckx, M, The vehicle routing problem with backhauls: properties and solution algorithms. Scholl of industrial and systems engineering - georgia institute of technology., Copyright 1992 – 1998. Disponible en: http://www.isye.gatech.edu/people/faculty/Marc_Goetschalckx/cali/Lineback/Vehicle%20Routing%20Problem%20with%20Backhauls.pdf, Consultado en Febrero de 2005.        [ Links ]

DETHLOFF, J., Relation between vehicle routing problems: an insertion heuristic for the vehicle routing problem with simultaneous delivery and pick-up applied to the vehicle routing problem with backhauls. En: Journal of the Operational Research Society, 2002, Consultado en: Febrero de 2005.        [ Links ]

FISHER, M. and JAIKUMUR, R. A generalized assignment heuristic for the vehicle routing problem., Networks 11, 1981, pp.109–124.        [ Links ]

FRANCIS, P. SMILOWITZ, K. and TZUR, M, The period vehicle routing problem with service choice; 2004., Disponible en: http://www.iems.northwestern.edu/ images/PDF/WP_04_005.pdf, Consultado en: Febrero 2005.        [ Links ]

GAMBARDELLA, L.M.; TAILLARD, E. and AGAZZI, G., MACS-VRPTW: A Multiple Ant Colony System for Vehicle Routing Problems with Time Windows, New Ideas in Optimization, Disponible en: http://www.idsia.ch. Consultado en: Febrero 2005, McGraw-Hill, 1999.        [ Links ]

GARCÍA, A. Optimización de rutas, seguridad en el transporte y sistemas GIS. 2000., Disponible en: http://www.imac.unavarra.es/SEMAOL/Ponencias/04_AlejandroG_del_Valle.pdf, Consultado en: Febrero de 2005.        [ Links ]

GENDREAU, M.; Hertz, A. and LAPORTE, G. A tabu search heuristic for the vehicle routing problem., Management Science 40, 1994, pp. 1276–1290.        [ Links ]

GENDREAU, M.; HERTZ, A., LAPORTE, G. and STAN, M.., A generalized insertion heuristic for the traveling salesman problem with time windows, Operations Research, Vol. 43, 1998, pp. 330–335.         [ Links ]

GUERRERO, L.; y OSORIO, L. La localización de instalaciones como decisión estratégica de la logística: un acercamiento al estado del arte., Trabajo de grado ingeniería industrial, Universidad Nacional de Colombia, Manizales, Marzo, 2003.        [ Links ]

HALSE, K. Modeling and Solving Complex Vehicle Routing Problems, PhD Thesis, Institute of Mathematical Statistics and Operations Research, Technical University of Denmark, Lyngby, 1992.        [ Links ]

HERMOSILLA, A. y BARÁN, B. Comparación de un sistema de colonias de hormigas y una estrategia evolutiva para un Problema Multiobjetivo de Ruteo de Vehículos con Ventanas de Tiempo. s/f., Disponible en: http://www.cnc.una.py/invest/paper2/augCLEI.pdf. Consultado en: Febrero de 2005.        [ Links ]

HURKENS, C. WHIZZKIDS '96 1997, Disponible en: http://www.win.tue.nl/whizzkids/1996/, Consultado en: Febrero de 2005.        [ Links ]

LAPORTE, G. The vehicle routing problem: an overview of exact and approximate algorithms. European Journal of Operational Research, 59, 1992, pp. 345–358.        [ Links ]

LEE, Ch.; EPELMAN, M.; CHELSEA, C. and BOZER Y., A Shortest Path Approach to the Multiple-Vehicle Routing Problem with Split Pick-Ups, 2002., Disponible en: http://osiris.tuwien.ac.at/~wgarn/VehicleRouting/mvrpsp_ts.pdf, Consultado en: Febrero de 2005.        [ Links ]

LIU, F. and SHEN, S. A Method for Vehicle Routing Problem with Multiple Vehicle Types and Time Windows. 1999, Disponible en: http://nr.stic.gov.tw/ejournal/ProceedingA/v23n4/526-536.pdf. Consultado en: Febrero de 2005, Proc. Natl. Sci. Counc. ROC(A), Vol. 23, No. 4, 1999.        [ Links ]

LOUIS, S., Yin, X. and YUAN, Z. Multiple Vehicle Routing with Time Windows Using Genetic Algorithms., Memorias de Proceedings of the Congress of Evolutionary Computation, 1999, pp. 1804-1808.        [ Links ]

MACHADO, P., TAVARES, J., Pereira, F. and Costa, E., Vehicle Routing Problem: Doing it the Evolutionary Way. 2000., Disponible en: http://osiris.tuwien.ac.at/~wgarn/VehicleRouting/GECCO02_VRPCoEvo.pdf. Consultado en: Febrero de 2005.        [ Links ]

MACHADO, P.; PEREIRA, F., TAVARES, J. and COSTA, E., GVR: a New Genetic Representation for the Vehicle Routing Problem. 2002., Disponible en: http://neo.lcc.uma.es/radi-aeb/WebVRP/data/articles/vrp_aics2002.pdf, Consultado en: Febrero de 2005.        [ Links ]

MACHADO, P.; TAVARES, J.; PEREIRA, F. and COSTA, E. Crossover and Diversity: A Study about GVR. 2003 (a)., Memorias de Proceedings of the Analysis and Design of Representations and Operators (ADoRo'2003) a bird-of-a-feather workshop at the 2003 Genetic and Evolutionary Computation Conference (GECCO-2003), Chicago, Illinois, USA, 12-16 July, 2003. Disponible en: http://cisuc.dei.uc.pt/ecos/view_pub.php?id_p=65. Consultado en: Febrero de 2005.        [ Links ]

MACHADO, P.; TAVARES, J.; PEREIRA, F. and COSTA, E. On the Influence of GVR in Vehicle Routing. 2003 (b)., Memorias de Proceedings of the 2003 ACM Symposium On Applied Computing (SAC 2003) - Evolutionary Computation And Optimization Track, pp.753-758, Melbourne, Florida, USA, 9-13 March, 2003. Disponible en: http://cisuc.dei.uc.pt/ecos/view_pub.php?id_p=65. Consultado en: Febrero de 2005         [ Links ]

MEDAGLIA, A. Combinatoria para Logística, Coloquio en Optimización Combinatoria Sesión Avanzada, Universidad de los Andes, marzo de 2005.        [ Links ]

MIN, H. The Multiple Vehicle Routing Problem with Simultaneous Delivery and Pick-Up Points., Transportation Research, vol.23-A, 1989, pp. 377-386.        [ Links ]

OLIVERA, A. Heurísticas para problemas de ruteo de vehículos, Instituto de Computación, Facultad de Ingeniería, Universidad de la República, Montevideo, Uruguay. 2004. Disponible en: http://www.fing.edu.uy/inco/pedeciba/bibliote/reptec/TR0408.pdf. Consultado en Febrero de 2005.        [ Links ]

REGO, C. Node Ejection Chains for the Vehicle Routing Problem: Sequential and Parallel Algorithms. (s/f)., Disponible en: http://hces.bus.olemiss.edu/reports/hces0201.pdf. Consultado en: febrero de 2005.        [ Links ]

RESTORI M. An Application of VRP Algorithms with Original Modifications. s/f., Disponible en: http://eal.asu.edu/Papers%5CMel_VRP_studentTech.pdf, Consultado en: Febrero 2005        [ Links ]

SHAW, P. Using Constraint Programming and Local Search Methods to Solve Vehicle Routing Problems., Memorias de Proceedings of the Fourth International Conference on Principles and Practice of Constraint Programming (CP '98), M. Maher and J. Puget (eds.), 1998, pp. 417-431.         [ Links ]

SCHROEDER, R. Administración de operaciones, toma de decisiones en la función de operaciones., 3ª. Ed., Editorial Mc Graw Hill, México, 1992.         [ Links ]

TANSINI, L.; URQUHART, M. and VIERA, O., Comparing assignment algorithms for the Multi-Depot VRP; s/f; Diponible en: http://www.fing.edu.uy/inco/pedeciba/bibliote/reptec/TR0108.pdf. Consultado en: Febrero 2005        [ Links ]

TAHA, H. Investigación de Operaciones, 5ª. Edición, alfaomega, Grupo Editor México, 1995.        [ Links ]

THOMPSON, P. and ORLIN, J. The Theory of Cyclic Transfers. 1989, Disponible en: http://web.mit.edu/jorlin/www/oldpapersfolder/cyclic_tranfers.pdf. Consultado en Febrero de 2005.        [ Links ]

TOTH, P. and VIGO, D. An Overview of Vehicle Routing Problems, Monographs on Discrete Mathematics and Applications. In: The Vehicle Routing Problem. SIAM, 2000, pp. 1–26.        [ Links ]

VACIC, V. Vehicle Routing Problem with Time Windows; 2002., Disponible en: http://www.bridgeport.edu/sed/projects/cs399/vacic/vrptw.html. Consultado en: Febrero 2005        [ Links ]

VOLKAN, A. A GA Based Meta-Heuristic for the Capacitated Vehicle Routing Problem with Simultaneous Pick-up and Deliveries. 2005., Disponible en: http://www.msie.sabanciuniv.edu/thesis/arif_volkan_vural.ppt. Consultado en: Marzo 2005        [ Links ]

WENDT, O. and KÖNIG, W., Cooperative Simulated Annealing: How much cooperation is enough?. s/f., Disponible en: http://www.wiiw.de/publikationen/Cooperative SimulatedAnnealing48.pdf. Consultado en: Febrero de 2005        [ Links ]

ZHU, K. A new genetic algorithm for VRPTW., Presented at IC-AI 2000, Las Vegas, USA, 2000.        [ Links ]

3 Disponible en: http://www.cs.tcd.ie/courses/baict/bass/4ict5/Networks2004.pdf. Recuperado el 1 de marzo de 2005.
4 Por recomendación de la empresa se considera la demanda constante, dado que la desviación estándar de la misma es poco representativa.
5 Se utilizará el término "minimización" con el ánimo de referirse al menor valor que pueda ser encontrado; sin embargo, los autores no garantizan la toma de la decisión basados en un valor óptimo.
6 Según notación de Ahuja et al., 1993.
7 La nomenclatura utilizada y las restricciones representadas por las expresiones 3 a 9 son tomadas de Ahuja et al., 1993.