I have two variables that are continuous: $$\gamma^n \geq 0\\ z_{ij}^n \geq 0$$ and one binary variable: $$y_{jk} = \{0,1\}$$
How can I linearize an equation that involves two kinds of variable products, specifically this? $$\sum_i z_{ij}^n(a_{jk} \varepsilon_{i}^n + b_{jk}) y_{jk} \leq \sum_i z_{ij}^nc_j + M(\gamma^n - 1)\sum_i z_{ij} ~~~ \forall j,k,n$$
$a_{jk}, \varepsilon_i^n, b_{jk}, c_j$ and $M$ are all constants.