[ABRIDGED] We aim to provide a holistic view on the typical size and kinematic evolution of massive early-type galaxies (ETGs), that encompasses their high-$z$ star-forming progenitors, their high-$z$ quiescent counterparts, and their configurations in the local Universe. Our investigation covers the main processes playing a relevant role in the cosmic evolution of ETGs. Specifically, their early fast evolution comprises: biased collapse of the low angular momentum gaseous baryons located in the inner regions of the host dark matter halo; cooling, fragmentation, and infall of the gas down to the radius set by the centrifugal barrier; further rapid compaction via clump/gas migration toward the galaxy center, where strong heavily dust-enshrouded star-formation takes place and most of the stellar mass is accumulated; ejection of substantial gas amount from the inner regions by feedback processes, which causes a dramatic puffing up of the stellar component. In the late slow evolution, passive aging of stellar populations and mass additions by dry merger events occur. We describe these processes relying on prescriptions inspired by basic physical arguments and by numerical simulations, to derive new analytical estimates of the relevant sizes, timescales, and kinematic properties for individual galaxies along their evolution. Then we obtain quantitative results as a function of galaxy mass and redshift, and compare them to recent observational constraints on half-light size $R_e$, on the ratio $v/\sigma$ between rotation velocity and velocity dispersion (for gas and stars) and on the specific angular momentum $j_\star$ of the stellar component; we find good consistency with the available multi-band data in average values and dispersion, both for local ETGs and for their $z\sim 1-2$ star-forming and quiescent progenitors.

[abridged] We investigate the coevolution of galaxies and hosted supermassive black holes throughout the history of the Universe by a statistical approach based on the continuity equation and the abundance matching technique. Specifically, we present analytical solutions of the continuity equation without source term to reconstruct the supermassive black hole (BH) mass function from the AGN luminosity functions. Such an approach includes physically-motivated AGN lightcurves tested on independent datasets, which describe the evolution of the Eddington ratio and radiative efficiency from slim- to thin-disc conditions. We nicely reproduce the local estimates of the BH mass function, the AGN duty cycle as a function of mass and redshift, along with the Eddington ratio function and the fraction of galaxies with given stellar mass hosting an AGN with given Eddington ratio. We exploit the same approach to reconstruct the observed stellar mass function at different redshift from the UV and far-IR luminosity functions associated to star formation in galaxies. These results imply that the buildup of stars and BHs in galaxies occurs via in-situ processes, with dry mergers playing a marginal role at least for stellar masses < 3 10^11 M_sun and BH masses < 10^9 M_sun, where the statistical data are more secure and less biased by systematic errors. In addition, we develop an improved abundance matching technique to link the stellar and BH content of galaxies to the gravitationally dominant dark matter component. The resulting relationships constitute a testbed for galaxy evolution models, highlighting the complementary role of stellar and AGN feedback in the star formation process. Finally, the clustering properties of BHs and galaxies are found to be in full agreement with current observations, so further validating our results from the continuity equation.