寒假Day16Dinic模板更新+优化

之前的模板好像比较简单，而且没进行cur优化，可能一些大数据过不去，并且建边也会有问题

貌似之前的Dinic模板适用范围有局限性

tot初始化为-1

 

Dinic模板：（邻接表add建边+dfs+bfs）

struct node
{
    // int u,v,flow;
    int nextt,v,flow;
} e[N*N];

void add(int u,int v,int flow)
{
    tot++;
    //nextt[tot]=head[u];
    e[tot].nextt=head[u];
    head[u]=tot;
    // e[tot].u=u;
    e[tot].v=v;
    e[tot].flow=flow;

    tot++;
    // nextt[tot]=head[v];
    e[tot].nextt=head[v];
    head[v]=tot;
    //e[tot].u=v;
    e[tot].v=u;
    e[tot].flow=0;
}

int dfs(int u,int flow)
{
    if(u==t)
        return flow;
    for(int &i=cur[u]; i!=-1; i=e[i].nextt) //注意这里的&符号，这样i增加的同时也能改变cur[u]的值，达到记录当前弧的目的
    {
        if((dep[e[i].v]==dep[u]+1)&&e[i].flow>0)
        {
            int di=dfs(e[i].v,min(flow,e[i].flow));
            if(di>0)
            {
                e[i].flow-=di;
                e[i^1].flow+=di;
                return di;
            }
        }
    }
    return -1;
}

bool bfs()
{
    if(s==t)return 0;
    queue<int>Q;
    while(!Q.empty())
        Q.pop();
    memset(dep,-1,sizeof(dep));
    dep[s]=1;
    Q.push(s);
    while(!Q.empty())
    {
        int u=Q.front();
        Q.pop();
        for (int i=head[u]; i!=-1; i=e[i].nextt)
        {
            if ((e[i].flow>0)&&(dep[e[i].v]==-1))
            {
                dep[e[i].v]=dep[u]+1;
                Q.push(e[i].v);
            }
        }
    }
    if(dep[t]!=-1)
        return 1;
    return 0;
}

int dinic()
{
    int sum=0;
    while(bfs())
    {
        for(int i=s; i<=t; i++)
            cur[i]=head[i];
        int di;
        while((di=dfs(s,inf))!=-1)
            sum+=di;
    }
    return sum;
}