/// <summary> /// 构建列车弧 /// </summary> /// <param name="arcCount"></param> private void create_Line_T_S_Arcs(ref int arcCount) { for (int i = 0; i < T_S_NodeList.Count - 1; i++) { for (int j = i + 1; j < i + 2 /*T_S_NodeList.Count*/; j++) { // 因为在创建并添加点的时候,同一Line的出发点和到达点是先后创建添加的,所以内层不用循环 T_S_Node startPoint = T_S_NodeList[i]; T_S_Node endPoint = T_S_NodeList[j]; if (startPoint.LineID == endPoint.LineID && startPoint.PointType == 1 && endPoint.PointType == 2 && ((startPoint.TimeCode < 1440 && endPoint.TimeCode < 1440) || (startPoint.TimeCode >= 1440 && endPoint.TimeCode >= 1440))) { T_S_Arc t_s_arc = new T_S_Arc(); t_s_arc.Init_db(startPoint, endPoint, arcCount, startPoint.LineID, startPoint.RountingID, startPoint.TrainCode, startPoint.ID, endPoint.ID, startPoint.StaCode, endPoint.StaCode, startPoint.TimeCode, endPoint.TimeCode, 1);//1——运行弧 T_S_ArcList.Add(t_s_arc); startPoint.Out_T_S_ArcList.Add(t_s_arc); arcCount++; } } } }
bool FeasibleDayCrewTime(T_S_S_Node curStateNode, T_S_Arc arc, ref Resource extendResource) { if (arc.ArcType == 31 || arc.ArcType == 33) { return(true); } int dayCrewTime_accumu = curStateNode.Resources.DayCrewTime_accumu; if (arc.ArcType != 32) { dayCrewTime_accumu += arc.Len; } if (dayCrewTime_accumu > BasicTimeRules.MaxDayCrewTime) { // 任何时候,dayCrewTime都不允许超过最大值 return(false); } else if (arc.ArcType == 32 && (dayCrewTime_accumu < BasicTimeRules.MinDayCrewTime || dayCrewTime_accumu > BasicTimeRules.MaxDayCrewTime)) { // 如果退乘时,dayCrewTime not in dayCrewTimeWindow [min, max],不可行 return(false); } extendResource.SetDayCrewTime(dayCrewTime_accumu); return(true); }
int getRealMealSlot(T_S_Arc arc, int mealwindow_begin, int mealwindow_end) { int slot_end = Math.Min(arc.EndTimeCode, mealwindow_end); int slot_begin = Math.Max(arc.StartTimeCode, mealwindow_begin); return(slot_end - slot_begin); }
// 处理用餐约束 void handleMealConstraint(T_S_S_Node curStateNode, T_S_Arc arc, ref T_S_S_Node extendStateNode) { if (withinMealWindow(arc, STLunch, ETLunch, minMealTime)) { if (/*curStateNode.Resources.LunchStatus < 2 &&*/ getRealMealSlot(arc, STLunch, ETLunch) < minMealTime) { // 还未用过午餐,且当前弧应该用餐但无法用餐 // 则施加惩罚 extendStateNode.PenaltyMealViolate = arcPenaltyMealViolate; extendStateNode.Resources.SetLunchStatus(curStateNode.Resources.LunchStatus); } else { extendStateNode.Resources.SetLunchStatus(2); } } else if (withinMealWindow(arc, STDinner, ETDinner, minMealTime)) { if (/*curStateNode.Resources.DinnerStatus < 2 &&*/ getRealMealSlot(arc, STDinner, ETDinner) < minMealTime) { // 还未用过晚餐,且当前弧应该用餐但无法用餐 // 则施加惩罚 extendStateNode.PenaltyMealViolate = arcPenaltyMealViolate; extendStateNode.Resources.SetDinnerStatus(curStateNode.Resources.DinnerStatus); } else { extendStateNode.Resources.SetDinnerStatus(2); } } else { extendStateNode.Resources.SetLunchStatus(curStateNode.Resources.LunchStatus); extendStateNode.Resources.SetDinnerStatus(curStateNode.Resources.DinnerStatus); } }
/// <summary> /// arc为运行弧时使用 /// 判断是否当前状态点(表示的部分路径)已经访问过该Line /// 避免同一Line被一条路径多次访问 /// </summary> /// <param name="curStateNode"></param> /// <param name="arc"></param> /// <returns></returns> bool visitedCurLine(T_S_S_Node curStateNode, T_S_Arc arc) { return(curStateNode.Resources.LinesVisitedArray[arc.LineID - 1] > 0); }
// end 处理用餐约束 // 更新时间属性值,不涉及用餐约束 bool FeasibleDriveTime(T_S_S_Node curStateNode, T_S_Arc arc, ref Resource extendResource) { if (arc.ArcType == 1) { //运行弧 if (!visitedCurLine(curStateNode, arc)) { int driveTime_accumu = curStateNode.Resources.DriveTime_accumu + arc.Len; if (driveTime_accumu <= BasicTimeRules.MaxDriveTime /*&& dayCrewTime_accumu <= BasicTimeRules.MaxDayCrewTime*/) { // NOTE:只检查驾驶时间,不检查dayCrewTime // 若执行该arc代表的Line后,驾驶时间 <= max,则可行 extendResource.SetDriveTime(driveTime_accumu); } else // 若驾驶时间大于最大允许值,则不可extend,cut掉 { return(false); } } } else if (arc.ArcType == 21) { //日间连接弧 int driveTime_accumu = curStateNode.Resources.DriveTime_accumu; if (BasicTimeRules.MinDriveTime <= driveTime_accumu && driveTime_accumu < BasicTimeRules.MaxDriveTime) { // 1.若驾驶时间 in [min. max),则可行,但要判断是否可以间休 // NOTE:这里将最小间休时间等于最小接续时间 // 那么下面这个判断始终是true,即只要是接续弧,就视为间休,驾驶时间清零 if (arc.Len >= BasicTimeRules.MinConnTime) { // 若弧长 >= 最小间休时长,视为在此间休 extendResource.SetDriveTime(0); } } else if (driveTime_accumu < BasicTimeRules.MinDriveTime) { // 2.若驾驶时间 < 最小驾驶时间,则不间休,视为接续 extendResource.SetDriveTime(driveTime_accumu + arc.Len); } else if (driveTime_accumu == BasicTimeRules.MaxDriveTime) { // 3.若驾驶时间 == 必须间休 if (arc.Len < BasicTimeRules.MinConnTime) { // 若弧长小于最小间休时长,则不可行 return(false); } // 否则,间休 extendResource.SetDriveTime(0); } } else if (arc.ArcType == 31 || arc.ArcType == 32 || arc.ArcType == 33) { //虚拟出乘弧, //虚拟退乘弧, //虚拟停驻弧 if (arc.ArcType != 32) //虚拟出乘,虚拟停驻,当前的状态点的驾驶时间都是0 { Debug.Assert(curStateNode.Resources.DriveTime_accumu == 0); Debug.Assert(curStateNode.Resources.DayCrewTime_accumu == 0); } extendResource.SetDriveTime(curStateNode.Resources.DriveTime_accumu); } else { Debug.Print("检查驾驶时间约束,出现未定义的T_S_Arc.ArcType:{0}\n", arc.ArcType); } return(true); }
bool withinMealWindow(T_S_Arc arc, int mealwindow_begin, int mealwindow_end, int minMealSpan) { return(arc.EndTimeCode - mealwindow_begin >= minMealSpan && mealwindow_end - arc.StartTimeCode >= minMealSpan); }
public void CreateT_S_S_Network_db(/*CrewRule crewRule*/) { Debug.Assert(BasicTimeRules.MinConnTime != -1); List <T_S_S_Node> allStateNodes = new List <T_S_S_Node>(); Queue <T_S_S_Node> queue = new Queue <T_S_S_Node>(); initVirStartNode(ref allStateNodes, ref queue); //List<T_S_S_Node> tempDNodeList = new List<T_S_S_Node>(); //Dictionary<int, T_S_S_Node> IDToNodeDict = new Dictionary<int, T_S_S_Node>(); //foreach (var startNode in allStateNodes) { // IDToNodeDict.Add(startNode.ID, startNode); //} int debug_count = 100000; int debug_iter = 0; int nodeCount = CrewBaseList.Count * 2; //其余时空状态点的编号从乘务基地虚拟起终点之后开始 T_S_S_Node top_node; while (queue.Count > 0) { //if (--debug_count < 0) { // debug_count = 100000; // debug_iter++; // Console.WriteLine("current Queue.Count is [{0}]", queue.Count); // Console.WriteLine("debug iter is [{0}]", debug_iter); //} top_node = queue.Dequeue(); //Console.WriteLine(Logger.stateNodeInfoToStr(top_node, OutputMode.console)); //if (debug_iter >= 10) { // if (debug_count < 3000) { // Console.WriteLine("cur debug count is [{0}]", debug_count); // } //} T_S_Node cur_super_point = top_node.SuperPoint; List <T_S_Arc> out_T_S_ArcList = cur_super_point.Out_T_S_ArcList; for (int i = 0; i < out_T_S_ArcList.Count; i++) { T_S_Arc cur_t_s_arc = out_T_S_ArcList[i]; // debug //if (top_node.PassLine.Contains(30) // && top_node.PassLine.Contains(143) // //&& top_node.PassLine.Contains(64) // //&& top_node.PassLine.Contains(157) // && cur_t_s_arc.EndTimeCode == 1029 // && cur_t_s_arc.ArcType == 21) { // Console.WriteLine("top_node.PenaltyMealViolate:" + top_node.PenaltyMealViolate); //} //end debug T_S_S_Node extendStateNode = new T_S_S_Node(); extendStateNode.Resources = new Resource(lineList.Count); if (!FeasibleDriveTime(top_node, cur_t_s_arc, ref extendStateNode.Resources)) { continue; } if (!FeasibleDayCrewTime(top_node, cur_t_s_arc, ref extendStateNode.Resources)) { continue; } // 交路退乘点必须与出乘点相同 (必须回基地的约束) if (cur_t_s_arc.ArcType == 32 && cur_t_s_arc.EndStaCode != top_node.OStation) { continue; } // 用餐软约束 // 必须是非运行弧,才可能可以用餐 if (cur_t_s_arc.ArcType != 1) { handleMealConstraint(top_node, cur_t_s_arc, ref extendStateNode); // NOTE:可以处理一下,若当前弧用餐成功,则标记为用餐弧 } else { extendStateNode.Resources.SetLunchStatus(top_node.Resources.LunchStatus); extendStateNode.Resources.SetDinnerStatus(top_node.Resources.DinnerStatus); } // extendStartNode基本属性赋值 extendStateNode.Init_db(nodeCount, top_node.ID, cur_t_s_arc.LineID, cur_t_s_arc.RountingID, cur_t_s_arc.TrainCode, cur_t_s_arc.EndStaCode, cur_t_s_arc.EndTimeCode, cur_t_s_arc.ArcType, top_node.OStation, cur_t_s_arc.EndPointID, cur_super_point.ID, top_node, cur_t_s_arc.EndPoint, cur_super_point); // 还有其他属性 // 1 途径运行线 extendStateNode.PassLine.AddRange(top_node.PassLine); if (cur_t_s_arc.ArcType == 1) { extendStateNode.PassLine.Add(cur_t_s_arc.LineID); extendStateNode.Resources.LinesVisitedArray[cur_t_s_arc.LineID - 1] += 1; } // 2 起终点的编号 if (cur_t_s_arc.ArcType == 31) { extendStateNode.Price = 1440; Debug.Assert(extendStateNode.PenaltyMealViolate == 0); } else if (cur_t_s_arc.ArcType == 32) { extendStateNode.ID = BaseStationToODIDDict[extendStateNode.StaCode][1]; extendStateNode.Price = 0; --nodeCount; // 终点有自己的编号,所以此次不用增加nodeCount; Debug.Assert(extendStateNode.PenaltyMealViolate == 0); } else if (cur_t_s_arc.ArcType == 33) { extendStateNode.ID = BaseStationToODIDDict[extendStateNode.StaCode][1]; extendStateNode.Price = virtualRoutingCost; --nodeCount; // 终点有自己的编号,所以此次不用增加nodeCount; Debug.Assert(extendStateNode.PenaltyMealViolate == 0); } else { extendStateNode.Price = cur_t_s_arc.Len;// + extendStateNode.PenaltyMealViolate; //后续对弧再赋值Penalty } // 比较标号优劣 //for (int stateIdx = cur_t_s_arc.EndPoint.StateNodeSet.Count; stateIdx >= 0; stateIdx--) { //} top_node.Out_T_S_S_NodeList.Add(extendStateNode); allStateNodes.Add(extendStateNode); queue.Enqueue(extendStateNode); ++nodeCount; //if (cur_t_s_arc.ArcType != 32 && cur_t_s_arc.ArcType != 33) { // IDToNodeDict.Add(extendStateNode.ID, extendStateNode); //} //else { // tempDNodeList.Add(extendStateNode); //} } // end of loop "out_T_S_ArcList" } // end of loop "while" // 找出所有可行的状态点(即可以到达虚拟终点的状态点) #region //基于回溯 //T_S_S_NodeList = new List<T_S_S_Node>(); //foreach (T_S_S_Node note in tempDNodeList) { // int d = note.PrevID; // do { // if (IDToNodeDict.ContainsKey(d)) { // IDToNodeDict[d].Isfeasible = 1; // T_S_S_NodeList.Add(IDToNodeDict[d]); // d = IDToNodeDict[d].PrevID; // } // } // while (d != -1); //} //T_S_S_NodeList.AddRange(tempDNodeList); //T_S_S_NodeList.Sort(T_S_S_Node_AscByTime.nodeAsc); //int virODNum = CrewBaseList.Count * 2; //int newIDCount = virODNum; //foreach (var node in T_S_S_NodeList) { // if (node.ID >= virODNum) { // node.ID = newIDCount; // ++newIDCount; // } //} //foreach (var node in T_S_S_NodeList) { // if (node.PrevID >= 0) { // node.PrevID = node.PrevNode.ID; // } //} #endregion keepAllFeasibleStateNode(allStateNodes); // 根据可行状态点建立时空状态弧 createT_S_S_Arcs(); } // end of function CreateT_S_S_Network_db
/// <summary> /// 构建出退乘弧,包括虚拟乘务组出退乘弧,即停驻弧 /// </summary> /// <param name="arcCount"></param> private void create_OnAndOff_T_S_Arcs(ref int arcCount) { #region 创建出退乘弧 // 只需要一次遍历,因为目的就是把位于基地所在站的stationNode与基地起终点相连 for (int i = 0; i < T_S_NodeList.Count; i++) { T_S_Node stationPoint = T_S_NodeList[i]; T_S_Node baseStartPoint = null; T_S_Node baseEndPoint = null; bool isBaseStation = false; foreach (var crewbase in CrewBaseList) { if (stationPoint.StaCode == crewbase.StaCode) { isBaseStation = true; baseStartPoint = T_S_NodeList[crewbase.VirStartPointID - 1]; baseEndPoint = T_S_NodeList[crewbase.VirEndPointID - 1]; break; } } if (!isBaseStation) { continue; } if (stationPoint.PointType == 1 && stationPoint.TimeCode < 1440) { T_S_Arc t_s_arc = new T_S_Arc(); t_s_arc.Init_db(baseStartPoint, stationPoint, arcCount, 0, stationPoint.RountingID, stationPoint.TrainCode, baseStartPoint.ID, stationPoint.ID, stationPoint.StaCode, stationPoint.StaCode, 0, stationPoint.TimeCode, 31);//31——出乘弧 T_S_ArcList.Add(t_s_arc); baseStartPoint.Out_T_S_ArcList.Add(t_s_arc); arcCount++; } if (stationPoint.PointType == 2) { T_S_Arc t_s_arc = new T_S_Arc(); t_s_arc.ArcType = 32;//32——退乘弧 t_s_arc.Init_db(stationPoint, baseEndPoint, arcCount, 0, 0, stationPoint.TrainCode, stationPoint.ID, baseEndPoint.ID, stationPoint.StaCode, stationPoint.StaCode, stationPoint.TimeCode, 1440, 32);//32——退乘弧 if (t_s_arc.StartTimeCode > 1440) { t_s_arc.EndTimeCode = 2880; } T_S_ArcList.Add(t_s_arc); stationPoint.Out_T_S_ArcList.Add(t_s_arc); arcCount++; } } // 创建虚拟停驻弧,即虚拟乘务组出乘弧 foreach (CrewBase crewbase in CrewBaseList) { T_S_Node startPoint = T_S_NodeList[crewbase.VirStartPointID - 1]; T_S_Node endPoint = T_S_NodeList[crewbase.VirEndPointID - 1]; T_S_Arc t_s_arc = new T_S_Arc(); t_s_arc.Init_db(startPoint, endPoint, arcCount, 0, 0, "", crewbase.VirStartPointID, crewbase.VirEndPointID, crewbase.StaCode, crewbase.StaCode, 0, 1440, //若是两天会有些变化 33); //33——虚拟乘务组出退乘弧 T_S_ArcList.Add(t_s_arc); startPoint.Out_T_S_ArcList.Add(t_s_arc); arcCount++; } #endregion }
/// <summary> /// 创建日间连接弧和夜间连接弧 /// </summary> /// <param name="arcCount"></param> /// <param name="minConnTime"></param> /// <param name="maxConnTime"></param> private void create_Connect_T_S_Arcs(ref int arcCount, int minConnTime, int maxConnTime) { Dictionary <string, List <T_S_Node> > StationToPointDict = new Dictionary <string, List <T_S_Node> >(); foreach (var point in T_S_NodeList) { if (point.PointType == 3 || point.PointType == 4) { continue; } if (!StationToPointDict.ContainsKey(point.StaCode)) { StationToPointDict.Add(point.StaCode, new List <T_S_Node>()); } StationToPointDict[point.StaCode].Add(point); } foreach (var pointset in StationToPointDict.Values) { pointset.Sort(T_S_Node_AscByTime.pointAsc); for (int i = 0; i < pointset.Count; i++) { T_S_Node startPoint = pointset[i]; if (startPoint.PointType == 1) { continue; // 连接弧的起点必须是到达点 } for (int j = i + 1; j < pointset.Count; j++) { // 因为排序了,所以后面的node j必然不早于前面的node i,所以可以直接 i + 1 T_S_Node endPoint = pointset[j]; if (endPoint.PointType == 2) { continue; // 连接弧的终点必须是出发点 } int len = endPoint.TimeCode - startPoint.TimeCode; if (len < minConnTime) { continue; } else if (len > maxConnTime) { break; } else { // 不考虑多天的情况,暂且仅考虑1天 // 故在此可以直接创建弧 T_S_Arc t_s_arc = new T_S_Arc(); t_s_arc.Init_db(startPoint, endPoint, arcCount, startPoint.LineID, startPoint.RountingID, startPoint.TrainCode, startPoint.ID, endPoint.ID, startPoint.StaCode, endPoint.StaCode, startPoint.TimeCode, endPoint.TimeCode, 21);//21——日间连接弧 T_S_ArcList.Add(t_s_arc); startPoint.Out_T_S_ArcList.Add(t_s_arc); arcCount++; } } } } }