We develop a solution theory for singular elliptic stochastic PDEs with fractional Laplacian, additive white noise and cubic nonlinearity. The method covers the whole subcritical regime. It is based on the Wilsonian renormalization group theory and the Polchinski flow equation.