Measurement of subjet multiplicities in neutral current deep inelastic scattering at HERA and determination of α_s