Calculation of centrality bias factors in $p$+A collisions based on a positive correlation of hard process yields with underlying event activity