HDU 1845

求最大匹配

用向量，数组会超内存

双向，除二

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;
}