一般图匹配(带花树)

int tol,head[maxn];
struct edge
{
    int u,v,next;
}es[maxm];
void addedge( int u , int v )
{
    es[tol].u = u;
    es[tol].v = v;
    es[tol].next = head[u];
    head[u] = tol++;
}
int n,m,f[maxn],match[maxn],id[maxn],pre[maxn];
int st,ed,q[maxn];
int tim,vis[maxn];
int Find( int x ){ return x==f[x]?f[x]:f[x]=Find(f[x]); }
int Lca ( int u , int v )
{
    tim++;
    while( vis[u]!=tim )
    {
        if ( u!=0 )
        {
            u = Find(u);
            if ( vis[u]==tim ) return u;
            vis[u] = tim;
            if ( match[u]!=0 ) u = Find(pre[match[u]]);
            else u = 0;
        }
        swap( u , v );
    }
    return u;
}
void change( int u , int v , int k )
{
    while ( Find(u)!=k )
    {
        pre[u] = v;
        int t = match[u];
        id[t] = 0;
        q[ed++] = t;
        if ( ed>=n-1 ) ed = 1;
        if ( Find(t)==t ) f[t] = k;
        if ( Find(u)==u ) f[u] = k;
        v = t; u = pre[v];
    }
}
void check( int X )
{
    for ( int i=1 ; i<=n ; i++ )
        f[i] = i,id[i] = -1;
    st = 1;
    ed = 2;
    q[st] = X,id[X] = 0;
    while ( st!=ed )
    {
        int u = q[st];
        for ( int i=head[u] ; i!=-1 ; i=es[i].next )
        {
            int v = es[i].v;
            if ( match[v]==0&&v!=X )
            {
                pre[v] = u;
                int last,t,now = v;
                while ( now!=0 )
                {
                    t = pre[now];
                    last = match[t];
                    match[t] = now;
                    match[now] = t;
                    now = last;
                }
                return;
            }
            if ( id[v]==-1 )
            {
                id[v] = 1;
                pre[v] = u;
                id[match[v]] = 0;
                q[ed++] = match[v];
                if ( ed>=n-1 ) ed = 1;
            }
            else if ( id[v]==0&&Find(u)!=Find(v) )
            {
                int lca = Lca( u , v );
                change( u , v , lca );
                change( v , u , lca );
            }
        }
        st++;
        if ( st>=n-1 ) st = 1;
    }
}
int solve()
{
    tim = 0; memset ( vis , 0 , sizeof(vis) );
    memset ( match , 0 , sizeof(match) );
    for ( int i=1 ; i<=n ; i++ )
        if ( match[i]==0 ) check(i);
    int res = 0;
    for ( int i=1 ; i<=n ; i++ )
        if ( match[i]!=0 ) res++;
    return res>>1;
}

 

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值