嘟嘟嘟


雖然我已經會網絡流了，但是還是學了一個匈牙利算法。
——就跟我會線段樹，但還是學了樹狀數組一樣。


其實匈牙利算法挺暴力的。簡單來說就是先貪心匹配，然後如果左部點\(i\)匹配不上了，就嘗試更改前面已經匹配好的點，騰出地給他匹配。
因此對於每一個點跑一遍匈牙利算法，如果這個點匹配成功，總匹配數就加1。
感覺沒啥好講的。
關於這道題怎麽做，看我這篇博客吧。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<vector>
#include<stack>
#include<queue>
using namespace std;
#define enter puts("") 
#define space putchar(‘ ‘)
#define Mem(a, x) memset(a, x, sizeof(a))
#define rg register
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-8;
const int maxn = 505;
const int maxe = 1e4 + 5;
inline ll read()
{
  ll ans = 0;
  char ch = getchar(), last = ‘ ‘;
  while(!isdigit(ch)) last = ch, ch = getchar();
  while(isdigit(ch)) ans = (ans << 1) + (ans << 3) + ch - ‘0‘, ch = getchar();
  if(last == ‘-‘) ans = -ans;
  return ans;
}
inline void write(ll x)
{
  if(x < 0) x = -x, putchar(‘-‘);
  if(x >= 10) write(x / 10);
  putchar(x % 10 + ‘0‘);
}

int n, k;
struct Edge
{
  int nxt, to;
}e[maxe];
int head[maxn], ecnt = -1;
void addEdge(int x, int y)
{
  e[++ecnt] = (Edge){head[x], y};
  head[x] = ecnt;
}

int fa[maxn], vis[maxn], vcnt = 0;
//fa:下標是右部點,表示右部點i和左部點fa[i]匹配上了
//vis:表示匹配點i的時候這個點是不是i要匹配的
//因此每一次dfs前應該清空,為了降低復雜度,改為累加標記
bool dfs(int now)
{
  for(int i = head[now], v; i != -1; i = e[i].nxt)
    {
      if(vis[v = e[i].to] != vcnt)
    {
      vis[v] = vcnt;
      if(!fa[v] || dfs(fa[v])) {fa[v] = now; return 1;} 
    }
    }
  return 0;
}


int main()
{
  Mem(head, -1);
  n = read(); k = read();
  for(int i = 1; i <= k; ++i)
    {
      int x = read(), y = read();
      addEdge(x, y);
    }
  int ans = 0;
  for(int i = 1; i <= n; ++i)
    {
      ++vcnt;
      if(dfs(i)) ans++;
    }
  write(ans), enter;
  return 0;
}
 

POJ3041 Asteroids（匈牙利算法）