\begin{equation}
\left\{
    \begin{array}{c}R_{xi\_yj}(t+1)=R_{xi\_yi}(t)+\Delta r_{xi\_yj}(t)\\
    \left\{
        \begin{array}{l} k_{xi\_yj}=1, \quad if\quad I_{xi\_yj}>i_t\\ 
        k_{xi\_yj}=0, \quad otherwise\\ 
        k=\sum_{i=1}^{m} \sum_{j=1}^{n} k_{xi\_yj} \\ 
        C=\sum_{i=1}^{m} \sum_{j=1}^{n} R_{xi\_yj}(t)\times k_{xi\_yj}\times L_r
        \end{array}
    \right.\\
    \Delta r_{xi\_yi}(t)=
    \left\{
        \begin{array}{l} -R_{xi\_yj}(t)\times L_r, \quad if\quad I_{xi\_yj}>I_t\\
        \frac{C}{m\times n-k} , \quad otherwise
        \end{array}
    \right.\\
    \sum_{i=1}^{m} \sum_{j=1}^{n} \Delta r_{xi\_yj}(t)=0
    \end{array}
\right.
\end{equation}
