Рекурсия при построении матрицы достигаемости

[[MI_nor]]

Дана матрица смежности некоторого графа из н узлов.
Нужно написать рекурсивную программу построения матрицы достижимостидля этого графа:

#include <stdio.h>
#include <stdlib.h>
#include <conio.h>

int matrs[50][50];
int matrd[50][50];

int i,j,k=0,l=0,n;

int add(int matr[50][50])
 {
  int a;
  printf("Vvedite razmernost matrici - ");
  scanf("%d",&n);
  printf("Zapolnit vruchnuiu(1) ili avtomaticheski(2)? - ");
  scanf("%d",&a);
  if (a!=1 && a!=2) {
		     while (a!=1 && a!=2){
					  printf("Incorrect value. Must be 2 or 1\n");
					  scanf("%d",&a);
					 }
		    }
  switch (a)
    {
     case 1:
	    {
	     printf("Vvedite matricy smegnosti - \n");
	     for (i=0;i<n;i++)
	      {
	       for (j=0;j<n;j++)
	       {
		scanf("%d",&matr[i][j]);
		if (matr[i][j]!=0 && matr[i][j]!=1) {
						     printf ("\nIncorrect element [%d,%d](not 0 or 1)\n",i,j);
						     while (matr[i][j]!=0 && matr[i][j]!=1)
						      {scanf(" %d",&matr[i][j]);}
						    }
	       }
	      }
	    }break;
     case 2: {
	      randomize();
	      for (i=0;i<n;i++)
	       {for (j=0;j<n;j++)
		 {
                  if (i==j) matr[i][j]==0; else
		  matr[i][j]=random(2);
		 }
	       }
	     } break;
    }
  return 1;
 }

int rekursiya(int a, int b)
 {
  for(i=a;i<n;i++)
   {
    for (j=b;j<n;j++)
     {
      if (matrs[i][j]==1) {
			   matrd[i][j]=1;
                                     if (matr[i][j]==1 && matr[i][j]==matr[k][l]) matr[k][j]=1;
			   k=i; l=j;
			   rekursiya(j,0);
			  }
     }
   }
 }

void main()
 {
  clrscr();
  add(matrs);
  clrscr();
  printf("Vasha matrica smegnosti - \n");
  for (i=0;i<n;i++)
   {
    for (j=0;j<n;j++)
     {printf("%d ",matrs[i][j]);}
    printf("\n");
   }
  rekursiya(0,0);
  printf("\nVasha matrica dostigaemosti - \n");
  for (i=0;i<n;i++)
   {
    for (j=0;j<n;j++)
     {printf("%d ",matrd[i][j]);}
    printf("\n");
   }

  getch();
 }

Собственно интересен алгоритм рекурсии) я так понимаю если в матр смежности 2 точки соединены напр дугой, то идем от второй дальше по возможным, параллельно ставя 1 в элементах первого, граничящего со 2м , итд...

Элемент матрицы смежности а(и,жи)=1, если существует дуга от вершины и к вершине жи, и
а(и,жи)=0 в противном случае. элемент матрицы достижимости р(и,жи)=1, если из вершины и можно попасть в вершину
жи, двигаясь по стрелкам, и р(и,жи)=0 в противном случае.