1 Although very successful in accounting for all the major cosmological observables, the ΛCDM model is plagued by fundamental This fact pushes cosmologists to explore alternative interpretations for the accelerated scenario. For instance, plausible landscapes include dynamical dark energy with evolving scalar fields [15][16][17] or evolving equation of state [18][19][20], attempts to unify dark matter and dark energy into a single fluid [21,22], and higher-dimensions braneworld models [23,24].An alternative view is to modify GR on cosmological scales. This turns out to explain the late-time acceleration without the need of dark energy. One of the most studied extensions of GR is represented by f (R) gravity [25][26][27], generalizing the Einstein-Hilbert action with higher-order curvature terms in the Lagrangian. Many theoretical studies carried out so far have focused on the cosmological viability of such theories [28][29][30][31]. Also, from the observational point of view, the viability of these models has been tested by means of several cosmological data surveys [32][33][34][35][36][37][38].Standard approaches toward the study of f (R) paradigms postulate the f (R) functions [39][40][41][42]. This approach is plagued by the fact that f (R) is assumed a priori without relying on first principles. For these reasons, we here implement the inverse procedure, i.e. we start issues, such as the coincidence problem [13] and the fine-tuning problem [14].