modeling/PublicTransit2.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Public Transit Play II \n",
    "\n",
    "Find \"The shortest route to iterate all subway lines in ShangHai.\"\n",
    "\n",
    "Answers can be found on Zhihu: https://www.zhihu.com/question/629922852/answer/54798222956\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import gurobipy as grb"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 185,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Station: 莘庄, Connections: {'外环路': 1303, '春申路': 1588}\n",
      "Station: 外环路, Connections: {'莘庄': 1303, '莲花路': 1459}\n",
      "Station: 莲花路, Connections: {'外环路': 1459, '锦江乐园': 1633}\n",
      "Station: 锦江乐园, Connections: {'莲花路': 1633, '上海南站': 2089}\n",
      "Station: 上海南站, Connections: {'锦江乐园': 2089, '漕宝路': 1661, '石龙路': 1314, '华东理工大学': 1623, '桂林公园': 1828}\n"
     ]
    }
   ],
   "source": [
    "def read_metro_data(file_path):\n",
    "    from collections import defaultdict\n",
    "    # 邻接表表示的地铁网络图\n",
    "    metro_network = defaultdict(dict)\n",
    "    line_routes = dict()\n",
    "    station_name_to_idx = dict()\n",
    "    idx_to_station_name = dict()\n",
    "    cnt_idx = 1\n",
    "    # transship = dict()\n",
    "    with open(file_path, 'r', encoding='utf-8') as file:\n",
    "        current_line = None\n",
    "        for line in file:\n",
    "            line = line.strip()\n",
    "            if line.startswith(\"Line\"):\n",
    "                # 读取线路名\n",
    "                current_line = line\n",
    "                line_routes[line] = list()\n",
    "            elif line:\n",
    "                # 读取站点和距离信息\n",
    "                station1, station2, distance = line.split(\",\")\n",
    "                if station1 not in station_name_to_idx:\n",
    "                    station_name_to_idx[station1] = cnt_idx\n",
    "                    idx_to_station_name[cnt_idx] = station1\n",
    "                    cnt_idx += 1\n",
    "                if station2 not in station_name_to_idx:\n",
    "                    station_name_to_idx[station2] = cnt_idx\n",
    "                    idx_to_station_name[cnt_idx] = station2\n",
    "                    cnt_idx += 1\n",
    "                distance = int(distance)  # 转换距离为整数\n",
    "                metro_network[station1][station2] = distance\n",
    "                metro_network[station2][station1] = distance\n",
    "                # metro_network[station1][station2] = 1\n",
    "                # metro_network[station2][station1] = 1\n",
    "                line_routes[current_line].append(f\"{station1},{station2}\")\n",
    "                \n",
    "\n",
    "    return metro_network, line_routes, station_name_to_idx, idx_to_station_name\n",
    "\n",
    "# 示例调用\n",
    "file_path = \"./data/PublicTransit/SH.txt\"  # 替换为你的实际文件路径\n",
    "metro_network, line_routes, station_name_to_idx, idx_to_station_name = read_metro_data(file_path)\n",
    "\n",
    "# 打印部分结果查看结构\n",
    "for station, connections in list(metro_network.items())[:5]:\n",
    "    print(f\"Station: {station}, Connections: {connections}\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 186,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Gurobi Optimizer version 11.0.2 build v11.0.2rc0 (mac64[arm] - Darwin 23.5.0 23F79)\n",
      "\n",
      "CPU model: Apple M3 Pro\n",
      "Thread count: 11 physical cores, 11 logical processors, using up to 11 threads\n",
      "\n",
      "Optimize a model with 2203 rows, 2586 columns and 9276 nonzeros\n",
      "Model fingerprint: 0xb9763a16\n",
      "Variable types: 0 continuous, 2586 integer (2181 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 7e+01]\n",
      "  Objective range  [6e+02, 1e+04]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 7e+01]\n",
      "Presolve time: 0.01s\n",
      "Presolved: 2203 rows, 2586 columns, 9276 nonzeros\n",
      "Variable types: 0 continuous, 2586 integer (2181 binary)\n",
      "\n",
      "Root relaxation: objective 1.541800e+04, 132 iterations, 0.00 seconds (0.00 work units)\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "     0     0 15418.0000    0   56          - 15418.0000      -     -    0s\n",
      "     0     0 16538.0000    0   60          - 16538.0000      -     -    0s\n",
      "     0     0 17443.4429    0   79          - 17443.4429      -     -    0s\n",
      "     0     0 18259.2500    0   93          - 18259.2500      -     -    0s\n",
      "     0     0 18703.7500    0   93          - 18703.7500      -     -    0s\n",
      "     0     0 18840.7500    0   97          - 18840.7500      -     -    0s\n",
      "     0     0 18941.7500    0   92          - 18941.7500      -     -    0s\n",
      "     0     0 19025.7500    0   95          - 19025.7500      -     -    0s\n",
      "     0     0 19102.7500    0  115          - 19102.7500      -     -    0s\n",
      "     0     0 19178.2500    0  102          - 19178.2500      -     -    0s\n",
      "     0     0 19221.6250    0  137          - 19221.6250      -     -    0s\n",
      "     0     0 19243.6667    0  145          - 19243.6667      -     -    0s\n",
      "     0     0 19273.2500    0  142          - 19273.2500      -     -    0s\n",
      "     0     0 19302.8542    0  162          - 19302.8542      -     -    0s\n",
      "     0     0 19333.2500    0  120          - 19333.2500      -     -    0s\n",
      "     0     0 19338.6181    0  157          - 19338.6181      -     -    0s\n",
      "     0     0 19360.5926    0  131          - 19360.5926      -     -    0s\n",
      "     0     0 19362.6852    0  132          - 19362.6852      -     -    0s\n",
      "     0     0 19362.9017    0  164          - 19362.9017      -     -    0s\n",
      "     0     0 21243.0000    0  150          - 21243.0000      -     -    0s\n",
      "     0     0 21777.4352    0  167          - 21777.4352      -     -    0s\n",
      "     0     0 22053.1378    0  199          - 22053.1378      -     -    0s\n",
      "     0     0 22348.1708    0  131          - 22348.1708      -     -    0s\n",
      "     0     0 22424.0613    0  206          - 22424.0613      -     -    0s\n",
      "     0     0 22482.7563    0  208          - 22482.7563      -     -    0s\n",
      "     0     0 22531.1968    0  207          - 22531.1968      -     -    0s\n",
      "     0     0 22544.7589    0  211          - 22544.7589      -     -    0s\n",
      "     0     0 22549.6220    0  234          - 22549.6220      -     -    0s\n",
      "     0     0 22573.5724    0  213          - 22573.5724      -     -    0s\n",
      "     0     0 22587.2985    0  242          - 22587.2985      -     -    0s\n",
      "     0     0 22601.0874    0  204          - 22601.0874      -     -    0s\n",
      "     0     0 22605.9471    0  241          - 22605.9471      -     -    0s\n",
      "     0     0 22611.4062    0  244          - 22611.4062      -     -    0s\n",
      "     0     0 22629.3305    0  227          - 22629.3305      -     -    0s\n",
      "     0     0 22630.5189    0  252          - 22630.5189      -     -    0s\n",
      "     0     0 24135.4422    0  259          - 24135.4422      -     -    0s\n",
      "     0     0 24618.8711    0  290          - 24618.8711      -     -    0s\n",
      "     0     0 25006.1518    0  324          - 25006.1518      -     -    0s\n",
      "     0     0 25041.1037    0  363          - 25041.1037      -     -    0s\n",
      "     0     0 25106.9477    0  368          - 25106.9477      -     -    0s\n",
      "     0     0 25115.1977    0  370          - 25115.1977      -     -    0s\n",
      "     0     0 25119.3908    0  386          - 25119.3908      -     -    0s\n",
      "     0     0 26222.9659    0  439          - 26222.9659      -     -    0s\n",
      "     0     0 26336.8880    0  355          - 26336.8880      -     -    0s\n",
      "     0     0 26422.1752    0  429          - 26422.1752      -     -    0s\n",
      "     0     0 26459.6391    0  455          - 26459.6391      -     -    0s\n",
      "     0     0 26462.2812    0  479          - 26462.2812      -     -    0s\n",
      "     0     0 26656.5322    0  570          - 26656.5322      -     -    0s\n",
      "     0     0 26672.9429    0  556          - 26672.9429      -     -    0s\n",
      "     0     0 26675.8522    0  550          - 26675.8522      -     -    0s\n",
      "     0     0 26865.7617    0  402          - 26865.7617      -     -    0s\n",
      "     0     0 26904.7921    0  385          - 26904.7921      -     -    0s\n",
      "     0     0 26910.7204    0  370          - 26910.7204      -     -    0s\n",
      "     0     0 27071.1648    0  433          - 27071.1648      -     -    0s\n",
      "     0     0 27102.0250    0  523          - 27102.0250      -     -    0s\n",
      "     0     0 27102.6794    0  514          - 27102.6794      -     -    0s\n",
      "     0     0 27238.4953    0  408          - 27238.4953      -     -    0s\n",
      "     0     0 27247.2433    0  437          - 27247.2433      -     -    0s\n",
      "     0     0 27280.5598    0  449          - 27280.5598      -     -    0s\n",
      "     0     0 27285.9967    0  449          - 27285.9967      -     -    0s\n",
      "     0     2 27286.5328    0  442          - 27286.5328      -     -    0s\n",
      "*  362   282              36    92157.000000 27315.0045  70.4%  81.6    1s\n",
      "H  406   304                    91522.000000 27315.0045  70.2%  76.9    1s\n",
      "H  406   304                    90657.000000 27315.0045  69.9%  76.9    1s\n",
      "H  412   304                    86860.000000 27315.0045  68.6%  75.9    1s\n",
      "H  476   347                    83517.000000 27315.0045  67.3%  71.4    1s\n",
      "H  545   360                    80843.000000 27315.0045  66.2%  69.5    1s\n",
      "H  585   384                    74306.000000 27315.0045  63.2%  72.1    1s\n",
      "H  649   410                    74190.000000 27315.0045  63.2%  78.7    2s\n",
      "H  651   410                    74164.000000 27315.0045  63.2%  78.8    2s\n",
      "H  720   419                    71158.000000 27319.5394  61.6%  85.9    2s\n",
      "H  741   425                    68797.000000 27319.5394  60.3%  88.2    2s\n",
      "H  800   435                    67836.000000 27319.5394  59.7%  90.2    2s\n",
      "H  820   435                    66123.000000 27319.5394  58.7%  89.5    2s\n",
      "H  869   449                    65877.000000 27319.5394  58.5%  92.2    2s\n",
      "  1550   559 59091.8333   23   52 65877.0000 27632.9772  58.1%  99.0    5s\n",
      "  5870   733 infeasible   36      65877.0000 54563.2079  17.2%  86.5   10s\n",
      " 18639  2710 59808.3803   45   58 65877.0000 58744.1604  10.8%  59.3   15s\n",
      " 34166  5922 60648.0126   39  449 65877.0000 60583.8768  8.03%  45.1   22s\n",
      " 34207  5949 60774.0598   47  662 65877.0000 60583.8768  8.03%  45.0   25s\n",
      " 34612  5969 60583.8768   43  106 65877.0000 60583.8768  8.03%  47.7   30s\n",
      " 38522  5458 61365.7949   49  213 65877.0000 60583.8768  8.03%  53.2   35s\n",
      " 46143  3943 62547.8333   49  125 65877.0000 62009.7105  5.87%  54.8   40s\n",
      " 57386  1698 65439.0166   52   97 65877.0000 63757.7077  3.22%  53.1   45s\n",
      "\n",
      "Cutting planes:\n",
      "  Learned: 187\n",
      "  Gomory: 67\n",
      "  Cover: 76\n",
      "  Implied bound: 203\n",
      "  Projected implied bound: 12\n",
      "  Clique: 156\n",
      "  MIR: 47\n",
      "  StrongCG: 8\n",
      "  Flow cover: 178\n",
      "  GUB cover: 2\n",
      "  Zero half: 138\n",
      "  RLT: 11\n",
      "  Relax-and-lift: 51\n",
      "\n",
      "Explored 63235 nodes (3233614 simplex iterations) in 46.33 seconds (87.57 work units)\n",
      "Thread count was 11 (of 11 available processors)\n",
      "\n",
      "Solution count 10: 65877 66123 67836 ... 83517\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 6.587700000000e+04, best bound 6.587700000000e+04, gap 0.0000%\n",
      "BEST VALUE:65877.0\n",
      "春申路,Dummy\n",
      "Dummy,国家会展中心L17\n",
      "莘庄,春申路\n",
      "外环路,莘庄\n",
      "莲花路,外环路\n",
      "锦江乐园,莲花路\n",
      "上海南站,锦江乐园\n",
      "漕宝路,桂林公园\n",
      "上海体育馆,漕宝路\n",
      "徐家汇,交通大学\n",
      "陕西南路,一大会址·黄陂南路\n",
      "一大会址·黄陂南路,大世界\n",
      "虹桥火车站,虹桥2号航站楼\n",
      "虹桥2号航站楼,虹桥1号航站楼\n",
      "陆家嘴,浦东南路\n",
      "浦东南路,世纪大道\n",
      "世纪大道,源深体育中心\n",
      "龙阳路,华夏中路\n",
      "宜山路,徐家汇\n",
      "虹桥路,宜山路\n",
      "东安路,上海体育场\n",
      "上海体育场,上海体育馆\n",
      "民生路,杨高中路\n",
      "源深体育中心,民生路\n",
      "东明路,成山路\n",
      "龙华中路,东安路\n",
      "后滩,龙华中路\n",
      "长清路,后滩\n",
      "耀华路,长清路\n",
      "大世界,豫园\n",
      "成山路,耀华路\n",
      "杨高中路,迎春路\n",
      "龙溪路,水城路\n",
      "虹桥1号航站楼,上海动物园\n",
      "上海动物园,龙溪路\n",
      "水城路,伊犁路\n",
      "伊犁路,宋园路\n",
      "宋园路,虹桥路\n",
      "交通大学,上海图书馆\n",
      "上海图书馆,陕西南路\n",
      "豫园,陆家嘴\n",
      "桂林公园,上海南站\n",
      "华鹏路,东明路\n",
      "下南路,华鹏路\n",
      "北蔡,下南路\n",
      "陈春路,北蔡\n",
      "莲溪路,陈春路\n",
      "华夏中路,莲溪路\n",
      "国家会展中心L17,虹桥火车站\n",
      "迎春路,龙阳路\n"
     ]
    }
   ],
   "source": [
    "ptn = grb.Model(\"SH_PTN\")\n",
    "idx_to_station_name[0] = \"Dummy\"\n",
    "station_name_to_idx[\"Dummy\"] = 0\n",
    "x = {}\n",
    "y = {}\n",
    "n = {}\n",
    "N = len(station_name_to_idx.keys())\n",
    "M = 70\n",
    "out_flow = dict()\n",
    "\n",
    "for i in range(1, N):\n",
    "    x[0, i] = ptn.addVar(vtype=grb.GRB.BINARY,obj = 0, name=f\"x_Dummy_{idx_to_station_name[i]}\")\n",
    "    x[i, 0] = ptn.addVar(vtype=grb.GRB.BINARY,obj = 0, name=f\"x_{idx_to_station_name[i]}_Dummy\")\n",
    "    y[i] = ptn.addVar(vtype=grb.GRB.BINARY,obj = 0, name=f\"y_{idx_to_station_name[i]}\")\n",
    "    n[i] = ptn.addVar(vtype=grb.GRB.INTEGER,obj = 0, lb = 0, name=f\"n_{idx_to_station_name[i]}\")\n",
    "    out_flow[i] = []\n",
    "\n",
    "for station, connections in list(metro_network.items()):\n",
    "    idx_1 = station_name_to_idx[station]\n",
    "    for next_station, dist in connections.items():\n",
    "        idx_2 = station_name_to_idx[next_station]\n",
    "        x[idx_1, idx_2] = ptn.addVar(vtype=grb.GRB.BINARY,obj = dist, name=f\"x_{station}_{next_station}\")\n",
    "        out_flow[idx_1].append(idx_2)\n",
    "# 针对换乘的情况, 给出带转运的解\n",
    "\n",
    "\n",
    "# 流出的\n",
    "for k, v in out_flow.items():\n",
    "\n",
    "    v_new = v + [0]\n",
    "    ptn.addConstr(grb.quicksum(x[k, j] for j in v_new) == y[k], name=f\"out_flow_{idx_to_station_name[k]}\")\n",
    "    ptn.addConstr(grb.quicksum(x[i, k] for i in v_new) == y[k], name=f\"in_flow_{idx_to_station_name[k]}\")\n",
    "\n",
    "\n",
    "ptn.addConstr(grb.quicksum(x[i, 0] for i in range(1, N)) == 1, name=f\"Dummy_balance1\")\n",
    "ptn.addConstr(grb.quicksum(x[0, i] for i in range(1, N)) == 1, name=f\"Dummy_balance2\")\n",
    "\n",
    "for i in range(1, N):\n",
    "    ptn.addConstr(n[i] <= M * y[i], name = f\"N_{i}\")\n",
    "\n",
    "for station, connections in list(metro_network.items()):\n",
    "    idx_1 = station_name_to_idx[station]\n",
    "    for next_station, dist in connections.items():\n",
    "        idx_2 = station_name_to_idx[next_station]\n",
    "        ptn.addConstr(n[idx_2] - n[idx_1] >= 1 + M * (x[idx_1, idx_2] - 1), name = f\"n_{idx_1}_{idx_2}\")\n",
    "\n",
    "for k, edges in line_routes.items():\n",
    "    left_edges = []\n",
    "    for edge in edges:\n",
    "        o_, d_ = edge.split(\",\")\n",
    "        o_idx, d_idx = station_name_to_idx[o_], station_name_to_idx[d_]\n",
    "        left_edges.append([o_idx, d_idx])\n",
    "        left_edges.append([d_idx, o_idx])\n",
    "    ptn.addConstr(grb.quicksum(x[i, j] for i, j in left_edges) >= 1, name = f\"line_{k}\")\n",
    "\n",
    "# 处理 Line 3 + Line 4 的情况\n",
    "edges_3, edges_4 = line_routes[\"Line 3\"], line_routes[\"Line 4\"]\n",
    "new_edges = list(set(edges_3 + edges_4))\n",
    "left_edges = []\n",
    "for edge in new_edges:\n",
    "    o_, d_ = edge.split(\",\")\n",
    "    o_idx, d_idx = station_name_to_idx[o_], station_name_to_idx[d_]\n",
    "    left_edges.append([o_idx, d_idx])\n",
    "    left_edges.append([d_idx, o_idx])\n",
    "ptn.addConstr(grb.quicksum(x[i, j] for i, j in left_edges) >= 2, name = f\"spec_line_3/4\")\n",
    "\n",
    "edges_2, edges_10 = line_routes[\"Line 2\"], line_routes[\"Line 10\"]\n",
    "new_edges = list(set(edges_2 + edges_10))\n",
    "left_edges = []\n",
    "for edge in new_edges:\n",
    "    o_, d_ = edge.split(\",\")\n",
    "    o_idx, d_idx = station_name_to_idx[o_], station_name_to_idx[d_]\n",
    "    left_edges.append([o_idx, d_idx])\n",
    "    left_edges.append([d_idx, o_idx])\n",
    "ptn.addConstr(grb.quicksum(x[i, j] for i, j in left_edges) >= 2, name = f\"spec_line_10/2\")\n",
    "\n",
    "\n",
    "ptn.modelSense = grb.GRB.MINIMIZE\n",
    "ptn.optimize()\n",
    "if (ptn.status == grb.GRB.status.OPTIMAL):\n",
    "    solType = 'IP_Optimal'\n",
    "    ofv = ptn.getObjective().getValue()\n",
    "    print(f\"BEST VALUE:{ofv}\")\n",
    "    for i, j in x:\n",
    "        if (x[i, j].x > 0.5):\n",
    "            print(f\"{idx_to_station_name[i]},{idx_to_station_name[j]}\")\n",
    "\n",
    "    gap = 0\n",
    "    lb = ofv\n",
    "    ub = ofv\n",
    "    runtime = ptn.Runtime"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 【所需边数最小的（44条边）】\n",
    "# 国家会展中心L17,虹桥火车站,2523\n",
    "# 虹桥火车站,虹桥2号航站楼,696\n",
    "# 虹桥2号航站楼,淞虹路,5921\n",
    "# 淞虹路,北新泾,1367\n",
    "# 北新泾,威宁路,1260\n",
    "# 威宁路,娄山关路,1661\n",
    "# 娄山关路,长风公园,1487\n",
    "# 长风公园,大渡河路,784\n",
    "# 大渡河路,金沙江路,1676\n",
    "# 金沙江路,曹杨路,903\n",
    "# 曹杨路,武宁,1050\n",
    "# 武宁,武定路,970\n",
    "# 武定路,静安寺,1030\n",
    "# 静安寺,常熟路,1091\n",
    "# 常熟路,陕西南路,939\n",
    "# 陕西南路,一大会址·新天地,1591\n",
    "# 一大会址·新天地,老西门,880\n",
    "# 老西门,陆家浜路,825\n",
    "# 陆家浜路,小南门,1475\n",
    "# 小南门,商城路,2398\n",
    "# 商城路,世纪大道,979\n",
    "# 世纪大道,源深体育中心,890\n",
    "# 源深体育中心,民生路,907\n",
    "# 民生路,杨高中路,928\n",
    "# 杨高中路,迎春路,633\n",
    "# 迎春路,龙阳路,2275\n",
    "# 龙阳路,华夏中路,4355\n",
    "# 华夏中路,莲溪路,1589\n",
    "# 莲溪路,御桥,1342\n",
    "# 御桥,浦三路,3221\n",
    "# 浦三路,三林东,1613\n",
    "# 三林东,三林,1126\n",
    "# 三林,东方体育中心,3516\n",
    "# 东方体育中心,龙耀路,2164\n",
    "# 龙耀路,云锦路,934\n",
    "# 云锦路,龙华,945\n",
    "# 龙华,龙漕路,1124\n",
    "# 龙漕路,石龙路,1484\n",
    "# 石龙路,上海南站,1314\n",
    "# 上海南站,锦江乐园,2089\n",
    "# 锦江乐园,莲花路,1633\n",
    "# 莲花路,外环路,1459\n",
    "# 外环路,莘庄,1303\n",
    "# 莘庄,春申路,1588\n",
    "# 69938 \n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "| 序号 | 线路 | 起点 | 终点| 距离 | \n",
    "|:---:| :---: |:----: | :---:| :---: |\n",
    "| 1 | L17 | 国家会展中心L17 | 虹桥火车站 | 2523 |\n",
    "| 2 | L2 | 虹桥火车站 | 虹桥2号航站楼 | 696 | \n",
    "| 3 | L10 | 虹桥2号航站楼 | 虹桥1号航站楼 | 2026 |\n",
    "|4 | L10  | 虹桥1号航站楼 | 上海动物园 | 2016 |\n",
    "| 5 | L10 | 上海动物园 | 龙溪路 | 1244| \n",
    "|6 | L10| 龙溪路 | 水城路 | 1303 |\n",
    "| 7| L10| 水城路 | 伊犁路 | 1167 |\n",
    "| 8| L10 | 伊犁路 | 宋园路 | 800 |\n",
    "| 9| L10| 宋园路 | 虹桥路 | 1053 |\n",
    "| 10 | L3 | 虹桥路 | 宜山路 | 1323 |\n",
    "| 11| L9 | 宜山路 | 徐家汇 | 1505 |\n",
    "| 12| L11 | 徐家汇 | 交通大学 | 905 |\n",
    "| 13| L11 | 交通大学 | 上海图书馆 | 1145 |\n",
    "| 14| L11 | 上海图书馆 | 陕西南路 | 1618 |\n",
    "| 15 | L1 | 陕西南路 | 一大会址·黄陂南路 | 1407 |\n",
    "|16 |  L14 | 一大会址·黄陂南路 | 大世界 | 710 | \n",
    "| 17| L14 | 大世界 | 豫园 | 740 |\n",
    "|8| L14 | 豫园 | 陆家嘴 | 1600 |\n",
    "|9| L14 | 陆家嘴 | 浦东南路 | 1010 |\n",
    "|20| L2 | 浦东南路 | 世纪大道 | 1245 |\n",
    "|21| L6 | 世纪大道 | 源深体育中心 | 890 |\n",
    "|22| L6| 源深体育中心 | 民生路 | 907 |\n",
    "|23| L18 | 民生路 | 杨高中路 | 928 |\n",
    "|24| L18| 杨高中路 | 迎春路 | 633 |\n",
    "|25| L18 | 迎春路 | 龙阳路 | 2275 |\n",
    "|26| L16 | 龙阳路 | 华夏中路 | 4355 |\n",
    "|27| L13 | 华夏中路 | 莲溪路 | 1589 |\n",
    "|28| L13| 莲溪路 | 陈春路 | 1239 |\n",
    "|29 |L13 | 陈春路 | 北蔡 | 980 |\n",
    "|30 | L13| 北蔡 | 下南路 | 1128 |\n",
    "|31|L13 | 下南路 | 华鹏路 | 1361 |\n",
    "|32| L13| 华鹏路 | 东明路 | 1464 |\n",
    "|33 | L13| 东明路 | 成山路 | 1397 |\n",
    "|34 |L8 | 成山路 | 耀华路 | 908 |\n",
    "|35| L7 | 耀华路 | 长清路 | 926 |\n",
    "|36| L7| 长清路 | 后滩 | 1263 |\n",
    "|37| L7| 后滩 | 龙华中路 | 2247 |\n",
    "|38| L7| 龙华中路 | 东安路 | 728 |\n",
    "|39| L4| 东安路 | 上海体育场 | 1220 |\n",
    "|40 | L4| 上海体育场 | 上海体育馆 | 685 |\n",
    "|41 |L1 | 上海体育馆 | 漕宝路 | 1534 |\n",
    "|42 |L12 | 漕宝路 | 桂林公园 | 1284 |\n",
    "|43 |L15 | 桂林公园 | 上海南站 | 1828 |\n",
    "|44 |L1 | 上海南站 | 锦江乐园 | 2089 |\n",
    "|45 |L1 | 锦江乐园 | 莲花路 | 1633 |\n",
    "|46 |L1 | 莲花路 | 外环路 | 1459 |\n",
    "|47 |L1 | 外环路 | 莘庄 | 1303 |\n",
    "|48 |L5 | 莘庄 | 春申路 | 1588 |\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "py310_arm",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.14"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}