为什么第二次Dinic会走反向边？

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
constexpr int maxn=505;
int n;
int arr[maxn];
int f[maxn];
int d[maxn*3];
int s;
int cur[maxn*3];
constexpr int S=0;
constexpr int T=2*maxn+50;
constexpr int INF=1e9;
constexpr int getright(int a)
{
    return a+maxn;
}
struct edge
{
    int u;
    int v;
    int cap;
    int flow;
};
vector<edge> eds;
vector<int> G[maxn*3];
void insed(int u,int v,int cap)
{
    //  static int de{};
    //  if(v==T)cout<<"\t"<<++de<<"\t"<<endl;
    eds.push_back({u,v,cap,0});
    eds.push_back({v,u,0,0});
    int m=eds.size();
    G[u].push_back(m-2);
    G[v].push_back(m-1);
}
bool bfs()
{
    static bool vis[maxn*3];
    memset(vis,0,sizeof vis);
    queue<int> q;
    q.push(S);
    d[S]=0;
    vis[S]=true;
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        for(int i=0;i<G[x].size();i++)
        {
            edge& e=eds[G[x][i]];
            if(!vis[e.v]&&e.cap>e.flow)
            {
                vis[e.v]=true;
                d[e.v]=d[x]+1;
                q.push(e.v);
            }
        }
    }
    return vis[T];
}
int dfs(int x,int a)
{
    //  static int clo{};
    //  cout<<"\t"<<++clo<<"\t"<<endl;
    if(x==T||a==0)return a;
    int flow=0,f;
    for(int& i=cur[x];i<G[x].size();i++)
    {
        edge& e=eds[G[x][i]];
        if(e.cap==0)
            continue;
        if(d[e.v]=d[x]+1&&(f=dfs(e.v,min(a,e.cap-e.flow)))>0)
        {
            flow+=f;
            eds[G[x][i]^1].flow-=f;
            a-=f;
            e.flow+=f;
            if(!a)break;
        }
    }
    return flow;
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        cin>>arr[i];
        f[i]=1;
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<i;j++)
            if(arr[j]<=arr[i])f[i]=max(f[i],f[j]+1);
    s=*max_element(f,f+n+1);
    cout<<s<<endl;
    if(s==1)
    {
        cout<<n<<endl<<n;
        return 0;
    }
    for(int i=1;i<=n;i++)
    {
        if(f[i]==1)insed(S,i,1);
        else if(f[i]==s)insed(getright(i),T,1);
        insed(i,getright(i),1);
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<i;j++)
            if(arr[j]<=arr[i]&&f[j]+1==f[i])insed(getright(j),i,1);
    int flow=0;
    while(bfs())
    {
        memset(cur,0,sizeof cur);
        flow+=dfs(S,INF);
    }
    cout<<flow<<endl;
    insed(S,1,INF);
    insed(1,getright(1),INF);
    if(f[n]==s)
    {
        insed(getright(n),T,INF);
        insed(n,getright(n),INF);
    }
    while(bfs())
    {
        memset(cur,0,sizeof cur);
        flow+=dfs(S,INF);
    }
    cout<<flow;
}

这是我的AC代码，但是虽然AC了，却是在爆栈无数次以后加了一行谁都没加的神奇代码

if(e.cap==0)
            continue;

A的。因为在我调试之后发现在第二次Dinic的时候在边(0,1)和边(1,0)之间一直处理。但在

if(d[e.v]=d[x]+1&&(f=dfs(e.v,min(a,e.cap-e.flow)))>0)
        {
            flow+=f;
            eds[G[x][i]^1].flow-=f;
            a-=f;
            e.flow+=f;
            if(!a)break;
        }

中加入反向边断言能成功通过。（i.e. if体中没有出现处理反向边情况）

有大佬能解释这个现象吗？

另外A题之后看了一下题解发现别人都是直接存储残量网络（我是个异端），是不是问题出在这里？