Obviously a necessary condition for [itex]\vec{B}[/itex] to be the curl of a vector field is
[tex]\vec{\nabla} \cdot \vec{B}=0.[/tex]
Then the vector potential is only determined up to the gradient of a scalar field. Thus one can impose one additional condition ("choice of a gauge"). In this case an axial gauge is chosen, i.e.,
[tex]\vec{\nabla} \times \vec{A}=\vec{B}, \quad A_x=0.[/tex]
Then we have in Cartesian coordinates
[tex]B_x=\partial_y A_z-\partial_z A_y, \quad B_y=-\partial_x A_z, \quad B_z=\partial_x A_y.[/tex]
From the last equation we get
[tex]A_y=\int_{x_0}^{x} \mathrm{d} x' B_z(x',y,z)+\tilde{A}_y(y,z)[/tex]
and from the 2nd equation
[tex]A_z=-\int_{x_0}^x \mathrm{d} x' B_y(x',y,z) + \tilde{A}_z(y,z).[/tex]
The first equation now reads
[tex]B_x(x,y,z)=-\int_{x_0}^{x} \mathrm{d} x' [\partial_y B_y(x',y,z)+\partial_z B_z(x',y,z)]+\partial_y \tilde{A}_z(y,z)-\partial_z \tilde{A}_y(y,z).[/tex]
Now from [itex]\vec{\nabla} \cdot \vec{B}=0[/itex] we find
[tex]B_x(x,y,z)=B_x(x,y,z)-B_x(x_0,y,z) +\partial_y \tilde{A}_z(y,z)-\partial_z \tilde{A}_y(y,z)[/tex]
or
[tex]B_x(x_0,y,z)=\partial_y \tilde{A}_z(y,z)-\partial_z \tilde{A}_y(y,z).[/tex]
This equation we can satisfy by setting
[tex]\tilde{A}_y=0, \quad \tilde{A}_z(y,z)=\int_{y_0}^{y} \mathrm{d}y' B_x(x_0,y',z).[/tex]
Putting all together we get
[tex]A_x=0, \quad A_y=\int_{x_0}^{x} \mathrm{d} x' B_z(x',y,z), \quad A_z=\int_{y_0}^{y} \mathrm{d}y' B_x(x_0,y',z)-\int_{x_0}^x \mathrm{d} x' B_y(x',y,z).[/tex]
This is what's claimed in the OP.