src/blocksorting/bs.c

   1
   2 #include "bs.h"
   3 #include <stdlib.h>
   4
   5 /* Block Sorting Optimizado en memoria! */
   6
   7 typedef struct _bs_decode_t_ {
   8         char c;
   9         Uint32 sig;
  10 } t_BlockSortDecode;
  11
  12 char es_menor(char *data, t_BlockSort *bs, int i, int j);
  13
  14 int _compare(const void *d1, const void *d2) {
  15         t_BlockSortDecode *s1, *s2;
  16
  17         s1 = (t_BlockSortDecode *)d1;
  18         s2 = (t_BlockSortDecode *)d2;
  19
  20         /* Hack si es el mismo caracter, dejo como mejor la de menor sig */
  21         if (s1->c == s2->c) {
  22                 return (s1->sig - s2->sig);
  23         }
  24
  25         return (s1->c - s2->c);
  26 }
  27
  28 int __compare(const void *d1, const void *d2) {
  29         t_BlockSortData *s1, *s2;
  30         char i;
  31
  32         s1 = (t_BlockSortData *)d1;
  33         s2 = (t_BlockSortData *)d2;
  34
  35         i = es_menor(s1->bs->data, s1->bs, s1->pos_inicial, s2->pos_inicial);
  36
  37 /*      if (i == 0) {*/
  38                 /* si ambos strings son iguales, fuerzo a que d1 quede primero */
  39                 /*return -1;
  40         }*/
  41         return i;
  42 }
  43
  44 char es_menor(char *data, t_BlockSort *bs, int i, int j)
  45 {
  46         Uint32 pi, pj, k;
  47
  48         for(k=0; k<bs->len; k++) {
  49                 pi = (i+k)%bs->len;
  50                 pj = (j+k)%bs->len;
  51
  52                 if (data[pi] > data[pj]) return 1; /* El primero es mayor */
  53                 if (data[pi] < data[pj]) return -1; /* El primero es menor */
  54         }
  55
  56         /* Son iguales! */
  57         /* Hack */
  58         return 0;
  59 }
  60
  61 void generar_array(char *data, t_BlockSort *bs)
  62 {
  63         int i;
  64         for(i=0; i<bs->len; i++) {
  65                 bs->array[i].pos_inicial = i;
  66                 bs->array[i].pos_final = (i+bs->len-1)%bs->len;
  67                 bs->array[i].ord = (i==0)?1:0;
  68                 bs->array[i].bs = bs;
  69         }
  70 }
  71
  72 void ordenar_array(char *data, t_BlockSort *bs)
  73 {
  74         qsort(bs->array, bs->len, sizeof(t_BlockSortData), __compare);
  75 }
  76
  77 void print_(char *data, Uint32 pos, Uint32 len)
  78 {
  79         Uint32 i;
  80
  81         for(i=0; i<len; i++)
  82                 fprintf(stderr, "%c", data[(pos+i)%len]);
  83         fprintf(stderr, "\n");
  84 }
  85
  86 int generar_salida(char *data, t_BlockSort *bs, char *salida)
  87 {
  88         Uint32 i, k;
  89         char *out;
  90
  91         /* Dejo lugar para guardar el k y el tamaño de este bloque */
  92         out = salida + sizeof(Uint32)*2;
  93
  94         k=-1;
  95         for(i=0; i<bs->len; i++) {
  96                 /* print_(data, bs->array[i].pos_inicial, bs->len); */
  97                 out[i] = data[bs->array[i].pos_final];
  98                 if (bs->array[i].ord == 1) k = i;
  99         }
 100         return k;
 101 }
 102
 103 void bs_solve(char *in, char *out, t_BlockSort *bs, Uint32 *k, Uint32 leido)
 104 {
 105         unsigned int l;
 106         l = bs->len;
 107         /* Hack para pedasos menores a la pagina */
 108         if (leido < bs->len) bs->len = leido;
 109         bs->data = in;
 110
 111         generar_array(in, bs);
 112         ordenar_array(in, bs);
 113         (*k) = generar_salida(in, bs, out);
 114
 115         /* Guardo el k y el tamaño en el array */
 116         memcpy(out, &leido, sizeof(Uint32));
 117         memcpy(out+sizeof(Uint32), k, sizeof(Uint32));
 118         bs->len = l;
 119 }
 120
 121 void bs_restore(char *dst, char *c, Uint32 k, Uint32 len)
 122 {
 123         Uint32 i, current;
 124         t_BlockSortDecode *in;
 125
 126         in = malloc(sizeof(t_BlockSortDecode)*len);
 127
 128         for(i=0; i<len; i++) {
 129                 in[i].c = c[i];
 130                 in[i].sig = i;
 131         }
 132
 133         qsort(in, len, sizeof(t_BlockSortDecode), _compare);
 134
 135         current = k;
 136         i=0;
 137         do {
 138                 dst[i++] = in[current].c;
 139                 current = in[current].sig;
 140         } while (current != k);
 141         free(in);
 142
 143         if (i < len) printf("Noo %ld != %ld\n", i, len);
 144 }
 145
 146 t_BlockSort *bs_create(Uint32 len)
 147 {
 148         t_BlockSort *tmp;
 149
 150         tmp = malloc(sizeof(t_BlockSort));
 151         tmp->array = malloc(sizeof(t_BlockSortData)*len);
 152         tmp->len = len;
 153
 154         return tmp;
 155 }
 156
 157
 158 void bs_destroy(t_BlockSort *bs)
 159 {
 160         free(bs->array);
 161         free(bs);
 162 }
 163