一般图匹配

最新推荐文章于 2024-08-28 18:40:20 发布
最新推荐文章于 2024-08-28 18:40:20 发布
阅读量177 收藏
点赞数
原文链接：http://www.cnblogs.com/zhujiangning/p/7157465.html
版权

tutte矩阵

一般图的最大匹配数为tutte矩阵的秩(必为偶数)除2

求出最大匹配后加入(n-2*最大匹配数)个点与所有点连边，新图一定存在最大匹配

与新点匹配的点在原图中不与任何点匹配

求一张图的完美匹配时可尝试删除一条边，两个点看是否仍存在完美匹配

有结论两点i,j可能在最大匹配中当且仅当

A的逆矩阵[j][i]不为0且i,j有边相连

删除一条边时可以以第i行的第j列为主元，第j行的第i列为主元消元维护逆矩阵

复杂度O(n^3)常数巨大。

#pragma GCC optimize(2)
#include <cstdio>
#include <iostream>
#include <cmath>
#define LL long long
using namespace std;

  const LL mo=1e9+7;
  LL sed=12939512;

  int n,m;
  int b[2001],fin[2001];
  int sid[124751][2],cnt,edge[2001][2001];
  LL ran(){
      sed*=21409124;sed+=1249134;sed^=1234904;sed%=mo;
      return(sed);
  }

  LL qpow(LL bas,int powe){
      LL ret=1;
      for (;powe;bas*=bas,bas%=mo){
        if (powe&1) ret*=bas,ret%=mo;
      powe>>=1;    
    }
    return(ret);
  }

  struct matrix{
    LL a[1201][1201],tmp[1201][1201];
    int n,m;
    
    int getrank(){
      for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) tmp[i][j]=a[i][j];
      for (int i=1;i<=n;i++){
          for (int j=i;j<=n;j++) if (a[j][i]) {swap(a[j],a[i]);break;}
          if (!a[i][i]) continue;
          LL inv=qpow(a[i][i],mo-2);
        for (int j=i+1;j<=n;j++){
          LL bas=a[j][i]*inv%mo;
          if (!bas) continue;
          for (int k=i;k<=n;k++)
            if (a[i][k])
              a[j][k]-=a[i][k]*bas%mo,a[j][k]%=mo,a[j][k]+=mo,a[j][k]%=mo;
          }    
      }
      LL ret=0;
      for (int i=1;i<=n;i++) if (a[i][i]) ret++;
      return(ret);
    }
    
    void getinv(){
      for (int i=1;i<=n;i++) a[i][n+i]=1;
      for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) a[i][j]=(a[i][j]%mo+mo)%mo;
      for (int i=1;i<=n;i++){
        int po;LL maxi=0;
        for (int j=i;j<=n;j++){
          if (abs(a[j][i])>maxi){
            maxi=abs(a[j][i]);po=j;
          }
        }
        for (int j=1;j<=2*n;j++){
          LL t=a[i][j];a[i][j]=a[po][j];a[po][j]=t;
        }
        if (abs(maxi)==0) continue;
        
        LL rev=qpow(a[i][i],mo-2);    
        for (int j=i+1;j<=n;j++){
          LL tim=a[j][i]*rev%mo;
          for (int k=i;k<=2*n;k++) a[j][k]-=a[i][k]*tim%mo,a[j][k]%=mo;
          }
        }
        
        for (int i=1;i<=n;i++) for (int j=1;j<=2*n;j++) a[i][j]=(a[i][j]%mo+mo)%mo;
        for (int i=n;i>=1;i--){
          for (int j=i+1;j<=n;j++){
            for (int k=n+1;k<=2*n;k++)
              a[i][k]-=a[i][j]*a[j][k]%mo,a[i][k]%=mo;
            a[i][j]=0;          
          }
          for (int j=n+1;j<=2*n;j++)
            a[i][j]*=qpow(a[i][i],mo-2),a[i][j]%=mo;
          a[i][i]=1;  
        }
        
        for (int i=1;i<=n;i++)
          for (int j=1;j<=n;j++)
            a[i][j]=(a[i][j+n]%mo+mo)%mo;
    }
    
    void elim(int x,int y){
      for (int i=1;i<=n;i++)
        if (!b[i]){
          LL bas=a[i][y]*qpow(a[x][y],mo-2)%mo;
          for (int j=1;j<=n;j++){
            a[i][j]-=a[x][j]*bas%mo;
            a[i][j]%=mo;a[i][j]+=mo;a[i][j]%=mo;
          }
        }
    }
  }a;

  int main(){      
      scanf("%d%d",&n,&m);
      a.n=a.m=n;
      for (int i=1;i<=m;i++){
        int t1,t2;
        scanf("%d%d",&t1,&t2);
        sid[++cnt][0]=t1;sid[cnt][1]=t2;
      LL num=ran();
      a.a[min(t1,t2)][max(t1,t2)]=num;
      a.a[max(t1,t2)][min(t1,t2)]=mo-num;
    }
    int ans=a.getrank()/2;
    printf("%d\n",ans);
    
    int addin=n-ans*2;
    a.n=n+addin;a.m=n+addin;
    for (int i=1;i<=a.n;i++) for (int j=1;j<=a.m;j++) a.a[i][j]=0;
    for (int i=1;i<=cnt;i++){
      LL num=ran();
      a.a[min(sid[i][0],sid[i][1])][max(sid[i][0],sid[i][1])]=num;
      a.a[max(sid[i][0],sid[i][1])][min(sid[i][0],sid[i][1])]=-num;
      edge[min(sid[i][0],sid[i][1])][max(sid[i][0],sid[i][1])]=1;
      edge[max(sid[i][0],sid[i][1])][min(sid[i][0],sid[i][1])]=1;
    }
    for (int i=1;i<=addin;i++) 
      for (int j=1;j<=n;j++){
          LL num=ran();
          a.a[j][n+i]=num;
          a.a[n+i][j]=-num;
          edge[j][n+i]=edge[n+i][j]=1;
      }
    a.getinv();
    
    for (int i=1;i<=n;i++) if (!b[i]){
      for (int j=i+1;j<=n+addin;j++)
        if (edge[i][j]&&a.a[j][i]!=0&&!b[j]){
          fin[i]=j;fin[j]=i;
          b[i]=1;b[j]=1;
          a.elim(i,j);
          a.elim(j,i);
          break;
        }
    }
    for (int i=1;i<=n;i++) printf("%d ",(fin[i]<=n ? fin[i]:0));
  }

 -------------------------------

以上方法可以解决必选点等问题。若只需求最大匹配,可发现最大匹配的点为tutte矩阵中最大线性无关组所对应的点。求最大匹配即可。

在求逆和消元中进行一定的常数优化。

#include <cstdio>
#include <iostream>
#include <cmath>
#define LL long long
using namespace std;

  const LL mo=1e9+7;
  LL sed=12939512;

  struct linknode{
      int pre,next;
  }row[2001],col[2001];

  int n,m;
  int b[2001],fin[2001];
  int sid[124751][2],cnt,edge[2001][2001],lis[2001],inlis[2001];
  LL ran(){
      sed*=21409124;sed+=1249134;sed^=1234904;sed%=mo;
      return(sed);
  }

  LL qpow(LL bas,int powe){
      LL ret=1;
      for (;powe;bas*=bas,bas%=mo){
        if (powe&1) ret*=bas,ret%=mo;
      powe>>=1;    
    }
    return(ret);
  }

  struct matrix{
    LL a[1201][1201],tmp[1201][1201];
    int n,m;
    
    int getrank(){
      for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) tmp[i][j]=a[i][j];
      for (int i=1;i<=n;i++){
          for (int j=i;j<=n;j++) if (a[j][i]) {swap(a[j],a[i]);break;}
          if (!a[i][i]) continue;
          LL inv=qpow(a[i][i],mo-2);
        for (int j=i+1;j<=n;j++){
          LL bas=a[j][i]*inv%mo;
          if (!bas) continue;
          for (int k=i;k<=n;k++)
            if (a[i][k])
              a[j][k]-=a[i][k]*bas%mo,a[j][k]%=mo,a[j][k]+=mo,a[j][k]%=mo;
          }    
      }
      LL ret=0;
      for (int i=1;i<=n;i++) if (a[i][i]) ret++;
      return(ret);
    }
    
    void getinv(){
      for (int i=1;i<=n;i++) a[i][n+i]=1;
      for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) a[i][j]=(a[i][j]%mo+mo)%mo;
      for (int i=1;i<=n;i++){
        int po;LL maxi=0;
        for (int j=i;j<=n;j++){
          if (abs(a[j][i])>maxi){
            maxi=abs(a[j][i]);po=j;
          }
        }
        for (int j=1;j<=2*n;j++){
          LL t=a[i][j];a[i][j]=a[po][j];a[po][j]=t;
        }
        if (abs(maxi)==0) continue;
        
        LL rev=qpow(a[i][i],mo-2);    
        for (int j=i+1;j<=n;j++){
          if (!a[j][i]) continue;
          LL tim=a[j][i]*rev%mo;
          for (int k=i;k<=2*n;k++) if (a[i][k]) a[j][k]-=a[i][k]*tim%mo,a[j][k]%=mo;
          }
        }
        
        for (int i=1;i<=n;i++) for (int j=1;j<=2*n;j++) a[i][j]=(a[i][j]%mo+mo)%mo;
        for (int i=n;i>=1;i--){
          for (int j=i+1;j<=n;j++) if (a[i][j]){
            for (int k=n+1;k<=2*n;k++) if (a[j][k])
              a[i][k]-=a[i][j]*a[j][k]%mo,a[i][k]%=mo;
            a[i][j]=0;          
          }
          LL rev=qpow(a[i][i],mo-2);
          for (int j=n+1;j<=2*n;j++)
            if (a[i][j])
              a[i][j]*=rev,a[i][j]%=mo;
          a[i][i]=1;  
        }
        
        for (int i=1;i<=n;i++)
          for (int j=1;j<=n;j++)
            a[i][j]=(a[i][j+n]%mo+mo)%mo;
    }
    
    void elim(int x,int y){
      for (int i=0;i<=n;i=row[i].next)
        if (!b[i]){
          LL bas=a[i][y]*qpow(a[x][y],mo-2)%mo;
          for (int j=col[0].next;j<=n;j=col[j].next){
            a[i][j]-=a[x][j]*bas%mo;
            a[i][j]%=mo;a[i][j]+=mo;a[i][j]%=mo;
          }
        }
    }
  }a;

  void del(int po){
      b[po]=1;
      row[row[po].pre].next=row[po].next;
      row[row[po].next].pre=row[po].pre;
      col[col[po].pre].next=col[po].next;
      col[col[po].next].pre=col[po].pre;
  }

  int main(){      
      scanf("%d%d",&n,&m);
      a.n=a.m=n;
      for (int i=1;i<=m;i++){
        int t1,t2;
        scanf("%d%d",&t1,&t2);
        sid[++cnt][0]=t1;sid[cnt][1]=t2;
      LL num=ran();
      a.a[min(t1,t2)][max(t1,t2)]=num;
      a.a[max(t1,t2)][min(t1,t2)]=mo-num;
    }
    int ans=a.getrank()/2;
    printf("%d\n",ans);
    
    m=0;
    for (int i=1;i<=n;i++) if (a.a[i][i]) lis[++m]=i,inlis[i]=m;
    for (int i=1;i<=m;i++){
      row[i].pre=i-1;row[i].next=i+1;
      col[i].pre=i-1;col[i].next=i+1;
    }
    row[0].next=col[0].next=1;
    for (int i=1;i<=a.n;i++) for (int j=1;j<=a.m;j++) a.a[i][j]=0;
    a.n=m;a.m=m;
    for (int i=1;i<=cnt;i++) if (inlis[sid[i][0]]&&inlis[sid[i][1]]){
      LL num=ran(),u=inlis[sid[i][0]],v=inlis[sid[i][1]];
      if (u>v) swap(u,v);
      a.a[u][v]=num;
      a.a[v][u]=-num;
      edge[u][v]=1;edge[v][u]=1;
    }
    a.getinv();
    
    for (int i=1;i<=m;i++) if (!b[i]){
      for (int j=i+1;j<=m;j++)
        if (edge[i][j]&&a.a[j][i]!=0&&!b[j]){
          fin[lis[i]]=lis[j];fin[lis[j]]=lis[i];
          del(i);del(j);
          a.elim(i,j);
          a.elim(j,i);
          break;
        }
    }
    for (int i=1;i<=n;i++) printf("%d ",fin[i]);
  }

 

转载于:https://www.cnblogs.com/zhujiangning/p/7157465.html

weixin_34240520
关注
一般图最大匹配
_Griefs
08-24 382
一、例题： Universal Online Judge #79. 一般图最大匹配 从前一个和谐的班级，所有人都是搞OI的。有 n 个是男生，有 0 个是女生。男生编号分别为 1,…,n。 现在老师想把他们分成若干个两人小组写动态仙人掌，一个人负责搬砖另一个人负责吐槽。每个人至多属于一个小组。 有若干个这样的条件：第 v 个男生和第 u 个男生愿意组成小组。 请问这个班级里最多产生多少个小组？ 输...
一般图匹配算法，带花树算法
05-23
【一般图匹配算法与带花树算法】 在图论中，匹配问题是一个经典的问题，它涉及到寻找图中的一组边，使得没有两条边共享相同的顶点。一般图匹配是指在一个无向图中寻找最大的匹配，即最多的独立边对。带花树算法...
一般图的匹配题集
夏天的风
09-17 5264
最大匹配(带花树) URAL 1099 ZOJ3316 Game HDU3446/TOJ 3557 daizhenyang's chess HDU3551 hard problem 最大权匹配(求解中国邮路问题) http://poj.org/problem?id=2404
一般图最大权匹配
最新发布
dingjiapeng1234的博客
08-28 937
匹配：对于图 G=(V,E)G=(V,E) 的一个匹配 MM 是边集 EE 的子集，其中两两不共用顶点。最大匹配：边数最多的匹配。完美匹配：包含原图所有点的匹配。最大权匹配：对于图 G(V,E)G(V,E)，边 e∈Ee∈E 的权重为 w(e)w(e)，最大权匹配 M=(V′,E′)M=(V′,E′) 为 ∑ew(e)∑e​w(e) 最大的匹配。最大权最大匹配：最大权最大匹配 M(V′,E′)M(V′,E′) 为 ∑ew(e)∑e​w(e) 最大的最大匹配。
URAL 1099. Work Scheduling 一般图匹配模板题
弱智就要多努力
02-26 514
题目：http://acm.timus.ru/problem.aspx?space=1&num=1099题意：给出一个无向图，求无向图的最大匹配思路：模板题#include <bits/stdc++.h> using namespace std;const int N = 310;struct edge { int to, next; }g[N*N*2]; int n; int cnt, h
一般图的匹配
soberman的专栏
03-09 662
YB的模板库-一般图匹配带花树算法.doc
05-29
YB的模板库-一般图匹配带花树算法.doc
matlab直方图匹配_直方图匹配_
10-02
直方图匹配是一种在图像处理领域中广泛应用的技术，主要用于调整图像的亮度和对比度，使得一张图像的色调分布与另一张图像相匹配。这个过程能够帮助我们统一不同图像的视觉效果，尤其在图像融合、图像增强或者在多源...
一般图中的加权最大匹配：计算一般无向图中的最大加权匹配。-matlab开发
05-30
计算一般无向图中的最大加权匹配。 可以选择只考虑最大基数匹配。 最初由 Joris van Rantwijk 用 Python 编写： http://jorisvr.nl/maximummatching.html 移植到 MATLAB，经许可（但未优化，例如模块化），Daniel ...
C# winform 图片定位 图片相似度匹配 查询图1在图2中的位置
05-06
2. 执行模板匹配：调用`MatchTemplate`，传入大图、小图以及匹配方法（如`CV_TM_CCOEFF_NORMED`）。 3. 使用`minMaxLoc`找出匹配度最高的位置。 4. 转换坐标：由于`minMaxLoc`返回的是匹配图像的坐标，需要根据大图...
Minimum-Cost-Perfect-Matching:用于在一般图中找到完美匹配的算法的C ++实现
05-10
一般图中最大基数匹配和最小成本完美匹配问题的算法 按照以下说明，我在2011年的博士学位期间实现了这些算法： Gerards，AMH（1995）。 匹配。 在M. Ball，T。Magnanti，C。Monma和G. Nemhauser中，网络模型的编辑，《运筹学与管理科学手册》第7卷，第3章，第135-224页。 爱思唯尔。 有关如何使用API​​的示例，请参见Example.cpp。 使用G ++进行编译： g++ -O3 Example.cpp BinaryHeap.cpp Matching.cpp Graph.cpp -o example 用作匹配求解器： ./example -f <filename> <--minweight | --max> --minweight实现最小重量的完美匹配 --max用于最大基数匹配 文件格式： 前两行分别给出n（顶点数）和m（边数）
一般图匹配（带花树算法）
weixin_45735431的博客
07-16 2172
一般图最大匹配 我们知道，二分图匹配有个著名的匈牙利算法（ntr算法 ），用于处理二分图的最大匹配关系。同样的，对于一般图，我们也有个常用的算法来处理一般图的匹配问题，这就是著名的带花树算法 一般图和二分图的区别 二分图是一个可以分为两个顶点集，并且这两个顶点集内部没有相互关系，只在一个集合和另一个集合之间存在关联 但一般图就是存在顶点集内部有关联的情况 二分图 一般图 这时候我们就不能用匈牙利算法，因为匈牙利算法便无法保证最优，这时带花树算法就基于匈牙利算法的基础上进行了处理环的情况 带花树算法 带花
一般图匹配（带花树）
maze_illusion的博客
09-06 466
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]...
[图论]------ 一般图最大匹配（带花树算法）
Berserker____的博客
03-13 1984
吐槽 这两天的训练赛接二连三地出这种又冷门又难的算法。 不得不硬刚。。。。 本博客默认阅读者对匈牙利算法有一定理解。 带花树算法 带花树算法，英文： Blossom algorithmBlossom \ algorithmBlossom algorithm，或Edmonds Matching AlgorithmEdmonds \ Matching \ AlgorithmEdmonds Matching Algorithm。是 Jack Edmond
【学习笔记】一般图最大匹配
七代目鸽影
08-09 791
尚有很多内容需补充，目前只是个空壳。
#79. 一般图最大匹配（带花树算法）
温姑娘的博客
10-15 573
题目链接：http://uoj.ac/problem/79 一般图的匹配算法，带花树算法。 今天上午无意看到了这个算法，出于好奇就像学学，结果感觉自己入了一个坑， 好难的算法，上午看到现在，还没怎么明白过来。水了一道模板题，暂且过一道 模板题吧。 AC代码： #include #include #include #include #include #i
一般图匹配--带花树算法
Free Cornel~~~跳动的犰狳
01-15 396
思路很简单，但是写起来很迷。 #include #include #include #include using namespace std; const int maxn=505,maxm=250005; int link[maxn],pre[maxn],fa[maxn],col[maxn],n; int Begin[maxn],Next[maxm],to[maxm],e; bool vis
一般图带花树匹配（URAL - 1099）模板
龙崎的专栏
08-26 1693
好久之前就想学一下什么是带花树，今天做了套题，终于看了看，不过目前的水平也就是理解它的大体思想，至于实现的细节，真的是无力。 下面推荐几个讲解的博客： 无向图匹配的带花树算法 带花树（一般图最大匹配） 在北京冬令营的时候，yby提到了“带花树开花”算法来解非二分图的最大匹配。 于是，我打算看看这是个什么玩意。其实之前，我已经对这个算法了解了个大概，但是。。。真的不敢去写。 有一个叫Gal
带花树算法实现一般图匹配
"本文主要介绍了一种用于一般图匹配的算法——带花树算法，它的时间复杂度为O(v^4)。虽然网络上存在声称时间复杂度为O(v^3)的算法，但根据描述，作者在实际使用中发现带花树算法更为可靠。带花树算法的核心在于处理...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值