import math

def solve_equation(a, b, c):
    if a == 0:
        if b == 0:
            return "No solution or infinite solutions."
        else:
            return f"Linear solution: x = {-c / b}"
    
    discriminant = b**2 - 4*a*c
    
    if discriminant > 0:
        x1 = (-b + math.sqrt(discriminant)) / (2*a)
        x2 = (-b - math.sqrt(discriminant)) / (2*a)
        return f"Two real roots: x1 = {x1}, x2 = {x2}"
    elif discriminant == 0:
        x = -b / (2*a)
        return f"One real root: x = {x}"
    else:
        return "No real roots."

print(f"a=1, b=1, c=1   => {solve_equation(1, 1, 1)}")