HDU3861 - The King’s Problem

数组开小了一直判WA太坑了。 SCC + 最小路径覆盖

#include <cstring>
#include <iostream>
#include <cstdlib>
#include <cctype>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#include <map>
#include <set>
#include <stack>
#include <queue>
#define INF (int)(1e9)
#define maxn 5050
using namespace std;
typedef long long ll;
int p = 0;
struct node {
    int to, w, next;
};
node edge[100010];
int head[maxn];
void add_edge(int u, int v, int w) {
    edge[p].to = v;
    edge[p].w = w;
    edge[p].next = head[u];
    head[u] = p;
    ++ p;
}
int pre[maxn], lowlink[maxn], sccno[maxn], dfs_clock, scc_cnt;
//sccno[i]为i所在scc编号
stack<int> S;
ll res, num, cnt;
void dfs1(int u) {
    pre[u] = lowlink[u] = ++dfs_clock;
    S.push(u);
    for (int k = head[u]; k != -1; k = edge[k].next) {
        int v = edge[k].to;
        if (!pre[v]) {
            dfs1(v);
            lowlink[u] = min(lowlink[u], lowlink[v]);
        } else if (!sccno[v]) {
            lowlink[u] = min(lowlink[u], pre[v]);
        }
    }
    if (lowlink[u] == pre[u]) {
        scc_cnt++;
        for (;;) {
            int x = S.top(); S.pop();
            sccno[x] = scc_cnt;
            if (x == u) break;
        }
    }
}
void find_scc(int n) {
    dfs_clock = scc_cnt = 0;
    memset(sccno, 0, sizeof(sccno));
    memset(pre, 0, sizeof(pre));
    for (int i = 0; i < n; ++ i) {
        if (!pre[i]) dfs1(i);
    }
}
int uN, vN;//u,v数目
int linker[maxn];
bool used[maxn];
bool dfs(int u) { //从左边开始找增广路径
    for (int k = head[u]; k != -1; k = edge[k].next) {
        int v = edge[k].to;
        if(!used[v]) {
            used[v] = true;
            if(linker[v] == -1 || dfs(linker[v])) {
                //找增广路,反向
                linker[v] = u;
                return true;
            }
        }
    }
    return false;
}
int max_match() {
    int res=0;
    int u;
    memset(linker,-1,sizeof(linker));
    for(u = 0; u < uN; u++) {
        memset(used, 0, sizeof(used));
        if(dfs(u)) res ++;
    }
    return res;
}
int no[100010][2];
int main() {
    int t;
    scanf("%d", &t);
    while (t --) {
        memset(head, -1, sizeof(head)), p = 0;
        int n, m, u, v;
        scanf("%d%d", &n, &m);
        for (int i = 0; i < m; ++ i) {
            scanf("%d%d", &u, &v);
            -- u, -- v;
            no[i][0] = u, no[i][1] = v;
            add_edge(u,v,0);
        }
        find_scc(n);
        memset(head, -1, sizeof(head)), p = 0;
        for (int i = 0; i < m; ++ i) {
            u = no[i][0], v = no[i][1];
            if (sccno[u] != sccno[v]) {
                add_edge(sccno[u]-1,sccno[v]-1,0);
            }
        }
        uN = vN = scc_cnt;
        cout << scc_cnt - max_match() << endl;
    }
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值