library(lpSolve)
coefiniciales <- c(28,18,42,36,40,28,54,48,36,24,52,40,38,30,54,46) #coeficientes de los costos de transporte
prueba<- rep(c(60,72,48,60),4)
coef<- coefiniciales+prueba #sumo los costos de producir de cada fabrica al costo de transporte
A1<-c(1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0)
A2<-c(0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0)
A3<-c(0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0)
A4<-c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1)
A5<-c(1,0,0,0,1,0,0,0,1,0,0,0,1,0,0,0)
A6<-c(0,1,0,0,0,1,0,0,0,1,0,0,0,1,0,0)
A7<-c(0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,0)
A8<-c(0,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1)
A <- rbind(A1,A2,A3,A4,A5,A6,A7,A8)#genero la matriz de los coeficientes de las restricciones, la disposición de 0's y 1's se hizo numerando del 1 al 16 cada una de las variables
#lo de arriba no es lo optimo pero no se me ocurrió /no encontré otra manera
print(A)
direccion<-c('<=','<=','<=','<=','=','=','=','=') #direccion de las restricciones
b<-c(140,260,360,220,180,280,150,200) #valor de las restricciones
solucion<-lp('min',coef,A,direccion,b)
print(solucion$solution)
print(solucion$objval)