Why would the universe care about dimensional analysis? Besides, the outside constant would do the unit conversion from whatever the right-hand side produces to units of force.
We're lucky (and maybe only in the linear domain, which is more accurately the case) that these ended up being the same.
But if it were sqrt(charge) then we'd define Culoumbs constant to have commensurate units to translate sqrt(charge) * sqrt(charge) / distance^2 to force units.
Partitioning a charge can't change the physics. If i have a bar of charge Q generating a force F on q which is proportional to sqrt Q then by partition i have n bars of charge Q_i = Q/n producing forces F_i
sqrt Q ~ F <> sum F_i = sum sqrt(Q_i) = sum sqrt(Q/n) = n * sqrt (Q/n) = sqrt (Qn)
Partitioning charge would lead to infinite forces.