# -*- coding: utf-8 -*-
from random import randrange

# Configuración inicial - MODIFICAR ESTO AL GUSTO
individuos_iniciales = 300000
dias_a_simular = 365


individuos_actuales = individuos_iniciales
suma_poblacion = individuos_actuales # para obtener la media poblacional
muertes = nacimientos = 0


for dia in range(1, dias_a_simular):
    evento = randrange(6) + 1 # rango 0 - 5
    
    # Actuar en base al evento aleatorio    
    # Muere uno
    if evento == 1:
        individuos_actuales += 1
        nacimientos += 1
    
    # Mantiene
    if evento == 2:
        pass        
    
	# Nacen dos
    if evento == 3:
        individuos_actuales += 2
        nacimientos += 2
    
	# Mueren dos
    if evento == 4:
        individuos_actuales -= 2    
        muertes += 2
    
	# Nace uno
    if evento == 5:
        individuos_actuales += 1
        nacimientos += 1

    # Muere uno
    if evento == 6:
        individuos_actuales -= 1
        muertes += 2        

    # Sumar la población cada día para obtener la media
    suma_poblacion += individuos_actuales
    
# Terminó la simulación
media_poblacional = suma_poblacion / dias_a_simular
tasa_mortalidad = (muertes / media_poblacional) * 1000
tasa_natalidad = (nacimientos / media_poblacional) * 1000

print("Individuos iniciales: %s" % individuos_iniciales)
print("Días simulados: %s " % dias_a_simular)
print("")
print("Tras simular los días indicados, este es el resultado:")
print("Individuos finales: %s" % individuos_actuales)
print("    - Ha habido una variación de %s individuos respecto a los iniciales" % str((individuos_actuales - individuos_iniciales)))
print("    - Se han producido %s defunciones" % muertes)
print("    - Se han producido %s nacimientos" % nacimientos)
print("    - La tasa de mortalidad ha sido de %s muertes por 1000 habitantes" % round(tasa_mortalidad, 1))
print("    - La tasa de natalidad ha sido de %s nacimientos por 1000 habitantes" % round(tasa_natalidad, 1))