Задается любое положительное число R.

#include <iostream>
#include <cmath>
#include <vector>
#include <iterator>
#include <numeric>
 
using namespace std;
 
// функция разбивает число с на составляющие так,
// что c = a*b = a+b
// возвращает false при невозможности
bool split(double c, double & a, double & b)
{
    if (c <= 4)
        return false;
    double d = sqrt(c*c-4*c);
    a = 0.5*(c-d);
    b = a/(a-1);
    return true;
}
 
int main()
{
    std::vector<double> v(1);
    cout << "Input value: ";
    cin >> v[0];
    // разбиение пока все составляющие не окажутся меньше 4
    for (size_t i = 0; i < v.size(); ++i)
    {
        if (v[i] < 4.0)
            ++i;
        else
        {
            double a, b;
            split(v[i], a, b);
            v[i] = a;
            v.push_back(b);
        }
    }
    double sum = 0;
    // вывод суммы
    copy(v.begin(), --v.end(), ostream_iterator<double>(cout,"+"));
    cout << v.back() << " = " 
        << accumulate(v.begin(), v.end(), 0.0)
        << endl;
    // вывод произведения
    copy(v.begin(), --v.end(), ostream_iterator<double>(cout,"*"));
    cout << v.back() << " = " 
        << accumulate(v.begin(), v.end(), 1.0, 
                multiplies<double>())
        << endl;
    return 0;
}

if (c < 4.0)