Dados dois naturais m e n determinar, entre todos os pares de números naturais
(x,y) tais que x < m e y < n, um par para o qual o valor da expressão
x*y - x**2 + y seja máximo e calcular também esse máximo.


Eu resolvi o problema assim:

    m = int(input("Digite m: "))
    
    n = int(input("Digite n: "))
    
    #m = 2
    #n =1
    x = 0
    y = 0
    maximo =0
    x_maximo = 0
    y_maximo = 0
    while x <m and y<n:
        for i in range(m):
            for j in range(n):
    
                soma = i*j -i**2 +j
                print (soma)
                if soma > maximo:
                    maximo = soma
                    x_maximo = i
                    y_maximo = j
        x +=1
        y +=1
    print ("O máximo é {} para o par (x,y) = ({},{})". format(maximo,x_maximo,y_maximo))

Aparentemente, o algoritmo está correto mas eu achei que a solução não está muito elegante. Alguma sugestão de como melhorá-lo? Alguma maneira mais simples?