Написать перегруженные шаблоны функций для нахождения корней линейного (a*x + b = 0) уравнения

template<typename A, typename B>
auto f(const A &a, const B&b) -> decltype(a + b) {
    return -b / a;
}

template<typename A, typename B, typename C>
auto f(const A &a, const B&b, const C &c) -> vector<complex<decltype(a + b + c)>> {
    using T = decltype(a + b + c);
    auto D_2 = static_cast<complex<T>>(sqrt(complex<double>(b * b - 4 * a * c)));
    return {(-T(b) - D_2) / T(2) / T(a), (-T(b) - D_2) / T(2) / T(a) };
}

#include <complex>
#include <vector>
#include <iostream>
 
using namespace std;
 
template<typename T, typename A, typename B, typename C>
vector<T> f(const A &a, const B&b, const C &c) {
    auto D_2 = static_cast<T>(sqrt(T(b * b - 4 * a * c)));
    return{(-T(b) - D_2) / T(2) / T(a), (-T(b) - D_2) / T(2) / T(a)};
}
 
int main() {
    auto res = f<complex<double>>(1, 1, 1);
 
    cout << "Roots: ";
    for (auto &it : res)
        cout << it.real() << (it.imag() >= 0 ? "+": "") << it.imag() << "*i; ";
    cout << endl;
 
    return 0;
}