import math  

def solve_quadratic(a, b, c):

    discriminant = (b**2) - (4 * a * c)

    if discriminant > 0:
        sol1 = (-b + math.sqrt(discriminant)) / (2 * a)
        sol2 = (-b - math.sqrt(discriminant)) / (2 * a)
        return f"Two real solutions: {sol1} and {sol2}"

    elif discriminant == 0:
        sol = -b / (2 * a)
        return f"One real solution: {sol}"

    else:
        return "No real solutions (discriminant < 0)"


print("Quadratic equation: ax^2 + bx + c = 0")

a = float(input("Enter a: "))
b = float(input("Enter b: "))
c = float(input("Enter c: "))

print(solve_quadratic(a, b, c))