#define H_PRIME   1000003
#define HASH_LEN  11

// my_m128_t mimics Intel __m128i type
typedef struct my_m128_t my_m128_t;
struct my_m128_t { short m128i_i16[8]; };

int inner(
   int startval, int endval,
   int m2, int d2, int c2, int l2, int x2, int v2, int i2,
   my_m128_t* ps)
{
  int m3, m4, m5, m6, m7, m9;
  int d3, d4, d5, d6, d7, d8, d9;
  int c3, c4, c5, c6, c7, c9;
  int l3, l4, l5, l6, l7, l9;
  int x3, x4, x5, x6, x7, x9;
  int q3, q4, q5, q6, q7, q8, q9;
  int dist;
  int iret = 0;

  for (q3 = startval; q3 < endval; ++q3) {
    if (q3 == 10 || q3 == 13) continue;
    m3 = (m2 ^ q3) * H_PRIME;
    d3 = (d2 ^ q3) * H_PRIME;
    c3 = (c2 ^ q3) * H_PRIME;
    l3 = (l2 ^ q3) * H_PRIME;
    x3 = (x2 ^ q3) * H_PRIME;
    for (q4 = 1; q4 < 128; ++q4) {
      if (q4 == 10 || q4 == 13) continue;
      m4 = (m3 ^ q4) * H_PRIME;
      d4 = (d3 ^ q4) * H_PRIME;
      c4 = (c3 ^ q4) * H_PRIME;
      l4 = (l3 ^ q4) * H_PRIME;
      x4 = (x3 ^ q4) * H_PRIME;
      for (q5 = 1; q5 < 128; ++q5) {
        if (q5 == 10 || q5 == 13) continue;
        m5 = (m4 ^ q5) * H_PRIME;
        d5 = (d4 ^ q5) * H_PRIME;
        c5 = (c4 ^ q5) * H_PRIME;
        l5 = (l4 ^ q5) * H_PRIME;
        x5 = (x4 ^ q5) * H_PRIME;
        for (q6 = 1; q6 < 128; ++q6) {
          if (q6 == 10 || q6 == 13) continue;
          m6 = (m5 ^ q6) * H_PRIME;
          d6 = (d5 ^ q6) * H_PRIME;
          c6 = (c5 ^ q6) * H_PRIME;
          l6 = (l5 ^ q6) * H_PRIME;
          x6 = (x5 ^ q6) * H_PRIME;
          for (q7 = 1; q7 < 128; ++q7) {
            if (q7 == 10 || q7 == 13) continue;
            m7 = (m6 ^ q7) * H_PRIME;
            d7 = (d6 ^ q7) * H_PRIME;
            c7 = (c6 ^ q7) * H_PRIME;
            l7 = (l6 ^ q7) * H_PRIME;
            x7 = (x6 ^ q7) * H_PRIME;
            for (q8 = 1; q8 < 128; ++q8) {
              if (q8 == 10 || q8 == 13) continue;

              d8 = (d7 ^ q8) * H_PRIME ^ HASH_LEN;
              d9 = d8 & ~127;
              d9 %= 1001; if (d9 < 0) d9 += 1001;
              if (d9 < 500 - 127) continue;
              if (d9 > 500 + 127) continue;
              dist = d9 < 500 ? 500 - d9 : d9 - 500;
              q9 = (d8 & 127) ^ dist;
              if (q9 == 0 || q9 == 10 || q9 == 13) continue;

              c9 = (c7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
              c9 %= 1001; if (c9 < 0) c9 += 1001;
              if (c9 != 100) continue;

              l9 = (l7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
              l9 %= 1001; if (l9 < 0) l9 += 1001;
              if (l9 != 50) continue;

              x9 = (x7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
              x9 %= 1001; if (x9 < 0) x9 += 1001;
              if (x9 != 10) continue;

              m9 = (m7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
              if (m9 == -1) m9 = -2;
              m9 %= 1001; if (m9 < 0) m9 += 1001;
              if (m9 != 1000) continue;

              ps[iret].m128i_i16[0] = q3;
              ps[iret].m128i_i16[1] = q4;
              ps[iret].m128i_i16[2] = q5;
              ps[iret].m128i_i16[3] = q6;
              ps[iret].m128i_i16[4] = q7;
              ps[iret].m128i_i16[5] = q8;
              ps[iret].m128i_i16[6] = q9;
              ++iret;
            }
          }
        }
      }
    }
  }
  return iret;
}

##</code><code>##

gcc -O3 -c inner2.c
g++ -o find1 find1.cpp inner2.o

##</code><code>##

cl /W3 /MD /O2 find1.cpp inner2.c

##</code><code>##

# find1 9  8 9  60 61  1 2
9: sv1=8 ev1=9 sv2=60 ev2=61 sv3=1 ev3=2:
9 8 60 1 94 6 51 38 78 100
(wall clock time:109 secs, cpu time:108.95 units)

##</code><code>##

109 * 125**4 = 26,611,328,125 seconds = 26611328125 / 60 / 60 / 24 / 365 = 844 years

##</code><code>##

m4 = (m3 ^ q4) * H_PRIME;
d4 = (d3 ^ q4) * H_PRIME;
c4 = (c3 ^ q4) * H_PRIME;
l4 = (l3 ^ q4) * H_PRIME;
x4 = (x3 ^ q4) * H_PRIME;

##</code><code>##

#include <emmintrin.h>
#include <smmintrin.h>

#define H_PRIME   1000003
#define HASH_LEN  11

// my_m128_t mimics Intel __m128i type
typedef struct my_m128_t my_m128_t;
struct my_m128_t { short m128i_i16[8]; };

int inner(
   int startval, int endval,
   int m2, int d2, int c2, int l2, int x2, int v2, int i2,
   my_m128_t* ps)
{
  __m128i s2 = _mm_set_epi32(x2, l2, c2, d2);
  __m128i hp = _mm_set1_epi32(H_PRIME);
  __m128i s3, s4, s5, s6, s7;

  int m3, m4, m5, m6, m7, m9;
  int d8, d9;
  int c9;
  int l9;
  int x9;
  int q3, q4, q5, q6, q7, q8, q9;
  int dist;
  int iret = 0;

  for (q3 = startval; q3 < endval; ++q3) {
    if (q3 == 10 || q3 == 13) continue;
    m3 = (m2 ^ q3) * H_PRIME;
    s3 = _mm_mullo_epi32(_mm_xor_si128(s2, _mm_set1_epi32(q3)), hp);
    for (q4 = 1; q4 < 128; ++q4) {
      if (q4 == 10 || q4 == 13) continue;
      m4 = (m3 ^ q4) * H_PRIME;
      s4 = _mm_mullo_epi32(_mm_xor_si128(s3, _mm_set1_epi32(q4)), hp);
      for (q5 = 1; q5 < 128; ++q5) {
        if (q5 == 10 || q5 == 13) continue;
        m5 = (m4 ^ q5) * H_PRIME;
        s5 = _mm_mullo_epi32(_mm_xor_si128(s4, _mm_set1_epi32(q5)), hp);
        for (q6 = 1; q6 < 128; ++q6) {
          if (q6 == 10 || q6 == 13) continue;
          m6 = (m5 ^ q6) * H_PRIME;
          s6 = _mm_mullo_epi32(_mm_xor_si128(s5, _mm_set1_epi32(q6)), hp);
          for (q7 = 1; q7 < 128; ++q7) {
            if (q7 == 10 || q7 == 13) continue;
            m7 = (m6 ^ q7) * H_PRIME;
            s7 = _mm_mullo_epi32(_mm_xor_si128(s6, _mm_set1_epi32(q7)), hp);
            for (q8 = 1; q8 < 128; ++q8) {
              if (q8 == 10 || q8 == 13) continue;

              d8 = (s7.m128i_i32[0] ^ q8) * H_PRIME ^ HASH_LEN;
              d9 = d8 & ~127;
              d9 %= 1001; if (d9 < 0) d9 += 1001;
              if (d9 < 500 - 127) continue;
              if (d9 > 500 + 127) continue;
              dist = d9 < 500 ? 500 - d9 : d9 - 500;
              q9 = (d8 & 127) ^ dist;
              if (q9 == 0 || q9 == 10 || q9 == 13) continue;

              c9 = (s7.m128i_i32[1] ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
              c9 %= 1001; if (c9 < 0) c9 += 1001;
              if (c9 != 100) continue;

              l9 = (s7.m128i_i32[2] ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
              l9 %= 1001; if (l9 < 0) l9 += 1001;
              if (l9 != 50) continue;

              x9 = (s7.m128i_i32[3] ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
              x9 %= 1001; if (x9 < 0) x9 += 1001;
              if (x9 != 10) continue;

              m9 = (m7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
              if (m9 == -1) m9 = -2;
              m9 %= 1001; if (m9 < 0) m9 += 1001;
              if (m9 != 1000) continue;

              ps[iret].m128i_i16[0] = q3;
              ps[iret].m128i_i16[1] = q4;
              ps[iret].m128i_i16[2] = q5;
              ps[iret].m128i_i16[3] = q6;
              ps[iret].m128i_i16[4] = q7;
              ps[iret].m128i_i16[5] = q8;
              ps[iret].m128i_i16[6] = q9;
              ++iret;
            }
          }
        }
      }
    }
  }
  return iret;
}

##</code><code>##

d8 = (d7 ^ q8) * H_PRIME ^ HASH_LEN;
d9 = d8 & ~127;
d9 %= 1001; if (d9 < 0) d9 += 1001;
if (d9 < 500 - 127) continue;
if (d9 > 500 + 127) continue;
dist = d9 < 500 ? 500 - d9 : d9 - 500;
q9 = (d8 & 127) ^ dist;
if (q9 == 0 || q9 == 10 || q9 == 13) continue;

c9 = (c7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
c9 %= 1001; if (c9 < 0) c9 += 1001;
if (c9 != 100) continue;

l9 = (l7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
l9 %= 1001; if (l9 < 0) l9 += 1001;
if (l9 != 50) continue;

x9 = (x7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
x9 %= 1001; if (x9 < 0) x9 += 1001;
if (x9 != 10) continue;

m9 = (m7 ^ q8) * H_PRIME ^ HASH_LEN ^ q9;
if (m9 == -1) m9 = -2;
m9 %= 1001; if (m9 < 0) m9 += 1001;
if (m9 != 1000) continue;

ps[iret].m128i_i16[0] = q3;
ps[iret].m128i_i16[1] = q4;
ps[iret].m128i_i16[2] = q5;
ps[iret].m128i_i16[3] = q6;
ps[iret].m128i_i16[4] = q7;
ps[iret].m128i_i16[5] = q8;
ps[iret].m128i_i16[6] = q9;
++iret;

##</code><code>##

q9 = bytevecM[(unsigned int)(m7 ^ q8)];
if (q9 != 0
&&  q9 == bytevecD[(unsigned int)(d7 ^ q8)]
&&  q9 == bytevecC[(unsigned int)(c7 ^ q8)]
&&  q9 == bytevecL[(unsigned int)(l7 ^ q8)]
&&  q9 == bytevecX[(unsigned int)(x7 ^ q8)]) {
   ps[iret].m128i_i16[0] = q3;
   ps[iret].m128i_i16[1] = q4;
   ps[iret].m128i_i16[2] = q5;
   ps[iret].m128i_i16[3] = q6;
   ps[iret].m128i_i16[4] = q7;
   ps[iret].m128i_i16[5] = q8;
   ps[iret].m128i_i16[6] = q9;
   ++iret;
}

##</code><code>##

# find3a 9  8 9  60 61  1 2
9: sv1=8 ev1=9 sv2=60 ev2=61 sv3=1 ev3=2:
9 8 60 1 94 6 51 38 78 100
(wall clock time:94 secs, cpu time:93.90 units)