为什么反向边容量设为0样例都过不了
查看原帖
为什么反向边容量设为0样例都过不了
556955
123wbl楼主2022/10/29 13:14
#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
struct edge
{
    int to,next,dis;
}e[6000001];
int n,m,head[1000001],cur[1000001],tot=1,anslow;
int dep[1000001],b[1000001],s=1,t;
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while (ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
	while (ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-48;ch=getchar();}
	return x*f;
}
inline void add(int u,int v,int w)
{
    e[++tot].to=v;
    e[tot].next=head[u];
    e[tot].dis=w;
    head[u]=tot;
}
inline bool bfs()
{
    queue<int> q;
    memset(dep,0x3f,sizeof(dep));
    memset(b,0,sizeof(b));
    memcpy(cur,head,sizeof(head));
    q.push(s);
    b[s]=1;
    dep[s]=0;
    while(!q.empty())
    {
        int u=q.front();
        b[u]=0;
        q.pop();
        for(int i=head[u];i;i=e[i].next)
        {
            int v=e[i].to;
            if(dep[v]>dep[u]+1&&e[i].dis)
            {
                dep[v]=dep[u]+1;
                if(!b[v])
                {
                    q.push(v);
                    b[v]=1;
                }
            }
        }
    }
    if(dep[t]!=0x3f3f3f3f) return true;
    else return false;
}
inline int dfs(int u,int flow)
{
    int used=0,rlow;
    if(u==t)
    {
        anslow+=flow;
        return flow;
    }
    for(int i=cur[u];i;i=e[i].next)
    {
        cur[u]=i;
        int v=e[i].to;
        if(e[i].dis&&dep[v]==dep[u]+1)
        {
            if(rlow=dfs(v,min(e[i].dis,flow-used)))
            {
                used+=rlow;
                e[i].dis-=rlow;
                e[i^1].dis+=rlow;
                if(used==flow)
                    break;
            }
        }
    }
    return used;
}
inline void dinic()
{
    while(bfs())
        dfs(s,0x7fffffff);
}
int main()
{
    n=read(),m=read();
    t=n*m;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m-1;j++)
        {
            int w;
            w=read();
            add((i-1)*m+j,(i-1)*m+j+1,w);
            add((i-1)*m+1+j,(i-1)*m+j,0);
        }
    for(int i=1;i<=n-1;i++)
        for(int j=1;j<=m;j++)
        {
            int w;
            w=read();
            add((i-1)*m+j,i*m+j,w);
            add(i*m+j,(i-1)*m+j,0);
        }
    for(int i=1;i<=n-1;i++)
        for(int j=1;j<=m-1;j++)
        {
            int w;
            w=read();
            add((i-1)*m+j,i*m+j+1,w);
            add(i*m+j+1,(i-1)*m+j,0);
        }
    dinic();
    cout<<anslow;
}
2022/10/29 13:14
加载中...