fork(1) download
copy 
  1. #include <iostream>
  2. #include <numeric>
  3. #include <cstdlib> 
  4. #include <cstddef>
  5. #include <stdint.h>
  6.  
  7. enum {
  8.   // Config
  9.   ZERO_POS = 2,
  10.   INDEX_SIZE = 5,
  11.   VALUE_LIMIT = -100,
  12.  
  13.   // Constants
  14.   NUM_ROWS = 37,
  15.   NUM_COLUMNS = 73,
  16.   VALUE_LIMIT_POS = -VALUE_LIMIT,
  17.   NUM_BIN_OUTLIERS = 8,
  18.   NUM_VALUES = NUM_ROWS * (NUM_COLUMNS - 1),
  19.   INDEX_PLUS = INDEX_SIZE - ZERO_POS - 1,
  20.   BIN4 = 0,
  21.   BIN12 = 1,
  22.   INDEX_FULL_SIZE = INDEX_SIZE + 2
  23. };
  24.  
  25. int16_t my_array[NUM_ROWS][NUM_COLUMNS] = {
  26. {150,145,140,135,130,125,120,115,110,105,100,95,90,85,80,75,70,65,60,55,50,45,40,35,30,25,20,15,10,5,0,-4,-9,-14,-19,-24,-29,-34,-39,-44,-49,-54,-59,-64,-69,-74,-79,-84,-89,-94,-99,104,109,114,119,124,129,134,139,144,149,154,159,164,169,174,179,175,170,165,160,155,150}, \
  27. {143,137,131,126,120,115,110,105,100,95,90,85,80,75,71,66,62,57,53,48,44,39,35,31,27,22,18,14,9,5,1,-3,-7,-11,-16,-20,-25,-29,-34,-38,-43,-47,-52,-57,-61,-66,-71,-76,-81,-86,-91,-96,101,107,112,117,123,128,134,140,146,151,157,163,169,175,178,172,166,160,154,148,143}, \
  28. {130,124,118,112,107,101,96,92,87,82,78,74,70,65,61,57,54,50,46,42,38,34,31,27,23,19,16,12,8,4,1,-2,-6,-10,-14,-18,-22,-26,-30,-34,-38,-43,-47,-51,-56,-61,-65,-70,-75,-79,-84,-89,-94,100,105,111,116,122,128,135,141,148,155,162,170,177,174,166,159,151,144,137,130}, \
  29. {111,104,99,94,89,85,81,77,73,70,66,63,60,56,53,50,46,43,40,36,33,30,26,23,20,16,13,10,6,3,0,-3,-6,-9,-13,-16,-20,-24,-28,-32,-36,-40,-44,-48,-52,-57,-61,-65,-70,-74,-79,-84,-88,-93,-98,103,109,115,121,128,135,143,152,162,172,176,165,154,144,134,125,118,111}, \
  30. {85,81,77,74,71,68,65,63,60,58,56,53,51,49,46,43,41,38,35,32,29,26,23,19,16,13,10,7,4,1,-1,-3,-6,-9,-13,-16,-19,-23,-26,-30,-34,-38,-42,-46,-50,-54,-58,-62,-66,-70,-74,-78,-83,-87,-91,-95,100,105,110,117,124,133,144,159,178,160,141,125,112,103,96,90,85}, \
  31. {62,60,58,57,55,54,52,51,50,48,47,46,44,42,41,39,36,34,31,28,25,22,19,16,13,10,7,4,2,0,-3,-5,-8,-10,-13,-16,-19,-22,-26,-29,-33,-37,-41,-45,-49,-53,-56,-60,-64,-67,-70,-74,-77,-80,-83,-86,-89,-91,-94,-97,101,105,111,130,109,84,77,74,71,68,66,64,62}, \
  32. {46,46,45,44,44,43,42,42,41,41,40,39,38,37,36,35,33,31,28,26,23,20,16,13,10,7,4,1,-1,-3,-5,-7,-9,-12,-14,-16,-19,-22,-26,-29,-33,-36,-40,-44,-48,-51,-55,-58,-61,-64,-66,-68,-71,-72,-74,-74,-75,-74,-72,-68,-61,-48,-25,2,22,33,40,43,45,46,47,46,46}, \
  33. {36,36,36,36,36,35,35,35,35,34,34,34,34,33,32,31,30,28,26,23,20,17,14,10,6,3,0,-2,-4,-7,-9,-10,-12,-14,-15,-17,-20,-23,-26,-29,-32,-36,-40,-43,-47,-50,-53,-56,-58,-60,-62,-63,-64,-64,-63,-62,-59,-55,-49,-41,-30,-17,-4,6,15,22,27,31,33,34,35,36,36}, \
  34. {30,30,30,30,30,30,30,29,29,29,29,29,29,29,29,28,27,26,24,21,18,15,11,7,3,0,-3,-6,-9,-11,-12,-14,-15,-16,-17,-19,-21,-23,-26,-29,-32,-35,-39,-42,-45,-48,-51,-53,-55,-56,-57,-57,-56,-55,-53,-49,-44,-38,-31,-23,-14,-6,0,7,13,17,21,24,26,27,29,29,30}, \
  35. {25,25,26,26,26,25,25,25,25,25,25,25,25,26,25,25,24,23,21,19,16,12,8,4,0,-3,-7,-10,-13,-15,-16,-17,-18,-19,-20,-21,-22,-23,-25,-28,-31,-34,-37,-40,-43,-46,-48,-49,-50,-51,-51,-50,-48,-45,-42,-37,-32,-26,-19,-13,-7,-1,3,7,11,14,17,19,21,23,24,25,25}, \
  36. {21,22,22,22,22,22,22,22,22,22,22,22,22,22,22,22,21,20,18,16,13,9,5,1,-3,-7,-11,-14,-17,-18,-20,-21,-21,-22,-22,-22,-23,-23,-25,-27,-29,-32,-35,-37,-40,-42,-44,-45,-45,-45,-44,-42,-40,-36,-32,-27,-22,-17,-12,-7,-3,0,3,7,9,12,14,16,18,19,20,21,21}, \
  37. {18,19,19,19,19,19,19,19,19,19,19,19,19,19,19,19,18,17,16,14,10,7,2,-1,-6,-10,-14,-17,-19,-21,-22,-23,-24,-24,-24,-24,-23,-23,-23,-24,-26,-28,-30,-33,-35,-37,-38,-39,-39,-38,-36,-34,-31,-28,-24,-19,-15,-10,-6,-3,0,1,4,6,8,10,12,14,15,16,17,18,18}, \
  38. {16,16,17,17,17,17,17,17,17,17,17,16,16,16,16,16,16,15,13,11,8,4,0,-4,-9,-13,-16,-19,-21,-23,-24,-25,-25,-25,-25,-24,-23,-21,-20,-20,-21,-22,-24,-26,-28,-30,-31,-32,-31,-30,-29,-27,-24,-21,-17,-13,-9,-6,-3,-1,0,2,4,5,7,9,10,12,13,14,15,16,16}, \
  39. {14,14,14,15,15,15,15,15,15,15,14,14,14,14,14,14,13,12,11,9,5,2,-2,-6,-11,-15,-18,-21,-23,-24,-25,-25,-25,-25,-24,-22,-21,-18,-16,-15,-15,-15,-17,-19,-21,-22,-24,-24,-24,-23,-22,-20,-18,-15,-12,-9,-5,-3,-1,0,1,2,4,5,6,8,9,10,11,12,13,14,14}, \
  40. {12,13,13,13,13,13,13,13,13,13,13,13,12,12,12,12,11,10,9,6,3,0,-4,-8,-12,-16,-19,-21,-23,-24,-24,-24,-24,-23,-22,-20,-17,-15,-12,-10,-9,-9,-10,-12,-13,-15,-17,-17,-18,-17,-16,-15,-13,-11,-8,-5,-3,-1,0,1,1,2,3,4,6,7,8,9,10,11,12,12,12}, \
  41. {11,11,11,11,11,12,12,12,12,12,11,11,11,11,11,10,10,9,7,5,2,-1,-5,-9,-13,-17,-20,-22,-23,-23,-23,-23,-22,-20,-18,-16,-14,-11,-9,-6,-5,-4,-5,-6,-8,-9,-11,-12,-12,-12,-12,-11,-9,-8,-6,-3,-1,0,0,1,1,2,3,4,5,6,7,8,9,10,11,11,11}, \
  42. {10,10,10,10,10,10,10,10,10,10,10,10,10,10,9,9,9,7,6,3,0,-3,-6,-10,-14,-17,-20,-21,-22,-22,-22,-21,-19,-17,-15,-13,-10,-8,-6,-4,-2,-2,-2,-2,-4,-5,-7,-8,-8,-9,-8,-8,-7,-5,-4,-2,0,0,1,1,1,2,2,3,4,5,6,7,8,9,10,10,10}, \
  43. {9,9,9,9,9,9,9,10,10,9,9,9,9,9,9,8,8,6,5,2,0,-4,-7,-11,-15,-17,-19,-21,-21,-21,-20,-18,-16,-14,-12,-10,-8,-6,-4,-2,-1,0,0,0,-1,-2,-4,-5,-5,-6,-6,-5,-5,-4,-3,-1,0,0,1,1,1,1,2,3,3,5,6,7,8,8,9,9,9}, \
  44. {9,9,9,9,9,9,9,9,9,9,9,9,8,8,8,8,7,5,4,1,-1,-5,-8,-12,-15,-17,-19,-20,-20,-19,-18,-16,-14,-11,-9,-7,-5,-4,-2,-1,0,0,1,1,0,0,-2,-3,-3,-4,-4,-4,-3,-3,-2,-1,0,0,0,0,0,1,1,2,3,4,5,6,7,8,8,9,9}, \
  45. {9,9,9,8,8,8,9,9,9,9,9,8,8,8,8,7,6,5,3,0,-2,-5,-9,-12,-15,-17,-18,-19,-19,-18,-16,-14,-12,-9,-7,-5,-4,-2,-1,0,0,1,1,1,1,0,0,-1,-2,-2,-3,-3,-2,-2,-1,-1,0,0,0,0,0,0,0,1,2,3,4,5,6,7,8,8,9}, \
  46. {8,8,8,8,8,8,9,9,9,9,9,9,8,8,8,7,6,4,2,0,-3,-6,-9,-12,-15,-17,-18,-18,-17,-16,-14,-12,-10,-8,-6,-4,-2,-1,0,0,1,2,2,2,2,1,0,0,-1,-1,-1,-2,-2,-1,-1,0,0,0,0,0,0,0,0,0,1,2,3,4,5,6,7,8,8}, \
  47. {8,8,8,8,9,9,9,9,9,9,9,9,9,8,8,7,5,3,1,-1,-4,-7,-10,-13,-15,-16,-17,-17,-16,-15,-13,-11,-9,-6,-5,-3,-2,0,0,0,1,2,2,2,2,1,1,0,0,0,-1,-1,-1,-1,-1,0,0,0,0,-1,-1,-1,-1,-1,0,0,1,3,4,5,7,7,8}, \
  48. {8,8,9,9,9,9,10,10,10,10,10,10,10,9,8,7,5,3,0,-2,-5,-8,-11,-13,-15,-16,-16,-16,-15,-13,-12,-10,-8,-6,-4,-2,-1,0,0,1,2,2,3,3,2,2,1,0,0,0,0,0,0,0,0,0,0,-1,-1,-2,-2,-2,-2,-2,-1,0,0,1,3,4,6,7,8}, \
  49. {7,8,9,9,9,10,10,11,11,11,11,11,10,10,9,7,5,3,0,-2,-6,-9,-11,-13,-15,-16,-16,-15,-14,-13,-11,-9,-7,-5,-3,-2,0,0,1,1,2,3,3,3,3,2,2,1,1,0,0,0,0,0,0,0,-1,-1,-2,-3,-3,-4,-4,-4,-3,-2,-1,0,1,3,5,6,7}, \
  50. {6,8,9,9,10,11,11,12,12,12,12,12,11,11,9,7,5,2,0,-3,-7,-10,-12,-14,-15,-16,-15,-15,-13,-12,-10,-8,-7,-5,-3,-1,0,0,1,2,2,3,3,4,3,3,3,2,2,1,1,1,0,0,0,0,-1,-2,-3,-4,-4,-5,-5,-5,-5,-4,-2,-1,0,2,3,5,6}, \
  51. {6,7,8,10,11,12,12,13,13,14,14,13,13,11,10,8,5,2,0,-4,-8,-11,-13,-15,-16,-16,-16,-15,-13,-12,-10,-8,-6,-5,-3,-1,0,0,1,2,3,3,4,4,4,4,4,3,3,3,2,2,1,1,0,0,-1,-2,-3,-5,-6,-7,-7,-7,-6,-5,-4,-3,-1,0,2,4,6}, \
  52. {5,7,8,10,11,12,13,14,15,15,15,14,14,12,11,8,5,2,-1,-5,-9,-12,-14,-16,-17,-17,-16,-15,-14,-12,-11,-9,-7,-5,-3,-1,0,0,1,2,3,4,4,5,5,5,5,5,5,4,4,3,3,2,1,0,-1,-2,-4,-6,-7,-8,-8,-8,-8,-7,-6,-4,-2,0,1,3,5}, \
  53. {4,6,8,10,12,13,14,15,16,16,16,16,15,13,11,9,5,2,-2,-6,-10,-13,-16,-17,-18,-18,-17,-16,-15,-13,-11,-9,-7,-5,-4,-2,0,0,1,3,3,4,5,6,6,7,7,7,7,7,6,5,4,3,2,0,-1,-3,-5,-7,-8,-9,-10,-10,-10,-9,-7,-5,-4,-1,0,2,4}, \
  54. {4,6,8,10,12,14,15,16,17,18,18,17,16,15,12,9,5,1,-3,-8,-12,-15,-18,-19,-20,-20,-19,-18,-16,-15,-13,-11,-8,-6,-4,-2,-1,0,1,3,4,5,6,7,8,9,9,9,9,9,9,8,7,5,3,1,-1,-3,-6,-8,-10,-11,-12,-12,-11,-10,-9,-7,-5,-2,0,1,4}, \
  55. {4,6,8,11,13,15,16,18,19,19,19,19,18,16,13,10,5,0,-5,-10,-15,-18,-21,-22,-23,-22,-22,-20,-18,-17,-14,-12,-10,-8,-5,-3,-1,0,1,3,5,6,8,9,10,11,12,12,13,12,12,11,9,7,5,2,0,-3,-6,-9,-11,-12,-13,-13,-12,-11,-10,-8,-6,-3,-1,1,4}, \
  56. {3,6,9,11,14,16,17,19,20,21,21,21,19,17,14,10,4,-1,-8,-14,-19,-22,-25,-26,-26,-26,-25,-23,-21,-19,-17,-14,-12,-9,-7,-4,-2,0,1,3,5,7,9,11,13,14,15,16,16,16,16,15,13,10,7,4,0,-3,-7,-10,-12,-14,-15,-14,-14,-12,-11,-9,-6,-4,-1,1,3}, \
  57. {4,6,9,12,14,17,19,21,22,23,23,23,21,19,15,9,2,-5,-13,-20,-25,-28,-30,-31,-31,-30,-29,-27,-25,-22,-20,-17,-14,-11,-9,-6,-3,0,1,4,6,9,11,13,15,17,19,20,21,21,21,20,18,15,11,6,2,-2,-7,-11,-13,-15,-16,-16,-15,-13,-11,-9,-7,-4,-1,1,4}, \
  58. {4,7,10,13,15,18,20,22,24,25,25,25,23,20,15,7,-2,-12,-22,-29,-34,-37,-38,-38,-37,-36,-34,-31,-29,-26,-23,-20,-17,-13,-10,-7,-4,-1,2,5,8,11,13,16,18,21,23,24,26,26,26,26,24,21,17,12,5,0,-6,-10,-14,-16,-16,-16,-15,-14,-12,-10,-7,-4,-1,1,4}, \
  59. {4,7,10,13,16,19,22,24,26,27,27,26,24,19,11,-1,-15,-28,-37,-43,-46,-47,-47,-45,-44,-41,-39,-36,-32,-29,-26,-22,-19,-15,-11,-8,-4,-1,2,5,9,12,15,19,22,24,27,29,31,33,33,33,32,30,26,21,14,6,0,-6,-11,-14,-15,-16,-15,-14,-12,-9,-7,-4,-1,1,4}, \
  60. {6,9,12,15,18,21,23,25,27,28,27,24,17,4,-14,-34,-49,-56,-60,-60,-60,-58,-56,-53,-50,-47,-43,-40,-36,-32,-28,-25,-21,-17,-13,-9,-5,-1,2,6,10,14,17,21,24,28,31,34,37,39,41,42,43,43,41,38,33,25,17,8,0,-4,-8,-10,-10,-10,-8,-7,-4,-2,0,3,6}, \
  61. {22,24,26,28,30,32,33,31,23,-18,-81,-96,-99,-98,-95,-93,-89,-86,-82,-78,-74,-70,-66,-62,-57,-53,-49,-44,-40,-36,-32,-27,-23,-19,-14,-10,-6,-1,2,6,10,15,19,23,27,31,35,38,42,45,49,52,55,57,60,61,63,63,62,61,57,53,47,40,33,28,23,21,19,19,19,20,22}, \
  62. {168,173,178,176,171,166,161,156,151,146,141,136,131,126,121,116,111,106,101,-96,-91,-86,-81,-76,-71,-66,-61,-56,-51,-46,-41,-36,-31,-26,-21,-16,-11,-6,-1,3,8,13,18,23,28,33,38,43,48,53,58,63,68,73,78,83,88,93,98,103,108,113,118,123,128,133,138,143,148,153,158,163,168}};
  63.  
  64. using namespace std;
  65.  
  66. uint16_t compacted[] = {
  67.     0, 0x4b, 0x292, 0x776, 0x9ce, 0xa18, 0xa55, 0xa68, 
  68.     0x601f, 0xea09, 0x4e8c, 0xafff, 0xbfed, 0x6ff1, 0x10fe, 0xfdb1, 
  69.     0xa100, 0xd1fe, 0xd800, 0x129f, 0xf010, 0x730f, 0x9ec1, 0x80e9, 
  70.     0x9a9c, 0x8888, 0x9880, 0x1, 0x2ab9, 0x8105, 0x101, 0xdfb8, 
  71.     0xfeaf, 0xc030, 0x8100, 0x930, 0x8c1, 0xcaf4, 0, 0xfff6, 
  72.     0xefb9, 0x416, 0x78eb, 0xe5d2, 0xeaf7, 0x597e, 0x8a50, 0xf278, 
  73.     0x30c, 0xad5a, 0x7f2a, 0x96cd, 0x860a, 0xbad7, 0x2145, 0x1b7, 
  74.     0x807f, 0xba8f, 0x8b22, 0x40c9, 0xcfe5, 0x43c, 0x56a9, 0x75e4, 
  75.     0x6b43, 0x23e, 0x3db4, 0xac48, 0x71bf, 0xeaf9, 0x4a9a, 0xe927, 
  76.     0x855, 0xe54, 0xd56f, 0xe170, 0x1c88, 0x57c, 0x99b2, 0xfd79, 
  77.     0xa8c8, 0xf5f6, 0xcd89, 0xfd66, 0x9815, 0xc126, 0x32e1, 0x56f, 
  78.     0xfe7a, 0xd86d, 0x155b, 0x4a22, 0x62d2, 0x36b8, 0x877e, 0xbd70, 
  79.     0x2746, 0xa87f, 0x169d, 0x4af0, 0xe059, 0x6cfb, 0xb27a, 0x49a6, 
  80.     0xae32, 0xd0ab, 0x14b1, 0xa498, 0xb5e9, 0xae0e, 0x9427, 0x9b4f, 
  81.     0x4957, 0xebd5, 0xf1fa, 0xf288, 0xff72, 0xc29f, 0x475e, 0xf4ca, 
  82.     0xaa0e, 0xd455, 0x1b1f, 0xeaee, 0xf2ba, 0xaee, 0x96d2, 0xa942, 
  83.     0xaeed, 0x8a5a, 0xdfda, 0x3ce7, 0x31d2, 0x60cb, 0xe98d, 0x8162, 
  84.     0xb0c, 0x27af, 0x6e6c, 0xd2e5, 0x6022, 0xe421, 0xe7f7, 0x394a, 
  85.     0x173b, 0xed47, 0xaef7, 0x4b78, 0x283b, 0x9d78, 0x3e13, 0x4a6, 
  86.     0xa6e, 0xcda6, 0x3bb0, 0x62b7, 0x2b9a, 0xee3f, 0xda27, 0x88ee, 
  87.     0x8280, 0xe658, 0x1668, 0xb701, 0x4e70, 0x147, 0xe058, 0xa1a0, 
  88.     0x3455, 0x5994, 0xd390, 0x1fac, 0xb58c, 0x3680, 0x3d42, 0xf12, 
  89.     0xba9e, 0x949f, 0x51df, 0x9603, 0x52af, 0x74d8, 0xbdc7, 0xefe9, 
  90.     0x3b5e, 0xe273, 0xbae, 0x13c4, 0x5844, 0x3c6d, 0xf045, 0x52ba, 
  91.     0x6ef2, 0xa428, 0x67b7, 0x91f7, 0x31ee, 0xbcb9, 0xfe7f, 0xefd3, 
  92.     0x367e, 0x4308, 0x2b58, 0xbd30, 0xc4bd, 0xa28f, 0xca5d, 0x2230, 
  93.     0x8f65, 0x134e, 0xaaa8, 0x6499, 0x4a14, 0x15ec, 0xf83f, 0x8c92, 
  94.     0x3928, 0x2d35, 0x176f, 0x1080, 0xc982, 0x1bae, 0xfb7a, 0x616a, 
  95.     0xc069, 0x8177, 0xe7f3, 0x70e3, 0x5d99, 0x8dfa, 0xdf80, 0x364b, 
  96.     0x6b3f, 0x1aec, 0x66c8, 0x87cc, 0xd62f, 0x687, 0x9a3d, 0x3704, 
  97.     0xe486, 0xa47e, 0xa87a, 0x1bd, 0xd8e2, 0x9f6d, 0xdd2, 0xd150, 
  98.     0xcdf3, 0xb4e4, 0x9bd7, 0xabc9, 0x52e8, 0x553c, 0xa2ab, 0x18b8, 
  99.     0x70fb, 0xa594, 0xd34c, 0x9450, 0xd928, 0xd581, 0x54a, 0x3096, 
  100.     0x30c7, 0x6899, 0xcfd0, 0x5a37, 0xa664, 0xfe0, 0x4b1a, 0x616c, 
  101.     0x9eba, 0x591a, 0xfa80, 0x2911, 0x2c5b, 0xe310, 0x32e1, 0x3680, 
  102.     0x55a8, 0x7e, 0xd7fd, 0xe81d, 0x13f2, 0x919d, 0x972f, 0x5c52, 
  103.     0xdf86, 0x4576, 0xd6ea, 0x2afd, 0x1759, 0xcf24, 0x1b0e, 0xc73f, 
  104.     0x75aa, 0x7f2, 0x7432, 0x8691, 0x389c, 0xb9a4, 0xf1ad, 0xb21a, 
  105.     0xa644, 0x2b0f, 0x8e52, 0x8a81, 0x9558, 0xed0a, 0xcff4, 0x217, 
  106.     0x5e98, 0x1d53, 0x59f2, 0x9b70, 0xb4b1, 0x2bd9, 0xbbec, 0xfa03, 
  107.     0xbf55, 0x3651, 0x6d8, 0xcdc0, 0x132d, 0x9c88, 0x349b, 0x4ac9, 
  108.     0x89c6, 0x5cc1, 0x7826, 0xc599, 0x763f, 0x735f, 0x1e1, 0xc668, 
  109.     0x67b6, 0xc22, 0xbc5, 0x2324, 0x9c9c, 0xc8fe, 0xf55a, 0x5218, 
  110.     0x4406, 0x2fff, 0xfffb, 0xe0e6, 0x12b0, 0xc3e9, 0xc119, 0x883c, 
  111.     0xd877, 0x8973, 0xf93f, 0x3684, 0xfffb, 0x7eeb, 0x730c, 0x7cbb, 
  112.     0xcea, 0x856, 0x8000, 
  113. };
  114.  
  115.       struct Taylor
  116.       {
  117. 	int16_t d[2];
  118. 	Taylor()
  119. 	{
  120. 	  fill(d, d + 2, 0);
  121. 	}
  122. 	Taylor(int16_t val, const Taylor& other)
  123. 	{
  124. 	  d[0] = val;
  125. 	  update(other);
  126. 	}
  127. 	int16_t predict()
  128. 	{
  129. 	  return accumulate(d, d + 2, 0);
  130. 	}
  131. 	void copyFrom(const Taylor& other, size_t end = 2)
  132. 	{
  133. 	  copy(other.d, other.d + end, d);
  134. 	}
  135. 	void update(const Taylor& other)
  136. 	{
  137. 	  d[1] = d[0] - other.d[0];
  138. 	}
  139. 	void negate()
  140. 	{
  141. 	  d[0] = -d[0];
  142. 	  d[1] = -d[1];
  143. 	}
  144. 	void clear()
  145. 	{
  146. 	  fill(d, d + 2, 0);
  147. 	}
  148.       };
  149.  
  150.       class NibbleReader
  151.       {
  152. 	uint16_t* stream;
  153. 	uint16_t block;
  154. 	uint16_t pos;
  155.       public:
  156. 	NibbleReader(uint16_t* p)
  157. 	: stream(p)
  158. 	, block(0)
  159. 	, pos(0)
  160. 	{}
  161. 	uint16_t read4()
  162. 	{
  163. 	  if (pos == 0)
  164. 	    load();
  165.  
  166. 	  uint16_t result = block & 15;
  167. 	  block >>= 4;
  168. 	  pos -= 4;
  169. 	  return result;
  170. 	}
  171. 	uint16_t read12()
  172. 	{
  173. 	  return read4() | (read4() << 4) | (read4() << 8);
  174. 	}
  175. 	void load()
  176. 	{
  177. 	  block = *stream++;
  178. 	  pos = 16;
  179. 	}
  180.       };
  181.  
  182.       class ArReader
  183.       { // http://w...content-available-to-author-only...s.com/ac_arithmetic.html
  184. 	uint16_t* stream;
  185. 	uint16_t* table;
  186. 	uint16_t* tableEnd;
  187. 	uint16_t high;
  188. 	uint16_t low;
  189. 	uint16_t code;
  190. 	uint16_t source;
  191. 	uint16_t sourceUsed;
  192.       public:
  193. 	ArReader(uint16_t* p, uint16_t* t, uint16_t* e)
  194. 	: stream(p + 1)
  195. 	, table(t) // Zero value in *(t-1) expected
  196. 	, tableEnd(e)
  197. 	, high(0xFFFF)
  198. 	, low(0)
  199. 	, code(*p)
  200. 	, source(0)
  201. 	, sourceUsed(0)
  202. 	{}
  203. 	uint16_t read()
  204. 	{
  205. 	  uint32_t range = static_cast<uint32_t>(high) - low + 1;
  206.  
  207. 	  uint16_t temp =
  208. 	    ((static_cast<uint32_t>(code) - low + 1) * *(tableEnd - 1) - 1) /
  209. 	    range;
  210.  
  211. 	  uint16_t result = searchTable(temp);
  212. 	  high = low + ((range * table[result]) / *(tableEnd - 1)) - 1;
  213. 	  low = low + (range * table[result - 1]) / *(tableEnd - 1);
  214.  
  215. 	  while (true) {
  216. 	    if ((high ^ low) & 0x8000) // different msb
  217. 	    {
  218. 	      if ((high & 0x4000) < (low & 0x4000))
  219. 	      {
  220. 		code ^= 0x4000;
  221. 		low &= 0x3FFF;
  222. 		high |= 0x4000;
  223. 		shift();
  224. 	      }
  225. 	      else
  226. 	      {
  227. 		break;
  228. 	      }
  229. 	    }
  230. 	    else // equal msb
  231. 	    {
  232. 	      shift();
  233. 	    }
  234. 	  }
  235.  
  236. 	  return result;
  237. 	}
  238. 	uint16_t searchTable(uint16_t x)
  239. 	{
  240. 	  uint16_t* p = table;
  241. 	  while (x >= *p)
  242. 	    ++p;
  243. 	  return p - table;
  244. 	}
  245. 	void shift()
  246. 	{
  247. 	  low <<= 1;
  248. 	  high <<= 1;
  249. 	  high |= 1;
  250. 	  code <<= 1;
  251. 	  code |= readBit();
  252. 	}
  253.         uint16_t readBit()
  254. 	{
  255. 	  if (sourceUsed == 0)
  256. 	    load();
  257.  
  258. 	  return (source >> --sourceUsed) & 1;
  259. 	}
  260. 	void load()
  261. 	{
  262. 	  source = *stream++;
  263. 	  sourceUsed = 16;
  264. 	}
  265.       };
  266.  
  267.       int16_t get(size_t r, size_t c)
  268.       {
  269. 	if (c == NUM_COLUMNS - 1)
  270. 	  c = 0; // last column is equal to the first one
  271.  
  272. 	Taylor lastTaylor;
  273. 	NibbleReader nibbleReader(compacted + 1 + INDEX_FULL_SIZE);
  274.  
  275. 	ArReader
  276. 	  arReader(compacted + 1 + INDEX_FULL_SIZE + nibbleReader.read12(),
  277. 		   compacted + 1,
  278. 		   compacted + 1 + INDEX_FULL_SIZE);
  279.  
  280. 	for (size_t row = 0; row < NUM_ROWS; ++row)
  281. 	{
  282. 	  for (size_t col = 0; col < NUM_COLUMNS - 1; ++col)
  283. 	  { // last column is not encoded
  284. 	    uint16_t value = arReader.read();
  285. 	    int16_t fixup = static_cast<int16_t>(value) - ZERO_POS;
  286.  
  287. 	    // Value was encoded in nibbles
  288. 	    if (value == INDEX_SIZE + BIN4)
  289. 	    {
  290. 	      uint16_t x16 = nibbleReader.read4();
  291.  
  292. 	      if (x16 & 0x8)
  293. 		fixup = -ZERO_POS - 1 - (x16 & ~0x8);
  294. 	      else
  295. 		fixup = INDEX_PLUS + 1 + x16;
  296. 	    }
  297. 	    else if (value == INDEX_SIZE + BIN12)
  298. 	    {
  299. 	      fixup = static_cast<int16_t>(nibbleReader.read12());
  300. 	      fixup |= -(fixup & (1 << 11));
  301. 	    }
  302.  
  303. 	    // Restore original value (extrapolate and add a correction)
  304. 	    Taylor taylor(lastTaylor.predict() + fixup, lastTaylor);
  305.  
  306. 	    // Hack for values, jumping from -100 to 100
  307. 	    if (taylor.d[0] <= VALUE_LIMIT)
  308. 	    {
  309. 	      taylor.d[0] = -taylor.d[0];
  310. 	      lastTaylor.negate();
  311. 	      taylor.update(lastTaylor);
  312. 	    }
  313.  
  314. 	    // Got the result
  315. 	    if (row == r && col == c)
  316. 	      return taylor.d[0];
  317.  
  318. 	    // This will be used to predict the next value
  319. 	    // (derivative is not updated for the first column)
  320. 	    lastTaylor.copyFrom(taylor, col? 3: 1);
  321. 	  }
  322. 	}
  323.  
  324. 	return 0; // Out of range
  325.       }
  326.  
  327.       int main()
  328.       {
  329. 	for (int row = 0; row < NUM_ROWS; ++row)
  330. 	  for (int col = 0; col < NUM_COLUMNS; ++col)
  331. 	    if (get(row, col) != my_array[row][col])
  332. 	    {
  333. 	      cerr
  334. 		<< "Mismatch in row " << row << ", column " << col
  335. 		<< "; got " << get(row, col)
  336. 		<< ", must be " << my_array[row][col]
  337. 		<< "\n";
  338. 	      exit(1);
  339. 	    }
  340.  
  341. 	return 0;
  342.       }
  343.  
Success #stdin #stdout 0.42s 2724KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
Standard output is empty
https://ideone.com/xJ9QI
language:
C++ (gcc 8.3)
created:
12 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!