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

#define N 50000 

int main(int argc, char *argv[]) {
    float a[N];
    float b[N];

    int ind[N]; 
    int i, j;

    double begin_t, end_t;

//    srand(time(NULL));

    for (i=0; i<N; i++)
    a[i]= (float) rand() / RAND_MAX;

    begin_t=omp_get_wtime();

    #pragma omp parallel for shared(a, b, ind) private(i,j)
    for (i=0; i<N; i++) {
        ind[i]=0;
//    #pragma omp parallel for shared(a, b, ind, i) private(j)
        for (j=0; j<N; j++) 
            if ((a[i]>a[j]) || ((a[i]==a[j]) && (i>j)))
                ind[i]++;
        b[ind[i]]=a[i]; 
    }

    end_t=omp_get_wtime();

    printf("a:\n"); 
    for (i=0; i<10; i++)
    printf("%10.6f", a[i]);
    printf("\nb:\n"); 
    for (i=0; i<10; i++)
    printf("%10.6f", b[i]);
    printf("\n");

    printf("Czas obliczen: %lf s\n",end_t-begin_t);
    return 0;
}