Mass-radius relation for neutron stars in f(R)=R+αR2 gravity: A comparison between purely metric and torsion formulations