Поиск гамильтоновых циклов в графе

matrix[i][line]=1;

//---------------------------------------------------------------------------
 
#include <vcl.h>
#include <stdio.h>
#pragma hdrstop
 
#include "Unit1.h"
#include "about.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
 
int numev=0;    //счетчик, необходимый для смены состояний построения графа
int xy[2][50];  //массив для хранения координат вершин на экране
int kv = 0;     //количество вершин
int line;       //переменная, необходимая для формирования матрицы смежности
bool select=false;    // отметка что вершину нашли
char **matrix;          //выделение памяти для будущей матрицы смежности
AnsiString way; //строка для вывода результатов в файл и на экран
int *buf;       //указатель на буфер
bool mark=false;     // удаление вершины
int kvdel=0;
 
 
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
 Status->Caption="Статус:\n\n...отмечаем вершины...";      //начальный статус
}
//---------------------------------------------------------------------------
 
// Функция проверки на повторы вершин в цепи
int __fastcall Repeat(int i, int num)
{
        for (int j=0; j<num; j++)       //цикл поиска такой же вершины
        {
        if (buf[j]==i) return 1;        //если  такая найдена, то повтор есть
        }
    return 0;   //нет повтора
}
// Функция для поиска новых вершин в цепи
 
void __fastcall Ways(int ver, int num)
{       select=false;
    if (mark==true) //если вершины уже удалялись
    {
        kvdel=buf[num]; //начинать отсчет вершин с
    }
    else kvdel=0; //иначе обнулить этот счетчик
    mark=false;
    for (int i=kvdel; i<kv; i++)        //цикл поиска возможных вершин
    {
        if (matrix[ver][i]==0) continue; //если нет ребра, следующий
        if (Repeat(i,num)) continue;    //если вершина уже была, следующий
        if (i==buf[num]) continue;
        ver=i;
        select=true; //отметка, что вершину нашли
        break;
    }
 
    if (select==true)   //если вершина найдена
    {
        buf[num]=ver;   //добавить ее в цепь
                num++;  //увеличить счетчик
                Ways(ver, num); //рекурсивно вызвать с новыми параметрами
     }
    else        //если вершина не найдена
    {
        if (num==kv)    //если длина цепи равна количеству вершин
        {
                way=buf[0]+1;
            for (int i=1; i<kv; i++)
            {
                way+=',';
                        way+=buf[i]+1;
            }
            if (matrix[buf[kv-1]][0]!=0) {way+=",1";
            Form1->Results->Items->Add(way); }  //вывод результата на экран
 
 
        }
       way=buf[0]+1;
        for (int i=1; i<num; i++)
            {
                way+=',';
                way+=buf[i]+1;
 
            }
            way+='\n';
        num--;  //уменьшаем счетчик
        mark=false;
        if (num!=0)     //если счетчик вершины не равен 0
        {
if (num!=kv-1) buf[num+1]=NULL; //и не равен количеству вершин, то обнулить
// следующую ячейку буфера
            mark=true;  //отметка о удалении вершины
            Ways(buf[num-1], num);      //поиск следующей вершины
        }
 
    }
 
}
 
 //----------------------------------------------------------------------------
//Функция поиска гамильтоновых циклов
void __fastcall Gamilton(void)
{
    buf=new int[kv]; //создание буфера
    buf[0]=0;
    for (int i=0; i<kv; i++) //обнуление буфера
    {
        buf[i]=NULL;
    }
    Ways(0, 1); //вызов функции для поиска новых вершин в цепи
 
}
//---------------------------------------------------------------------------
 
//---------------------------------------------------------------------------
 
//---------------------------------------------------------------------------
void __fastcall TForm1::DeskMouseDown(TObject *Sender, TMouseButton Button,
      TShiftState Shift, int X, int Y)
{               Label2->Visible=false;
        if (numev==0)   //если нажатий на кнопку «Далее» не производилось, то ставятся вершины
        {
 
                Desk->Canvas->Pen->Color=clGreen;       //цвет вершины = зеленный
                Desk->Canvas->Pen->Width=10;    //толщина пера=10 пикселей
                Desk->Canvas->MoveTo(X,Y);      //ставится точка
                Desk->Canvas->LineTo(X,Y);
                xy[0][kv]=X; //координаты вершины заносятся в массив
                xy[1][kv]=Y;
                Desk->Canvas->TextOutA(X+10, Y+10, AnsiString(kv=kv+1));
                //рядом выводится номер вершины
        }
        if (numev==1) //построение ребер
        {
                for (int i=0; i<kv; i++)
            {
                if ((X>(xy[0][i]-5))&&(X<(xy[0][i]+5)))
                {
                        if ((Y>(xy[1][i]-5))&&(Y<(xy[1][i]+5)))
                    {
                        if (select==false) //если вершина отмечена первой
                        {
                            Desk->Canvas->Pen->Color=clRed; //отметить ее красным цветом
                            Desk->Canvas->Pen->Width=7;
                            Desk->Canvas->MoveTo(xy[0][i],xy[1][i]);
                            Desk->Canvas->LineTo(xy[0][i],xy[1][i]);
                            line=i; //занести в переменную номер вершины
                            select=true; //поставить отметку, что одна вершина уже выделена
                        }
                        else //иначе
                        {
                            Desk->Canvas->Pen->Color=clBlack;
                            Desk->Canvas->Pen->Width=10;
                            Desk->Canvas->MoveTo(xy[0][line],xy[1][line]);
                            Desk->Canvas->LineTo(xy[0][line],xy[1][line]);
                            Desk->Canvas->Pen->Width=3;
                            Desk->Canvas->LineTo(xy[0][i],xy[1][i]); //провести линию
  /*постоянно тут ошибка*/    matrix[line][i]=1; //поставить отметку о ребре в матрице смежности
                            matrix[i][line]=1;
                            select=false; //ребро построено, выделенных вершин нет
                        }
                        break;
                    }
                }
            }
        }
}
//---------------------------------------------------------------------------
void __fastcall TForm1::BNextClick(TObject *Sender)
{
if (numev==0) //если нажатие производится в первый раз
{
    matrix=new char* [kv];      //создать массив для матрицы смежности
    for (int i=0; i<kv; i++) //обнулить этот массив
    {
        matrix[i]=new char [kv];
        for (int j=0; j<kv; j++)
        {
                matrix[i][j]=0;
        }
    }
        Status->Caption="Статус:\n\n...отмечаем ребра..."; //сменить статус
 
}
if (numev==1) //повторное нажатие на кнопку
{
    Status->Caption="Статус:\n\n Готово!"; //смена статуса
    BNext->Enabled=false;
    Gamilton(); //вывов функции Gamilton
            //если нет циклов
      if (Form1->Results->Items->Count == 0)   Form1->Results->Items->Add(" циклов нет");
}
numev++; //увеличение счетчика событий на единицу
 
 
}
//---------------------------------------------------------------------------
// новый граф
void __fastcall TForm1::Button1Click(TObject *Sender)
{     //очистка компонента
 Desk->Picture->Bitmap->FreeImage();
 Desk->Picture->Bitmap = NULL;
 Label2->Visible=true;
 // нулевые значения
 numev=0; // кол нажатий кнопки нехт
 BNext->Enabled=true; // видимость кнопки
 Status->Caption="Статус:\n\n...отмечаем вершины...";      //начальный статус
kv=0;
select=false;
char **matrix;          //выделение памяти для будущей матрицы смежности
AnsiString way; //строка для вывода результатов в файл и на экран
int *buf;       //указатель на буфер
mark=false;
kvdel=0;
Results->Items->Clear();
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button2Click(TObject *Sender)
{
// окно о программе
Form2->Show();
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button3Click(TObject *Sender)
{
// кнопка ДОБАВИТЬ ВЕРШИНУ
numev=0;
        Status->Caption="Статус:\n\n...отмечаем вершины..."; //сменить статус
 
        matrix=new char* [kv];      //создать массив для матрицы смежности
    for (int i=0; i<kv; i++) //обнулить этот массив
    {
        matrix[i]=new char [kv];
        for (int j=0; j<kv; j++)
        {
                matrix[i][j]=0;
        }
    }
 
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button5Click(TObject *Sender)
{
// кнопка ДОБАВИТЬ РЕБРО
numev=1;
Status->Caption="Статус:\n\n ...отмечаем ребра...!"; //смена статуса
 
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button7Click(TObject *Sender)
{
// кнопка ОК
    Status->Caption="Статус:\n\n Готово!"; //смена статуса
     Results->Items->Clear();// очищаем старый список
    Gamilton(); //вывов функции Gamilton
     if (Form1->Results->Items->Count == 0)   Form1->Results->Items->Add(" циклов нет");
}

First chance exception at $7742B727. Exception class EAccessViolation with message 'Access violation at address 00402B85 in module 'Project1.exe'. Write of address 01010100'. Process Project1.exe (15080)

matrix[i][line]=1;

void __fastcall TForm1::BNextClick(TObject *Sender)
{
if (numev==0) //если нажатие производится в первый раз
{
    matrix=new char* [kv];      //создать массив для матрицы смежности
    for (int i=0; i<kv; i++) //обнулить этот массив
    {
        matrix[i]=new char [kv];
        for (int j=0; j<kv; j++)
        {
                matrix[i][j]=0;
        }
    }
        Status->Caption="Статус:\n\n...отмечаем ребра..."; //сменить статус
 
}
if (numev==1) //повторное нажатие на кнопку
{
    Status->Caption="Статус:\n\n Готово!"; //смена статуса
    BNext->Enabled=false;
    Gamilton(); //вывов функции Gamilton
            //если нет циклов
      if (Form1->Results->Items->Count == 0)   Form1->Results->Items->Add(" циклов нет");
}
numev++; //увеличение счетчика событий на единицу

    if (matrix[ver][i]==0) continue; //если нет ребра, следующий

char matrix[1000][1000]; //заранее известный размер

char **matrix;          //выделение памяти для будущей матрицы смежности

matrix=new char* [kv];