二分匹配加速算法（简单，就是先把能匹配的先找出来，不用Hopcroft_Karp算法）

#include
using namespace std;
int n,m;
int INF=0x3f3f3f;
int dis;
int cx[3003];
int dy[3002];
int cy[3003];
int dx[3003];
int used[3002];
int a[3003];
int b[3002];
int c[3002];
int d[3002];
int e[3002];
int mp[3002][3002];
bool searchpath()
{
    queueQ;
  dis=INF;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=0; idis)
            break;
        for(int i=0; i>t;
    int k;
    for(k=1; k<=t; k++)
    {
        int rain;
        scanf("%d",&rain);
        memset(mp,0,sizeof(mp));
        scanf("%d",&m);
        int i,j;
        for(i=0; i

----HDU 2389测试和hp算法时间复杂都差不多

主要函数部分

這裡介紹一個加速的手段：一開始就把一些明顯可以配對的點給配對起來，這樣就不必替每一個未匹配點建立交錯樹了。這個手段在圖很龐大的時候可以發揮作用。

 

1. int greedy_matching()
2. {
3.     int c = 0;
4.     for (int x=0; x
5.         if (mx[x] == -1)
6.             for (int y=0; y
7.                 if (my[y] == -1)
8.                     if (adj[x][y])
9.                     {
10.                         mx[x] = y; my[y] = x;
11.                         c++;
12.                         break;
13.                     }
14.     return c;
15. }
16.  
17. int bipartite_matching()
18. {
19.     memset(mx, -1, sizeof(mx));
20.     memset(my, -1, sizeof(my));
21.  
22.     int c = greedy_matching();  // 能連的先連一連
23.  
24.     for (int x=0; x
25.         if (mx[x] == -1)    // 這行記得補上來
26.         {
27.             memset(vy, false, sizeof(vy));
28.             if (DFS(i)) c++;
29.         }
30.     return c;
31. }

它的時間複雜度僅為一次 Graph Traversal 的時間，不太會影響原演算法的運行效率。