Untitled

// -*- mode: c++ -*-// $Id: heapsort.g++,v 1.6 2001/08/01 23:41:52 doug Exp $// http://www.bagley.org/~doug/shootout///** This copy courtesy of http://www.kano.net/javabench///**  Added Timer() code#include <iostream>#include <math.h>#include <cstdlib>#include <ctime>using namespace std;#define IM 139968#define IA   3877#define IC  29573doublegen_random(double max) {    static long last = 42;    return( max * (last = (last * IA + IC) % IM) / IM );}voidheapsort(int n, double *ra) {    int i, j;    int ir = n;    int l = (n >> 1) + 1;    double rra;    for (;;) {	if (l > 1) {	    rra = ra[--l];	} else {	    rra = ra[ir];	    ra[ir] = ra[1];	    if (--ir == 1) {		ra[1] = rra;		return;	    }	}	i = l;	j = l << 1;	while (j <= ir) {	    if (j < ir && ra[j] < ra[j+1]) { ++j; }	    if (rra < ra[j]) {		ra[i] = ra[j];		j += (i = j);	    } else {		j = ir + 1;	    }	}	ra[i] = rra;    }}intmain(int argc, char *argv[]) {    int N = ((argc == 2) ? atoi(argv[1]) : 1);    clock_t start=clock();     double *ary;    int i;        /* create an array of N random doubles */    ary = (double *)malloc((N+1) * sizeof(double));    for (i=1; i<=N; i++) {	ary[i] = gen_random(1);    }    heapsort(N, ary);    printf("%.10f\n", ary[N]);    free(ary);    clock_t endt=clock(); 	std::cout <<"Time: " << double(endt-start)/CLOCKS_PER_SEC * 1000 << " ms\n";     return(0);}