максимальное сопоставление с Эдмондом Карпом

#c #matching #bipartite #max-flow

#c #сопоставление #двудольный #максимальный поток

Вопрос:

Я хотел бы спросить, может ли кто-нибудь сказать мне, что не так с моим кодом. Я пытаюсь решить проблему BerSU Ball на codeforces с Эдмондом-Карпом (я знаю о жадном солне)

Проблема: https://codeforces.com/problemset/problem/489/B

Я получаю неправильный ответ на тест 37. Я не могу найти свою ошибку. Вот моя реализация (код чистый, адаптированный из cp-алгоритмов)

 #include<vector>
#include<iostream>
#include<queue>
using namespace std;
typedef vector<int> vi;
typedef vector<vi> vii;
const int INF=1e9 5;
int maxflow(int source, int sink , const viiamp; graph , viiamp; capacity);
int main(){
    int n;
    scanf("%d",amp;n);
    vi boys(n); 

    for(int i = 0 ; i <n ;i  ){
        scanf("%d",amp;boys[i]);
    }
    int m;

    scanf("%d",amp;m);
    vi girls(m);
    for(int i = 0 ; i < m ;i  ){
        scanf("%d",amp;girls[i]);
    }
    int len = n m 2;
    int source = len-2;
    int sink = len-1;
    vii graph(len);
    vii capacity(len , vi(len,0));
    for(int i = 0 ; i < n ; i  ){
        for(int j  = 0 ; j < m ; j  ){
            if(abs(boys[i]-girls[j]) <= 1){
                graph[i].push_back(j n);
                graph[n j].push_back(i);
                capacity[i][j n]=1;
                capacity[j n][i]=1;
            }
        }
    }
    for(int i = 0 ; i < n ;i  ){
        graph[i].push_back(source);
        graph[source].push_back(i);
        capacity[i][source]=1;
        capacity[source][i]=1;
    }


    for(int i = 0 ; i < m ; i  ){
        graph[sink].push_back(n i);
        graph[n i].push_back(sink);
        capacity[sink][n i]=1;
        capacity[n i][sink]=1;
    }
    // for(int i = 0 ; i < len; i  ){
    //      cout << i << " : ";
    //      for(int j = 0 ; j < len ; j  ){
    //          cout << capacity[i][j] << " ";
    //      }
    //      cout << endl;
    //  }

    cout << maxflow(source,sink,graph,capacity) << endl;

    return 0;
}


int bfsEdmond(int source , int sink , const viiamp; graph , viiamp; capacity , viamp; parent){
    int len = (int)graph.size();
    parent.assign(len,-1);
    parent[source]=-2;
    queue<pair<int,int>>q;
    q.push({source,INF});
    while(!q.empty()){
        int cnode = q.front().first;
        int cflow=q.front().second;
        q.pop();
        //cout << "exploring " << cnode << endl;
        for(int next: graph[cnode]){
            //cout << next << " " << endl;
            if(parent[next] == - 1 amp;amp; capacity[cnode][next]){
                // explore adjacency list
                parent[next]=cnode;
                int newflow = min(cflow,capacity[cnode][next]);
                if(next == sink){
                    return newflow;
                }
                q.push({next,newflow}); 
            }   
        }
    }
    return 0;
}

int maxflow(int source, int sink , const viiamp; graph , viiamp; capacity){
    int len = (int)graph.size();
    vi parent(len);
    int maxx=0;
    int cflow=0;
    while( (cflow = bfsEdmond(source,sink,graph,capacity,parent))){
        maxx =cflow;
        int cnode = sink;
        while(cnode!=source){
            int prev = parent[cnode];
            capacity[prev][cnode]-=cflow;
            capacity[cnode][prev] =cflow;
            cnode=prev;
            //cout << "dine " << endl;
        }
    }
    return maxx;
}