最大流Dinic

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f3f3f3f3f
#define ll long long
#define MAXN 10005    //这里为点的最大个数
using namespace std;
int n,m;             //点数、边数
int sp,tp;           //原点、汇点
struct node
{
    int v,next;
    ll cap;
} mp[MAXN*10];             //MAXN*10表示边的最大可能数（这里开的比较大）
int pre[MAXN],dis[MAXN],cur[MAXN];//cur为当前弧优化，dis存储分层图中每个点的层数（即到原点的最短距离），pre建邻接表
int cnt=0;
void init()//不要忘记初始化
{
    cnt=0;
    memset(pre,-1,sizeof(pre));
}
void add(int u,int v,int w)//加边
{
    mp[cnt].v=v;
    mp[cnt].cap=w;
    mp[cnt].next=pre[u];
    pre[u]=cnt++;
}
bool bfs()//建分层图
{
    memset(dis,-1,sizeof(dis));
    queue<int>q;
    while(!q.empty())
        q.pop();
    q.push(sp);
    dis[sp]=0;
    int u,v;
    while(!q.empty())
    {
        u=q.front();
        q.pop();
        for(int i=pre[u]; i!=-1; i=mp[i].next)
        {
            v=mp[i].v;
            if(dis[v]==-1&&mp[i].cap>0)
            {
                dis[v]=dis[u]+1;
                q.push(v);
                if(v==tp)
                    break;
            }
        }
    }
    return dis[tp]!=-1;
}
ll dfs(int u,ll cap)//寻找增广路
{
    if(u==tp||cap==0)
        return cap;
    ll res=0,f;
    for(int &i=cur[u]; i!=-1; i=mp[i].next)
    {
        int v=mp[i].v;
        if(dis[v]==dis[u]+1&&(f=dfs(v,min(cap-res,mp[i].cap)))>0)
        {
            mp[i].cap-=f;
            mp[i^1].cap+=f;
            res+=f;
            if(res==cap)
                return cap;
        }
    }
    if(!res)
        dis[u]=-1;
    return res;
}
ll dinic()
{
    ll ans=0;
    while(bfs())
    {                                             **************关键问题*************
        for(int i=1; i<=n; i++)         注意这里的所有点第一个下标为1，如果有超级源点下标为0需要更改
                                                    因为这里为n的循环，所以主函数里面n的值需要给出一个尽可能小的n
            cur[i]=pre[i];                    n为网络图中最大可能的点的个数
        ans+=dfs(sp,inf);
    }
    return ans;
}
int main()
{
    while(~scanf("%d%d",&m,&n))    输入边的条数m和点的个数n ，如果没有输入需要手动赋值，因为dinic中有关于n的循环
                                                         
    {
        int u,v,w;
       init();
        n=**;  //点的最大个数
        sp=**;  tp=**;    //赋值起点的位置和终点的位置
       
        for(int i=1;i<=m;i++)          ****自定义----->构建网络图连边
        {   scanf("%d%d%d",&u,&v,&w);      //依次输入m条边的起点，终点和容量
            add(u,v,w);
            add(v,u,0);                  //记得加容量为0的反向边
        }    
                          
        ll ans=dinic();
        printf("%lld\n",ans);
        return 0;
    }

***输出二分图匹配边的每组匹配点的情况

for(int i=0;i<=cnt;i=i+2)
    {
        if(mp[i].v!=sp&&mp[i^1].v!=sp)
        if(mp[i].v!=tp&&mp[i^1].v!=tp)
        if(mp[i^1].cap>0)
        {
            printf("%d %d\n",mp[i^1].v,mp[i].v);
        }
    }