apis/examples_python/tsp_8py_source.html

fromitertools가져오기제품

가져오기os

가져오기sys

from게임가져오기GamsWorkspace, GamsModifier, SymbolUpdateType


GAMS_MODEL ="""

$tspdata가 설정되지 않은 경우 $abort 'tspdata가 설정되지 않음'


설정

ii '도시'

i(ii) '도시의 하위 집합';


별칭 (ii,jj), (i,j,k);


매개변수 c(ii,jj) '거리 행렬';


$gdxIn %tspdata%

$load ii c

$gdxIn


$nrCities가 설정되지 않은 경우 $set nrCities 20

i(ii)$(ord(ii) < %nrCities%) = 예;


변수

x(ii,jj) '결정 변수 - 여행 구간'

z '객관적 변수';


바이너리 변수 x;

x.fx(ii,ii) = 0;


수식

목표 '총 비용'

rowsum(ii) '각 도시를 한 번만 떠나세요'

colsum(jj) '각 도시에 한 번만 도착';


* 할당 문제는 슬롯 사이트 추천의 완화입니다.

목적.. z =e= sum((i,j), c(i,j)*x(i,j));


rowsum(i).. sum(j, x(i,j)) =e= 1;


colsum(j)..sum(i, x(i,j)) =e= 1;


$cmax가 설정되지 않은 경우 $set cmax 2

컷 설정 /c0*c%cmax%/;


매개변수

acut(슬롯 사이트 추천ii,jj) '제약조건 행렬 자르기'

rhscut(슬롯 사이트 추천) '제한조건 rhs 자르기';


수식 sscut(슬롯 사이트 추천) 'sub_tour 제거 컷';

sscut(슬롯 사이트 추천).. sum((i,j), acut(슬롯 사이트 추천i,j)*x(i,j)) =l= rhscut(슬롯 사이트 추천);


cc(슬롯 사이트 추천) '이전 컷' 설정;

cc(잘라내기) = 아니요;


$if set cutdataexecute_load '%cutdata%', cc, acut, rhscut;


acut(슬롯 사이트 추천i,j)$(cc(슬롯 사이트 추천) 아님) = eps;


rhscut(컷)$(cc(컷) 아님) = 카드(ii);


모델 할당 /all/;


옵션 optcr=0;

"""


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)


    # 모델 인스턴스에 추가할 수 있는 컷 수

cuts_per_round = 10


    # 현재 컷

cur_cut = 0


    # 현재 모델 인스턴스에 대한 절단 제한(cmax = cur_cut + cuts_per_round)

cmax = 0


    # 생성된 모든 컷을 수집하는 데 사용되는 데이터베이스

cut_data = ws.add_database()

cc = cut_data.add_set("cc", 1)

acut = cut_data.add_parameter("어컷", 3)

rhscut = cut_data.add_parameter("rhscut", 1)


n = []# 도시 목록(i1, i2, i3, ...)

    그동안 참:

        # 절단 한계에 도달하면 새 모델 인스턴스를 생성

        ifcur_cut >= cmax:

cmax = cur_cut + cuts_per_round

cut_data.export()


            # 체크포인트 생성

작업 = ws.add_job_from_string(GAMS_MODEL)

cp = ws.add_checkpoint()

opt = ws.add_options()

opt.defines["nrcities"] ="20"

opt.defines["cmax"] = str(cmax - 1)

opt.defines["컷데이터"] = cut_data.name

opt.defines["tspdata"] = (

                '"'

+ os.path.join(ws.system_directory,"api파일", "데이터", "슬롯 사이트 추천gdx")

                + '"'

            )

job.run(gams_options=opt, checkpoint=cp)


            # 도시 목록을 한 번만 채웁니다.

            if 아님n:

                fori에job.out_db["나"]:

n += i.keys


            # 수정자 mi_acut 및 mi_rhscut을 사용하여 모델 인스턴스를 인스턴스화합니다.

mi = cp.add_modelinstance()

mi_acut = mi.sync_db.add_parameter("어컷", 3)

mi_rhscut = mi.sync_db.add_parameter("rhscut", 1)

mi.instantiate(

                "mip min z 사용 할당", [GamsModifier(mi_acut), GamsModifier(mi_rhscut)]

            )


        # 업데이트 유형 누적을 사용하여 모델 인스턴스를 해결하고 나중에 acut 및 rhscut 지우기

mi.solve(SymbolUpdateType.Accumulate)

mi.sync_db["어컷"].clear()

mi.sync_db["rhscut"].clear()


        # 솔루션에서 그래프 정보 수집

그래프 =

        fori, j에제품(n, n):

            ifmi.sync_db["x"].find_record([i, j]).level > 0.5:

그래프[i] = j


        # 모든 하위 투어를 찾고 acut 및 rhscut을 수정하여 필요한 컷을 추가

not_visited = 목록(n)# 복사하기

        그동안방문하지 않음:

i = 방문하지 않음[0]

sub_tour = [i]

            그동안그래프[i] != 방문하지 않음[0]:

i = 그래프[i]

sub_tour.append(i)

not_visited = 목록(필터(람다x:x아님 insub_tour, not_visited))


            # 두 데이터베이스에 컷 추가(cut_data, mi.sync_db)

            fori, j에제품(sub_tour, sub_tour):

acut.add_record([f"ccur_cut", i, j]).값 = 1

mi_acut.add_record([f"ccur_cut", i, j]).값 = 1

rhscut.add_record(f"ccur_cut").value = len(sub_tour) - 0.5

mi_rhscut.add_record([f"ccur_cut"]).value = len(sub_tour) - 0.5

cc.add_record(f"ccur_cut")

cur_cut += 1


        # 하위 투어에 모든 도시가 포함되어 있으면 루프를 종료하고 솔루션을 인쇄합니다.

        iflen(sub_tour) == len(n):

인쇄(f"z=mi.sync_db['z'].first_record().level")

인쇄("sub_tour: ")

            foriinsub_tour:

인쇄(f"i -> ", 끝="")

인쇄(sub_tour[0])

            휴식