HDU - 4687(带花树求最大匹配的变形)

本题目设定n<40 , m < 123

n个点m个点对匹配，求找出所有的冗余对（它不在任何一个最大匹配中）。

分析：

就按照题意，先求个最大值，然后尝试删除每条边，做一次带花树匹配看是不是最大值-1，是的话，那么改边保留，反之为冗余边。第一开始读错了,题意也是醉了。

还有板子是抄别人的，用自己的一直wa。

#include<iostream>
#include<cstdlib>
#include<stdio.h>
#include<algorithm>
#include<memory.h>
#include<cstring>
#include <deque>
using namespace std;
typedef long long ll;
const int N = 1e5+100;

#define MAXN 300
#define SET(a,b) memset(a,b,sizeof(a))
deque<int> Q;
//g[i][j]存放关系图：i,j是否有边,match[i]存放i所匹配的点
bool g[MAXN][MAXN],inque[MAXN],inblossom[MAXN];
int match[MAXN],pre[MAXN],base[MAXN];

//找公共祖先
int findancestor(int u,int v)
{
    bool inpath[MAXN]={false};
    while(1)
    {
        u=base[u];
        inpath[u]=true;
        if(match[u]==-1)break;
        u=pre[match[u]];
    }
    while(1)
    {
        v=base[v];
        if(inpath[v])return v;
        v=pre[match[v]];
    }
}

//压缩花
void reset(int u,int anc)
{
    while(u!=anc)
    {
        int v=match[u];
        inblossom[base[u]]=1;
        inblossom[base[v]]=1;
        v=pre[v];
        if(base[v]!=anc)pre[v]=match[u];
        u=v;
    }
}

void contract(int u,int v,int n)
{
    int anc=findancestor(u,v);
    SET(inblossom,0);
    reset(u,anc);reset(v,anc);
    if(base[u]!=anc)pre[u]=v;
    if(base[v]!=anc)pre[v]=u;
    for(int i=1;i<=n;i++)
        if(inblossom[base[i]])
        {
            base[i]=anc;
            if(!inque[i])
            {
                Q.push_back(i);
                inque[i]=1;
            }
        }
}

bool dfs(int S,int n)
{
    for(int i=0;i<=n;i++)pre[i]=-1,inque[i]=0,base[i]=i;
    Q.clear();Q.push_back(S);inque[S]=1;
    while(!Q.empty())
    {
        int u=Q.front();Q.pop_front();
        for(int v=1;v<=n;v++)
        {
            if(g[u][v]&&base[v]!=base[u]&&match[u]!=v)
            {
                if(v==S||(match[v]!=-1&&pre[match[v]]!=-1))contract(u,v,n);
                else if(pre[v]==-1)
                {
                    pre[v]=u;
                    if(match[v]!=-1)Q.push_back(match[v]),inque[match[v]]=1;
                    else
                    {
                        u=v;
                        while(u!=-1)
                        {
                            v=pre[u];
                            int w=match[v];
                            match[u]=v;
                            match[v]=u;
                            u=w;
                        }
                        return true;
                    }
                }
            }
        }
    }
    return false;
}

int solve(int n)
{
    SET(match,-1);
    int ans=0;
    for(int i=1;i<=n;i++)
        if(match[i]==-1&&dfs(i,n))
            ans++;
    return ans;
}
typedef pair<int,int> pii;
pii st[MAXN];
int vis[MAXN],n,m;
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        memset(g,false,sizeof(g));
        memset(vis,0,sizeof(vis));
        for(int i=1;i<=m;i++){
              int u,v;
              scanf("%d %d",&u,&v);
              if(u > v) swap(u,v);
              st[i] = pii(u,v);
              g[u][v] = g[v][u] = 1;
        }
        int ans = solve(n);
        int res = 0;
        for(int i=1;i<=m;i++){
            int x = st[i].first;
            int y = st[i].second;
            memset(g,false,sizeof(g));
            for(int j=1;j<=m;j++){
                 if(st[j].first ==  x) continue;
                 if(st[j].second == x) continue;
                 if(st[j].first ==  y) continue;
                 if(st[j].second == y) continue;
                 g[st[j].first][st[j].second]=g[st[j].second][st[j].first]=true;
            }
            int ttt = solve(n);
            if(ttt == ans-1){
                 vis[i] = 1;
            }
            else res++;
        }
        cout<<res<<endl;
        int tt = 0;
        for(int i=1;i<=m;i++){
          if(!vis[i]){
               if(tt++) printf(" ");
               printf("%d",i);
          }
        }
        printf("\n");
    }
}