https://blog.sengxian.com/solutions/uva-10735 感谢思路和测试数据
存在重边!!!
如上面博客的第二个测试数据
#include <bits/stdc++.h>
//https://blog.sengxian.com/solutions/uva-10735 感谢思路和测试数据
using namespace std;
const int maxn=1e5+5;
const int inf=0x3f3f3f3f;
struct edge
{
int to,cap,rev;
int flag;//添加标记1是正向边 0是反向边
edge(int a=0,int b=0,int c=0,int d=0)
{
to=a;
cap=b;
rev=c;
flag=d;
}
};
vector<edge>G[maxn];
int level[maxn];///顶点到源点的距离标号
int iter [maxn];///当前弧,在其之前的便已经没用了
void add_edge(int from,int to,int cap)
{
G[from].push_back(edge(to,cap,G[to].size(),1));
///因为G[from]添加了一个新的 所以在下面要 减1
G[to].push_back(edge(from,0,G[from].size()-1,0));
}
///初始化 不要忘记在主函数里面添加!!!
void init(int z)
{
for(int i=0; i<z; i++)
G[i].clear();
}
void bfs(int s)
{
memset(level,-1,sizeof level);
queue<int>q;
level[s]=0;
q.push(s);
while(!q.empty())
{
int v=q.front();
q.pop();
for(int i=0; i<G[v].size(); i++)
{
edge &e=G[v][i];
if(e.cap>0&&level[e.to]<0)
{
///确定了远近先后关系
level[e.to]=level[v]+1;
q.push(e.to);
}
}
}
}
int dfs(int v,int t,int f)
{
///到达汇点 返回值
if(v==t||f==0)
return f;
///用到了 引用 改变iter
for(int &i=iter[v]; i<G[v].size(); i++)
{
edge &e=G[v][i];
///正向 且还有流量
if(e.cap>0&&level[v]<level[e.to])
{
int d=dfs(e.to,t,min(f,e.cap));
if(d>0)
{
e.cap-=d;
G[e.to][e.rev].cap+=d;
return d;
}
}
}
return 0;
}
int max_flow(int s,int t)
{
int flow=0;
while(1)
{
bfs(s);
///没有路径能够到达源点了
if(level[t]<0)
return flow;
memset(iter,0,sizeof iter);
int f;
///一遍又一遍的压榨 注意不是等于关系 是等于号!!
while(f=dfs(s,t,inf))
{
flow+=f;
}
}
}
int n,m;
struct oula
{
int from;
int to;
bool vis;
};
vector<oula>mp[105];
stack<int>path;
int vis[105][105];
int in[105],out[105];
void init()
{
for(int i=0; i<105; i++)
mp[i].clear();
while(!path.empty())
path.pop();
memset(vis,0,sizeof vis);
memset(in,0,sizeof in);
memset(out,0,sizeof out);
}
void get_path(int st)
{
for(int i=0; i<mp[st].size(); i++)
{
oula &e=mp[st][i];
if(!e.vis)///通过判断vis[st][e.to]==0 是行不通的!!!
{
e.vis=1;
get_path(e.to);
path.push(e.to);
}
}
}
void print_ans()
{
cout<<1<<" ";
while(!path.empty())
{
cout<<path.top()<<" ";
path.pop();
}
cout<<endl<<endl;
}
int main()
{
int t;
cin>>t;
while(t--)
{
init(105);
init();
cin>>n>>m;
for(int i=0; i<m; i++)
{
int x;
int y;
string c;
cin>>x>>y>>c;
if(c=="U")
add_edge(x,y,1);
else
mp[x].push_back({x,y,0});
in[y]++;
out[x]++;
}
int sum=0;
int ok=true;
for(int i=1; i<=n; i++)
{
int diff=out[i]-in[i];
if(diff%2!=0)
ok=false;
if(diff>0)
{
//cout<<i<<endl;
add_edge(0,i,diff/2);
sum+=diff/2;
}
if(diff<0)
add_edge(i,n+1,-1*diff/2);
}
if(!ok||sum!=max_flow(0,n+1))
{
cout<<"No euler circuit exist"<<endl<<endl;
continue;
}
for(int i=1; i<=n; i++)
{
for(int j=0; j<G[i].size(); j++)
{
edge e=G[i][j];
if(e.to<=n&&e.flag)//不能涉及到汇点 且只统计正向边
if(e.cap==0)
{
//cout<<"st: "<<e.to<<"en: "<<i<<endl;
mp[e.to].push_back({e.to,i,0});
}
else
{
//cout<<"st: "<<i<<"en: "<<e.to<<endl;
mp[i].push_back({i,e.to,0});
}
}
}
get_path(1);
print_ans();
}
return 0;
}