# J03
# The set of nodes in the network
set N := 1 2 3 4 5 6 7 8 9;

# The set of modular sizes for structures and couplers
set W := 4 8 16 20 40 80;

# The set of OD pairs and 
# the number wavelengths required for each o-d pair
param: D: r :=
1	3	9
1	6	6
1	7	17
2	3	12
2	5	8
2	6	12
3	5	17
3	6	8
4	8	14
5	7	15
6	7	11
6	8	4
7	8	6
7	9	8
8	9	7
;

# The set of available structures
set S := 1 2 3 4 5;

# The set of available couplers
set C := 1 2 3 4 5 6;

# The set of routing paths
set P := 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17;

# The set of pairs of routing paths that intersect
set H :=
(1,2)     (4,17)    (10,12)   (2,14)    (7,13)    (2,17)    (9,14)    (4,15)
(2,9)     (7,9)     (13,17)   (4,8)     (9,11)    (4,11)    (12,16)   (6,15)
(3,17)    (8,13)    (1,13)    (5,10)    (11,13)   (6,10)    (2,7)     (8,10)
(4,16)    (10,11)   (2,13)    (7,12)    (2,3)     (7,16)    (3,13)    (9,16)
(7,8)     (13,16)   (4,7)     (9,10)    (2,16)    (9,13)    (4,13)    (13,14)
(8,11)    (1,9)     (5,7)     (10,16)   (4,10)    (12,13)   (6,13)
(9,17)    (2,12)    (7,11)    (16,17)   (6,8)     (2,6)     (8,9)
(13,15)   (4,6)     (8,17)    (1,17)    (7,15)    (3,9)     (9,15)
(1,3)     (5,6)     (10,15)   (2,15)    (9,12)    (4,12)    (12,17)
(2,11)    (7,10)    (14,17)   (4,9)     (11,15)   (6,11)    (2,8)
(4,5)     (8,15)    (1,14)    (6,7)     (2,4)     (7,17)    (3,14);

# Sets of routing paths serving each OD pair
set J[1,3] := 1;
set J[1,6] := 2;
set J[1,7] := 3;
set J[2,3] := 4;
set J[2,5] := 5;
set J[2,6] := 6;
set J[3,5] := 7;
set J[3,6] := 8;
set J[4,8] := 9;
set J[5,7] := 10;
set J[6,7] := 11 12;
set J[6,8] := 13;
set J[7,8] := 14;
set J[7,9] := 15 16;
set J[8,9] := 17;

# Sets of routing paths using each structure
set K[1] := 4 5 6 7 10;
set K[2] := 2 4 6 7 8 9 10 11 13 15 17;
set K[3] := 1 2 3 9 13 14 17;
set K[4] := 4 6 7 10;
set K[5] := 2 4 7 9 10 12 13 16 17;

# Sets of routing paths using each coupler
set L[1] := 4 6;
set L[2] := 4 6 7 10;
set L[3] := 2 9 13 17;
set L[4] := 4 7 10;
set L[5] := 4 7 10;
set L[6] := 2 9 13 17;

# The structure costs
param a:
	4	8	16	20	40	80	 :=
1	14	25	42	56	112	210
2	26	47	78	104	208	390
3	37	67	111	148	296	555
4	43	77	129	172	344	645
5	32	58	96	128	256	480
;

# The coupler costs
param f:
	4	8	16	20	40	80	 :=
1	5	9	15	20	40	74
2	6	10	17	23	47	89
3	7	12	20	27	55	104
4	10	17	29	39	79	149
5	7	12	20	27	55	104
6	8	14	23	31	63	119
;