using JuMP, Gurobi

m, n = 10, 10
w1, w2 = 1, 5
P = 1 # 1 pour modèle linéaire, 2 pour quadratique
L = 3
g = 1.26157

# Matrice t à importer de tes données
t = [[84,68,97,98,64,89,82,71,74,76],
   [87,83,98,75,60,90,78,67,92,94],
   [84,68,70,81,67,61,73,92,86,90],
   [79,62,86,79,73,84,76,98,84,90],
   [62,72,66,72,92,80,71,91,87,70],
   [85,77,63,93,90,94,76,81,99,98],
   [76,63,66,84,94,93,72,92,79,65],
   [76,63,92,69,60,88,79,93,66,73],
   [92,82,77,72,77,81,89,95,80,80],
   [88,89,83,86,69,78,91,64,94,92]
]

# Création du tenseur 2D : t_val[i, j]
t_val = [t[i][j] for i in 1:m, j in 1:m]

# Calcul des voisins (N[i,j])
voisins(i, j) = [(k, l) for (k, l) in [(i-1,j), (i+1,j), (i,j-1), (i,j+1)] if 1 <= k <= m && 1 <= l <= n]
N_voisins = [length(voisins(i, j)) for i in 1:m, j in 1:n]