Не могу грамотно написать цикл поиска (зацикливается)

#include <iostream>
#include <conio.h>
#include <math.h>
#include <fstream>
#include <vector>
using namespace std;
 
int main()
{
    int d = 0,y = 0,max = 0;
    int v = 0,e = 0,k,f = 0;
    int raz = 4;
    bool gh=true, dach=false;
    bool matrix[20][20];
    bool dad[10];
    bool Matrix[30][30];
    setlocale(0,"Russian");
    cout << "\t\t***Поиск вершинного покрытия***\n\n";
    cout << "Ввод из файла matrix.txt\n\n";
    
    ifstream file;
    file.open("matrix.txt");
 
    if(!file.is_open())
    {
        cout << "Не удалось открыть файл matrix.txt\n";
        _getch();
        return -1;
    }
 
    while (!file.eof())
    {
        for(int i = 0;i < raz;i++)
        {   
            for(int j = 0;j < raz;j++)
            {
                file >> v;
                cout<<v<<"  ";
                 matrix[i][j] = v;
            }
        cout<<endl;
        }
        cout <<endl;
    }
 
    cout << "Граф заданный матрицей инцедентности:\n"<<endl;
 
        for(int i = 0;i < raz; i++)
        {   
            for(int j = i + 1;j < raz; j++)
            {
                if(matrix[i][j] == true)
                    e = e + 1;
                
            }
        }
 
        for(int i = 0;i < raz;i++)
        {   
            for(int j = 0;j < e;j++)
            {
                Matrix[i][j] = false;
            }
        }
 
        for(int i = 0;i < raz; i++)
        {   
            for(int j = i + 1;j < raz; j++)
            {
                
                if(matrix[i][j] == true)
                {
                    Matrix[i][f] = true;
                    Matrix[j][f] = true;
                    f++;
                }
            }
        }
 
 
 
        for(int i = 0;i < raz;i++)
        {   
            for(int j = 0;j < e;j++)
            {
 
                cout<<Matrix[i][j]<<"  ";
            }
        cout<<endl;
        }
 
 
while(gh)
{
 
 
    for(int i = 0;i < raz;i++)
    {   
        for(int j = 0;j < e;j++)
        {
            if(Matrix[i][j] == true)
            {
                d++;
            }
        }
        if(y <= d)
        {
            y = d;
            max = i;
        }
        d = 0;
    }
    //cout<<y<<endl;
    //cout<<max<<endl;
    
    dad[max]=true;
 
 
 
    for(int i = 0;i < raz;i++)
    {   
        for(int j = 0;j < e;j++)
        {
            if(Matrix[max][j]==true)
            {
                Matrix[max][j]=false;
                for(int i = 0;i < raz;i++)
                {   
                    Matrix[i][max]=false;                       
                }
            }
        }
    }
 
    for(int i = 0;i < raz;i++)
    {   
        for(int j = 0;j < e;j++)
        {
            if(Matrix[i][j]==true)
                dach=true;
        }
    }
    if(dach==false)
    {
        gh = false;
    }
    dach=false;
}
    
    for(int i=0;i<raz;i++)
    {
        if(dad[i]==true)
            cout<<i+1;
    }
    file.close();
    _getch();
    return 0;
}