The requested generalization of Jacobi's two-square theorem is a remarkably recent result: N. Bagis and M.L Glasser, On the Number of Representations of Integers by various Quadratic and Higher Forms (2015):
http://ilorentz.org/beenakker/MO/sumofsquares_2.pngwhere $r(n)$ is Jacobi's formula for the number of representations of $n=x^2+y^2$.
I understand from the comments that this (general but complicated) formula is not what the OP was looking for. Berkovich and Yesilyurt in Ramanujan's Identities and Representation of Integers by Certain Binary and Quaternary Quadratic Forms (2006) give rather simple expressions for specific cases of $n=x^2+ay^2$, I reproduce below the result for $a=5$, and there are others (6,15,27,...).
http://ilorentz.org/beenakker/MO/arXiv_0611300.png