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

int *populate ( int n )
{
   int i, *A;
   
   A = (int *)malloc(n * sizeof(int));
   for (i=0; i<n; ++i) A[i] = rand() % n;
   return A;
}

#ifdef __VERBOSE
void print ( int *A, int n )
{
   int i;

   for (i=0; i<n; ++i) {
      printf("%7d", A[i]);
      if (i % 10 == 9) printf("\n");
   }
   if (i % 10) printf("\n");
}
#endif

void cycle ( int *A, int n, int s )
{
   int *B, k, r;

   B = (int *)malloc(n * sizeof(int));
   for (k=0; k<n; ++k) B[k] = 0;

   r = s;
   while (B[r] == 0) {
      B[r] = 1;
      r = A[r];
   }

   printf("Initial part:\n"); k = 0;
   while (s != r) {
      printf("%7d", s); ++k; if (k % 10 == 0) printf("\n");
      s = A[s];
   }
   if (k % 10) printf("\n");

   printf("Cyclic part:\n");
   printf("%7d", s); k = 1;
   r = A[s];
   while (r != s) {
      printf("%7d", r); ++k; if (k % 10 == 0) printf("\n");
      r = A[r];
   }
   if (k % 10) printf("\n");
   printf("Length of the cycle = %d\n", k);

   free(B);
}

int main ()
{
   int n, *A, s;

   srand((unsigned int)time(NULL));
   printf("n = "); scanf("%d", &n);
   if (n < 10) n = 1000;
   A = populate(n);
#ifdef __VERBOSE
   print(A,n);
#endif
   printf("s = "); scanf("%d", &s);
   cycle(A,n,s);
   exit(0);
}