hdu 5727(匈牙利算法)

这里

因为是环状,固定一个阴石,其他阴石枚举全排列,然后枚举每个阳石放的位置,建二分图(如果阳石放的位置不会引起爆炸,就连边),跑匈牙利求出最大匹配。n - 最大匹配就是 最少的引起爆炸的阳石数量,再取min

#include <bits/stdc++.h>

#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;

#define LL long long
#define pii pair<int,int>
#define MP make_pair
#define ls i << 1
#define rs ls | 1
#define md (ll + rr >> 1)
#define lson ll, md, ls
#define rson md + 1, rr, rs
#define Pi acos(-1.0)
#define mod 1000000007
#define eps 1e-10
#define inf 0x3f3f3f3f
#define N 10
#define M 800020

bool g[N][N], vis[N];
int fst[N], vv[M], nxt[M], e;
int n, m, p[N], a[N], match[N];
void init(){
    memset(fst, -1, sizeof fst); e = 0;
}
void add(int u, int v){
    vv[e] = v, nxt[e] = fst[u], fst[u] = e++;
}
bool hungry(int u){
    for(int i = fst[u]; ~i; i = nxt[i]){
        int v = vv[i];
        if(!vis[v]){
            vis[v] = 1;
            if(match[v] == 0 || hungry(match[v])){
                match[v] = u;
                return 1;
            }
        }
    }
    return 0;
}
int calc(){
    init();
    p[1] = 1;
    for(int i = 2; i <= n; ++i) p[i] = a[i-1] + 1;
    for(int i = 1; i <= n; ++i){
        for(int j = 1; j <= n; ++j){
            int l = p[j], r = p[j+1];
            if(j == n) r = p[1];
            if(!g[i][l] && !g[i][r]){
                add(i, j);
            }
        }
    }
    int ret = 0;
    memset(match, 0, sizeof match);
    for(int i = 1; i <= n; ++i){
        if(!match[i]){
            memset(vis, 0, sizeof vis);
            if(hungry(i))
                ret++;
        }
    }
    return n - ret;
}
int main(){
    while(scanf("%d%d", &n, &m) != EOF){
        memset(g, 0, sizeof g);
        for(int i = 1; i <= m; ++i){
            int u, v;
            scanf("%d%d", &u, &v);
            g[u][v] = 1;
        }
        if(n == 0 || m == 0){
            puts("0"); continue;
        }
        if(n == 1){
            printf("%d\n", m); continue;
        }
        for(int i = 1; i < n; ++i) a[i] = i;
        int ans = n;
        do{
            ans = min(ans, calc());
        }while(next_permutation(a + 1, a + n));
        printf("%d\n", ans);
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值