apis/examples_python/transport10_8py_source.html

가져오기os

가져오기sys

from게임가져오기GamsWorkspace

에서openpyxl가져오기load_workbook


GAMS_MODEL ="""

설정

i '통조림 식물'

j '시장';


매개변수

a(i) '경우에 따른 식물 i의 용량'

b(j) '경우에 따라 시장 j의 수요'

d(i,j) '거리(천 마일)';


스칼라 f '1,000마일당 케이스당 운임(달러)' / 90 /;


$if 설정되지 않음 gdxincname $abort '제공된 데이터 파일에 대한 포함 파일 이름 없음'

$gdxIn %gdxincname%

$load i j a b d

$gdxIn


매개변수 c(i,j) '케이스당 수천 달러의 운송 비용';

c(i,j) = f*d(i,j)/1000;


변수

x(i,j) '케이스 내 배송 수량'

z '총 운송 비용(수천 달러)';


양수 변수 x;


수식

비용 '목적 함수 정의'

supply(i) '공장 i의 공급 제한을 준수하세요'

demand(j) '시장 j의 수요를 충족';


비용.. z =e= sum((i,j), c(i,j)*x(i,j));


공급(i).. 합계(j, x(i,j)) =l= a(i);


수요(j).. 합계(i, x(i,j)) =g= b(j);


모델 수송 /all/;


z를 최소화하는 lp를 사용하여 전송을 해결합니다.


x.l, x.m 표시;

"""


if__name__ =="__main__":

sys_dir = sys.argv[1]iflen(sys.argv) > 1그밖에 없음

work_dir = sys.argv[2]iflen(sys.argv) > 2그밖에 없음

ws = GamsWorkspace(system_directory=sys_dir,working_directory=work_dir)


wb = load_workbook(

os.path.join(ws.system_directory,"api파일", "데이터", "transport.xlsx")

    )

용량 = wb["용량"]

수요 = wb["수요"]

거리 = wb["거리"]


    #개의 시장/공장 수는 모든 스프레드시트에서 동일해야 합니다

    if (

distance.max_column - 1 != 수요.max_column

        또는distance.max_row - 1 != 용량.max_row

):

        올림예외("스프레드시트의 크기가 일치하지 않습니다")


db = ws.add_database()

i = db.add_set("나", 1,"식물")

j = db.add_set("j", 1,"시장")

capacity_param = db.add_parameter_dc("a", [i],"용량")

demand_param = db.add_parameter_dc("b", [j],"수요")

distance_param = db.add_parameter_dc("d", [i, j],"거리")


    forcincapacity.iter_cols():

키 = c[0].값

i.add_record(키)

capacity_param.add_record(key).value = c[1].value


    forcindemand.iter_cols():

키 = c[0].값

j.add_record(키)

demand_param.add_record(key).value = c[1].value


    forcin범위(2, distance.max_column + 1):

        forr에범위(2, distance.max_row + 1):

키 = (

거리.셀(행=r, 열=1).값,

거리.셀(행=1, 열=c).값,

            )

v = distance.cell(행=r, 열=c).값

distance_param.add_record(keys).value = v


    # GamsJob 생성 및 실행

작업 = ws.add_job_from_string(GAMS_MODEL)

opt = ws.add_options()

opt.defines["gdxincname"] = db.name

opt.all_model_types ="익스프레스"

job.run(opt, 데이터베이스=db)

    forrec에job.out_db["x"]:

인쇄(

f"x(rec.key(0),rec.key(1)): 레벨=rec.level 한계=rec.marginal"

        )