Balanced numbers n such that n mod 12 = 1. All seem to have balance quotient sigma(n)/phi(n) equal to 2. List has been checked out to 0<=m<=500000 with n=12*m+1.
with(numtheory); BNM1:=[]: for z from 1 to 1 do for m from 0 to 500000 do n:=12*m+1; if sigma(n) mod phi(n) = 0 then BNM1:=[op(BNM1), n] fi; od; od; BNM1;
