URAL 1099 一般图带花树匹配

点击打开链接

题意:给n个人和许多关系,问能够匹配的人数和匹配的人

思路:据说是带花树的模版题,来练一练,原理什么的根本不是本弱可以理解的,权当是一个模版直接用算了,就像最大流一样,而且据本弱做的一些弱题发现这种题目考验的是建图的能力,也就是上边的所有的模版几乎不会动,然后直接建好图就行了,说白了就是太弱看不懂原理直接当模版使了

#include <queue>
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3fll;
const int maxn=230;
int f[maxn],spouse[maxn],Next[maxn],Q[maxn],mark[maxn],vis[maxn];
int n,bot;
vector<int>G[maxn];
int find1(int x){
    if(x!=f[x]) f[x]=find1(f[x]);
    return f[x];
}
void unite(int a,int b){
    int aa=find1(a);
    int bb=find1(b);
    if(aa==bb) return ;
    f[aa]=bb;
}
int findLCA(int x,int y){
    static int t=0;
    t++;
    while(1){
        if(x!=-1){
            x=find1(x);
            if(vis[x]==t) return x;
            vis[x]=t;
            if(spouse[x]!=-1) x=Next[spouse[x]];
            else x=-1;
        }
        swap(x,y);
    }
}
void goup(int a,int p){
    while(a!=p){
        int b=spouse[a],c=Next[b];
        if(find1(c)!=p) Next[c]=b;
        if(mark[b]==2) mark[Q[bot++]=b]=1;
        if(mark[c]==2) mark[Q[bot++]=c]=1;
        unite(a,b);unite(b,c);a=c;
    }
}
void findaugment(int s){
    for(int i=0;i<n;i++) Next[i]=-1,f[i]=i,mark[i]=0,vis[i]=-1;
    Q[0]=s;bot=1;mark[s]=1;
    for(int head=0;spouse[s]==-1&&head<bot;head++){
        int x=Q[head];
        for(unsigned int i=0;i<G[x].size();i++){
            int y=G[x][i];
            if(spouse[x]!=y&&find1(x)!=find1(y)&&mark[y]!=2){
                if(mark[y]==1){
                    int p=findLCA(x,y);
                    if(find1(x)!=p) Next[x]=y;
                    if(find1(y)!=p) Next[y]=x;
                    goup(x,p);goup(y,p);
                }else if(spouse[y]==-1){
                    Next[y]=x;
                    for(int j=y;j!=-1;){
                        int k=Next[j];
                        int l=spouse[k];
                        spouse[j]=k;spouse[k]=j;j=l;
                    }
                    break;
                }else{
                    Next[y]=x;
                    mark[Q[bot++]=spouse[y]]=1;
                    mark[y]=2;
                }
            }
        }
    }
}
int num[maxn][maxn];
int main(){
    scanf("%d",&n);
    int x,y;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++) num[i][j]=0;
    }
    while(scanf("%d%d",&x,&y)!=EOF){
        x--;y--;
        if(x!=y&&num[x][y]==0){
            num[x][y]=num[y][x]=1;
            G[x].push_back(y);
            G[y].push_back(x);
        }
    }
    for(int i=0;i<n;i++) spouse[i]=-1;
    for(int i=0;i<n;i++){
        if(spouse[i]==-1) findaugment(i);
    }
    int ans=0;
    for(int i=0;i<n;i++) if(spouse[i]!=-1) ans++;
    printf("%d\n",ans);
    for(int i=0;i<n;i++){
        if(spouse[i]!=-1&&spouse[i]>i) printf("%d %d\n",i+1,spouse[i]+1);
    }
    return 0;
}
时间: 2024-10-07 03:15:30

URAL 1099 一般图带花树匹配的相关文章

URAL 1099 Work scheduling 一般图的最大匹配 带花树算法(模板)

R - Work scheduling Time Limit:500MS     Memory Limit:65536KB     64bit IO Format:%I64d & %I64u Submit Status Practice URAL 1099 Description There is certain amount of night guards that are available to protect the local junkyard from possible junk r

hdu 2255 二分图带权匹配 模板题

模板+注解在 http://blog.csdn.net/u011026968/article/details/38276945 hdu 2255 代码: //KM×î´ó×îСƥÅä #include <cstdio> #include <cstring> #include <algorithm> #include <iostream> using namespace std; #define INF 0x0fffffff const int MAXN

二分图带权匹配(推箱子问题的思考)

转载请附上原文链接: http://blog.csdn.net/u013351484/article/details/51598270 二分图带权匹配(也叫二分图最优(最佳)匹配,Kuhn-Munkres 算法) 预备知识:二分图最大匹配,二分图完备(完美)匹配 二分图带权匹配,可以把它看作集合 X 的每个顶点到集合 Y 的每个顶点均有边的二分图(设权重均为正数,则原来分别在集合 X 与集合 Y 中没边的顶点 xi,yi 就可以添加一条权重为 0 的边,即 w[xi,yi] = 0):这样的二分

URAL 1099. Work Scheduling 一般图匹配带花树

一般图匹配带花树模版题: 将奇环缩成圈(Blossom),然后找增广路..... 1099. Work Scheduling Time limit: 0.5 second Memory limit: 64 MB There is certain amount of night guards that are available to protect the local junkyard from possible junk robberies. These guards need to sche

Work Scheduling URAL - 1099 一般图最大匹配带花树

#include <stdio.h> #include <string.h> #include <iostream> #include <algorithm> #include <string> using namespace std; const int MAXN = 250; int N; //点的个数,点的编号从1到N bool Graph[MAXN][MAXN]; int Match[MAXN]; bool InQueue[MAXN],I

URAL 1099 Work Scheduling

一般图的最大匹配  带花树开花算法 有两个模板,一个kuangbin大神的,另一个不知道谁写的. #include<stdio.h> #include<string.h> #include<string.h> #include<iostream> #include<queue> #include<algorithm> using namespace std; const int MAXN = 250; int N; //点的个数,点的

poj 2195 二分图带权匹配

题意:有一个矩阵,某些格有人,某些格有房子,每个人可以上下左右移动,问给每个人进一个房子,所有人需要走的距离之和最小是多少. 貌似以前见过很多这样类似的题,都不会,现在知道是用KM算法做了 KM算法目前还没弄懂,先套模板做 Sample Input 2 2 .m H. 5 5 HH..m ..... ..... ..... mm..H 7 8 ...H.... ...H.... ...H.... mmmHmmmm ...H.... ...H.... ...H.... 0 0 Sample Out

二分图带权匹配

带权匹配是指在最大匹配的基础上,使匹配边的边权和最大.一般有两种写法,一个是KM算法(只针对可以完备匹配的二分图),一个是费用流. KM算法在稠密图上比费用流更优秀一些,不过应用范围太小,所以还是鼓励大家用费用流.当然啦,作为一种算法KM也是我们需要了解的(况且我不会网络流??). KM算法有个流程: 1.存图 2.用贪心算法初始化标杆 3.运用匈牙利算法找到完备匹配 4.如果找不到完备匹配,通过修改标杆增加一些边 5.重复做3和4两个步骤,直到找到完备匹配. 其他的话,我表述的也不是很好,理解

二分图带权匹配-Kuhn-Munkres算法模板 [二分图带权匹配]

尴尬...理解不太好T T 1 #include<cstdio> 2 #include<cstring> 3 #include<iostream> 4 #include<algorithm> 5 using namespace std; 6 #define inf 0x3f3f3f3f 7 8 const int maxn=1005; 9 10 int n; 11 //标杆序号 12 int lx[maxn],ly[maxn]; 13 //是否被搜索过 14