Performance-guaranteed regularization in maximum likelihood method: Gauge symmetry in Kullback -- Leibler divergence