#include <stdio.h>
#include <math.h>
#define true 1
#define false 0
#define and &&
#define or ||
#define boolean int

#define sqr(a) ((a)*(a))

#define INFINITY_LONG 10000000L

static r_seed = 12345678L;


double uniform ()
{
  static long m = 2147483647;
  static long a = 48271;
  static long q = 44488;
  static long r = 3399;
  long t, lo, hi;

  hi = r_seed / q;
  lo = r_seed - q * hi;
  t = a * lo - r * hi;
  if (t > 0)
    r_seed = t;
  else
    r_seed = t + m;
  return ( (double) r_seed / (double) m );

}

double* makeRandomArray (int length)
{
  int i;

  double *A = (double*) malloc (sizeof(double) * length);
  for (i=0; i<length; i++) {
    A[i] = uniform ();
  }
  return A;
}

void algorithm1 (double *A, int length)
{
  double max = fabs (A[0] - A[1]);
  double x = A[0];
  double y = A[1];
  int i, j;
  double diff;

  for (i=0; i<length-1; i++){
    for (j=i+1; j<length; j++) {
      diff = fabs (A[i] - A[j]);
      if (diff > max) {
        max = diff;
        x = A[i];
        y = A[j];
      }
    }
  }
  printf ("Algorithm 1: the numbers: %lf,%lf  difference=%lf\n", x, y, max);
}

int main ()
{
  double *A;
  int i;

  A = makeRandomArray (10000);
  algorithm1 (A, 10000);
}