Consistent Second-Order Conic Integer Programming for Learning Bayesian Networks