complete_randomization_with_forced_balance_plus_one_min_one function

Implements the balanced complete randomization design (BCRD)