POJ 1637 Sightseeing tour (SAP | Dinic 混合欧拉图的判断)
SAP 判断 poj 混合 欧拉
2023-09-27 14:22:23 时间
Sightseeing tour
Description The city executive board in Lund wants to construct a sightseeing tour by bus in Lund, so that tourists can see every corner of the beautiful city. They want to construct the tour so that every street in the city is visited exactly once. The bus should also start and end at the same junction. As in any city, the streets are either one-way or two-way, traffic rules that must be obeyed by the tour bus. Help the executive board and determine if it's possible to construct a sightseeing tour under these constraints.
Input On the first line of the input is a single positive integer n, telling the number of test scenarios to follow. Each scenario begins with a line containing two positive integers m and s, 1 <= m <= 200,1 <= s <= 1000 being the number of junctions and streets, respectively. The following s lines contain the streets. Each street is described with three integers, xi, yi, and di, 1 <= xi,yi <= m, 0 <= di <= 1, where xi and yi are the junctions connected by a street. If di=1, then the street is a one-way street (going from xi to yi), otherwise it's a two-way street. You may assume that there exists a junction from where all other junctions can be reached.
Output For each scenario, output one line containing the text "possible" or "impossible", whether or not it's possible to construct a sightseeing tour.
Sample Input 4 5 8 2 1 0 1 3 0 4 1 1 1 5 0 5 4 1 3 4 0 4 2 1 2 2 0 4 4 1 2 1 2 3 0 3 4 0 1 4 1 3 3 1 2 0 2 3 0 3 2 0 3 4 1 2 0 2 3 1 1 2 0 3 2 0 Sample Output possible impossible impossible possible Source |
给出一张混合图(有有向边,也有无向边),判断是否存在欧拉回路。
首先是对图中的无向边随意定一个方向,然后统计每个点的入度(indeg)和出度(outdeg),如果(indeg - outdeg)是奇数的话,一定不存在欧拉回路;
如果所有点的入度和出度之差都是偶数,那么就开始网络流构图:
1,对于有向边,舍弃;对于无向边,就按照最开始指定的方向建权值为 1 的边;
2,对于入度小于出度的点,从源点连一条到它的边,权值为(outdeg - indeg)/2;出度小于入度的点,连一条它到汇点的权值为(indeg - outdeg)/2 的边;
构图完成,如果满流(求出的最大流值 == 和汇点所有连边的权值之和),那么存在欧拉回路,否则不存在。
另附一个讲解欧拉图不错的博客:http://www.cnblogs.com/destinydesigner/archive/2009/09/28/1575674.html
SAP果然快,0ms:
#include<iostream> #include<cstdio> #include<cstring> #include<queue> using namespace std; const int VM=210; const int EM=1010; const int INF=0x3f3f3f3f; struct Edge{ int to,nxt; int cap; }edge[EM<<1]; int n,m,cnt,head[VM]; int src,des,tot,sum,indeg[VM],outdeg[VM]; int dep[VM],gap[VM],cur[VM],aug[VM],pre[VM]; void addedge(int cu,int cv,int cw){ edge[cnt].to=cv; edge[cnt].cap=cw; edge[cnt].nxt=head[cu]; head[cu]=cnt++; edge[cnt].to=cu; edge[cnt].cap=0; edge[cnt].nxt=head[cv]; head[cv]=cnt++; } int SAP(int n){ int max_flow=0,u=src,v; int id,mindep; aug[src]=INF; pre[src]=-1; memset(dep,0,sizeof(dep)); memset(gap,0,sizeof(gap)); gap[0]=n; for(int i=0;i<=n;i++) cur[i]=head[i]; // 初始化当前弧为第一条弧 while(dep[src]<n){ int flag=0; if(u==des){ max_flow+=aug[des]; for(v=pre[des];v!=-1;v=pre[v]){ // 路径回溯更新残留网络 id=cur[v]; edge[id].cap-=aug[des]; edge[id^1].cap+=aug[des]; aug[v]-=aug[des]; // 修改可增广量,以后会用到 if(edge[id].cap==0) // 不回退到源点,仅回退到容量为0的弧的弧尾 u=v; } } for(int i=cur[u];i!=-1;i=edge[i].nxt){ v=edge[i].to; // 从当前弧开始查找允许弧 if(edge[i].cap>0 && dep[u]==dep[v]+1){ // 找到允许弧 flag=1; pre[v]=u; cur[u]=i; aug[v]=min(aug[u],edge[i].cap); u=v; break; } } if(!flag){ if(--gap[dep[u]]==0) /* gap优化,层次树出现断层则结束算法 */ break; mindep=n; cur[u]=head[u]; for(int i=head[u];i!=-1;i=edge[i].nxt){ v=edge[i].to; if(edge[i].cap>0 && dep[v]<mindep){ mindep=dep[v]; cur[u]=i; // 修改标号的同时修改当前弧 } } dep[u]=mindep+1; gap[dep[u]]++; if(u!=src) // 回溯继续寻找允许弧 u=pre[u]; } } return max_flow; } void Init(){ cnt=0; memset(head,-1,sizeof(head)); memset(indeg,0,sizeof(indeg)); memset(outdeg,0,sizeof(outdeg)); } int main(){ //freopen("input.txt","r",stdin); int t; scanf("%d",&t); while(t--){ Init(); scanf("%d%d",&n,&m); int u,v,c; for(int i=0;i<m;i++){ scanf("%d%d%d",&u,&v,&c); indeg[v]++; outdeg[u]++; if(c==0) addedge(u,v,1); } int flag=1; for(int i=1;i<=n;i++) if((indeg[i]-outdeg[i])%2==1){ flag=0; break; } if(!flag) puts("impossible"); else{ sum=0; src=0, des=n+1; for(int i=1;i<=n;i++){ //无向边建图,有向边舍弃 if(indeg[i]<outdeg[i]) addedge(src,i,(outdeg[i]-indeg[i])/2); else if(indeg[i]>outdeg[i]){ addedge(i,des,(indeg[i]-outdeg[i])/2); sum+=(indeg[i]-outdeg[i])/2; } } int ans=SAP(des+1); if(sum==ans) puts("possible"); else puts("impossible"); } } return 0; }
#include<iostream> #include<cstdio> #include<cstring> #include<queue> using namespace std; const int VM=210; const int EM=1010; const int INF=0x3f3f3f3f; struct Edge{ int u,v,nxt; int cap; }edge[EM<<1]; int n,m,cnt,head[VM]; int src,des,tot,sum,dep[VM],indeg[VM],outdeg[VM]; void addedge(int cu,int cv,int cw){ edge[cnt].u=cu; edge[cnt].v=cv; edge[cnt].cap=cw; edge[cnt].nxt=head[cu]; head[cu]=cnt++; edge[cnt].u=cv; edge[cnt].v=cu; edge[cnt].cap=0; edge[cnt].nxt=head[cv]; head[cv]=cnt++; } void Init(){ cnt=0; memset(head,-1,sizeof(head)); memset(indeg,0,sizeof(indeg)); memset(outdeg,0,sizeof(outdeg)); } int BFS(){ queue<int> q; while(!q.empty()) q.pop(); memset(dep,-1,sizeof(dep)); dep[src]=0; q.push(src); while(!q.empty()){ int u=q.front(); q.pop(); for(int i=head[u];i!=-1;i=edge[i].nxt){ int v=edge[i].v; if(edge[i].cap>0 && dep[v]==-1){ //没有标记,且可行流大于0 dep[v]=dep[u]+1; q.push(v); } } } return dep[des]!=-1; //汇点是否成功标号,也就是说是否找到增广路 } int DFS(int u,int minx){ if(u==des) return minx; int tmp; for(int i=head[u];i!=-1;i=edge[i].nxt){ int v=edge[i].v; if(edge[i].cap>0 && dep[v]==dep[u]+1 && (tmp=DFS(v,min(minx,edge[i].cap)))){ edge[i].cap-=tmp; edge[i^1].cap+=tmp; return tmp; } } dep[u]=-1; return 0; } int Dinic(){ int ans=0,tmp; while(BFS()){ while(1){ tmp=DFS(src,INF); if(tmp==0) break; ans+=tmp; } } return ans; } int main(){ //freopen("input.txt","r",stdin); int t; scanf("%d",&t); while(t--){ Init(); scanf("%d%d",&n,&m); int u,v,c; for(int i=0;i<m;i++){ scanf("%d%d%d",&u,&v,&c); indeg[v]++; outdeg[u]++; if(c==0) addedge(u,v,1); } int flag=1; for(int i=1;i<=n;i++) if((indeg[i]-outdeg[i])%2==1){ flag=0; break; } if(!flag) puts("impossible"); else{ sum=0; src=0, des=n+1; for(int i=1;i<=n;i++){ //无向边建图,有向边舍弃 if(indeg[i]<outdeg[i]) addedge(src,i,(outdeg[i]-indeg[i])/2); else if(indeg[i]>outdeg[i]){ addedge(i,des,(indeg[i]-outdeg[i])/2); sum+=(indeg[i]-outdeg[i])/2; } } int ans=Dinic(); if(sum==ans) puts("possible"); else puts("impossible"); } } return 0; }
相关文章
- SAP HANA 详细介绍
- SAP系统在企业内部控制方面有哪些独到之处?
- SAP WM高阶之上架策略B (Bulk Storage)
- SAP WM高阶之同一个货架上相同物料不同批次的库存可以有同一个SU号码?
- SAP 盘盈盘亏移动类型701&702 Vs 711&712
- SAP MM 公司间STO的BILLING输出报错 - Inbound partner profile does not exist –
- SAP MM 移动类型101与103+105组合混用?
- 2020年SAP项目艰辛曲折的开工历程 III
- SAP S4HANA TR传输之操作
- SAP QA32试图做UD,系统报错-工厂 BTYC中的 QM 基选设置需要维护
- 「SAP技术」SAP不够严谨?
- SAP MM 有了采购订单历史的PO行项目里的采购附加费不允许再改了?
- SAP 如何得到交货单上的序列号清单?
- SAP MM MI01事务代码里的批次确定
- 《SAP CRM管理与实施指南》一一1.3 常用名词解释
- 以人致胜 SAP开启HR云端变革之旅
- SAP ABAP编程 字符串加密-MD5_CALCULATE_HASH_FOR_CHAR
- 没有永远的对手只有永远的利益,浪潮与SAP合作了!