Tarea 3: Búsqueda del subarreglo máximo

import sys

def encuentra_max_subarreglo_cruzando(A, bajo, medio, alto):
    """
    Encuentra el subarreglo máximo que cruza el punto medio.

    Args:
        A: Lista de números
        bajo: Índice inferior
        medio: Índice medio
        alto: Índice superior

    Returns:
        Tupla (max_izq, max_der, suma_total)
    """
    suma_izq = float('-inf')
    suma = 0
    max_izq = medio

    # Buscar hacia la izquierda desde el medio
    for i in range(medio, bajo - 1, -1):
        suma += A[i]
        if suma > suma_izq:
            suma_izq = suma
            max_izq = i

    suma_der = float('-inf')
    suma = 0
    max_der = medio + 1

    # Buscar hacia la derecha desde el medio
    for j in range(medio + 1, alto + 1):
        suma += A[j]
        if suma > suma_der:
            suma_der = suma
            max_der = j

    return (max_izq, max_der, suma_izq + suma_der)


def encuentra_subarreglo_maximo(A, bajo, alto):
    """
    Encuentra el subarreglo máximo usando divide y vencerás.

    Args:
        A: Lista de números
        bajo: Índice inferior
        alto: Índice superior

    Returns:
        Tupla (inicio, fin, suma)
    """
    # Caso base: un solo elemento
    if bajo == alto:
        return (bajo, alto, A[bajo])
    else:
        medio = (bajo + alto) // 2

        # Encontrar máximo en la mitad izquierda
        (bajo_izq, alto_izq, suma_izq) = encuentra_subarreglo_maximo(A, bajo, medio)

        # Encontrar máximo en la mitad derecha
        (bajo_der, alto_der, suma_der) = encuentra_subarreglo_maximo(A, medio + 1, alto)

        # Encontrar máximo que cruza el medio
        (bajo_cruza, alto_cruza, suma_cruza) = encuentra_max_subarreglo_cruzando(A, bajo, medio, alto)

        # Retornar el máximo de los tres
        if suma_izq >= suma_der and suma_izq >= suma_cruza:
            return (bajo_izq, alto_izq, suma_izq)
        elif suma_der >= suma_izq and suma_der >= suma_cruza:
            return (bajo_der, alto_der, suma_der)
        else:
            return (bajo_cruza, alto_cruza, suma_cruza)


def leer_datos_archivo(nombre_archivo):
    """
    Lee datos de un archivo y retorna una lista de precios.

    Args:
        nombre_archivo: Ruta del archivo

    Returns:
        Lista de precios
    """
    precios = []
    try:
        with open(nombre_archivo, 'r') as archivo:
            for linea in archivo:
                linea = linea.strip()
                if linea and not linea.startswith('#'):
                    # Dividir y tomar el segundo valor (precio)
                    partes = linea.split()
                    if len(partes) >= 2:
                        try:
                            precio = float(partes[1])
                            precios.append(precio)
                        except ValueError:
                            continue
    except FileNotFoundError:
        print(f"Error: No se pudo encontrar el archivo {nombre_archivo}")
        return []
    except Exception as e:
        print(f"Error al leer el archivo {nombre_archivo}: {e}")
        return []

    return precios


def calcular_cambios_diarios(precios):
    """
    Calcula los cambios diarios (diferencias) entre precios consecutivos.

    Args:
        precios: Lista de precios

    Returns:
        Lista de cambios diarios
    """
    cambios = []
    for i in range(1, len(precios)):
        cambio = precios[i] - precios[i-1]
        cambios.append(cambio)
    return cambios


def main():
    """
    Función principal que ejecuta el algoritmo de análisis financiero.
    Uso: python max_subarray.py <archivo>
    """
    if len(sys.argv) != 2:
        print("Uso: python max_subarray.py <archivo>")
        print("Ejemplo: python max_subarray.py datos1.txt")
        sys.exit(1)

    archivo = sys.argv[1]

    print("=" * 60)
    print(f"ANALISIS FINANCIERO - ARCHIVO: {archivo}")
    print("=" * 60)

    # Leer precios del archivo
    precios = leer_datos_archivo(archivo)

    if not precios or len(precios) < 2:
        print(f"Error: No se pudieron leer suficientes datos de {archivo}")
        sys.exit(1)

    print(f"\nPRECIOS ORIGINALES ({len(precios)} dias):")
    print("-" * 40)
    for i, precio in enumerate(precios):
        print(f"Dia {i:2d}: ${precio:6.1f}")

    # Calcular cambios diarios
    cambios = calcular_cambios_diarios(precios)
    print(f"\nCAMBIOS DIARIOS ({len(cambios)} periodos):")
    print("-" * 40)
    for i, cambio in enumerate(cambios):
        signo = "+" if cambio >= 0 else ""
        print(f"Dia {i+1:2d}->{i+2:2d}: {signo}{cambio:6.1f}")

    # Ejecutar el algoritmo en los cambios diarios
    inicio, fin, suma_maxima = encuentra_subarreglo_maximo(cambios, 0, len(cambios) - 1)

    print(f"\nESTRATEGIA OPTIMA DE INVERSION:")
    print("=" * 60)

    # Mostrar el subarreglo máximo de forma destacada
    subarreglo = cambios[inicio:fin+1]
    print(f"SUBARREGLO MAXIMO ENCONTRADO:")
    print("-" * 40)
    print(f"Valores: {subarreglo}")
    print(f"Indices: [{inicio} - {fin}]")
    print(f"Suma total: {suma_maxima:.1f}")

    print(f"\nRESUMEN FINANCIERO:")
    print("-" * 40)
    print(f"GANANCIA MAXIMA: ${suma_maxima:.1f}")
    print(f"COMPRAR: Dia {inicio + 1} (Precio: ${precios[inicio + 1]:.1f})")
    print(f"VENDER:  Dia {fin + 2} (Precio: ${precios[fin + 2]:.1f})")
    print(f"DURACION: {fin - inicio + 1} dias")

    print(f"\nDETALLES DEL PERIODO OPTIMO:")
    print("-" * 40)
    for i in range(inicio, fin + 1):
        dia_inicio = i + 1
        dia_fin = i + 2
        cambio = cambios[i]
        signo = "+" if cambio >= 0 else ""
        print(f"Dia {dia_inicio:2d}->{dia_fin:2d}: {signo}{cambio:6.1f}")

    # Verificación con precios reales
    precio_compra = precios[inicio + 1]
    precio_venta = precios[fin + 2]
    ganancia_real = precio_venta - precio_compra

    print(f"\nVERIFICACION:")
    print("-" * 40)
    print(f"Precio de compra: ${precio_compra:.1f}")
    print(f"Precio de venta:  ${precio_venta:.1f}")
    print(f"Ganancia real:    ${ganancia_real:.1f}")
    print(f"Suma de cambios:  ${sum(cambios[inicio:fin+1]):.1f}")

    print("\n" + "=" * 60)


if __name__ == "__main__":
    main()