auxilin-prediction/src/config.py

'''Set configurations for all the raw and processed datasets
Takes in path to raw data, tracks, processed df, and then later interim data
'''
import os
from os.path import join as oj
import numpy as np


# data paths
DIR_REPO = os.path.dirname(os.path.realpath(__file__)) # directory of the config file
DIR_TRACKS = oj(DIR_REPO, '..', 'data', 'tracks')
DIR_PROCESSED = oj(DIR_REPO, '..', 'data', 'processed')
DIR_INTERIM = oj(DIR_REPO, '..', 'data', 'interim')
DIR_FIGS = oj(DIR_REPO, '..', 'auxilin-paper', 'figs')
DIR_RESULTS = '/scratch/users/vision/chandan/abc'
FINAL_MODEL = oj(DIR_REPO, '..', 'models/dnn_full_long_normalized_across_track_1_feat_dynamin.pkl')
FINAL_MODEL_VPS = oj(DIR_REPO, '..', 'models/dnn_vps_fit_extended_lifetimes>3.pkl')

# raw data ################################################################################
data_dir_orig = '/scratch/users/vision/data/abc_data/auxilin_data_tracked'
d_new = '/scratch/users/vision/data/abc_data/auxilin_new_data/AI_Clathrin_molecularPrediction'
clath_aux_folder = 'CLTA-TagRFP+-+ EGFP-Aux1-A7D2+-+ EGFP-GAK-F6+-+  TIRF data'
dynamin_folder = 'CLTA-TagRFP EGFP-Aux1-GAK-F6 Dyn2-Halo-E1-JF646'
vps_folder = '/scratch/users/vision/data/abc_data/vps_snf'


# data splitting
s_orig = 'A7D2'
s_orig_gak = 'EGFP-GAK F6'
s_clath_aux_no_a7d2 = 'CLTA-TagRFP EGFP-Aux1 EGFP-GAK F6' # goes [1-11]
s_clath_aux = 'CLTA-TagRFP EGFP-Aux1-A7D2 EGFP-GAK-F6' # goes [1, 2] + [5-12] - 13 is there but missing some tifs
s_a8 = 'CLTA-TagRFP EGFP-GAK A8'
s_dynamin = '488-1.5mW 561-1.5mW 647-1.5mW Exp100ms Int1.5s'
s_clath_pi4p = 'cme'
s_ap2_pi4p = 'cme'

DSETS = {
    ############################################################################vps/snf
    
    # the most-relevant best dataset to use (make sure to use the key with wt - wild type and not mt - mutant)
    'vps4_snf7': {
        'data_dir': oj(vps_folder, 
                       '2016-07-10_mergedTracks_vps4_snf7.mat'),        
        'feature_selection': None,
        'train': ['train'],
        'test': ['test'],
    },
    
    # the second part of this key says to save the mutant key from this same file (in load_tracking.py)
    'vps4_snf7___key=mt': {
        'data_dir': oj(vps_folder, 
                       '2016-07-10_mergedTracks_vps4_snf7.mat'),        
        'feature_selection': None,
        'train': ['train'],
        'test': ['test'],
    },    
    
    # more complicated...
    'AnchorAway_vps4_snf7': {
        'data_dir': oj(vps_folder, 
                       '2016-09-01_mergedTracks_AnchorAway_vps4_snf7_RapExposure-004.mat'),        
        'feature_selection': None,
        'train': ['train'],
        'test': ['test'],
    },
    
    # only has 1 channel
    'vps24': {
        'data_dir': oj(vps_folder, 
                       '2016-07-10_mergedTracks_vps24.mat'),        
        'feature_selection': None,
        'train': ['train'],
        'test': ['test'],
    },    
    
    # only has mutant key
    'vps4_snf7vps2': {
        'data_dir': oj(vps_folder, 
                       '2016-10-26_mergedTracks_vps4_snf7vps2.mat'),        
        'feature_selection': None,
        'train': ['train'],
        'test': ['test'],
    },    
    
    ## mutant (EQ)
    'vps24_VPS4EQ': {
        'data_dir': oj(vps_folder, 
                       '2016-09-06_mergedTracks_vps24_VPS4EQ.mat'),        
        'feature_selection': None,
        'train': ['train'],
        'test': ['test'],
    },
    
    ## mutant (EQ)
    'vps4EQ_snf7-eGFP': {
        'data_dir': oj(vps_folder, 
                       '2016-09-06_mergedTracks_vps4EQ_snf7-eGFP.mat'),        
        'feature_selection': None,
        'train': ['train'],
        'test': ['test'],
    },

    # ignore for now
    # '2017-05-22_newExp_tracks_vps4_snf7_dpep-002.mat' 

    # clath/aux############################################################################
    'clath_aux+gak_a7d2': {  # this was the original data
        'data_dir': data_dir_orig,
        'feature_selection': np.array([f'{s_orig}/1']),
        'train': np.array([f'{s_orig}/{i}' for i in [1, 2, 3, 4, 5, 6]]),
        'test': np.array([f'{s_orig}/{i}' for i in [7, 8]])
    },
    'clath_aux+gak': {
        'data_dir': data_dir_orig,
        'feature_selection': np.array([f'{s_orig_gak}/1']),
        'train': np.array([f'{s_orig_gak}/{i}' for i in [1, 2, 3, 4, 5, 6]]),
        'test': np.array([f'{s_orig_gak}/{i}' for i in [7, 8]])
    },
    'clath_aux+gak_a7d2_new': {
        'data_dir': oj(d_new, clath_aux_folder),        
        'feature_selection': np.array([f'{s_clath_aux}/1']),
        'train': np.array([f'{s_clath_aux}/{i}' for i in [1, 2, 5, 6, 7, 8, 9, 10]]),
        'test': np.array([f'{s_clath_aux}/{i}' for i in [11, 12]])
    },
    'clath_aux_dynamin': {  # this is 3-channel data
        'data_dir': oj(d_new, dynamin_folder),        
        'feature_selection': None,
        'train': np.array([f'{s_dynamin}_{i}/1_1.5s' for i in [4, 5, 6]] + 
                          [f'{s_dynamin}_4_Pos{i}/1_1.5s' for i in np.arange(13) if i != 9]),
        'test': np.array([f'{s_dynamin}_7_Pos{i}/1_1.5s' for i in np.arange(6)]),
    },
    'clath_aux+gak_new': {
        'data_dir': oj(d_new, clath_aux_folder),
        'feature_selection': np.array([f'{s_clath_aux_no_a7d2}/1']),
        'train': np.array([f'{s_clath_aux_no_a7d2}/{i}' for i in [1, 2, 3, 4, 5, 6, 7, 9]]),
        'test': np.array([f'{s_clath_aux_no_a7d2}/{i}' for i in [10, 11]])
    }, 
    'clath_gak': {
        'data_dir': oj(d_new, 'CLTA-TagRFP+-+ EGFP-GAK-A8+-+  TIRF data'),
        'feature_selection': np.array([f'{s_a8}/1']),
        'train': np.array([f'{s_a8}/{i}' for i in [1, 2, 3, 4, 5, 6, 7]]),
        'test': np.array([f'{s_a8}/{i}' for i in [8, 9]])    
    },
    'clath_pi4p_notreatment': {
        'data_dir': oj(d_new, 'CLTA-TagRFP+-+ EGFP-DrrA-Aux1 PI4P probe TIRF data', 'notreatment'), 
        'feature_selection': np.array([f'{s_clath_pi4p}/2']),
        'train': np.array([f'{s_clath_pi4p}/{i}' for i in [2, 5]]),
        'test': np.array([f'{s_clath_pi4p}/{i}' for i in [9, 10]])  
    },
    'ap2_pi4p': {
        'data_dir': oj(d_new, 'AP2-TagRFP+-+ EGFP-DrrA-Aux1 PI4P probe TIRF data'),
        'feature_selection': np.array([f'{s_ap2_pi4p}/1']),
        'train': np.array([f'{s_ap2_pi4p}/{i}' for i in [1, 2, 4, 5, 6, 7, 8]]),
        'test': np.array([f'{s_ap2_pi4p}/{i}' for i in [9, 11]])  
    },
}


# gt labels for original data (by pid) ################################################################################

'''Note these labels constitute corrections to the y_consec_thresh defn (as of 02/09/2020)
'''
LABELS_ORIG = {
    'neg': [6725, 909, 5926, 983, 8224, 3363] +
           [221, 248, 255, 274, 291, 294, 350, 358, 364, 384, 434, 481, 518,
            535, 588, 590, 602, 604, 638, 668, 675, 708, 752, 778, 808, 827,
            852, 867, 884, 888, 897, 905, 952, 953, 967, 984, 987, 997, 1002,
            1003, 1015, 1030, 1045, 1051, 1055, 1070, 1091, 1100, 1144, 1178,
            1226, 1246, 1258, 1260, 1276, 1307, 1317, 1370, 1394, 1399, 1400,
            1423, 1432, 1451, 1613, 2115, 2143, 2286, 2330, 2369, 2381, 2399,
            2416, 2435, 2612, 2613, 2631, 2667, 2679, 2705, 2714, 2750, 2756,
            2857, 2962, 3002, 3017, 3032, 3061, 3128, 3181, 3216, 3801, 3819,
            3916, 3945, 3954, 3967, 3970, 4042, 4108, 4120, 4141, 4142, 4155,
            4162, 4165, 4172, 4221, 4279, 4290, 4297, 4328, 4333, 4341, 4384,
            4422, 4428, 4468, 4547, 4558, 4573, 4586, 4616, 4621, 4674, 4875,
            4936, 4941, 4947, 5151, 5391, 5602, 5693, 5706, 5707, 5946, 5973,
            6007, 6072, 6147, 6157, 6253, 6489, 6618, 7340, 7363, 7413, 7423,
            7424, 7455, 7483, 7485, 7512, 7565, 7568, 7641, 7647, 7690, 7728,
            7872, 7877, 7883, 7892, 7942, 7969, 7985, 8005, 8018, 8020, 8024,
            8111, 8120, 8166, 8255, 8264, 8375, 8462, 8484, 8512, 8532, 8829,
            8833, 8918, 8944] +
           [3964, 10054, 735, 846, 1362, 3823, 5389, 8834, 9565, 10882, 11166, 718] +
           [3964, 718, 735, 846, 1362, 3823, 5389, 8834, 9565, 10882, 11166, 10054],
    # note: first list were flipped, second list were already correct
    'pos': [3982, 8243, 777, 3940, 7559, 2455, 4748, 633, 2177, 1205, 603, 7972, 8458, 3041, 924, 8786, 4116, 885,
            6298, 4658, 7889, 982, 829, 1210, 3054, 504, 1164, 347, 627, 1470, 2662, 2813, 422, 8400, 7474, 1273,
            6365, 1559, 4348, 1156, 6250, 4864, 639, 930, 5424, 7818, 8463, 4358, 7656, 843, 890, 4373, 2737, 7524,
            2590, 3804, 7667, 2148, 8585, 2919, 5712, 4440, 1440, 4699, 1089, 3004, 3126, 2869, 4183, 7335, 3166,
            8461, 2180, 849, 6458, 4575, 4091, 3966, 4725, 2514, 7626, 3055, 4200, 6429, 1220, 4472, 8559, 412, 903,
            5440, 1084, 2136, 6833, 1189, 7521, 8141, 7939, 8421, 944, 1264, 298, 6600, 1309, 3043, 243, 4161, 6813,
            5464] +
           [238, 251, 389, 524, 556, 758, 759, 830, 1213, 1255, 1290, 1422, 1463, 1484, 1504, 2016, 2046,
            2061, 2077, 2083, 2106, 2112, 2116, 2124, 2129, 2135, 2195, 2209, 2217, 2234, 2273, 2291, 2313,
            2338, 2353, 2402, 2460, 2651, 2658, 2703, 2791, 2805, 2848, 3101, 3138, 3142, 3150, 3170, 3343,
            3367, 3918, 3946, 4223, 4386, 4430, 4583, 4608, 4620, 4724, 5366, 5383, 5384, 5407, 5412, 5415,
            5422, 5442, 5447, 5449, 5478, 5492, 5503, 5506, 5516, 5548, 5550, 5558, 5589, 5634, 5694, 5708,
            5728, 5760, 5780, 5787, 5788, 5800, 5811, 5813, 5814, 5879, 5882, 5885, 5888, 5891, 5899, 5911,
            5912, 5950, 5951, 5953, 5957, 5960, 5986, 5988, 6000, 6011, 6012, 6020, 6021, 6032, 6049, 6053,
            6065, 6096, 6106, 6113, 6118, 6123, 6152, 6155, 6202, 6237, 6246, 6248, 6263, 6266, 6272, 6273,
            6302, 6305, 6321, 6325, 6327, 6363, 6368, 6398, 6407, 6410, 6423, 6424, 6431, 6444, 6449, 6461,
            6462, 6482, 6490, 6517, 6518, 6526, 6568, 6586, 6594, 6601, 6608, 6640, 6656, 6662, 6683, 6684,
            6693, 6703, 6771, 6774, 6801, 6802, 6823, 6851, 7348, 7352, 7448, 7470, 7496, 7511, 7596, 7720,
            7787, 7805, 7819, 7826, 7885, 7900, 7908, 7926, 7930, 7951, 7965, 8000, 8072, 8109, 8122, 8123,
            8143, 8159, 8211, 8242, 8248, 8257, 8259, 8265, 8286, 8321, 8330, 8357, 8368, 8372, 8385, 8407,
            8430, 8436, 8444, 8448, 8454, 8490, 8507, 8513, 8556, 8604, 8639, 8750, 8751, 8755, 8764, 8777,
            8822, 8852, 8863, 8911, 8981, 2668, 889, 6066, 9529, 9676, 9990, 10157, 10183, 10243,
            10434] +
           [2321, 11032, 4484, 4750, 8084, 6770, 6624, 2749, 6378, 7833, 4399, 9547, 2253] +
           [2321, 10996, 4454, 10431, 11032, 10057, 4484, 8084, 10754, 2382, 938, 2228, 10887, 6770, 10895,
            10863, 6624, 10333, 4069, 10113, 4849, 9719, 10116, 245, 10077, 9547, 9557, 10457, 10037, 9900,
            10146, 5507, 10517, 2749, 9563, 6378, 2014, 9714, 1353, 10117, 7504, 9724, 3141, 5797, 10508, 10374,
            5593, 9932, 4399, 10632, 1039, 9904, 9930, 8505, 429, 10331, 5470, 8557, 7773, 10830, 10749, 2031,
            3822, 7833, 5791, 10602, 2203, 542, 10843, 7759, 10483, 4827, 225, 7679, 9617, 2378, 5409, 10142,
            9975, 10264, 918, 10148, 10066, 9917, 9485, 6400, 5961, 10023, 10418, 231, 10695, 3065, 6420, 7865,
            9813, 10765, 6290, 2270, 729, 2626, 8424, 10199, 2200, 2854, 2253],
    # note: first list were flipped, second list were already correct
    'hotspots': [6510, 6606, 2373, 6135, 6023, 7730, 2193, 8307, 5626, 4109, 2921, 4614, 2573, 7490, 6097, 7836,
                 1011, 6493, 5779, 8660, 6232, 6009, 2579, 929, 3824, 357, 6162, 477, 5640, 6467, 244, 2922, 4288,
                 2926, 1480, 4441, 4683, 8239, 9749, 9826, 9844, 10945, 11037] +
                [10069, 5485, 3146, 5560, 5600, 5937, 7688, 6055, 5670, 10235,
                 5583, 6151, 5720, 2553, 6040, 292, 5456, 2437, 5966, 5499, 10043,
                 10232, 5434, 6224, 5785, 6210, 2761, 6359, 6438, 5423, 5774,
                 7556, 5766, 7882, 7732, 5798, 2711, 2562, 5939, 2214, 2881,
                 2588, 10123, 6527, 10309, 2038, 2683, 5617, 2146, 4117, 10821,
                 2538, 5408, 5527, 6079, 7499, 6641, 2930, 5683, 6353, 5958, 2154, 5835] +
                [10015, 6339, 3168, 7481, 7779, 646, 4117, 7891, 4324, 2146, 10007, 1162,
                 10330, 285, 5527, 7616, 5617, 4196, 7771, 2085, 1104, 5512, 8303, 4409,
                 7343, 2538, 7570, 3977, 2683, 2038, 10250, 2494, 10309, 8423, 2417, 6353,
                 2564, 685, 6471, 6527, 10123, 2173, 2588, 205, 2881, 5863, 5958, 5408,
                 2607, 2214, 5939, 2562, 6079, 8208, 2154, 2799, 3909, 5798, 7663, 2574,
                 7732, 5835, 7882, 5766, 7556, 6438, 6641, 6359, 2761, 6210, 5434, 10232,
                 10043, 6109, 2323, 9550, 2437, 5456, 2930, 5683, 413, 2553, 7499, 6151,
                 5583, 10235, 6055, 7688, 5600, 5560, 3146, 5485, 10069, 5423, 5499, 6224, 5937, 292]
}

# neg should actually be neg, pos should actually be pos
LABELS_DYNAMIN = {
    'neg': [204777, 255597, 257856, 254414, 210090, 228965, 223120, 244787,
       210219, 205670, 207817, 208985, 243014, 244540, 243101, 205523,
       246994, 243813, 247387, 223499, 226986, 207715, 207461, 228022,
       246466, 246049, 249456, 249260, 243906, 248469, 248677, 246461,
       216470, 246446, 224584, 248959, 218106, 212688, 216575, 221288,
       246639, 213189, 250127, 216917, 239495, 218845, 249898, 248365,
       213853, 238928, 234087, 245726, 246329, 242735, 247350, 233576,
       243179, 242756, 247149, 249987, 248686, 245505, 245172, 237980,
       235687, 235962, 257613, 218102, 246146, 245821, 238930, 248793,
       243616, 243001, 242160, 229700, 222705, 224733, 214902, 205387,
       232491, 236402, 243552, 215134, 246401, 250057, 234118, 217468,
       243037, 239471, 217535, 249669, 235951, 238382, 237019, 208402,
       249106, 239288, 245737, 232361, 219713, 248476, 223001, 232697,
       231401, 215084, 252052, 214974, 205220, 257383, 219267, 212179,
       233113, 205274, 239585, 252200, 235908, 252770, 228385, 226748,
       218512, 242516, 256045, 252589, 238125, 222314] + \
       [258864, 228390, 207373, 205484, 207847, 258914], # from 40 worst errs
    'pos': [212256, 214854, 235817, 252489, 229125, 208935, 224734, 225681,
       257459, 220275, 237926, 207089, 237666, 215606, 252622, 227713,
       244095, 249564, 218127, 207014, 228064, 218584, 234681, 244427,
       205577, 221330, 222762, 259677, 217438, 227043, 239332, 225555,
       218725, 229340, 232956, 249714, 253432, 216647, 229976, 218328,
       243715, 205099, 252075, 216391, 218423, 215674, 223778, 238083,
       217819, 256932, 224605, 224576, 219263, 209148, 213457, 244520,
       219232, 208333, 228467, 228204, 255908, 249674, 213107, 221335,
       227820, 219138, 213488, 226788, 217890, 221283, 253437, 223868,
       246066, 238332, 256658, 210055, 230096, 204197, 212192, 205678,
       226214, 206266, 246347, 228949, 244897, 231391, 224139, 252028,
       219621, 253661, 245841, 214386, 222897, 225659, 247270, 204155,
       205157, 246161, 252953, 215964, 239835, 226266, 257059, 225921,
       211909, 228291, 211968, 224783, 246821, 249302, 224306, 247821,
       249017, 206983, 226609, 227327, 227044, 216729, 207840, 212872,
       225000, 228672, 234528, 225612, 234261, 214450, 212401, 236587,
       205938, 223768, 242149, 216730, 249207, 216074, 224790, 210343,
       212138, 253250, 213593, 235697, 235482, 237720, 217202, 233687,
       214668, 213157, 213751, 252837, 212942, 214825, 253023, 246366,
       234697] + \
       [238377, 216355, 212351, 244573, ], # from 40 worst errs
    'hotspots': [205537, 215233, 216903, 216241, 236608, 216941, 208591, 223088, # from tp
       217357, 214308, 237087, 247979, 207586, 216892, 217105, 213597,
       225263, 214879, 222906, 245501, 247460, 217602, 207269, 215196,
       216337, 245702, 224723, 245320, 213849, 209150, 226611, 248664,
       255361, 249998, 215481, 204719, 214674, 214432, 205126, 205146,
       218024, 206783, 225708, 247813, 255434, 208411, 212022, 212829,
       243765, 224915, 214699, 243846, 219117, 214503, 214622, 212061,
       208072, 243264, 213813, 228133, 224472, 247999, 213403, 208150,
       242500, 222445, 217523, 207271, 206599, 204323, 245793, 213013,
       218184, 216085, 236403, 208438, 246518, 225456, 216652, 245564,
       234101, 222997, 245241, 215456, 234109, 207612, 246110, 236069,
       204970, 225392, 235906, 212937, 234456, 246400, 225547, 213433,
       241137, 224283, 215557, 218920, 245648, 215118, 231402, 227704,
       235268, 242653, 254434, 244044, 229557, 246020, 204639, 223538,
       255904, 235393, 208086, 207861, 214203, 213260, 218675, 216595,
       237485, 204597, 239637, 212874, 246936, 243518, 235349, 254114,
       235877, 227335, 214540, 232583, 216441, 215495, 244976, 259687,
       204475, 224961, 233204, 218681, 242455, 239226, 223669, 214875,
       236654, 216784, 206824, 242479, 206267, 217824, 251958, 226253,
       214332] + \
       [242828, 254428, 257377, 207581, 235660, 205824, 225114, 213307, # from fn
       209487, 206981, 239467, 252051, 205115, 206939, 226274, 242758,
       228777, 256040, 255719, 235465, 209012, 222720, 256464, 223698,
       254836, 223231, 254747, 255906, 253835, 254133, 204951, 254969] + \
       [215974, 224389, 207496, 212337, 215555, 246347, 244847, 233344,
        244324, 242252], # from 40 worst errs
}


LABELS = {data_dir: None for data_dir in DSETS.keys()}
LABELS['clath_aux+gak_a7d2'] = LABELS_ORIG
LABELS['clath_aux_dynamin'] = LABELS_DYNAMIN



LABELS_DYNAMIN_NEW = {
        'hotspots': [172868, 77374, 164533, 173308, 76013, 214830, 115275, 86766, 176277, 139509, 218580, 97279, 130184, 83200, 244850, 176241, 243233, 108341, 147163, 191970, 168032, 76785, 180315, 147732, 94289, 153579, 224739, 244326, 108770, 103430, 183172, 185366, 112627, 191468, 76164, 242576, 126272, 147430, 124255, 175312, 137799, 147237, 173860, 186038, 147764, 124678, 146161, 74685, 216992, 93914, 92189, 218179, 182612, 123187, 96054, 243671, 223192, 97534, 216832, 181858, 117229, 165534, 79522, 186204, 148968, 194481, 95770, 232757, 136997, 192046, 69346, 125579, 181534, 130633, 147538, 75868, 187927, 244127, 254203, 92488, 103684, 162291, 75118, 96993, 146527, 154818, 184729, 68905, 255364, 150068, 84924, 125150, 166460, 132829, 187010, 183607, 78896, 255814, 96870, 246264, 181464, 249701, 249161, 107385, 107287, 148087, 183597, 154229, 152942, 186369, 191563, 178048, 173433, 190905, 85832, 165772, 98544, 182742, 140938, 164643, 212316, 106303, 205390, 185324, 209675, 149368, 152904, 76422, 248797, 183027, 162434, 97579, 182524, 254812, 252289, 133231, 222823, 172582, 180307, 147291, 209280, 227087, 188342, 236954, 244798, 130274, 223188, 243718, 147058, 212717, 102274, 192222, 114093, 155922, 166332, 87302, 136111, 138144, 115099, 212480, 191805, 145530, 145582, 175466, 114967, 92276, 77767, 175470, 192498, 121269, 93132, 92172, 183575, 138708, 83755, 113261, 94914, 153725, 225190, 159803, 175368, 174557, 92962, 157631, 104817, 191940, 191916, 115175, 95375, 113570, 148276, 115750, 175401, 86910, 68789, 147991, 207430, 205560, 146309, 115981, 148031, 132978, 83305, 156023, 118008, 131759, 94725, 191664, 223332, 110954, 94821, 93913, 138156, 252166, 104496, 116753, 138991, 69340, 114679, 139789, 147662],
        'pos': [142117, 144527, 225634, 97901, 249850, 236228, 233588, 125715, 96242, 158505, 257382, 127169, 133540, 118404, 134882, 123899, 135231, 97213, 138249, 84243, 176592, 104323, 225498, 131338, 209033, 148040, 126268, 218314, 75956, 123548, 77685, 133283, 148798, 172630, 194904, 109032, 164976, 190867, 147479, 209511, 122980, 113164, 121481, 125555, 122018, 213648, 121256, 76555, 132973, 245669, 121388, 190718, 122979, 123685, 124759, 109580, 78776, 183782, 252108, 126438, 121651, 126768, 244340, 125063, 92061, 94819, 174116, 245280, 195507, 187969, 166027, 122298, 182878, 235674, 112441, 123890, 207924, 125319, 242240, 84713, 124026, 108079, 150464, 181764, 130304, 141661, 194263, 186217, 125141, 121536, 138377, 257162, 74961, 158955, 137674, 172050, 254289, 85750, 94144, 80111, 131850, 79811, 122413, 82502, 185538, 172815, 182183, 218955, 234444, 80718, 116492, 84199, 258232, 181950, 96006, 237379, 220387, 101651, 80711, 151900, 142284, 168443, 81580, 237545, 94285, 137095, 159347, 104523, 216449, 171765, 118487, 125560, 89533, 92486, 80625, 229964, 224943, 212140, 249301, 119292, 108747, 240273, 235503, 217340, 75715, 89535, 215440, 149008, 89059, 123283, 112882, 137656, 176142, 86993, 108931, 237414, 157102, 187521, 208106, 234010, 245158, 186466, 209408, 208964, 252988, 75413, 166469, 124232, 172233, 87985, 181502, 145697, 122101, 89181, 90053, 118157, 134781, 244578, 228182, 184702, 166156, 123629, 130748, 186421, 205492, 212843, 224582, 247355, 109927, 184255, 124035, 232320, 177438, 126145, 157952, 94049, 155994, 127086, 174370, 145522, 146247, 149918, 229962, 244520, 147669, 106262, 75691, 173174, 164327, 204824, 147434, 93380, 160595, 172621, 159948, 251023, 166886, 134043, 126219, 181782, 126929, 208301, 178108, 230097, 139636, 97032, 257961, 149116, 193476, 114288, 114096, 125686, 116887, 126734, 157430, 79623, 216044, 103205, 91917, 75313, 101533, 258093, 223181, 92827, 83568, 74504, 83486, 154065, 176904, 174177, 186715, 239573, 233372, 154646, 137418, 246522, 117056, 129965, 141523, 109675, 96159, 191538, 78905, 156615, 235859, 218105, 253803, 104099, 105598, 101911, 123190, 212774, 148247, 216392, 133403, 230510, 121524, 84582, 192151, 151302, 206932, 92394, 146165, 178975, 244128, 228337, 90844, 218921, 132720, 75217, 103349, 209072, 244459, 156870, 191725, 186936, 125104, 126271, 134745, 244499, 88630, 127032, 246108, 253816, 107767, 206722, 88042, 128178, 89478, 129522, 124468, 219440, 212306, 254473, 192719, 121847, 138535, 134207, 141220, 104374, 190866, 84592, 158460, 106125, 126001, 75883, 207125, 217640, 154488, 143050, 89327, 79937, 124854, 76992, 188848, 218079, 214045, 127010, 162853, 79565, 252447, 228269, 156807, 192343, 101442, 184120, 171840, 174650, 93743, 95444, 182704, 105649, 253147, 117630, 153368, 129520, 182423, 188032, 95088, 126257, 256818, 115587, 148568, 157314, 159027, 174214, 85971, 106727, 148173, 89521, 227568, 101247, 211917, 134817, 160063, 83637, 242764, 143116, 137336, 224736, 255477, 76888, 137035, 123824, 123483, 166066, 115549, 108866, 160132, 171888, 242483, 222998, 103621, 187273, 91736, 123375, 137414, 162978, 228977, 116195, 177654, 257029, 236947, 138537, 176412, 138297, 256684, 188669, 147795, 183558, 191253, 172184, 249249, 92917, 238479, 181689, 125963, 176837, 177248, 68717, 80595, 184532, 156058, 192356, 159755],
        'neg': [166613, 146641, 216624, 184913, 249095, 175142, 118630, 133842, 82524, 248981, 150032, 155124, 156864, 150322, 141959, 236619, 118581, 237255, 92893, 166291, 122713, 143115, 75134, 76085, 159986, 98254, 149329, 117374, 142536, 215951, 133330, 234061, 121538, 134194, 139369, 194718, 225005, 248148, 121984, 144531, 206217, 114600, 181859, 104132, 183783, 170107, 137875, 145936, 132066, 177011, 129717, 192858, 108472, 162882, 78844, 142939, 184566, 115156, 78945, 75516, 98245, 97083, 205873, 176416, 82404, 173709, 227427, 245028, 178204, 209045, 118976, 115401, 131935, 129966, 244596, 225118, 77403, 245927, 256700, 238403, 208664, 132007, 95599, 117751, 69850, 97140, 186132, 167254, 163573, 117766, 133820, 193806, 135174, 134353, 146231, 122926, 150555, 131279, 74390, 208023, 146259, 132828, 125653, 131292, 134309, 129191, 129259, 149695, 76177, 188290, 130716, 209477, 94743, 132573, 146664, 75358, 121349, 77278, 75507, 253593, 125589, 251932, 209870, 96933, 222420, 252157, 132753, 149122, 226080, 206526, 258741, 206530, 131452, 133427, 80025, 79441, 147746, 258717, 130913, 222370, 151419, 130079, 130803, 133518, 223585, 129759, 93685, 227923, 226935, 121961, 222993, 182072, 253141, 148009, 76302, 230257, 193829, 129607, 208194, 105883, 223819, 205003, 204621, 133246, 130752, 258711, 228246, 122228, 148103, 113600, 129209, 103649, 94650, 189012, 101367, 94412, 121458, 133651, 118881, 98191, 75164, 146517, 132576, 186777, 130857, 76698, 132628, 148288, 78553, 106455, 74763, 225525, 80318, 130953, 105750, 135163, 77131, 78222, 148361, 97897, 75247, 130849, 98042, 196415, 187131, 129111, 204642, 148615, 116020, 149473, 95601, 74325, 224638, 92139, 133145, 129144, 106403, 252327, 226560, 204102, 140642, 126966, 218133, 169088, 102715, 252621, 226005, 98680, 108166, 79083, 168734, 149475, 146248, 97816, 254475, 195545, 95000, 254024, 76515, 188033, 127487, 192599, 150199, 91778, 132389, 155432, 163487, 186581, 121739, 107093, 147071, 226534, 184825, 75775, 101729, 188459, 249914, 147968, 253172, 77185, 109284, 75949, 129498, 252700, 129493, 130674, 223595, 108670, 154434, 102403, 191046, 255453, 106703, 258772, 226429, 105582, 155092, 167269, 77005, 146910, 78611, 105504, 167319, 153431, 163762, 136150, 107648, 251961, 91966, 96162, 130885, 88493, 134993, 230037, 110029, 204664, 144495, 126538, 164573, 93061, 166290, 225782, 209903, 229456, 103820, 155555, 109462, 168611, 124822, 167289, 131299, 207123, 78831, 147737, 129525, 208745, 75156, 95780, 167209, 113676, 78683, 96331, 153439, 103879, 177208, 132776, 80716, 129611, 243548, 217071, 179023, 178134, 98429, 222619, 150382, 132594, 132990, 257925, 174469, 126197, 173602, 89610, 186234, 185241, 160174, 188783, 255859, 130815, 132794, 132494, 181527, 150057, 130387, 134004, 121168, 167810, 88737, 110966, 165687, 132745, 133236, 256515, 164555, 148676, 190923, 157318, 79690, 86791, 130723, 165620, 76558, 79532, 122891, 151926, 171595, 81594, 167772, 250055, 182584, 130666, 129438, 144497, 134210, 146328, 208038, 102597, 175092, 147384, 96088, 166784, 134035, 165491, 79998, 246490, 87832, 136110, 104436, 131970, 176905, 150832, 171634, 103200, 84296, 150776, 215152, 99863, 131204, 215561, 131824, 104041, 76014, 109494, 125886, 183803, 103087, 158255, 125117, 130174, 179075, 215881, 131258, 101365, 124441, 105998, 168342, 130430, 150381, 117206, 148385, 173521, 215375, 104825, 109169, 85976, 133075, 163368, 252699, 235241, 146831, 129299, 130299, 86664, 257433, 169275, 155337, 219808, 114028, 98251, 132134, 153199, 158999, 148124, 159870, 74231, 242738, 252093, 157139, 204624, 86192, 237067, 93496, 123475, 149733, 215322, 145807, 134994, 155314, 78348, 168682, 185870, 174952, 215062, 176230, 253535, 115338, 253791, 145876, 228996, 134383, 167363, 79606, 162775, 190756, 96299, 134738, 112085, 219265, 109260, 192400, 249563, 79100, 101176, 182040, 176525, 143533, 146391, 95032, 257387, 140187, 78163, 243187, 183293, 256465, 133401, 254874, 148074, 132422, 164228, 255279, 126270, 82885, 85861, 133171, 105989, 182330, 209346, 155874, 108463, 243555, 129151, 155754, 149589, 102807, 130831, 187235, 183618, 232490, 121467, 163992, 154725, 154673, 204638, 153498, 164913, 76989, 157953, 133237, 257278, 77914, 134680, 158046, 74989, 244556, 155455, 76480, 244010, 177720, 147490, 179519, 156086, 227999, 114580, 109176, 213019, 158742, 132668, 181389, 258130, 254672, 108523, 102803, 164974, 114555, 155038, 80022, 181663, 129173, 129066, 135396, 253838, 114471, 248712, 185202, 141209, 204325, 145911, 158476, 157662, 99864, 134661, 124849, 75919, 121312, 83937, 122068, 156352, 206223, 252523, 133526, 124526, 186673, 89330, 88797, 154103, 129067, 132305, 180329, 176331, 234590, 74960, 92530, 205900, 256642, 213209, 193773, 131434, 254028, 239590, 222262, 103368, 132852, 242729, 71002, 96438, 123984, 94486, 95113, 88995, 134889, 80824, 175661, 233305, 125142, 147595, 135004, 101097, 177940, 166204, 167775, 253768, 104820, 82943, 222374, 112494, 165399, 243267, 115274, 214818, 132936, 113935, 156340, 167554, 248969, 147871, 103025, 122604, 91822, 252149, 107441, 172905, 168808, 153617, 105405, 109890, 145736, 102614, 218828, 95836, 76965, 209531, 149691, 69820, 76926, 123691, 238912, 181639, 246550, 205369, 248134, 141921, 152957, 182074, 230381, 159404, 226497, 208084, 184951, 208661, 167579, 164623, 123323, 104415, 234097, 258487, 118809, 158551, 106473, 131130, 129192, 148471, 104802, 140137, 93658, 205808, 86075, 86947, 235256, 170166, 139150]
}