Normal Masuda Method shiny odds: 1/512 or 0.1953125%
(1÷512)×(3÷5)×(1÷25)×(1÷2) = 0.0000234375 or 0.00234375% for one Lapras to have each of those specific features you listed, including HA pass down chance, (60/100 or 3/5,) Brave Nature, (1/25,) and Male, (1/2).
Using binomial probability, we can find the odds of finding 2 Lapras with these features in 31 attempts. 0.0000234375 being the chance of an event occurring, 31 being the number of attempts, and 2 being the number of successes, we get 0.00000025525857 or a 0.000025525857% chance. That's 25525857/100 trillion.
I don't really want to do these calculations myself because I know they don't make up a perfect model for Pokemon's RNG...
This number isn't exactly right because of Pokemon's complicated RNG, but it paints a pretty close picture as to how rare this occurrence is.