Bilinear optimal control for the fractional Laplacian: analysis and discretization