coreutils/test/factor.rs

/*
* This file is part of the uutils coreutils package.
*
* (c) kwantam <kwantam@gmail.com>
*
* For the full copyright and license information, please view the LICENSE file
* that was distributed with this source code.
*/

extern crate libc;
extern crate rand;

use rand::{weak_rng, Rng};
use rand::distributions::{IndependentSample, Range};
use sieve::Sieve;
use std::io::Write;
use std::process::{Command, Stdio};

#[path="../src/factor/sieve.rs"]
mod sieve;

const NUM_PRIMES: usize = 10000;
const LOG_PRIMES: f64 = 14.0;   // ceil(log2(NUM_PRIMES))

const NUM_TESTS: usize = 1000;
const PROGNAME: &'static str = "./factor";

#[test]
fn test_random() {
    let primes = Sieve::primes().take(NUM_PRIMES).collect::<Vec<u64>>();

    let mut rng = weak_rng();
    let mut rand_gt = move |min: u64| {
        let mut product = 1u64;
        let mut factors = Vec::new();
        while product < min {
            // log distribution---higher probability for lower numbers
            let mut factor;
            loop {
                let next = rng.gen_range(0f64, LOG_PRIMES).exp2().floor() as usize;
                if next < NUM_PRIMES {
                    factor = primes[next];
                    break;
                }
            }
            let factor = factor;

            match product.checked_mul(factor) {
                Some(p) => {
                    product = p;
                    factors.push(factor);
                },
                None => break,
            };
        }

        factors.sort();
        (product, factors)
    };

    // build an input and expected output string from factor
    let mut instring = String::new();
    let mut outstring = String::new();
    for _ in 0..NUM_TESTS {
        let (product, factors) = rand_gt(1 << 63);
        instring.push_str(&(format!("{} ", product))[..]);

        outstring.push_str(&(format!("{}:", product))[..]);
        for factor in factors {
            outstring.push_str(&(format!(" {}", factor))[..]);
        }
        outstring.push_str("\n");
    }

    run(instring.as_bytes(), outstring.as_bytes());
}

#[test]
fn test_random_big() {
    let mut rng = weak_rng();
    let bitrange_1 = Range::new(14usize, 51);
    let mut rand_64 = move || {
        // first, choose a random number of bits for the first factor
        let f_bit_1 = bitrange_1.ind_sample(&mut rng);
        // how many more bits do we need?
        let rem = 64 - f_bit_1;

        // we will have a number of additional factors equal to nfacts + 1
        // where nfacts is in [0, floor(rem/14) )  NOTE half-open interval
        // Each prime factor is at least 14 bits, hence floor(rem/14)
        let nfacts = Range::new(0usize, rem / 14).ind_sample(&mut rng);
        // we have to distribute extrabits among the (nfacts + 1) values
        let extrabits = rem - (nfacts + 1) * 14;
        // (remember, a Range is a half-open interval)
        let extrarange = Range::new(0usize, extrabits + 1);

        // to generate an even split of this range, generate n-1 random elements
        // in the range, add the desired total value to the end, sort this list,
        // and then compute the sequential differences.
        let mut f_bits = Vec::new();
        for _ in 0..nfacts {
            f_bits.push(extrarange.ind_sample(&mut rng));
        }
        f_bits.push(extrabits);
        f_bits.sort();

        // compute sequential differences here. We leave off the +14 bits
        // so we can just index PRIMES_BY_BITS
        let mut f_bits = f_bits.iter().scan(0, |st,&x| {
            let ret = x - *st; // + 14 would give actual number of bits
            *st = x;
            Some(ret)
        }).collect::<Vec<usize>>();
        // finally, add f_bit_1 in there
        f_bits.push(f_bit_1 - 14);  // index of f_bit_1 in PRIMES_BY_BITS
        let f_bits = f_bits;

        let mut nbits = 0;
        let mut product = 1u64;
        let mut factors = Vec::new();
        for bit in f_bits {
            assert!(bit < 37);
            nbits += 14 + bit;
            let elm = Range::new(0, PRIMES_BY_BITS[bit].len()).ind_sample(&mut rng);
            let factor = PRIMES_BY_BITS[bit][elm];
            factors.push(factor);
            product *= factor;
        }
        assert_eq!(nbits, 64);

        factors.sort();
        (product, factors)
    };

    let mut instring = String::new();
    let mut outstring = String::new();
    for _ in 0..NUM_TESTS {
        let (product, factors) = rand_64();
        instring.push_str(&(format!("{} ", product))[..]);

        outstring.push_str(&(format!("{}:", product))[..]);
        for factor in factors {
            outstring.push_str(&(format!(" {}", factor))[..]);
        }
        outstring.push_str("\n");
    }

    run(instring.as_bytes(), outstring.as_bytes());
}

#[test]
fn test_big_primes() {
    let mut instring = String::new();
    let mut outstring = String::new();
    for prime in PRIMES64 {
        instring.push_str(&(format!("{} ", prime))[..]);
        outstring.push_str(&(format!("{0}: {0}\n", prime))[..]);
    }

    run(instring.as_bytes(), outstring.as_bytes());
}

fn run(instring: &[u8], outstring: &[u8]) {
    // now run factor
    let mut process = Command::new(PROGNAME)
                                   .stdin(Stdio::piped())
                                   .stdout(Stdio::piped())
                                   .spawn()
                                   .unwrap_or_else(|e| panic!("{}", e));

    process.stdin.take().unwrap_or_else(|| panic!("Could not take child process stdin"))
        .write_all(instring).unwrap_or_else(|e| panic!("{}", e));

    let output = process.wait_with_output().unwrap_or_else(|e| panic!("{}", e));
    assert_eq!(&output.stdout[..], outstring);
}

const PRIMES_BY_BITS: &'static [&'static [u64]] = &[PRIMES14, PRIMES15, PRIMES16, PRIMES17,
    PRIMES18, PRIMES19, PRIMES20, PRIMES21, PRIMES22, PRIMES23, PRIMES24, PRIMES25, PRIMES26,
    PRIMES27, PRIMES28, PRIMES29, PRIMES30, PRIMES31, PRIMES32, PRIMES33, PRIMES34, PRIMES35,
    PRIMES36, PRIMES37, PRIMES38, PRIMES39, PRIMES40, PRIMES41, PRIMES42, PRIMES43, PRIMES44,
    PRIMES45, PRIMES46, PRIMES47, PRIMES48, PRIMES49, PRIMES50,
];

const PRIMES64: &'static [u64] = &[18446744073709551557, 18446744073709551533,
    18446744073709551521, 18446744073709551437, 18446744073709551427, 18446744073709551359,
    18446744073709551337, 18446744073709551293, 18446744073709551263, 18446744073709551253,
    18446744073709551191, 18446744073709551163, 18446744073709551113, 18446744073709550873,
    18446744073709550791, 18446744073709550773, 18446744073709550771, 18446744073709550719,
    18446744073709550717, 18446744073709550681, 18446744073709550671, 18446744073709550593,
    18446744073709550591, 18446744073709550539, 18446744073709550537, 18446744073709550381,
    18446744073709550341, 18446744073709550293, 18446744073709550237, 18446744073709550147,
    18446744073709550141, 18446744073709550129, 18446744073709550111, 18446744073709550099,
    18446744073709550047, 18446744073709550033, 18446744073709550009, 18446744073709549951,
    18446744073709549861, 18446744073709549817, 18446744073709549811, 18446744073709549777,
    18446744073709549757, 18446744073709549733, 18446744073709549667, 18446744073709549621,
    18446744073709549613, 18446744073709549583, 18446744073709549571,
];

const PRIMES14: &'static [u64] = &[16381, 16369, 16363, 16361, 16349, 16339, 16333, 16319,
    16301, 16273, 16267, 16253, 16249, 16231, 16229, 16223, 16217, 16193, 16189, 16187, 16183,
    16141, 16139, 16127, 16111, 16103, 16097, 16091, 16087, 16073, 16069, 16067, 16063, 16061,
    16057, 16033, 16007, 16001, 15991, 15973, 15971, 15959, 15937, 15923, 15919, 15913, 15907,
    15901, 15889, 15887, 15881, 15877, 15859, 15823, 15817, 15809, 15803, 15797, 15791, 15787,
    15773, 15767, 15761, 15749, 15739, 15737, 15733, 15731, 15727, 15683, 15679, 15671, 15667,
    15661, 15649, 15647, 15643, 15641, 15629, 15619, 15607, 15601, 15583, 15581, 15569, 15559,
    15551, 15541, 15527, 15511, 15497, 15493, 15473, 15467, 15461, 15451, 15443, 15439, 15427,
    15413, 15401, 15391, 15383, 15377, 15373,
];

const PRIMES15: &'static [u64] = &[32749, 32719, 32717, 32713, 32707, 32693, 32687, 32653,
    32647, 32633, 32621, 32611, 32609, 32603, 32587, 32579, 32573, 32569, 32563, 32561, 32537,
    32533, 32531, 32507, 32503, 32497, 32491, 32479, 32467, 32443, 32441, 32429, 32423, 32413,
    32411, 32401, 32381, 32377, 32371, 32369, 32363, 32359, 32353, 32341, 32327, 32323, 32321,
    32309, 32303, 32299, 32297, 32261, 32257, 32251, 32237, 32233, 32213, 32203, 32191, 32189,
    32183, 32173, 32159, 32143, 32141, 32119, 32117, 32099, 32089, 32083, 32077, 32069, 32063,
    32059, 32057, 32051, 32029, 32027, 32009, 32003, 31991, 31981, 31973, 31963, 31957, 31907,
    31891, 31883, 31873, 31859, 31849, 31847, 31817, 31799, 31793, 31771, 31769, 31751,
];

const PRIMES16: &'static [u64] = &[65521, 65519, 65497, 65479, 65449, 65447, 65437, 65423,
    65419, 65413, 65407, 65393, 65381, 65371, 65357, 65353, 65327, 65323, 65309, 65293, 65287,
    65269, 65267, 65257, 65239, 65213, 65203, 65183, 65179, 65173, 65171, 65167, 65147, 65141,
    65129, 65123, 65119, 65111, 65101, 65099, 65089, 65071, 65063, 65053, 65033, 65029, 65027,
    65011, 65003, 64997, 64969, 64951, 64937, 64927, 64921, 64919, 64901, 64891, 64879, 64877,
    64871, 64853, 64849, 64817, 64811, 64793, 64783, 64781, 64763, 64747, 64717, 64709, 64693,
    64679, 64667, 64663, 64661, 64633, 64627, 64621, 64613, 64609, 64601, 64591, 64579, 64577,
    64567, 64553,
];

const PRIMES17: &'static [u64] = &[131071, 131063, 131059, 131041, 131023, 131011, 131009,
    130987, 130981, 130973, 130969, 130957, 130927, 130873, 130859, 130843, 130841, 130829,
    130817, 130811, 130807, 130787, 130783, 130769, 130729, 130699, 130693, 130687, 130681,
    130657, 130651, 130649, 130643, 130639, 130633, 130631, 130621, 130619, 130589, 130579,
    130553, 130547, 130531, 130523, 130517, 130513, 130489, 130483, 130477, 130469, 130457,
    130447, 130439, 130423, 130411, 130409, 130399, 130379, 130369, 130367, 130363, 130349,
    130343, 130337, 130307, 130303, 130279, 130267, 130261, 130259, 130253, 130241, 130223,
    130211, 130201, 130199, 130183, 130171, 130147, 130127, 130121, 130099, 130087, 130079,
    130073, 130069, 130057, 130051,
];

const PRIMES18: &'static [u64] = &[262139, 262133, 262127, 262121, 262111, 262109, 262103,
    262079, 262069, 262051, 262049, 262027, 262007, 261983, 261977, 261973, 261971, 261959,
    261917, 261887, 261881, 261847, 261823, 261799, 261791, 261787, 261773, 261761, 261757,
    261739, 261721, 261713, 261707, 261697, 261673, 261643, 261641, 261637, 261631, 261619,
    261601, 261593, 261587, 261581, 261577, 261563, 261557, 261529, 261523, 261509, 261467,
    261463, 261451, 261439, 261433, 261431, 261427, 261407, 261389, 261379, 261353, 261347,
    261337, 261329, 261323, 261301, 261281, 261271, 261251, 261241, 261229, 261223, 261169,
    261167, 261127,
];

const PRIMES19: &'static [u64] = &[524287, 524269, 524261, 524257, 524243, 524231, 524221,
    524219, 524203, 524201, 524197, 524189, 524171, 524149, 524123, 524119, 524113, 524099,
    524087, 524081, 524071, 524063, 524057, 524053, 524047, 523997, 523987, 523969, 523949,
    523937, 523927, 523907, 523903, 523877, 523867, 523847, 523829, 523801, 523793, 523777,
    523771, 523763, 523759, 523741, 523729, 523717, 523681, 523673, 523669, 523667, 523657,
    523639, 523637, 523631, 523603, 523597, 523577, 523573, 523571, 523553, 523543, 523541,
    523519, 523511, 523493, 523489, 523487, 523463, 523459, 523433, 523427, 523417, 523403,
    523387, 523357, 523351, 523349, 523333, 523307, 523297,
];

const PRIMES20: &'static [u64] = &[1048573, 1048571, 1048559, 1048549, 1048517, 1048507,
    1048447, 1048433, 1048423, 1048391, 1048387, 1048367, 1048361, 1048357, 1048343, 1048309,
    1048291, 1048273, 1048261, 1048219, 1048217, 1048213, 1048193, 1048189, 1048139, 1048129,
    1048127, 1048123, 1048063, 1048051, 1048049, 1048043, 1048027, 1048013, 1048009, 1048007,
    1047997, 1047989, 1047979, 1047971, 1047961, 1047941, 1047929, 1047923, 1047887, 1047883,
    1047881, 1047859, 1047841, 1047833, 1047821, 1047779, 1047773, 1047763, 1047751, 1047737,
    1047721, 1047713, 1047703, 1047701, 1047691, 1047689, 1047671, 1047667, 1047653, 1047649,
    1047647, 1047589, 1047587, 1047559,
];

const PRIMES21: &'static [u64] = &[2097143, 2097133, 2097131, 2097097, 2097091, 2097083,
    2097047, 2097041, 2097031, 2097023, 2097013, 2096993, 2096987, 2096971, 2096959, 2096957,
    2096947, 2096923, 2096911, 2096909, 2096893, 2096881, 2096873, 2096867, 2096851, 2096837,
    2096807, 2096791, 2096789, 2096777, 2096761, 2096741, 2096737, 2096713, 2096693, 2096687,
    2096681, 2096639, 2096629, 2096621, 2096599, 2096597, 2096569, 2096539, 2096533, 2096483,
    2096449, 2096431, 2096429, 2096411, 2096407, 2096401, 2096399, 2096377, 2096357, 2096291,
    2096273, 2096261, 2096233, 2096231, 2096221, 2096209, 2096191, 2096183, 2096147,
];

const PRIMES22: &'static [u64] = &[4194301, 4194287, 4194277, 4194271, 4194247, 4194217,
    4194199, 4194191, 4194187, 4194181, 4194173, 4194167, 4194143, 4194137, 4194131, 4194107,
    4194103, 4194023, 4194011, 4194007, 4193977, 4193971, 4193963, 4193957, 4193939, 4193929,
    4193909, 4193869, 4193807, 4193803, 4193801, 4193789, 4193759, 4193753, 4193743, 4193701,
    4193663, 4193633, 4193573, 4193569, 4193551, 4193549, 4193531, 4193513, 4193507, 4193459,
    4193447, 4193443, 4193417, 4193411, 4193393, 4193389, 4193381, 4193377, 4193369, 4193359,
    4193353, 4193327, 4193309, 4193303, 4193297,
];

const PRIMES23: &'static [u64] = &[8388593, 8388587, 8388581, 8388571, 8388547, 8388539,
    8388473, 8388461, 8388451, 8388449, 8388439, 8388427, 8388421, 8388409, 8388377, 8388371,
    8388319, 8388301, 8388287, 8388283, 8388277, 8388239, 8388209, 8388187, 8388113, 8388109,
    8388091, 8388071, 8388059, 8388019, 8388013, 8387999, 8387993, 8387959, 8387957, 8387947,
    8387933, 8387921, 8387917, 8387891, 8387879, 8387867, 8387861, 8387857, 8387839, 8387831,
    8387809, 8387807, 8387741, 8387737, 8387723, 8387707, 8387671, 8387611, 8387609, 8387591,
];

const PRIMES24: &'static [u64] = &[16777213, 16777199, 16777183, 16777153, 16777141, 16777139,
    16777127, 16777121, 16777099, 16777049, 16777027, 16776989, 16776973, 16776971, 16776967,
    16776961, 16776941, 16776937, 16776931, 16776919, 16776901, 16776899, 16776869, 16776857,
    16776839, 16776833, 16776817, 16776763, 16776731, 16776719, 16776713, 16776691, 16776689,
    16776679, 16776659, 16776631, 16776623, 16776619, 16776607, 16776593, 16776581, 16776547,
    16776521, 16776491, 16776481, 16776469, 16776451, 16776401, 16776391, 16776379, 16776371,
    16776367, 16776343, 16776337, 16776317, 16776313, 16776289, 16776217, 16776211,
];

const PRIMES25: &'static [u64] = &[33554393, 33554383, 33554371, 33554347, 33554341, 33554317,
    33554291, 33554273, 33554267, 33554249, 33554239, 33554221, 33554201, 33554167, 33554159,
    33554137, 33554123, 33554093, 33554083, 33554077, 33554051, 33554021, 33554011, 33554009,
    33553999, 33553991, 33553969, 33553967, 33553909, 33553901, 33553879, 33553837, 33553799,
    33553787, 33553771, 33553769, 33553759, 33553747, 33553739, 33553727, 33553697, 33553693,
    33553679, 33553661, 33553657, 33553651, 33553649, 33553633, 33553613, 33553607, 33553577,
    33553549, 33553547, 33553537, 33553519, 33553517, 33553511, 33553489, 33553463, 33553451,
    33553417,
];

const PRIMES26: &'static [u64] = &[67108859, 67108837, 67108819, 67108777, 67108763,
    67108757, 67108753, 67108747, 67108739, 67108729, 67108721, 67108709, 67108693, 67108669,
    67108667, 67108661, 67108649, 67108633, 67108597, 67108579, 67108529, 67108511, 67108507,
    67108493, 67108471, 67108463, 67108453, 67108439, 67108387, 67108373, 67108369, 67108351,
    67108331, 67108313, 67108303, 67108289, 67108271, 67108219, 67108207, 67108201, 67108199,
    67108187, 67108183, 67108177, 67108127, 67108109, 67108081, 67108049, 67108039, 67108037,
    67108033, 67108009, 67108007, 67108003, 67107983, 67107977, 67107967, 67107941, 67107919,
    67107913, 67107883, 67107881, 67107871, 67107863,
];

const PRIMES27: &'static [u64] = &[134217689, 134217649, 134217617, 134217613, 134217593,
    134217541, 134217529, 134217509, 134217497, 134217493, 134217487, 134217467, 134217439,
    134217437, 134217409, 134217403, 134217401, 134217367, 134217361, 134217353, 134217323,
    134217301, 134217277, 134217257, 134217247, 134217221, 134217199, 134217173, 134217163,
    134217157, 134217131, 134217103, 134217089, 134217079, 134217049, 134217047, 134217043,
    134217001, 134216987, 134216947, 134216939, 134216933, 134216911, 134216899, 134216881,
    134216869, 134216867, 134216861, 134216837, 134216827, 134216807, 134216801, 134216791,
    134216783, 134216777, 134216759, 134216737, 134216729,
];

const PRIMES28: &'static [u64] = &[268435399, 268435367, 268435361, 268435337, 268435331,
    268435313, 268435291, 268435273, 268435243, 268435183, 268435171, 268435157, 268435147,
    268435133, 268435129, 268435121, 268435109, 268435091, 268435067, 268435043, 268435039,
    268435033, 268435019, 268435009, 268435007, 268434997, 268434979, 268434977, 268434961,
    268434949, 268434941, 268434937, 268434857, 268434841, 268434827, 268434821, 268434787,
    268434781, 268434779, 268434773, 268434731, 268434721, 268434713, 268434707, 268434703,
    268434697, 268434659, 268434623, 268434619, 268434581, 268434577, 268434563, 268434557,
    268434547, 268434511, 268434499, 268434479, 268434461,
];

const PRIMES29: &'static [u64] = &[536870909, 536870879, 536870869, 536870849, 536870839,
    536870837, 536870819, 536870813, 536870791, 536870779, 536870767, 536870743, 536870729,
    536870723, 536870717, 536870701, 536870683, 536870657, 536870641, 536870627, 536870611,
    536870603, 536870599, 536870573, 536870569, 536870563, 536870561, 536870513, 536870501,
    536870497, 536870473, 536870401, 536870363, 536870317, 536870303, 536870297, 536870273,
    536870267, 536870239, 536870233, 536870219, 536870171, 536870167, 536870153, 536870123,
    536870063, 536870057, 536870041, 536870027, 536869999, 536869951, 536869943, 536869937,
    536869919, 536869901, 536869891,
];

const PRIMES30: &'static [u64] = &[1073741789, 1073741783, 1073741741, 1073741723, 1073741719,
    1073741717, 1073741689, 1073741671, 1073741663, 1073741651, 1073741621, 1073741567, 1073741561,
    1073741527, 1073741503, 1073741477, 1073741467, 1073741441, 1073741419, 1073741399, 1073741387,
    1073741381, 1073741371, 1073741329, 1073741311, 1073741309, 1073741287, 1073741237, 1073741213,
    1073741197, 1073741189, 1073741173, 1073741101, 1073741077, 1073741047, 1073740963, 1073740951,
    1073740933, 1073740909, 1073740879, 1073740853, 1073740847, 1073740819, 1073740807,
];

const PRIMES31: &'static [u64] = &[2147483647, 2147483629, 2147483587, 2147483579, 2147483563,
    2147483549, 2147483543, 2147483497, 2147483489, 2147483477, 2147483423, 2147483399, 2147483353,
    2147483323, 2147483269, 2147483249, 2147483237, 2147483179, 2147483171, 2147483137, 2147483123,
    2147483077, 2147483069, 2147483059, 2147483053, 2147483033, 2147483029, 2147482951, 2147482949,
    2147482943, 2147482937, 2147482921, 2147482877, 2147482873, 2147482867, 2147482859, 2147482819,
    2147482817, 2147482811, 2147482801, 2147482763, 2147482739, 2147482697, 2147482693, 2147482681,
    2147482663, 2147482661,
];

const PRIMES32: &'static [u64] = &[4294967291, 4294967279, 4294967231, 4294967197, 4294967189,
    4294967161, 4294967143, 4294967111, 4294967087, 4294967029, 4294966997, 4294966981, 4294966943,
    4294966927, 4294966909, 4294966877, 4294966829, 4294966813, 4294966769, 4294966667, 4294966661,
    4294966657, 4294966651, 4294966639, 4294966619, 4294966591, 4294966583, 4294966553, 4294966477,
    4294966447, 4294966441, 4294966427, 4294966373, 4294966367, 4294966337, 4294966297,
];

const PRIMES33: &'static [u64] = &[8589934583, 8589934567, 8589934543, 8589934513, 8589934487,
    8589934307, 8589934291, 8589934289, 8589934271, 8589934237, 8589934211, 8589934207, 8589934201,
    8589934187, 8589934151, 8589934141, 8589934139, 8589934117, 8589934103, 8589934099, 8589934091,
    8589934069, 8589934049, 8589934027, 8589934007, 8589933973, 8589933971, 8589933967, 8589933931,
    8589933917, 8589933907, 8589933853, 8589933827, 8589933823, 8589933787, 8589933773, 8589933733,
    8589933731, 8589933721, 8589933683, 8589933647, 8589933641, 8589933637, 8589933631, 8589933629,
    8589933619, 8589933601, 8589933581,
];

const PRIMES34: &'static [u64] = &[17179869143, 17179869107, 17179869071, 17179869053,
    17179869041, 17179869019, 17179868999, 17179868977, 17179868957, 17179868903, 17179868899,
    17179868887, 17179868879, 17179868873, 17179868869, 17179868861, 17179868843, 17179868833,
    17179868809, 17179868807, 17179868777, 17179868759, 17179868729, 17179868711, 17179868683,
    17179868681, 17179868597, 17179868549, 17179868543, 17179868521, 17179868513, 17179868479,
    17179868443, 17179868437, 17179868429, 17179868383, 17179868369, 17179868357, 17179868353,
    17179868351, 17179868333, 17179868317, 17179868309, 17179868297, 17179868287, 17179868249,
    17179868243, 17179868183,
];

const PRIMES35: &'static [u64] = &[34359738337, 34359738319, 34359738307, 34359738299,
    34359738289, 34359738247, 34359738227, 34359738121, 34359738059, 34359738043, 34359738011,
    34359737917, 34359737869, 34359737849, 34359737837, 34359737821, 34359737813, 34359737791,
    34359737777, 34359737771, 34359737717, 34359737591, 34359737567, 34359737549, 34359737519,
    34359737497, 34359737479, 34359737407, 34359737393, 34359737371,
];

const PRIMES36: &'static [u64] = &[68719476731, 68719476719, 68719476713, 68719476671,
    68719476619, 68719476599, 68719476577, 68719476563, 68719476547, 68719476503, 68719476493,
    68719476479, 68719476433, 68719476407, 68719476391, 68719476389, 68719476377, 68719476361,
    68719476323, 68719476307, 68719476281, 68719476271, 68719476257, 68719476247, 68719476209,
    68719476197, 68719476181, 68719476169, 68719476157, 68719476149, 68719476109, 68719476053,
    68719476047, 68719476019, 68719475977, 68719475947, 68719475933, 68719475911, 68719475893,
    68719475879, 68719475837, 68719475827, 68719475809, 68719475791, 68719475779, 68719475771,
    68719475767, 68719475731, 68719475729,
];

const PRIMES37: &'static [u64] = &[137438953447, 137438953441, 137438953427, 137438953403,
    137438953349, 137438953331, 137438953273, 137438953271, 137438953121, 137438953097,
    137438953037, 137438953009, 137438952953, 137438952901, 137438952887, 137438952869,
    137438952853, 137438952731, 137438952683, 137438952611, 137438952529, 137438952503,
    137438952491,
];

const PRIMES38: &'static [u64] = &[274877906899, 274877906857, 274877906837, 274877906813,
    274877906791, 274877906759, 274877906753, 274877906717, 274877906713, 274877906687,
    274877906647, 274877906629, 274877906627, 274877906573, 274877906543, 274877906491,
    274877906477, 274877906473, 274877906431, 274877906419, 274877906341, 274877906333,
    274877906327, 274877906321, 274877906309, 274877906267, 274877906243, 274877906213,
    274877906209, 274877906203, 274877906179, 274877906167, 274877906119, 274877906063,
    274877906053, 274877906021, 274877905931,
];

const PRIMES39: &'static [u64] = &[549755813881, 549755813869, 549755813821, 549755813797,
    549755813753, 549755813723, 549755813669, 549755813657, 549755813647, 549755813587,
    549755813561, 549755813513, 549755813507, 549755813461, 549755813417, 549755813401,
    549755813371, 549755813359, 549755813357, 549755813351, 549755813339, 549755813317,
    549755813311, 549755813281, 549755813239, 549755813231, 549755813213, 549755813207,
    549755813197, 549755813183, 549755813161, 549755813149, 549755813147, 549755813143,
    549755813141, 549755813059, 549755813027, 549755813003, 549755812951, 549755812937,
    549755812933, 549755812889, 549755812867,
];

const PRIMES40: &'static [u64] = &[1099511627689, 1099511627609, 1099511627581, 1099511627573,
    1099511627563, 1099511627491, 1099511627483, 1099511627477, 1099511627387, 1099511627339,
    1099511627321, 1099511627309, 1099511627297, 1099511627293, 1099511627261, 1099511627213,
    1099511627191, 1099511627177, 1099511627173, 1099511627143, 1099511627089, 1099511626987,
    1099511626949, 1099511626937, 1099511626793, 1099511626781, 1099511626771,
];

const PRIMES41: &'static [u64] = &[2199023255531, 2199023255521, 2199023255497, 2199023255489,
    2199023255479, 2199023255477, 2199023255461, 2199023255441, 2199023255419, 2199023255413,
    2199023255357, 2199023255327, 2199023255291, 2199023255279, 2199023255267, 2199023255243,
    2199023255203, 2199023255171, 2199023255137, 2199023255101, 2199023255087, 2199023255081,
    2199023255069, 2199023255027, 2199023255021, 2199023254979, 2199023254933, 2199023254913,
    2199023254907, 2199023254903, 2199023254843, 2199023254787, 2199023254699, 2199023254693,
    2199023254657, 2199023254567,
];

const PRIMES42: &'static [u64] = &[4398046511093, 4398046511087, 4398046511071, 4398046511051,
    4398046511039, 4398046510961, 4398046510943, 4398046510939, 4398046510889, 4398046510877,
    4398046510829, 4398046510787, 4398046510771, 4398046510751, 4398046510733, 4398046510721,
    4398046510643, 4398046510639, 4398046510597, 4398046510577, 4398046510547, 4398046510531,
    4398046510463, 4398046510397, 4398046510391, 4398046510379, 4398046510357, 4398046510331,
    4398046510327, 4398046510313, 4398046510283, 4398046510279, 4398046510217, 4398046510141,
    4398046510133, 4398046510103, 4398046510093,
];

const PRIMES43: &'static [u64] = &[8796093022151, 8796093022141, 8796093022091, 8796093022033,
    8796093021953, 8796093021941, 8796093021917, 8796093021899, 8796093021889, 8796093021839,
    8796093021803, 8796093021791, 8796093021769, 8796093021763, 8796093021743, 8796093021671,
    8796093021607, 8796093021587, 8796093021533, 8796093021523, 8796093021517, 8796093021493,
    8796093021467, 8796093021461, 8796093021449, 8796093021409, 8796093021407, 8796093021371,
    8796093021347, 8796093021337, 8796093021281, 8796093021269,
];

const PRIMES44: &'static [u64] = &[17592186044399, 17592186044299, 17592186044297, 17592186044287,
    17592186044273, 17592186044267, 17592186044129, 17592186044089, 17592186044057, 17592186044039,
    17592186043987, 17592186043921, 17592186043889, 17592186043877, 17592186043841, 17592186043829,
    17592186043819, 17592186043813, 17592186043807, 17592186043741, 17592186043693, 17592186043667,
    17592186043631, 17592186043591, 17592186043577, 17592186043547, 17592186043483, 17592186043451,
    17592186043433, 17592186043409,
];

const PRIMES45: &'static [u64] = &[35184372088777, 35184372088763, 35184372088751, 35184372088739,
    35184372088711, 35184372088699, 35184372088693, 35184372088673, 35184372088639, 35184372088603,
    35184372088571, 35184372088517, 35184372088493, 35184372088471, 35184372088403, 35184372088391,
    35184372088379, 35184372088363, 35184372088321, 35184372088319, 35184372088279, 35184372088259,
    35184372088249, 35184372088241, 35184372088223, 35184372088183, 35184372088097, 35184372088081,
    35184372088079, 35184372088051, 35184372088043, 35184372088039, 35184372087937, 35184372087929,
    35184372087923, 35184372087881, 35184372087877, 35184372087869,
];

const PRIMES46: &'static [u64] = &[70368744177643, 70368744177607, 70368744177601, 70368744177587,
    70368744177497, 70368744177467, 70368744177427, 70368744177377, 70368744177359, 70368744177353,
    70368744177331, 70368744177289, 70368744177283, 70368744177271, 70368744177257, 70368744177227,
    70368744177167, 70368744177113, 70368744177029, 70368744176959, 70368744176921, 70368744176909,
    70368744176879, 70368744176867, 70368744176833, 70368744176827, 70368744176807, 70368744176779,
    70368744176777, 70368744176729, 70368744176719, 70368744176711,
];

const PRIMES47: &'static [u64] = &[140737488355213, 140737488355201, 140737488355181,
    140737488355049, 140737488355031, 140737488354989, 140737488354893, 140737488354787,
    140737488354709, 140737488354679, 140737488354613, 140737488354557, 140737488354511,
    140737488354431, 140737488354413, 140737488354409, 140737488354373, 140737488354347,
    140737488354329,
];

const PRIMES48: &'static [u64] = &[281474976710597, 281474976710591, 281474976710567,
    281474976710563, 281474976710509, 281474976710491, 281474976710467, 281474976710423,
    281474976710413, 281474976710399, 281474976710339, 281474976710327, 281474976710287,
    281474976710197, 281474976710143, 281474976710131, 281474976710129, 281474976710107,
    281474976710089, 281474976710087, 281474976710029, 281474976709987, 281474976709891,
    281474976709859, 281474976709831, 281474976709757, 281474976709741, 281474976709711,
    281474976709649, 281474976709637,
];

const PRIMES49: &'static [u64] = &[562949953421231, 562949953421201, 562949953421189,
    562949953421173, 562949953421131, 562949953421111, 562949953421099, 562949953421047,
    562949953421029, 562949953420973, 562949953420871, 562949953420867, 562949953420837,
    562949953420793, 562949953420747, 562949953420741, 562949953420733, 562949953420727,
    562949953420609, 562949953420571, 562949953420559, 562949953420553, 562949953420523,
    562949953420507, 562949953420457, 562949953420403, 562949953420373, 562949953420369,
    562949953420343, 562949953420303, 562949953420297,
];

const PRIMES50: &'static [u64] = &[1125899906842597, 1125899906842589, 1125899906842573,
    1125899906842553, 1125899906842511, 1125899906842507, 1125899906842493, 1125899906842463,
    1125899906842429, 1125899906842391, 1125899906842357, 1125899906842283, 1125899906842273,
    1125899906842247, 1125899906842201, 1125899906842177, 1125899906842079, 1125899906842033,
    1125899906842021, 1125899906842013, 1125899906841973, 1125899906841971, 1125899906841959,
    1125899906841949, 1125899906841943, 1125899906841917, 1125899906841901, 1125899906841883,
    1125899906841859, 1125899906841811, 1125899906841803, 1125899906841751, 1125899906841713,
    1125899906841673, 1125899906841653, 1125899906841623, 1125899906841613,
];