Cod sursa(job #1981664)

Utilizator MaligMamaliga cu smantana Malig Data 16 mai 2017 13:38:42
Problema Radix Sort Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.31 kb
#include <iostream>
#include <fstream>
#include <cstdlib>
#include <time.h>

using namespace std;
ifstream in("radixsort.in");
ofstream out("radixsort.out");

#define ll long long
#define pb push_back
const int NMax = 1e7 + 5;
const int inf = 1e9 + 5;
const int Radix = 0b11111111;
const int nrBytes = 4;
const int bits = 8;

ll N,A,B,C;
int v1[NMax],v2[NMax];

void countingSort(int*,int*,int);

int main() {
    in>>N>>A>>B>>C;

    //cout<<Radix<<'\n';

    v1[1] = B;
    A %= C;
    B %= C;
    for (int i=2;i <= N;++i) {
        v1[i] = (A * v1[i-1] + B) % C;
    }

    for (int i=0;i < nrBytes;++i) {
        if (i&1) {
            countingSort(v2,v1,i);
        }
        else {
            countingSort(v1,v2,i);
        }
    }

    for (int i=1;i <= N;i+=10) {
        out<<v1[i]<<' ';
    }
    out<<'\n';

    in.close();out.close();
    return 0;
}

void countingSort(int* one,int* two,int byte) {
    int occ[(1<<bits) + 5] = {};

    for (int i=1;i <= N;++i) {
        ++occ[(one[i]>>(bits*byte)) & Radix];
    }

    int total = 0;
    for (int i=0;i <= Radix;++i) {
        int old = occ[i];
        occ[i] = total + 1;
        total += old;
    }

    for (int i=1;i <= N;++i) {
        two[occ[(one[i]>>(bits*byte)) & Radix]++] = one[i];
    }
}