// src/mt19937_64.rs
//
// Copyright 2015 David Creswick
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

use std::cmp::max;
use std::default::Default;
use std::mem;
use std::num::Wrapping;

use rand::{Rng, SeedableRng};

const NN: usize = 312;
const MM: usize = 156;
const MATRIX_A: Wrapping<u64> = Wrapping(0xB5026F5AA96619E9);
const UM: Wrapping<u64> =  Wrapping(0xFFFFFFFF80000000); // Most significant 33 bits
const LM: Wrapping<u64> =  Wrapping(        0x7FFFFFFF); // Least significant 31 bits

/// The 64-bit flavor of the Mersenne Twister pseudorandom number
/// generator.
#[allow(non_camel_case_types)]
#[derive(Copy)]
pub struct MT19937_64 {
    idx: usize,
    state: [Wrapping<u64>; NN],
}

// FIXME: This manual implementation will probably become unnecessary someday.
impl Clone for MT19937_64 {
    #[inline]
    fn clone(&self) -> MT19937_64 {
        use std::ptr;
        unsafe {
            let mut mt = mem::uninitialized();
            ptr::copy_nonoverlapping(self, &mut mt, 1);
            mt
        }
    }
}

impl SeedableRng<u64> for MT19937_64 {
    #[inline]
    fn from_seed(seed: u64) -> MT19937_64 {
        let mut mt: MT19937_64 = unsafe { mem::uninitialized() };
        mt.reseed(seed);
        mt
    }

    fn reseed(&mut self, seed: u64) {
        self.idx = NN;
        self.state[0] = Wrapping(seed);
        for i in 1..NN {
            self.state[i] = Wrapping(6364136223846793005) * (self.state[i-1] ^ (self.state[i-1] >> 62)) + Wrapping(i as u64);
        }
    }
}

impl<'a> SeedableRng<&'a [u64]> for MT19937_64 {
    #[inline]
    fn from_seed(seed: &[u64]) -> MT19937_64 {
        let mut mt: MT19937_64 = unsafe { mem::uninitialized() };
        mt.reseed(seed);
        mt
    }

    fn reseed(&mut self, key: &[u64]) {
        self.reseed(19650218u64);
        let mut i = 1;
        let mut j = 0;
        for _ in 0 .. max(NN, key.len()) {
            self.state[i] = (self.state[i] ^ ((self.state[i-1] ^ (self.state[i-1]>>62)) * Wrapping(3935559000370003845))) + Wrapping(key[j]) + Wrapping(j as u64);
            i += 1;
            if i >= NN {
                self.state[0] = self.state[NN-1];
                i = 1;
            }
            j += 1;
            if j >= key.len() {
                j = 0;
            }
        }
        for _ in 0 .. NN-1 {
            self.state[i] = (self.state[i] ^ ((self.state[i-1] ^ (self.state[i-1]>>62)) * Wrapping(2862933555777941757))) - Wrapping(i as u64);
            i += 1;
            if i >= NN {
                self.state[0] = self.state[NN-1];
                i = 1;
            }
        }
        self.state[0] = Wrapping(1 << 63);
    }
}

impl Rng for MT19937_64 {
    #[inline]
    fn next_u32(&mut self) -> u32 {
        self.next_u64() as u32
    }

    #[inline]
    fn next_u64(&mut self) -> u64 {
        if self.idx >= NN {
            self.fill_next_state();
        }
        let Wrapping(x) = self.state[self.idx];
        self.idx += 1;
        temper(x)
    }
}

impl MT19937_64 {
    /// Create a new Mersenne Twister random number generator using
    /// the default fixed seed.
    #[inline]
    pub fn new_unseeded() -> MT19937_64 {
        SeedableRng::from_seed(5489u64)
    }

    fn fill_next_state(&mut self) {
        static MAG01: [Wrapping<u64>; 2] = [Wrapping(0), MATRIX_A];
        for i in 0 .. NN-MM {
            let x = (self.state[i]&UM) | (self.state[i+1]&LM);
            self.state[i] = self.state[i+MM] ^ (x>>1) ^ MAG01[(x.0&1) as usize];
        }
        for i in NN-MM .. NN-1 {
            let x = (self.state[i]&UM) | (self.state[i+1]&LM);
            self.state[i] = self.state[i+MM-NN] ^ (x>>1) ^ MAG01[(x.0&1) as usize];
        }
        let x = (self.state[NN-1]&UM) | (self.state[0]&LM);
        self.state[NN-1] = self.state[MM-1] ^ (x>>1) ^ MAG01[(x.0&1) as usize];
        self.idx = 0;
    }

    /// Recover the internal state of a Mersenne Twister instance
    /// from 312 consecutive outputs of the algorithm.
    ///
    /// The returned `MT19937_64` is guaranteed to identically
    /// reproduce subsequent outputs of the RNG that was sampled. This
    /// is why one does not use Mersenne Twister in cryptography.
    ///
    /// Fails if the length of the slice is not exactly 312.
    pub fn recover(samples: &[u64]) -> MT19937_64 {
        assert!(samples.len() == NN);
        let mut mt: MT19937_64 = unsafe { mem::uninitialized() };
        for (in_, out) in Iterator::zip(samples.iter(), mt.state.iter_mut()) {
            *out = Wrapping(untemper(*in_));
        }
        mt.idx = NN;
        mt
    }
}

#[inline]
fn temper(mut x: u64) -> u64 {
    x ^= (x >> 29) & 0x5555555555555555;
    x ^= (x << 17) & 0x71D67FFFEDA60000;
    x ^= (x << 37) & 0xFFF7EEE000000000;
    x ^=  x >> 43;
    x
}

#[inline]
fn untemper(mut x: u64) -> u64 {
    // reverse "x ^=  x >> 43;"
    x ^=  x >> 43;

    // reverse "x ^= (x << 37) & 0xFFF7EEE000000000;"
    x ^= (x << 37) & 0xFFF7EEE000000000;

    // reverse "x ^= (x << 17) & 0x71D67FFFEDA60000;"
    x ^= (x << 17) & 0x00000003eda60000;
    x ^= (x << 17) & 0x00067ffc00000000;
    x ^= (x << 17) & 0x71d0000000000000;

    // reverse "x ^= (x >> 29) & 0x5555555555555555;"
    x ^= (x >> 29) & 0x0000000555555540;
    x ^= (x >> 29) & 0x0000000000000015;

    x
}

impl Default for MT19937_64 {
    #[inline]
    fn default() -> MT19937_64 {
        MT19937_64::new_unseeded()
    }
}

#[test]
fn test_64bit_seeded() {
    let mt: MT19937_64 = SeedableRng::from_seed(0x123456789abcdefu64);
    for (&Wrapping(x), &y) in mt.state.iter().zip(STATE_SEEDED_BY_U64.iter()) {
        assert!(x == y);
    }
}

#[test]
fn test_64bit_slice_seeded() {
    let mt: MT19937_64 = SeedableRng::from_seed(&[0x12345u64, 0x23456u64,
                                                  0x34567u64, 0x45678u64][..]);
    for (&Wrapping(x), &y) in mt.state.iter().zip(STATE_SEEDED_BY_SLICE.iter()) {
        assert!(x == y);
    }
}

#[test]
fn test_64bit_output() {
    let mut mt: MT19937_64 = SeedableRng::from_seed(&[0x12345u64, 0x23456u64,
                                                      0x34567u64, 0x45678u64][..]);
    for x in TEST_OUTPUT.iter() {
        assert!(mt.next_u64() == *x);
    }
}

#[cfg(test)]
static STATE_SEEDED_BY_U64: [u64; NN] = [
    81985529216486895,  1791783545744114692,   604013923869356214,
    296213109303876097,  6923710694692163889, 13205119584621876597,
    11033091707158564593,  6240080577330458046,  6391091230317693851,
    3142493543514746907, 12293512030652881353,  2062181304702107066,
    4654451180172173246, 11441927625268696480, 14365762244911147912,
    18029629756502187134,  7331740749253838089,  6610823624717752953,
    18246827935145324330, 10680693382789357896, 14638587269153607702,
    8633864969252503238,  9566721819554316817,  6143884850529323630,
    10906314227572090523, 10907689704186884094, 14802274577639469926,
    3401801631809463260, 12344530338405822920,  4956444630876681887,
    14793227046740343780,  2037099132120965306,  6441991539629590226,
    12964898550824984632, 12259055633860243540, 14928425551652759873,
    6183427467171036862,  5400792687527135416,  6408845498294423467,
    14659391812871673609,  4904880327708328682, 16004025791281387640,
    15740202225231696585, 10802686250382853037,  7329769865301751535,
    13158972625497029123,  1462611033182269787, 14970257040982151726,
    1785898624104826137,  4018756387164900502,  4658545214302937232,
    12929223690873008,  2189306884804589348,  7407145675715022217,
    17926372328502417694, 10393381392936406352,  1519619275753850274,
    3949038713401990067, 11368953568800795569, 14644212788164115378,
    3789816819491622233, 11192801356205178850,  4208023585462672030,
    1023784985917378565, 10420842316752161825,  2134329768317386600,
    10009382106317716874,  1061284023672113451, 12181229230564241875,
    2507813008873369858,   491377794161482656,  1753720522338319207,
    13405595244999805539,    64138164073787990, 12320498326396355432,
    11779553027707203565,  1112643061962194511, 10649695524108772656,
    3606504838176593944,  9766617496074413191, 14913709254545614513,
    17786322741581487003,  5869942273187627530,  3835450227478597954,
    13916038151171491566,  9003731200188680958, 13696512428165460009,
    12769326485216300262, 12033565608077516908, 15827437595453457839,
    9081822962023519126, 13729953554469296870, 14683166359047427696,
    15953207475273187220, 11421240586025801449,  6102838545700560046,
    8293613563303025699, 18329916755578359899,  3729358436700194778,
    13602102386364221365, 18211263000630633615,  2749239093744591105,
    12313868987179764883, 17343761002274333924,  2968322519436596739,
    5919354133661866992,  5667644606957959367, 16954289952443151673,
    6384564465401602462,  6843583206213995104,  3776716724405842555,
    2493814834786559246, 11695884210276578790,  5277431840939578757,
    10770107586796995750, 17730900742717009223,  8605167815254274408,
    1834669206241314538, 13386897725877112216,  8421817962304670345,
    3844006852001387104, 17749959801471998809,  2889900833057928524,
    12405269285066822359,  3981740344450333421,  3599739118559850278,
    8389005402870866988,  8056933947569598312, 10816660871222865141,
    11505725300589228268,  3985627793803701336, 15054413190071596027,
    11203097170853309212, 15982971736834602699,  6491081411731397038,
    5032971170848265546,  6707451110701085623, 16660695088837779847,
    1292203148674769086,  6831779758890451953,  1273082991379300796,
    4816977719034172057,  1227740039402122566, 11774987201411346397,
    11000284286582049731, 17607239656417890430,  3905066968366541195,
    4536740449007338242,  2456916080186164718, 13537108370327412075,
    13061598491450163979, 11333202480645029932,   417941510859004590,
    9359271011890692911, 15516319166549451395, 13173397503737614619,
    9674015854587717889, 14319184628944941860, 11461157460156307321,
    8470403413024138814, 15360905514977934003, 12434310806890127761,
    5668365606955471993, 18003197684936171963,  1713810861376697852,
    15749877514364305905,  4878375384529970736, 10666006561241761348,
    14715617963938206966,  1521473769572924090,  2657449321132939612,
    2734687391377742295,  9191020263043986551,  3643264948674757483,
    3578760580494634877, 14476122573168434344,  5845171853656694719,
    2788922203841648407,  8033902436969920701, 15882240616055029183,
    10128649212328844992,  1006259818353534671, 10264281777791075353,
    14847900867811057270,  4583748711890019913, 12309703074219471502,
    6941586872827696982,  3745802569518482438,  2894905614301835210,
    7335498349179313471, 17387470308519765668,  5433932443266606362,
    14371289029466467199, 10374778237833390477,  4776388970536189925,
    1816258018112460759,  7276434563887405711, 16169106984301305275,
    2254368830142317086,   862540640621281805, 10871008795911557137,
    5418842661241340192,  8974351759208567447, 11487152290323754793,
    6225694545187332187,  8855197444497053343, 18275914506735836308,
    15923775749857325146,  9368682834059667317, 11886054335274482108,
    7248432174334290232, 12948783033451316952, 16413608400513391398,
    5395767780443283798, 18177623543956751905,    79375280521917137,
    88402592242998933, 12424141698864532490,  3506566252331882050,
    15975266567245850741,  1432897815408109466, 16349251380638817519,
    6935098887364586074, 18134822759325374430, 13289363508364170681,
    17834238808526510016, 17849560451691129641,  9625769816461176133,
    11863959125305356127, 10159009876047017790,  2270215115277990514,
    10335825866811872497, 10689024433546771615, 15402424147055348610,
    15938065922691302807,  2186961644917581807, 17743571035416418031,
    18189648471277689961,  6166355766821279120, 10383219040510737004,
    11833675146901362758,  5015772504555552997,  9466405520041559302,
    3890913253695856295,  4199704116067417423,  9536371304943548632,
    2803365044720657736,  8294255653172709023,  4949919456077810878,
    13179261481143762828, 17169909754051256048,  4131432844970114738,
    4710999410297164870,    79362067765594744, 16470327280077266966,
    9045473539587198896, 12919627169836945949, 13687654190330302836,
    10291072957815184576, 12669641865365818653,  7897093283919933559,
    8774498508392009667, 17733475113271198752,  7563755234646661166,
    14400591327379984971,  3324264236052574129, 17443813642144444967,
    14253897042581220703,  2880238993097547320, 15284741228796762341,
    7340174240767294332, 12625292554856200456,  5758675782428868050,
    9228068019606031400, 10429450678612884084, 10614459586086506961,
    8139473603542569003, 10814373026410622583,  4011828072370763431,
    5540639547564078450,  8161872968472326991, 14300024824953534415,
    2135593418061136886,   628742652522331481,  6231251893857464513,
    4921212502846493405,  3364411375653074378, 11642530140974623649,
    3585719091449989575, 13407410611541963292,   818686227520123496,
    14065717611881608555, 12577192992972897388, 11463784563875803771,
    13136698746312533867, 12860699992103060892, 14670774660241317870,
    17890504795855548626,  9605805324699469031,  6412165298585204844,
    678652354371239765,  1676294081870876190,  8995213228755328628,
    16153912129483269824,  4748371455456104055, 15348364511313581551,
    10147917971634761646,  8545060171396429423,  1413465926917694602,
    14925734419407102839, 12220731818240280730, 16069120541665094639];

#[cfg(test)]
static STATE_SEEDED_BY_SLICE: [u64; NN] = [
    9223372036854775808, 11906848963885322876,  2366922794640711022,
    1039674428528475923,  3532693117775135603,  4008869120475271563,
    17747978335774359060, 14972770412725955390,  3858352546733417902,
    5345612019402505742,  7758276406156335750,  6001919363687313191,
    2547684760428594376, 10392438263683240802,  8968713240322292543,
    14406031157106643586,  3964892302884586199,  4062779262181776338,
    16580562407461089692, 12589582373071700978, 15750410628907513263,
    11310313695732438417, 15633947903704351334, 18222915079249959102,
    14276842509123192698,  6930516820607322720, 15365297350207276392,
    8213498164667572384,  2152779499883018824, 13966133250965870373,
    10240771306433844368, 14307168365450430565, 11040552514370714895,
    2533639886572881353,  9185674266382191104, 14346522715072559330,
    14953130301049970203, 15180875257012729653,  1964218244082486892,
    6405341251447791135, 11761351367919151087,  3167849151490256233,
    1558463723413910204,  4047863190093408546, 18058286074055427032,
    7317716653786836408, 17652746742586216402,  6760726663271410936,
    17327162216619411332, 14378483978609909724,  2131305364409541273,
    16726816105924499124,  9210182618446396614,  6445236935855240210,
    2449852682157161554,  6263322181339295035,  6914873478768120792,
    5860130632837132104,   186277607674026464,  8693717914902489281,
    7758418146359524665, 15036826958618606534,  9020793635699266615,
    16681251223473385742,    28790294402912016, 16206435735609334231,
    4276151532523641810, 17083626852979795873, 15898385609762430430,
    7413408202278515221,  6817627683834618267, 14308485276827746090,
    18090035704918431532, 17439280401835765253,  9920708817472459377,
    17596494468204334260,  8987881797008539223,  5918674948541707799,
    1543803321776487042,  4821966630366015761,  6984206624831694542,
    18035097434221633402, 16518401100072381614,  5428007001392054535,
    9346296170248113790,  8189955722678051590,  8166058216851039534,
    5495705342415866955, 17808354540446325627,   239520990467221973,
    6774849326041393260, 10588207904624720754, 10508799966524468092,
    7314664134865101995, 13330428261370692994,  7315653225976003522,
    14752942780607268272,  7386201019337306472, 11074689804650560421,
    13251660986032530719,  5804947702644260322,  4725780123182197186,
    14723414396390934648, 17218644267506543502,  3765965502701684261,
    3604066792039542938, 17362622132271438638, 16875585620489891286,
    4514873179377317837,  9308332464342002048,   188567056399553729,
    2705891390840254977, 13168211036505445211,  6269486280874873673,
    7485168387351336182, 18163642631999775128,  1152810103088079572,
    10314639940536358822,  3898255220108939804, 12461785437811001117,
    14080290123823101293,  8186155435419807960, 15303633033386294539,
    16368372373799902572, 16908740919128790793,  7790227374552432870,
    3370942488781431421,  5321468336044485163,  9099464595185099272,
    12686959047356864605, 15044140022848320070, 17853451835866976714,
    15090981758503849092,  7259421845759239401, 17006897347539967307,
    11864671541689053919,  3589888654233765802, 17133675018297505478,
    10350789726188603215,  5805562075215454138,  6078313175004902257,
    15159259062264755988,  4606458898055115768,  8595623819596928458,
    10874665851480906795,  9541918022094056642,  7910492575696069731,
    14253990448274524591, 11587237029817781157,  8898150318182787055,
    9129838589367321414, 13328388735954292593, 17796364720220167616,
    8597915663953875317,  2018393986500621684, 13175088550176860755,
    1094714902621613987,  8807072083132332364, 14066712852747710486,
    8999937541229388765,  9112579609010886366,  3356356224211045166,
    4853220202257544572, 10444170003278214107,  1052947048789639998,
    1827363901483126686, 15601796185052218891, 11983393206205800981,
    13121259426767810380,  9553182408890238788,  3524616048466714316,
    6892649984567035338, 12306438715782265106,  6929080643604729797,
    10488981816167940982,  3526412707713863908, 13925166661639043721,
    3120081366166455125, 12074161212633821290,  8925556738888186426,
    8406789658316452783,  9694538352161485523, 10289430848224100722,
    18339763433413173518,   920701930237947644,  2818534229746382484,
    15586266845147401960, 15899701243238871280, 16635647107881635975,
    9436756036899151216, 15869514808861706046,  5377602059434346830,
    881770225321498601,  4972312278793239748, 13470351019631907494,
    15034318809458171738,  5731070042043014264,  8320132766717278470,
    10458189336915616638,  1667073997662158604, 16517673739089400781,
    9780691302643590844,  1888289618321818453,  9933541508716288366,
    7342610457445606218,  7803853065241579933, 16080562130327610447,
    1972475373062043809,  2445190066574935778, 16665102149891804629,
    3729422737555637347, 15627314164145375682,  5165153149317749684,
    1851884315110708417, 17883927180228563680, 18254359763092694099,
    8051365941650618695,  9309548635890345699,   263899551166810219,
    8033066166620504604,  1803538698706947772, 15690816866062867851,
    2223071627096777838, 16579019743023052535, 13162141123527406604,
    303902486500873444, 16218914908155140970,  7664023458078374736,
    4211782740976596951,  9852041272241461170,  6869872483205843001,
    17682149354114760050,  2193325669751198964, 18183657083348212654,
    17956232649388583703, 17885330999340609820, 15086254682221937578,
    2661870345357193374,   824878816789277867, 12195389379115398039,
    8118221630987421587, 18090490544414672070,  3156155883929578921,
    12762751134459068851,  5164889741276341708,  7716579664961922616,
    6242366333202544756, 12664253830362087207,  9768406626100492062,
    18276764521683238961,  1981044956697052467, 10832499469197748318,
    5822985766637060664,  6589889711659930475,  9415661998487594717,
    10372021223737917970,  9008931203668510507, 12048148024357412075,
    3686532183913588371,  5820540035396407131,  2469949645283598917,
    7565651341821095299,  5197273755353226573,  5880387491135838713,
    6730520764971726984,  1627471928992184271,  3164256474773483571,
    10891562408149301222, 16127745993321845961,  6089459304242537067,
    14350402415135140161, 18114662877700924330,  9359478022685723144,
    13171912360297387648, 11237048723429932701,  7800636006979166579,
    9289540047600185720,   998919556795606739,  2392987300775238592,
    12059856439917591208, 10857265173257971923, 15122926213571157845,
    12840523308034798986,  1895411170536695373,  4057105479876795998,
    6310373902706727113,  8171018879503955734,  1863802676851052784,
    14129048472806811984,  5840360500395533120,  6408991116533430303,
    17554903680613470616,  7251872690709372415,  3092878727914765317,
    11961922532714192260,  6989054962715140021,  7236766257488832035,
    5794159225877267393, 15592220508405117453, 13266267152869308514,
    12964887929723524270, 15463019828589740680,  5488263871426190118,
    18181645947417462531,  3280221468542061885,  4905596419059094060,
    3548266682460467135,  1326982712296517686, 12179707347278731003,
    14375460083628818157, 12930964925920180049,  9341456266356607437];

#[cfg(test)]
static TEST_OUTPUT: [u64; 1000] = [
    7266447313870364031,  4946485549665804864, 16945909448695747420,
    16394063075524226720,  4873882236456199058, 14877448043947020171,
    6740343660852211943, 13857871200353263164,  5249110015610582907,
    10205081126064480383,  1235879089597390050, 17320312680810499042,
    16489141110565194782,  8942268601720066061, 13520575722002588570,
    14226945236717732373,  9383926873555417063, 15690281668532552105,
    11510704754157191257, 15864264574919463609,  6489677788245343319,
    5112602299894754389, 10828930062652518694, 15942305434158995996,
    15445717675088218264,  4764500002345775851, 14673753115101942098,
    236502320419669032, 13670483975188204088, 14931360615268175698,
    8904234204977263924, 12836915408046564963, 12120302420213647524,
    15755110976537356441,  5405758943702519480, 10951858968426898805,
    17251681303478610375,  4144140664012008120, 18286145806977825275,
    13075804672185204371, 10831805955733617705,  6172975950399619139,
    12837097014497293886, 12903857913610213846,   560691676108914154,
    1074659097419704618, 14266121283820281686, 11696403736022963346,
    13383246710985227247,  7132746073714321322, 10608108217231874211,
    9027884570906061560, 12893913769120703138, 15675160838921962454,
    2511068401785704737, 14483183001716371453,  3774730664208216065,
    5083371700846102796,  9583498264570933637, 17119870085051257224,
    5217910858257235075, 10612176809475689857,  1924700483125896976,
    7171619684536160599, 10949279256701751503, 15596196964072664893,
    14097948002655599357,   615821766635933047,  5636498760852923045,
    17618792803942051220,   580805356741162327,   425267967796817241,
    8381470634608387938, 13212228678420887626, 16993060308636741960,
    957923366004347591,  6210242862396777185,  1012818702180800310,
    15299383925974515757, 17501832009465945633, 17453794942891241229,
    15807805462076484491,  8407189590930420827,   974125122787311712,
    1861591264068118966,   997568339582634050, 18046771844467391493,
    17981867688435687790,  3809841506498447207,  9460108917638135678,
    16172980638639374310,   958022432077424298,  4393365126459778813,
    13408683141069553686, 13900005529547645957, 15773550354402817866,
    16475327524349230602,  6260298154874769264, 12224576659776460914,
    6405294864092763507,  7585484664713203306,  5187641382818981381,
    12435998400285353380, 13554353441017344755,   646091557254529188,
    11393747116974949255, 16797249248413342857, 15713519023537495495,
    12823504709579858843,  4738086532119935073,  4429068783387643752,
    585582692562183870,  1048280754023674130,  6788940719869959076,
    11670856244972073775,  2488756775360218862,  2061695363573180185,
    6884655301895085032,  3566345954323888697, 12784319933059041817,
    4772468691551857254,  6864898938209826895,  7198730565322227090,
    2452224231472687253, 13424792606032445807, 10827695224855383989,
    11016608897122070904, 14683280565151378358,  7077866519618824360,
    17487079941198422333,  3956319990205097495,  5804870313319323478,
    8017203611194497730,  3310931575584983808,  5009341981771541845,
    6930001938490791874, 14415278059151389495, 11001114762641844083,
    6715939435439735925,   411419160297131328,  4522402260441335284,
    3381955501804126859, 15935778656111987797,  4345051260540166684,
    13978444093099579683,  9219789505504949817,  9245142924137529075,
    11628184459157386459,  7242398879359936370,  8511401943157540109,
    11948130810477009827,  6865450671488705049, 13965005347172621081,
    15956599226522058336,  7737868921014130584,  2107342503741411693,
    15818996300425101108, 16399939197527488760, 13971145494081508107,
    3910681448359868691,  4249175367970221090,  9735751321242454020,
    12418107929362160460,   241792245481991138,  5806488997649497146,
    10724207982663648949,  1121862814449214435,  1326996977123564236,
    4902706567834759475, 12782714623891689967,  7306216312942796257,
    15681656478863766664,   957364844878149318,  5651946387216554503,
    8197027112357634782,  6302075516351125977, 13454588464089597862,
    15638309200463515550, 10116604639722073476, 12052913535387714920,
    2889379661594013754, 15383926144832314187,  7841953313015471731,
    17310575136995821873,  9820021961316981626, 15319619724109527290,
    15349724127275899898, 10511508162402504492,  6289553862380300393,
    15046218882019267110, 11772020174577005930,  3537640779967351792,
    6801855569284252424, 17687268231192623388, 12968358613633237218,
    1429775571144180123, 10427377732172208413, 12155566091986788996,
    16465954421598296115, 12710429690464359999,  9547226351541565595,
    12156624891403410342,  2985938688676214686, 18066917785985010959,
    5975570403614438776, 11541343163022500560, 11115388652389704592,
    9499328389494710074,  9247163036769651820,  3688303938005101774,
    2210483654336887556, 15458161910089693228,  6558785204455557683,
    1288373156735958118, 18433986059948829624,  3435082195390932486,
    16822351800343061990,  3120532877336962310, 16681785111062885568,
    7835551710041302304,  2612798015018627203, 15083279177152657491,
    6591467229462292195, 10592706450534565444,  7438147750787157163,
    323186165595851698,  7444710627467609883,  8473714411329896576,
    2782675857700189492,  3383567662400128329,  3200233909833521327,
    12897601280285604448,  3612068790453735040,  8324209243736219497,
    15789570356497723463,  1083312926512215996,  4797349136059339390,
    5556729349871544986, 18266943104929747076,  1620389818516182276,
    172225355691600141,  3034352936522087096,  1266779576738385285,
    3906668377244742888,  6961783143042492788, 17159706887321247572,
    4676208075243319061, 10315634697142985816, 13435140047933251189,
    716076639492622016, 13847954035438697558,  7195811275139178570,
    10815312636510328870,  6214164734784158515, 16412194511839921544,
    3862249798930641332,  1005482699535576005,  4644542796609371301,
    17600091057367987283,  4209958422564632034,  5419285945389823940,
    11453701547564354601,  9951588026679380114,  7425168333159839689,
    8436306210125134906, 11216615872596820107,  3681345096403933680,
    5770016989916553752, 11102855936150871733, 11187980892339693935,
    396336430216428875,  6384853777489155236,  7551613839184151117,
    16527062023276943109, 13429850429024956898,  9901753960477271766,
    9731501992702612259,  5217575797614661659, 10311708346636548706,
    15111747519735330483,  4353415295139137513,  1845293119018433391,
    11952006873430493561,  3531972641585683893, 16852246477648409827,
    15956854822143321380, 12314609993579474774, 16763911684844598963,
    16392145690385382634,  1545507136970403756, 17771199061862790062,
    12121348462972638971, 12613068545148305776,   954203144844315208,
    1257976447679270605,  3664184785462160180,  2747964788443845091,
    15895917007470512307, 15552935765724302120, 16366915862261682626,
    8385468783684865323, 10745343827145102946,  2485742734157099909,
    916246281077683950, 15214206653637466707, 12895483149474345798,
    1079510114301747843, 10718876134480663664,  1259990987526807294,
    8326303777037206221, 14104661172014248293, 15531278677382192198,
    3874303698666230242,  3611366553819264523,  1358753803061653874,
    1552102816982246938, 14492630642488100979, 15001394966632908727,
    2273140352787320862, 17843678642369606172,  2903980458593894032,
    16971437123015263604, 12969653681729206264,  3593636458822318001,
    9719758956915223015,  7437601263394568346,  3327758049015164431,
    17851524109089292731, 14769614194455139039,  8017093497335662337,
    12026985381690317404,   739616144640253634, 15535375191850690266,
    2418267053891303448, 15314073759564095878, 10333316143274529509,
    16565481511572123421, 16317667579273275294, 13991958187675987741,
    3753596784796798785,  9078249094693663275,  8459506356724650587,
    12579909555010529099,  7827737296967050903,  5489801927693999341,
    10995988997350541459, 14721747867313883304,  7915884580303296560,
    4105766302083365910, 12455549072515054554, 13602111324515032467,
    5205971628932290989,  5034622965420036444,  9134927878875794005,
    11319873529597990213, 14815445109496752058,  2266601052460299470,
    5696993487088103383,  6540200741841280242,  6631495948031875490,
    5328340585170897740, 17897267040961463930,  9030000260502624168,
    14285709137129830926, 12854071997824681544, 15408328651008978682,
    1063314403033437073, 13765209628446252802,   242013711116865605,
    4772374239432528212,  2515855479965038648,  5872624715703151235,
    14237704570091006662,   678604024776645862, 12329607334079533339,
    17570877682732917020,  2695443415284373666,  4312672841405514468,
    6454343485137106900,  8425658828390111343, 16335501385875554899,
    5551095603809016713, 11781094401885925035,  9395557946368382509,
    9765123360948816956, 18107191819981188154, 16049267500594757404,
    16349966108299794199,  1040405303135858246,  2366386386131378192,
    223761048139910454, 15375217587047847934, 15231693398695187454,
    12916726640254571028,  8878036960829635584,  1626201782473074365,
    5758998126998248293, 18077917959300292758, 10585588923088536745,
    15072345664541731497,  3559348759319842667, 12744591691872202375,
    2388494115860283059,  6414691845696331748,  3069528498807764495,
    8737958486926519702, 18059264986425101074,  3139684427605102737,
    12378931902986734693,   410666675039477949, 12139894855769838924,
    5780722552400398675,  7039346665375142557,  3020733445712569008,
    2612305843503943561, 13651771214166527665, 16478681918975800939,
    566088527565499576,  4715785502295754870,  6957318344287196220,
    11645756868405128885, 13139951104358618000, 17650948583490040612,
    18168787973649736637,  5486282999836125542,  6122201977153895166,
    17324241605502052782, 10063523107521105867, 17537430712468011382,
    10828407533637104262, 10294139354198325113, 12557151830240236401,
    16673044307512640231, 10918020421896090419, 11077531235278014145,
    5499571814940871256,  2334252435740638702, 18177461912527387031,
    2000007376901262542,  7968425560071444214,  1472650787501520648,
    3115849849651526279,  7980970700139577536, 12153253535907642097,
    8109716914843248719,  3154976533165008908,  5553369513523832559,
    10345792701798576501,  3677445364544507875, 10637177623943913351,
    7380255087060498096, 14479400372337014801, 15381362583330700960,
    204531043189704802, 13699106540959723942,  3817903465872254783,
    10972364467110284934,  2701394334530963810,  2931625600749229147,
    16428252083632828910, 11873166501966812913,  5566810080537233762,
    7840617383807795056, 10699413880206684652, 18259119259617231436,
    10332714341486317526, 10137911902863059694,   669146221352346842,
    8373571610024623455, 10620002450820868661, 12220730820779815970,
    5902974968095412898,  7931010481705150841, 16413777368097063650,
    11273457888324769727, 13719113891065284171,  8327795098009702553,
    10333342364827584837,  6202832891413866653,  9137034567886143162,
    14514450826524340059,   473610156015331016,   813689571029117640,
    13776316799690285717, 10429708855338427756,  8995290140880620858,
    2320123852041754384,  8082864073645003641,  6961777411740398590,
    10008644283003991179,  3239064015890722333, 16762634970725218787,
    16467281536733948427, 10563290046315192938,  5108560603794851559,
    15121667220761532906, 14155440077372845941, 10050536352394623377,
    15474881667376037792,  3448088038819200619,  3692020001240358871,
    6444847992258394902,  8687650838094264665,  3028124591188972359,
    16945232313401161629, 15547830510283682816,  3982930188609442149,
    14270781928849894661, 13768475593433447867, 13815150225221307677,
    8502397232429564693,   718377350715476994,  7459266877697905475,
    8353375565171101521,  7807281661994435472, 16924127046922196149,
    10157812396471387805,  2519858716882670232,  7384148884750265792,
    8077153156180046901,  3499231286164597752,  2700106282881469611,
    14679824700835879737, 14188324938219126828,  3016120398601032793,
    10858152824243889420,  9412371965669250534,  4857522662584941069,
    984331743838900386,  4094160040294753142,  2368635764350388458,
    15101240511397838657, 15584415763303953578,  7831857200208015446,
    1952643641639729063,  4184323302594028609, 16795120381104846695,
    3541559381538365280, 15408472870896842474,  5628362450757896366,
    16277348886873708846, 12437047172652330846, 10172715019035948149,
    1999700669649752791,  6217957085626135027, 11220551167830336823,
    16478747645632411810,  5437280487207382147, 11382378739613087836,
    15866932785489521505,  5502694314775516684, 16440179278067648435,
    15510104554374162846, 15722061259110909195, 10760687291786964354,
    10736868329920212671,  4166148127664495614, 14303518358120527892,
    9122250801678898571, 10028508179936801946,   216630713752669403,
    10655207865433859491,  4041437116174699233,  6280982262534375348,
    297501356638818866, 13976146806363377485, 13752396481560145603,
    11472199956603637419, 16393728429143900496, 14752844047515986640,
    1524477318846038424,  6596889774254235440,  1591982099532234960,
    8065146456116391065,  3964696017750868345, 17040425970526664920,
    11511165586176539991,  3443401252003315103, 16314977947073778249,
    16860120454903458341,  5370503221561340846, 15362920279125264094,
    2822458124714999779, 14575378304387898337,  9689406052675046032,
    2872149351415175149, 13019620945255883050, 14929026760148695825,
    8503417349692327218,  9677798905341573754,   828949921821462483,
    16110482368362750196, 15794218816553655671, 14942910774764855088,
    12026350906243760195, 13610867176871462505, 18324536557697872582,
    2658962269666727629,   327225403251576027,  9207535177029277544,
    8744129291351887858,  6129603385168921503, 18385497655031085907,
    13024478718952333892, 14547683159720717167,  5932119629366981711,
    325385464632594563,  3559879386019806291,  6629264948665231298,
    14358245326238118181, 15662449672706340765, 13975503159145803297,
    3609534220891499022,  4224273587485638227,  9274084767162416370,
    13156843921244091998, 18284750575626858789, 14664767920489118779,
    11292057742031803221, 13919998707305829132, 14473305049457001422,
    9696877879685767807,  1406758246007973837,  2429517644459056881,
    14361215588101587430, 11386164476149757528, 10474116023593331839,
    2921165656527786564, 15604610369733358953, 12955027028676000544,
    10314281035410779907,  3167047178514709947,  1088721329408346700,
    17930425515478182741,  7466411836095405617, 15534027454610690575,
    10879629128927506091, 11502219301371200635, 13915106894453889418,
    4226784327815861027, 12335222183627106346,  3648499746356007767,
    18441388887898023393, 18117929843327093625,  4237736098094830438,
    14229123019768296655,  3930112058127932690, 12663879236019645778,
    9281161952002617309,  4978473890680876319,   845759387067546611,
    1386164484606776333,  8008554770639925512, 11159581016793288971,
    18065390393740782906, 17647985458967631018,  9092379465737744314,
    2914678236848656327,  4376066698447630270, 16057186499919087528,
    3031333261848790078,  2926746602873431597,  7931945763526885287,
    147649915388326849, 15801792398814946230,  5265900391686545347,
    16173686275871890830,  7562781050481886043,  5853506575839330404,
    14957980734704564792, 10944286556353523404,  1783009880614150597,
    9529762028588888983,   822992871011696119,  2130074274744257510,
    8000279549284809219,  3514744284158856431,   128770032569293263,
    3737367602618100572, 16364836605077998543,   783266423471782696,
    4569418252658970391, 11093950688157406886, 14888808512267628166,
    4217786261273670948, 17047486076688645713, 14133826721458860485,
    17539744882220127106, 12394675039129853905,  5757634999463277090,
    9621947619435861331,  1182210208559436772, 14603391040490913939,
    17481976703660945893, 14063388816234683976,  2046622692581829572,
    8294969799792017441,  5293778434844788058, 17976364049306763808,
    399482430848083948, 16495545010129798933, 15241340958282367519,
    989828753826900814, 17616558773874893537,  2471817920909589004,
    11764082277667899978,  9618755269550400950,  1240014743757147125,
    1887649378641563002,  1842982574728131416, 13243531042427194002,
    7688268125537013927,  3080422097287486736,  2562894809975407783,
    12428984115620094788,  1355581933694478148,  9895969242586224966,
    8628445623963160889,  4298916726468199239, 12773165416305557280,
    5240726258301567487,  4975412836403427561,  1842172398579595303,
    7812151462958058676, 17974510987263071769, 14980707022065991200,
    18294903201142729875, 12911672684850242753,  8979482998667235743,
    16808468362384462073,  5981317232108359798, 12373702800369335100,
    16119707581920094765,  2782738549717633602, 15454155188515389391,
    16495638000603654629, 16348757069342790497,  7769562861984504567,
    17504300515449231559,  5557710032938318996, 11846125204788401203,
    13957316349928882624,  2738350683717432043, 15738068448047700954,
    6224714837294524999,  6081930777706411111, 11366312928059597928,
    4355315799925031482, 12393324728734964015, 15277140291994338591,
    1406052433297386355, 15859448364509213398,  1672805458341158435,
    2926095111610982994, 11056431822276774455, 12083767323511977430,
    3296968762229741153, 12312076899982286460, 17769284994682227273,
    15349428916826953443,  1056147296359223910, 18305757538706977431,
    6214378374180465222, 14279648441175008454, 17791306410319136644,
    956593013486324072,  2921235772936241950, 10002890515925652606,
    10399654693663712506,  6446247931049971441,  6380465770144534958,
    11439178472613251620, 10131486500045494660,  3692642123868351947,
    10972816599561388940,  4931112976348785580,  8213967169213816566,
    15336469859637867841, 15026830342847689383,  7524668622380765825,
    17309937346758783807,   372780684412666438,  5642417144539399955,
    18303842993081194577, 11085303253831702827, 15658163165983586950,
    8517521928922081563, 16091186344159989860, 17614656488010863910,
    4736067146481515156, 13449945221374241354, 17755469346196579408,
    13300502638545717375,  6611828134763118043, 14177591906740276597,
    9340430243077460347,  7499765399826404087,  3409518087967832469,
    9013253864026602045,  4444307427984430192,  3729283608700519712,
    13642048880719588383, 16486557958022946240,  2996465014991157904,
    10020049344596426576, 12302485648009883778,  8492591321344423126,
    17407986443716172520, 10530482934957373052, 15740662350540828750,
    1790629986901049436,  6305948377669917188, 15092985352503125323,
    928505047232899787, 14404651977039851607,  7564177565277805597,
    3411236815351677870,  7752718145953236134, 12315979971311483798,
    12477729506691004724, 14654956300924793305,  6689803038918974388,
    1540738812233000153, 13508351811701989957, 15864432023192136053,
    7990997967273843917,  7424300239290765161,    39585249496300263,
    3877436595063283319, 10710642254398044448,  4653804418844456375,
    1232267496410380283,  3690525514009038824, 15459770765077428485,
    13240346522153894145,  5674964360688390624, 16973644653010587289,
    15924280764204855206, 15196708627253442662, 17596174821341373274,
    16196745023027393691,  6980050627399795351, 17582264380857746637,
    18170372407506856324, 12108126025631005514, 15687749089493373169,
    5814107289258228434,  9381977959648494876, 15895601183088112734,
    16267869075651604263, 15228381979765852785, 11949618678312581999,
    4545324791131029438,   582725409406225185, 15282520250746126790,
    14758446535973412711,  7605613563088071833,  1111140641057375915,
    5364843095234852245,   218335432181198977,  4891472444796201742,
    4564628942836375772, 15500501278323817088,  4913946328556108657,
    2684786251736694229, 12090498456116310122,  5310885782157038567,
    5032788439854011923, 12627401038822728242, 11869662610126430929,
    17650156853043540226, 12126672500118808436, 10437658933435653256,
    13133995470637873311,  4601324715591152820,  1874350460376708372,
    5808688626286061164, 13777088437302430376,  5018451954762213522,
    2588296738534474754,  5503414509154170711,  5230497186769951796,
    13261090710400573914,  8515217303152165705, 11074538219737365303,
    15481562385740613213, 12705484409881007350, 14221931471178549498,
    12905633420087112297, 17337759164357146506, 14081997515778175224,
    17384320185513122939,  7131793076779216692, 17483217190312403109,
    900692047897995877, 14723287313048560400,  6132094372965340305,
    7572797575350925726, 12725160700431903514,   380860122911632449,
    1900504978569024571,  8423729759529914138,  7305587201606052334,
    12446871355267313320,  4615812356515386206,  3361817115406652303,
    17690418922000878428, 14632214537567910559,  2709702289926174775,
    3459675155951086144,  7788364399926538150, 16043992474431955950,
    15830963823784930267,  4216893617835797954,   538159724689093771,
    16029152738918251363, 14444848757576686696, 12941757045272633696,
    10900480525147953314, 12547307449905859302, 16001571796892398181,
    407942194622690676, 13873235372903944444, 18071603799493008777,
    1015646077646778622,  9387605808959554815, 11566702442022019410,
    7061722181092883183,  2629032108249254109,  5271820053177594520,
    12640880742139693547, 10098688629735675775,  5716304472850923064,
    3312674502353063071,  7295926377425759633,   833281439103466115,
    16316743519466861667,  9912050326606348167, 11651133878100804242,
    18026798122431692459,  6157758321723692663,  4856021830695749349,
    7074321707293278978, 10748097797809573561,  2949954440753264783,
    9813922580940661152,  9949237950172138336, 15643982711269455885,
    16078663425810239127, 12508044395364228880, 12920301578340189344,
    15368071871011048915,  1610400750626363239, 11994736084146033126,
    6042574085746186088,  4154587549267685807, 15915752367312946034,
    1191196620621769193,   467437822242538360,  2836463788873877488,
    10476401302029164984,  1716169985450737419,  5327734953288310341,
    3994170067185955262,   884431883768190063, 11019001754831208284,
    14322807384384895215,   161011537360955545,  1466223959660131656,
    5227048585229497539, 12410731857504225031,  2142243279080761103,
    17682826799106851430,  1792612570704179953, 14727410295243056025,
    1459567192481221274,  5669760721687603135, 17507918443756456845,
    10354471145847018200, 10362475129248202288, 13143844410150939443,
    6861184673150072028, 18396524361124732580,   543906666394301875,
    12476817828199026728, 11853496871128122868, 12747674713108891748,
    7986179867749890282,  9158195177777627533,  2217320706811118570,
    8631389005200569973,  5538133061362648855,  3369942850878700758,
    7813559982698427184,   509051590411815948, 10197035660403006684,
    13004818533162292132,  9831652587047067687,  7619315254749630976,
    994412663058993407];


#[bench]
fn benchmark_seeding(b: &mut ::test::Bencher) {
    b.iter(|| MT19937_64::new_unseeded());
}

#[bench]
fn benchmark_fill_next_state(b: &mut ::test::Bencher) {
    let mut mt = MT19937_64::new_unseeded();
    b.iter(|| mt.fill_next_state());
}

#[test]
fn test_untemper() {
    let x = ::rand::thread_rng().gen::<u64>();
    assert_eq!(x, untemper(temper(x)));
}

#[test]
fn test_recovery() {
    let seed = ::rand::thread_rng().gen::<u64>();
    let mut orig_mt: MT19937_64 = SeedableRng::from_seed(seed);
    // skip some samples so the RNG is in an intermediate state
    let to_skip = ::rand::thread_rng().gen_range(1, NN);
    for _ in 0..to_skip {
        orig_mt.next_u64();
    }
    let samples = orig_mt.gen_iter::<u64>().take(NN).collect::<Vec<_>>();
    let mut recovered_mt = MT19937_64::recover(&samples[..]);
    for _ in 0..NN*2 {
        assert!(orig_mt.next_u64() == recovered_mt.next_u64());
    }
}