0

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))
3
  • 念のため、書籍の正誤表 も確認してみてください。
    – cubick
    Commented 2022年1月13日 12:29
  • コメントありがとうございます。 正誤表にも特に記載はありません。 Commented 2022年1月13日 13:17
  • 私の手元にある「python 実践データ分析100本ノック第1版第1刷」では1296が結果として記載されています。kunifさんが回答で示されている値と一緒です。 Commented 2022年1月14日 1:23

1 件の回答 1

0

配布されているsample_100knocks.zip中の該当answerと比べると、23行目末尾df_demand.iloc[0][i]のインデックスが間違っているようです。

正しくはiではなくjであり、df_demand.iloc[0][j]となっています。

ちなみにAnaconda環境ではなくjupyter notebookでもない、Windows10の素のPython 3.10.1 では総輸送コストが1296というまた少し違った値になったので、何か他のモジュールの版数とかあるいはCSVデータとかに違いがあるのかもしれません。


@ほーりー さんのコメントとか以下記事に総輸送コストは1296と書かれていて特に注釈とか無いため1296が正しいようなので、質問に「書籍の1298」とあるのはtypoのようですね。

『Python実践データ分析100本ノック』ノック61
Python実践データ分析100本ノックを読む

1
  • コメントありがとうございます。ご回答の通りに修正したところ正しい数値を出力できました。ご指導ありがとうございました Commented 2022年1月14日 11:45

この質問に回答するには、ログインする必要があります。

求めていた回答ではありませんか? のタグが付いた他の質問を参照する。