• [Ural1099]工作安排 带花树


    题目大意

    一般图最大匹配.

    题解:

    求解一般图最大匹配.
    直接使用带花树即可.
    (带花树也是刚学)...

    马上写带花树的ppt,会很通俗易懂.
    (充分证明了本苣智商不够,写不出高深的课件)

    如果有想要的依然可以直接私信我 QAQ

    带花树模板

    #include <cstdio>
    #include <cstring>
    #include <algorithm>
    typedef long long ll;
    inline void read(int &x){
        x=0;char ch;bool flag = false;
        while(ch=getchar(),ch<'!');if(ch == '-') ch=getchar(),flag = true;
        while(x=10*x+ch-'0',ch=getchar(),ch>'!');if(flag) x=-x;
    }
    #define rg register int
    #define rep(i,a,b) for(rg i=(a);i<=(b);++i)
    #define per(i,a,b) for(rg i=(a);i>=(b);--i)
    const int maxn = 256;
    int n;
    bool map[maxn][maxn],inq[maxn],inBlo[maxn],mark[maxn];
    int link[maxn],fa[maxn],belong[maxn],q[maxn],l,r;
    inline void bloom(int x,int y){
        rep(i,0,n) mark[i] = inBlo[i] = false;
    #define nx(x) fa[link[(x)]]
        int lca;
        for(rg i=x;i;i=nx(i)) mark[i = belong[i]] = true;
        for(rg i=y;i;i=nx(i)){
            i = belong[i];
            if(mark[i]){lca = i;break;}
        }
        for(rg i=x;belong[i] != lca;i=nx(i)){
            if(belong[nx(i)] != lca) fa[nx(i)] = link[i];
            inBlo[belong[i]] = inBlo[belong[link[i]]] = true;
        }
        for(rg i=y;belong[i] != lca;i=nx(i)){
            if(belong[nx(i)] != lca) fa[nx(i)] = link[i];
            inBlo[belong[i]] = inBlo[belong[link[i]]] = true;
        }
    #undef nx
        if(belong[x] != lca) fa[x] = y;
        if(belong[y] != lca) fa[y] = x;
        rep(i,1,n) if(inBlo[belong[i]]){
            belong[i] = lca;
            if(!inq[i]) inq[q[++r] = i] = true;
        }
    }
    inline void modify(int s){
        int x,y,z = s;
        while(z){
            y = fa[z];x = link[y];
            link[y] = z;link[z] = y;
            z = x;
        }
    }
    inline void bfs(int s){
        rep(i,0,n){
            fa[i] = inq[i] = 0;
            belong[i] = i;
        }
        l = 0;r = -1;q[++r] = s;
        inq[s] = true;
        while(l <= r){
            int u = q[l++];
            rep(v,1,n){
                if(map[u][v] == 0 || belong[u] == belong[v] || link[u] == v) continue;
                if(s == v || link[v] && fa[link[v]]) bloom(u,v);
                else if(fa[v] == 0){
                    fa[v] = u;
                    if(link[v] != 0){
                        q[++r] = link[v];
                        inq[link[v]] = true;
                    }else{
                        modify(v);
                        return ;
                    }
                }
            }
        }return ;
    }
    int main(){
        read(n);int u,v;
        while(scanf("%d%d",&u,&v) != EOF){
            map[u][v] = map[v][u] = true;
        }
        rep(i,1,n) if(link[i] == 0) bfs(i);
        int ans = 0;
        rep(i,1,n) if(link[i] != 0) ++ ans;
        printf("%d
    ",ans);
        rep(i,1,n){
            if(link[i] > i) printf("%d %d
    ",i,link[i]);
        }
        return 0;
    }
    
  • 相关阅读:
    返回顶部按钮效果实现
    WebAPI Angularjs 上传文件
    C# 单元测试
    C# 如何获取Url的host以及是否是http
    Dapper批量操作实体
    易优CMS:type的基础用法
    易优CMS:arcview基础用法
    易优CMS:channel的基础用法
    易优CMS:arclist 文档列表
    c语言必背代码
  • 原文地址:https://www.cnblogs.com/Skyminer/p/6815608.html
Copyright © 2020-2023  润新知