Линейное программирование (lpsolve) дает другое решение при добавлении константы в функцию

Question

Я явно делаю что-то не так, но не могу понять, что. У меня есть матрица функции для максимизации названного здесь ранга (если интересует ожидаемая производительность потомства от спаривания самца и самки) с некоторыми ограничениями - все значения функции положительны. Я запускаю его и получаю решение. Затем я добавляю константу (здесь 10) ко всем значениям функции и получаю другое решение. Я не эксперт по LP, но это не имеет смысла - конечно, добавление константы ко всему не должно изменить вещи?

library(lpSolve)
rank <- matrix(c(1, 2, 3, 4, 5,6, 7, 8), nrow=2, ncol=4,byrow=T)
rank<- rank+10  # switch on and off
rank
exclusions<-  matrix(c(0,0,0,0,0,0,0,0), nrow=2, ncol=4,byrow=T)
sires <- matrix(c(2,2), nrow=1, ncol=2,byrow=T)
m <- NROW(rank) # n of males
f <- NCOL(rank) # n of females
obj <- as.numeric(rank)    # objective function
nMalePerFemaleRhs <- rep(1,f)    # one mating per female 
nMalePerFemaleSign <- rep("<=",f)   # every female must gets no more than 1 mating 
nMalePerFemaleConstr <- matrix(0,nrow=f,ncol=m*f)

for(i in 1:f){
  nMalePerFemaleConstr[i,(i-1)*m+(1:m)] <- 1    # cows are rows and the number of columns are bulls times cows
}

nFemalePerMaleRhs <- sires[1,] # number of matings per sire
nFemalePerMaleSign <- rep("<=",m)  
nFemalePerMaleConstr <- matrix(0,nrow=m,ncol=m*f)

for(i in 1:m){
  nFemalePerMaleConstr[i,seq(from=i,length.out=f,by=m)] <- 1
}


######################################################
#         Sibling exclusions constraint              #
######################################################

siblingConstr <- t(as.numeric(exclusions))
siblingRhs <- 0
siblingSign <- '='


######################################################
#           Solve the linear program                 #
######################################################

res <- lp(direction='max',
          objective.in=obj,
          const.mat = rbind
(nMalePerFemaleConstr,nFemalePerMaleConstr,siblingConstr),
          const.dir = c(nMalePerFemaleSign,nFemalePerMaleSign,siblingSign),
          const.rhs = c(nMalePerFemaleRhs,nFemalePerMaleRhs,siblingRhs),
          all.int = TRUE
          )

solutionLP <- matrix(res$solution,nrow=m)

solutionLP