<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def hmean_int(a, a_min=5778, a_max=1149851):
 """ Harmonic mean of an array, returns the closest int """

    from scipy.stats import hmean
    return int(round(hmean(np.clip(a, a_min, a_max))))