You use the same formula. However, unlike in a QFT generating functional, where we can describe general correlation functions of off-shell objects, here we can only calculate on-shell bulk amplitudes. On the boundary, we have correlators of off-shell operators, but the sources for the bulk fields are the on-shell boundary values of the fields. Explicitly, we can say something like for an ##n##-point function in the bulk
$$ \int d\mathbf{x}'_1 dr_1 \cdots d\mathbf{x}_n dr_n A_{1\cdots n}(\mathbf{x}'_1,r_1;\ldots \mathbf{x}'_n,r_n) K(\mathbf{x}'_1,r_1; \mathbf{x}_1) \cdots K(\mathbf{x}'_n,r_n; \mathbf{x}_n) =\langle \mathcal{O}_1(\mathbf{x}_1) \cdots \mathcal{O}_n(\mathbf{x}_n)\rangle,$$
where ##\mathbf{x}_I## are the boundary coordinates, ##r_I## is the radial variable and ##K(\mathbf{x}'_I,r_I; \mathbf{x}_I) ## are the boundary to bulk propagators. It is not expected that one can manipulate these expressions to obtain the bare ##A_{1\cdots n}##, but the amplitudes for all combinations of physical states are in principle determined from the boundary correlators.