poj1637 Sightseeing tour 混合图欧拉回路判定

　　传送门

第一次做这种题， 尽管ac了但是完全不知道为什么这么做。

题目就是给一些边， 有向边与无向边混合， 问你是否存在欧拉回路。

做法是先对每个点求入度和出度， 如果一条边是无向边， 就随便指定一个方向， 然后连一条边， 权值为1。 最后统计入度出度， 如果一个点的（入度-出度）%2==1， 就说明不存在欧拉回路。 如果全都满足， 就判断每个点的入度出度的大小关系， 入度>出度， 就向汇点连一条边， 权值为（入度-出度）/2， 相反的话就向源点连边。

跑一遍最大流， 看是否满流， 如果满流就说明存在。

完全不理解.....还是太弱。

 #include<bits/stdc++.h>
 using namespace std;
 #define pb(x) push_back(x)
 #define ll long long
 #define mk(x, y) make_pair(x, y)
 #define lson l, m, rt<<1
 #define mem(a) memset(a, 0, sizeof(a))
 #define rson m+1, r, rt<<1|1
 #define mem1(a) memset(a, -1, sizeof(a))
 #define mem2(a) memset(a, 0x3f, sizeof(a))
 #define rep(i, a, n) for(int i = a; i<n; i++)
 #define ull unsigned long long
 typedef pair<int, int> pll;
 const double PI = acos(-1.0);
 const double eps = 1e-;
 const int mod = 1e9+;
 const int inf = ;
 const int dir[][] = { {-, }, {, }, {, -}, {, } };
 const int maxn = ;
 const int maxE = ;
 int indeg[maxn], outdeg[maxn];
 int head[maxE], s, t, num, q[maxE], dis[maxn];
 struct node
 {
     int to, nextt, c;
 }e[maxE];
 void init() {
     mem1(head);
     num = ;
     mem(indeg);
     mem(outdeg);
 }
 void add(int u, int v, int c) {
     e[num].to = v; e[num].nextt = head[u]; e[num].c = c; head[u] = num++;
     e[num].to = u; e[num].nextt = head[v]; e[num].c = ; head[v] = num++;
 }
 int bfs() {
     mem(dis);
     int st = , ed = ;
     q[ed++] = s;
     dis[s] = ;
     while(st<ed) {
         int u = q[st++];
         for(int i = head[u]; ~i; i = e[i].nextt) {
             int v = e[i].to;
             if(e[i].c&&!dis[v]) {
                 dis[v] = dis[u]+;
                 if(v == t)
                     return ;
                 q[ed++] = v;
             }
         }
     }
     return ;
 }
 int dfs(int u, int limit) {
     if(u == t)
         return limit;
     int cost = ;
     for(int i = head[u]; ~i; i = e[i].nextt) {
         int v = e[i].to;
         if(e[i].c&&dis[v] == dis[u]+) {
             int tmp = dfs(v, min(limit-cost, e[i].c));
             if(tmp>) {
                 e[i].c -= tmp;
                 e[i^].c += tmp;
                 cost += tmp;
                 if(limit == cost)
                     break;
             } else {
                 dis[v] = -;
             }
         }
     }
     return cost;
 }
 int dinic() {
     int ans = ;
     while(bfs()) {
         ans += dfs(s, inf);
     }
     return ans;
 }
 int main()
 {
     int T, n, m, x, y, z;
     cin>>T;
     while(T--) {
         scanf("%d%d", &n, &m);
         init();
         s = , t = n+;
         while(m--) {
             scanf("%d%d%d", &x, &y, &z);
             indeg[y]++;
             outdeg[x]++;
             if(z == ) {
                 add(x, y, );
             }
         }
         int flag = ;
         for(int i = ; i<=n; i++) {
             if(abs(indeg[i]-outdeg[i])%==) {
                 flag = ;
                 break;
             }
         }
         if(flag) {
             cout<<"impossible"<<endl;
             continue;
         }
         int sum = ;
         for(int i = ; i<=n; i++) {
             if(indeg[i]<outdeg[i]) {
                 add(s, i, (outdeg[i]-indeg[i])/);
             } else {
                 add(i, t, (indeg[i]-outdeg[i])/);
                 sum += (indeg[i]-outdeg[i])/;
             }
         }
         int ans = dinic();
         if(ans == sum) {
             cout<<"possible"<<endl;
         } else {
             cout<<"impossible"<<endl;
         }
     }
 }