【图论】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;
    }
}