【图论】hdu6073 Matching In Multiplication

比赛的时候，一眼看到二分图，哇！掏出了dinic模板。之后发现，好像完全不是那么回事。隐隐约约有一点感觉，但就是没有清楚地理清思路。先找出一度节点就没有想到，非联通块更是不知道怎么处理。

本期官方题解说的特别清楚，强烈好评~~
先把一度节点挑出来，得到每个二分图都确定一样的边。
然后剩下一堆度数全部是二的联通快，每个联通快有两种匹配方式，间隔取边，分别得权值积num1,num2;

最后 ans = ans （num1+num2）(num1,num2)… …所有联通快的num1,num2;

有几个注意：
1. 去一度节点，可以不用拓扑那么高端（瑟瑟发抖，差点因为这个放弃了），用了一个队列。
2. 最后的ans部分，我本来又写了个DFS，看了别人代码才发现原来这么简单。
3. 我的程序中点和边都标了号，应该还能省不少空间，丑代码就留给自己看吧QWQ;

#include<iostream>#include<cstdio>#include<cstring>#include<vector>#include<queue>using namespace std;typedef long long ll;const int maxn = 300010;const int mod = 998244353;struct edge{    int to; ll cost;    edge(int t = 0, ll c = 0){        to = t; cost = c;    }};struct Ans{    ll a1, a2;};int d[maxn * 2] = { 0 };int vis[maxn * 2] = { 0 };int CNT = 0;Ans ans2[maxn];vector<int> G[maxn * 2];edge es[maxn * 4];int vise[maxn * 4] = { 0 };queue<int>q;int t, n;ll bfs(){    for (int i = n + 1; i <= 2 * n; i++){        if (d[i] == 1){            q.push(i);        }    }    ll ans = 1;    while (!q.empty()){        int now = q.front(); q.pop();        vis[now] = 1;        edge to;        for (int i = 0; i < G[now].size(); i++){            if (!vis[es[G[now][i]].to])                to = es[G[now][i]];        }        vis[to.to] = 1;        ans = (ans* to.cost) % mod;        for (int i = 0; i < G[to.to].size(); i++){            d[es[G[to.to][i]].to]--; d[to.to]--;            if (d[es[G[to.to][i]].to] == 1) {                //printf("!%d %d %d", now, to.to, es[G[to.to][i]].to);                q.push(es[G[to.to][i]].to);            }        }    }    return ans;}void dfs(int u, int s, ll sum1, ll sum2, int swi){    if (vis[u] && u == s) {        ans2[CNT].a1 = sum1; ans2[CNT++].a2 = sum2;    }    if (vis[u]) return; vis[u] = 1;    for (int i = 0; i < G[u].size(); i++){        edge &e = es[G[u][i]];        int v = e.to;        if (!vise[G[u][i]]){            vise[G[u][i]] = 1;            if (G[u][i] & 1) vise[G[u][i] - 1] = 1;            else vise[G[u][i] + 1] = 1;            if (swi) dfs(v, s, (sum1*e.cost) % mod, sum2, swi ^ 1);            else dfs(v, s, sum1, (sum2*e.cost) % mod, swi ^ 1);        }    }}int main(){    scanf("%d", &t);    while (t--){        scanf("%d", &n);        int cnt = 0;        for (int i = 1; i <= n; i++){            int v1, v2; ll c1, c2;            scanf("%d%lld%d%lld", &v1, &c1, &v2, &c2);            es[cnt].to = n + v1;            es[cnt].cost = c1;            G[i].push_back(cnt++);            es[cnt].to = i;            es[cnt].cost = c1;            G[n + v1].push_back(cnt++);            es[cnt].to = n + v2;            es[cnt].cost = c2;            G[i].push_back(cnt++);            es[cnt].to = i;            es[cnt].cost = c2;            G[v2 + n].push_back(cnt++);            d[v1 + n]++; d[v2 + n]++; d[i] = 2;        }        //去除一度点        ll ans = bfs();        for (int i = 1; i <= n; i++){            if (!vis[i]) dfs(i, i, 1, 1, 1);        }        for (int i = 0; i < CNT; i++){            ans = (ans* (ans2[i].a1 + ans2[i].a2) % mod) % mod;        }        printf("%lld\n", ans);        for (int i = 0; i <= 2 * n + 2; i++) {            vis[i] = 0;`````````            d[i] = 0;            G[i].clear();        }        for (int i = 0; i < cnt; i++){            vise[i] = 0;        }        CNT = 0; cnt = 0;    }}