SPOJ 4206 Fast Maximum Matching(HK算法)

n到5w了，显然要用邻接链表存图了，把邻接矩阵的改了改，就AC了。

#include<iostream>
#include<cstring>
#include<queue>
#include<cmath>
#include<cstdio>
using namespace std;
const int maxn=50000+5;
const int inf=(1<<30);
int head[maxn];
int cx[maxn],cy[maxn];
int nx,ny;
int dx[maxn],dy[maxn];
int dis;
bool used[maxn];
int cnt;
struct edge{int u,v,next;}edges[maxn*3];
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(cx[i]==-1) {q.push(i);dx[i]=0;}//对于未遍历的点入队
    }
    //准备分层
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        if(dx[u]>dis) break;//如果目前的层次大于最小增广长度，那么退出
        for(int j=head[u];j!=-1;j=edges[j].next)//对于一切可能的点遍历
        {
            int v=edges[j].v;
            if(dy[v]==-1){//只对未分层的点遍历
                dy[v]=dx[u]+1;
                if(cy[v]==-1) dis=dy[v];
                else{
                    dx[cy[v]]=dy[v]+1;
                    q.push(cy[v]);
                }
            }
        }
    }
    return dis!=inf;
}
//寻找路径配对
bool findpath(int x)
{
    for(int j=head[x];j!=-1;j=edges[j].next)
    {
        int v=edges[j].v;
        if(!used[v]&&dy[v]==dx[x]+1)//符合继续搜索的条件有三个：未访问过，图上联通和层次符合
        {
            used[v]=1;
            if(cy[v]!=-1&&dis==dy[v]) continue;//如果下一个点还是匹配点且目前已经到达增广最小层次，不需要扩展了
            if(cy[v]==-1||findpath(cy[v]))
            {
                cy[v]=x;cx[x]=v;
                return true;
            }
        }
    }
    return false;
}
int HK()
{
    int ans=0;
    memset(cx,-1,sizeof(cx));
    memset(cy,-1,sizeof(cy));
    while(searchP())
    {
        memset(used,0,sizeof(used));
        for(int i=1;i<=nx;i++){
            if(cx[i]==-1)
            {
                if(findpath(i)) ans++;
            }
        }
    }
    return ans;
}
void addedge(int u,int v)
{
    edges[cnt].v=v;
    edges[cnt].next=head[u];
    head[u]=cnt++;
}
int main()
{
    int p;
    while(~scanf("%d%d%d",&nx,&ny,&p))
    {
        memset(head,-1,sizeof(head));
        cnt=0;
        for(int i=0;i<p;i++)
        {
            scanf("%d%d",&edges[i].u,&edges[i].v);
            addedge(edges[i].u,edges[i].v);
        }
        printf("%d\n",HK());
    }
    return 0;
}