python 実践データ分析100本ノック ノック61に関する質問です。
上記のコードを実行すると、総輸送コストが1450になり、書籍の1298とは違った結果になってしまいます。
実行環境はjupiter notebook でpulpとortoolpyはどちらもconda promptでインストール済みです。
コードに問題があれば教えて頂きたいです。
import numpy as np
import pandas as pd
from itertools import product
from pulp import LpVariable,lpSum,value
from ortoolpy import model_min,addvars,addvals
df_tc = pd.read_csv("trans_cost.csv",index_col="工場")
df_demand = pd.read_csv("demand.csv")
df_supply= pd.read_csv("supply.csv")
np.random.seed(1)
nw=len(df_tc.index)
nf=len(df_tc.columns)
pr=list(product(range(nw),range(nf)))
m1=model_min()
v1={(i,j):LpVariable("v%d_%d"%(i,j), lowBound=0) for i,j in pr}
m1 +=lpSum(df_tc.iloc[i][j]*v1[i,j] for i,j in pr)
for i in range(nw):
m1 += lpSum(v1[i,j] for j in range(nf)) <= df_supply.iloc[0][i]
for j in range(nf):
m1 += lpSum(v1[i,j] for i in range(nw)) >= df_demand.iloc[0][i]
m1.solve()
df_tr_sol = df_tc.copy()
total_cost=0
for k,x in v1.items():
i,j = k[0],k[1]
df_tr_sol.iloc[i][j]=value(x)
total_cost +=df_tc.iloc[i][j]*value(x)
print(df_tr_sol)
print("総輸送コスト:"+str(total_cost))