Counting extreme U₁ matrices and characterizing quadratic doubly stochastic operators