Гамильтонов цикл.

[SrgGld]

Нужно было написать программу для нахождения гамильтонова цикла в графе из n вершин.
Я пока пытался писать поиск гамильтонова пути в графе, в котором заведомо есть гамильтонов цикл.
Но для n=1 и n=2 программа выдаёт не то что нужно. А для n=3 вообще иногда с переполнением стека вылетает.

Алгоритм пытался реализовать такой:
Запрашивается количество вершин.
Запрашивается матрица смежности.
В функцию передаётся первая точка.
Печатается сформированный функцией путь в массиве "c".
В функции:

если точка уже используется (содержится в пути), ретёрн фэлс

если точка - тупик (не имеет смежных не использованных) и не последняя оставшаяся - ретёрн фэлс

если точка - тупик (не имеет смежных не использованных) и последняя - добавляем её в путь и ретёрн труЪ

добавляем точку в путь

перебираем циклом смежные с ней вершины:

для каждой из них вызываем эту же функцию:

если она возвращает труЪ то ретёрн труЪ

ну, а нет так нет -передаём ей следующую подходящую точку

// Гамильтонов цикл.cpp: определяет точку входа для консольного приложения.
//
#include "stdafx.h"
#include <iostream>
using namespace std;

bool find(int *c, int &n, int &x)
{
	for (int i=0; i<n; i++) 
		if (c[i]==x) return true;
	return false;
}
bool deadlock(int *aj, int *c, int n)
{
	for (int i=0; i<n; i++) {
		if (aj[i] && !find(c,n,i))//при обнаружении первой неиспользованной смежной точки
			 return false;
	}
	return true;
}
bool path (int **a, int n, int i, int *c)//не правильно
{
	if (find(c,n,i))
		return false;	// точка уже используется
	if (deadlock(a[i],c,n)) //точка не имеет смежных неиспользованных точек
		if (find(c,n,n)!=n-1) //путь ещё не завершается
			return false;
		else {
			c[n-1]=i;
			return true;
		}
	c[find(c,n,n)]=i;//данная точка помещается в конец уже найденного пути
	for (int j=0; j<n; j++) {
		if (a[i][j]!=1) continue;
		if	(path(a, n, j, c))
			return true;
		//else
		//	c[find(c,n,i)+1]=n;
	}
	return false; //если не осталось смежных неиспользованных точек у данной возвращаемся назад для следующей итерации
}

int _tmain(int argc, _TCHAR* argv[])
{
	int n;
	cout<<"The program is correct for correct graph."<<endl;
	cout<<"Enter quantity of points!"<<endl;
	cin>>n;
	cout<<"Enter contiguity matrix."<<endl;
	int **a=new int * [n];
	for (int i=0; i<n; i++) {
		a[i]=new int [n];
		for (int j=0;j<n; j++)
			cin>>a[i][j];
	}
	int *c= new int [n];//массив для хранения использованных вершин в порядке обхода
	for (int i=1; i<n; i++)
		c[i]=n;//в неиспользованной части массива будут невозможные значения
	bool f=path(a, n, 0, c) ;
		for (int i=0; i<n; i++)
			cout<<c[i]<<' ';
	return 0;
}

Может, что посоветуете?