Let $n\geq3$ and $A\subset\mathbb{R}^n$ be a discrete subset. How to characterize all the positive harmonic function on $\mathbb{R}^n\setminus A$ ?
Rmk: This is from Exercise 17 in Chapter 3 of the textbook named 'Harmonic Function Theory' (GTM 137) by Axler, Bourdon, and Ramey. If $A = \{p_1,\dotsc,p_m\}$ is finite, It's easy to show the positive harmonic functions $V$ can only be of the form $$V = \sum_{i = 1}^m \frac{a_i}{|x - p_i|^{n-2}} + C$$ where $C>0$ is constant and $a_i >0$ (You could find one proof at Sheldon Axler's answer to Characterize a positive harmonic functions with multi-singularities). As Noam commented, if the sum $$V = \sum_{i = 1}^\infty \frac{a_i}{|x - p_i|^{n-2}} + C $$ converges for $a_i\geq0$, then $V$ is a solution. But it's harder to show they are the only solutions with singularities at $A$. Moreover, since Bôcher's theorem classifies positive harmonic function on punctured ball, I also want to know is there any possibility to classify positive harmonic function on any open set in $\mathbb{R}^n$ ?