求最大匹配
用向量,数组会超内存
双向,除二
hungary算法
#include <iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
const int MAXN=5001;
vector<int> g[MAXN];
int uN,vN; //u,v数目
int linker[MAXN];
bool used[MAXN];
bool dfs(int u)
{
int v;
for(int i=0;i<g[u].size();i++)
{v=g[u][i];
if(!used[v])
{
used[v]=true;
if(linker[v]==-1||dfs(linker[v]))
{
linker[v]=u;
return true;
}
}
}
return false;
}
int hungary()
{
int res=0;
int u;
memset(linker,-1,sizeof(linker));
for(u=1;u<=uN;u++)
{
memset(used,0,sizeof(used));
if(dfs(u)) res++;
}
return res;
}
int main()
{
int t,n,i,x,y;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(i=0;i<=n;i++)
g[i].clear();
for(i=0;i<n*3/2;i++)
{
scanf("%d%d",&x,&y);
g[x].push_back(y);
g[y].push_back(x);
}
uN=n;
printf("%d\n",hungary()/2);
}
return 0;
}
Hopcroft-Carf算法
#include <iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
const int MAXN=5001;
const int INF=1<<28;
vector<int> g[MAXN];
int Mx[MAXN],My[MAXN],Nx,Ny;
int dx[MAXN],dy[MAXN],dis;
bool vst[MAXN];
bool searchP()
{
queue<int> Q;
dis=INF;
memset(dx,-1,sizeof(dx));
memset(dy,-1,sizeof(dy));
for(int i=1;i<=Nx;i++)
if(Mx[i]==-1)
{
Q.push(i);
dx[i]=0;
}
while(!Q.empty())
{
int u=Q.front();
Q.pop();
if(dx[u]>dis)
break;
for(int i=0;i<g[u].size();i++)
if(dy[g[u][i]]==-1)
{
int v=g[u][i];
dy[v]=dx[u]+1;
if(My[v]==-1)
dis=dy[v];
else
{
dx[My[v]]=dy[v]+1;
Q.push(My[v]);
}
}
}
return dis!=INF;
}
bool DFS(int u)
{
for(int i=0;i<g[u].size();i++)
if(!vst[g[u][i]]&&dy[g[u][i]]==dx[u]+1)
{
int v=g[u][i];
vst[v]=1;
if(My[v]!=-1&&dy[v]==dis)
continue;
if(My[v]==-1||DFS(My[v]))
{
My[v]=u;
Mx[u]=v;
return 1;
}
}
return 0;
}
int MaxMatch()
{
int res=0;
memset(Mx,-1,sizeof(Mx));
memset(My,-1,sizeof(My));
while(searchP())
{
memset(vst,0,sizeof(vst));
for(int i=1;i<=Nx;i++)
if(Mx[i]==-1&&DFS(i))
res++;
}
return res;
}
int main()
{
int t,n,i,x,y;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(i=1;i<=n;i++)
g[i].clear();
for(i=0;i<n*3/2;i++)
{
scanf("%d%d",&x,&y);
g[x].push_back(y);
g[y].push_back(x);
}
Nx=n;Ny=n;
printf("%d\n",MaxMatch()/2);
}
return 0;
}
神算法
#include<stdio.h>
int main()
{
int t,n,i,x,y;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(i=0;i<n*3/2;i++)
scanf("%d%d",&x,&y);
printf("%d\n",n/2);
}
return 0;
}