#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <time.h>

#define H_PRIME 1000003

// Most unlikely to get this many hits in a single call of the inner loop function.
#define MAX_HITS 64

// my_m128_t mimics Intel __m128i type (which will be used later)
typedef struct my_m128_t my_m128_t;
struct my_m128_t { short m128i_i16[8]; };

// the inner loop function (lives in a separate file)
extern "C" int inner(
   int startval, int endval,
   int m2, int d2, int c2, int l2, int x2, int v2, int i2,
   my_m128_t* ps
);

void do_one(int q0, int sv1, int ev1, int sv2, int ev2, int sv3, int ev3)
{
  time_t  tstart = time(NULL);
  clock_t cstart = clock();
  time_t  tend;
  clock_t cend;
  my_m128_t soln[MAX_HITS];    
  int m1, m2;
  int m_char = 'M';
  int m0 = ( ( H_PRIME * (m_char << 7) ) ^ m_char ) * H_PRIME;
  int d1, d2;
  int d_char = 'D';
  int d0 = ( ( H_PRIME * (d_char << 7) ) ^ d_char ) * H_PRIME;
  int c1, c2;
  int c_char = 'C';
  int c0 = ( ( H_PRIME * (c_char << 7) ) ^ c_char ) * H_PRIME;
  int l1, l2;
  int l_char = 'L';
  int l0 = ( ( H_PRIME * (l_char << 7) ) ^ l_char ) * H_PRIME;
  int x1, x2;
  int x_char = 'X';
  int x0 = ( ( H_PRIME * (x_char << 7) ) ^ x_char ) * H_PRIME;
  int v1, v2;
  int v_char = 'V';
  int v0 = ( ( H_PRIME * (v_char << 7) ) ^ v_char ) * H_PRIME;
  int i1, i2;
  int i_char = 'I';
  int i0 = ( ( H_PRIME * (i_char << 7) ) ^ i_char ) * H_PRIME;
  int q1, q2;
  int mm0 = (m0 ^ q0) * H_PRIME;
  int dd0 = (d0 ^ q0) * H_PRIME;
  int cc0 = (c0 ^ q0) * H_PRIME;
  int ll0 = (l0 ^ q0) * H_PRIME;
  int xx0 = (x0 ^ q0) * H_PRIME;
  int vv0 = (v0 ^ q0) * H_PRIME;
  int ii0 = (i0 ^ q0) * H_PRIME;
  fprintf(stderr, "%d: sv1=%d ev1=%d sv2=%d ev2=%d sv3=%d ev3=%d:\n",
     q0, sv1, ev1, sv2, ev2, sv3, ev3);
  // We ignore 0, 10, 13 because these three require escaping in Python.
  for (q1 = sv1; q1 < ev1; ++q1) {
    if (q1 == 10 || q1 == 13) continue;
    m1 = (mm0 ^ q1) * H_PRIME;
    d1 = (dd0 ^ q1) * H_PRIME;
    c1 = (cc0 ^ q1) * H_PRIME;
    l1 = (ll0 ^ q1) * H_PRIME;
    x1 = (xx0 ^ q1) * H_PRIME;
    v1 = (vv0 ^ q1) * H_PRIME;
    i1 = (ii0 ^ q1) * H_PRIME;
    for (q2 = sv2; q2 < ev2; ++q2) {
      if (q2 == 10 || q2 == 13) continue;
      m2 = (m1 ^ q2) * H_PRIME;
      d2 = (d1 ^ q2) * H_PRIME;
      c2 = (c1 ^ q2) * H_PRIME;
      l2 = (l1 ^ q2) * H_PRIME;
      x2 = (x1 ^ q2) * H_PRIME;
      v2 = (v1 ^ q2) * H_PRIME;
      i2 = (i1 ^ q2) * H_PRIME;
      int isoln = inner(sv3, ev3, m2, d2, c2, l2, x2, v2, i2, soln);
      if (isoln > 0) {
        int k;
        for (k = 0; k < isoln; ++k) {
          fprintf(stderr, "%d %d %d %hd %hd %hd %hd %hd %hd %hd\n",
          q0, q1, q2,           soln[k].m128i_i16[0],
          soln[k].m128i_i16[1], soln[k].m128i_i16[2], soln[k].m128i_i16[3],
          soln[k].m128i_i16[4], soln[k].m128i_i16[5], soln[k].m128i_i16[6]);
        }
      }
    }
  }
  tend = time(NULL);
  cend = clock();
  fprintf(stderr,  "(wall clock time:%ld secs, cpu time:%.2f units)\n",
    (long)   (difftime(tend, tstart)+0.5),
    (double) (cend-cstart) / (double)CLOCKS_PER_SEC);
}

int main(int argc, char* argv[])
{
   int sv0, sv1, sv2, sv3, ev1, ev2, ev3;
   if (argc != 8) { fprintf(stderr, "usage: prog sv0 sv1 ev1 sv2 ev2 sv3 ev3\n"); return 1; }
   sv0 = atoi(argv[1]);
   sv1 = atoi(argv[2]);
   ev1 = atoi(argv[3]);
   sv2 = atoi(argv[4]);
   ev2 = atoi(argv[5]);
   sv3 = atoi(argv[6]);
   ev3 = atoi(argv[7]);
   do_one(sv0, sv1, ev1, sv2, ev2, sv3, ev3);
   return 0;
}