densityfold_lin.cpp

#include <iostream>
#include <string>
#include <cstring>
#include <vector>
#include <fstream>
using namespace std;
#define INF 10000

#define STACK 101
#define BULGE 102
#define MULTILOOP 103
#define MISC 104
#define HAIRPIN 105

#define MULTI1 1001
#define MULTI2 1002
#define MULTI3 1003
#define MULTI4 1004
#define MULTIS 1005

struct data
{	
	int type;
	int i;
	int j;
	int i_prime;
	int j_prime;
};

// --------------------------------
//params for multibranch loop
float a = 3.4;
float b = 0;
float c = 0.4;
float sigma = 0;
// --------------------------------

////////////////
// GLOBALS
///////////////
//"-GGGGGUAUAGCUCAGGGGUAGAGCAUUUGACUGCAGAUCAAGAGGUCCCUGGUUCAAAUCCAGGUGCCCCCU"; 
// string str = "GCCGCGACAAGCGGUCCGGGCGCCCUAGGGGCCCGCGACGCCUGCUCGUACCCAUCUUGCUCCUUGGAGGAUUUGGCUAUGAGGA";
//char* str = "-AAAAAAGGGUUUUUU"; // GGGAACAAAGCUGAAGUACUUACCC";
// char* str = "-AAAAAAGGGUUUUUU";//"-AAAAAAAAAAAAGGGUUUUUUUUUUUU";
char* str = "-GCCGCGACAAGCGGUCCGGGCGCCCUAGGGGCCCGGCGGAGACGGGCGCCG"; // GAGGUGUCCGACGCCUGCUCGUACCCAUCUUGCUCCUUGGAGGAUUUGGCUAUGAGGA
// char* str = "-AAAAAAAAAAAAGGGUUUUUUUUUUUU"; 
int n = strlen(str) - 1;

int isPairing(char A, char B);
void print_matrix(const vector<vector<float> > &A);
void printStructure(const std::vector<std::pair<int, int> >& result);
void initializeTables(int length);

void tbELC(int j, std::vector<std::pair<int, int> >* result, float sigma);
void tbELCs(int i, int j, std::vector<std::pair<int, int> >* result, float sigma);
void reportloop(int i,int j,int k,int h, struct data ***solM, std::vector<std::pair<int, int> >* result, float sigma);

vector<std::pair<int, int> > result;

float* ELC;
float** ELCs;
float** ELCbi;
float** ELCm;
float** Es;
float** Em;
float** Ebi;

data** solM;
float** ELCm_dash;

// tracebacks 
data* traceback_ELC;
data** traceback_multi;
data** traceback_bulge;
data** traceback_ELCs;
int** hairpin_number;


////////////////////////////////////////////////////
/*Special cases for internal loop 1x2.
  Y 
  ------------------ 
  (X)  A    C    G    U   
  ------------------ 
  5' --> 3' 
  X
  P1 P2 
  YZ 
  3' <-- 5' 
  (A)  .    .    .    .    
  (C)  .    .    .    .    
  (G)  .    .    .    .    
  (U)  .    .    .    .    
 */
float sint1_2[] = {
	3.40,  3.40,  2.20,  3.90,  2.90,  2.90,  1.70,  3.40,  2.90,  2.90,  1.70,  3.40,  3.40,  3.40,  2.20,  3.90,  3.40,  3.40,  2.20,  3.90,  2.90,  2.90,  1.70,  3.40, 
	3.40,  3.40,  4.20,  3.40,  2.90,  2.90,  3.70,  2.90,  2.90,  2.90,  3.70,  2.90,  3.40,  3.40,  4.20,  3.40,  3.40,  3.40,  4.20,  3.40,  2.90,  2.90,  3.70,  2.90, 
	1.00,  2.70,  2.20,  2.70,  0.50,  2.20,  1.70,  2.20,  0.50,  2.20,  1.70,  2.20,  1.00,  2.70,  2.20,  2.70,  1.00,  2.70,  2.20,  2.70,  0.50,  2.20,  1.70,  2.20, 
	4.70,  4.00,  4.70,  3.00,  3.80,  3.10,  3.80,  2.10,  3.20,  2.50,  3.20,  1.50,  4.60,  3.90,  4.60,  2.90,  4.70,  4.00,  4.70,  3.00,  3.80,  3.10,  3.80,  2.10, 
	3.40,  3.40,  2.20,  3.90,  2.90,  2.90,  1.70,  3.40,  2.90,  2.90,  1.70,  3.40,  3.40,  3.40,  2.20,  3.90,  3.40,  3.40,  2.20,  3.90,  2.90,  2.90,  1.70,  3.40, 
	3.40,  3.40,  4.20,  3.40,  2.90,  2.90,  3.70,  2.90,  2.90,  2.90,  3.70,  2.90,  3.40,  3.40,  4.20,  3.40,  3.40,  3.40,  4.20,  3.40,  2.90,  2.90,  3.70,  2.90, 
	3.40,  5.10,  4.60,  5.10,  2.60,  4.30,  3.80,  4.30,  2.20,  3.90,  3.40,  3.90,  3.00,  4.70,  4.20,  4.70,  3.40,  5.10,  4.60,  5.10,  2.60,  4.30,  3.80,  4.30, 
	4.10,  3.40,  4.10,  2.40,  3.60,  2.90,  3.60,  1.90,  3.60,  2.90,  3.60,  1.90,  4.10,  3.40,  4.10,  2.40,  4.10,  3.40,  4.10,  2.40,  3.60,  2.90,  3.60,  1.90, 
	2.20,  2.20,  1.00,  2.70,  1.70,  1.70,  0.50,  2.20,  1.70,  1.70,  0.50,  2.20,  2.20,  2.20,  1.00,  2.70,  2.20,  2.20,  1.00,  2.70,  1.70,  1.70,  0.50,  2.20, 
	4.00,  4.00,  4.80,  4.00,  2.90,  2.90,  3.70,  2.90,  2.50,  2.50,  3.30,  2.50,  3.90,  3.90,  4.70,  3.90,  4.00,  4.00,  4.80,  4.00,  2.90,  2.90,  3.70,  2.90, 
	2.20,  3.90,  3.40,  3.90,  1.70,  3.40,  2.90,  3.40,  1.70,  3.40,  2.90,  3.40,  2.20,  3.90,  3.40,  3.90,  2.20,  3.90,  3.40,  3.90,  1.70,  3.40,  2.90,  3.40, 
	4.70,  4.00,  4.70,  3.00,  3.60,  2.90,  3.60,  1.90,  3.20,  2.50,  3.20,  1.50,  4.60,  3.90,  4.60,  2.90,  4.70,  4.00,  4.70,  3.00,  3.60,  2.90,  3.60,  1.90, 
	4.60,  4.60,  3.40,  5.10,  3.60,  3.60,  2.40,  4.10,  3.00,  3.00,  1.80,  3.50,  4.10,  4.10,  2.90,  4.60,  4.60,  4.60,  3.40,  5.10,  3.60,  3.60,  2.40,  4.10, 
	3.40,  3.40,  4.20,  3.40,  2.90,  2.90,  3.70,  2.90,  2.90,  2.90,  3.70,  2.90,  3.40,  3.40,  4.20,  3.40,  3.40,  3.40,  4.20,  3.40,  2.90,  2.90,  3.70,  2.90, 
	3.40,  5.10,  4.60,  5.10,  2.40,  4.10,  3.60,  4.10,  1.80,  3.50,  3.00,  3.50,  2.90,  4.60,  4.10,  4.60,  3.40,  5.10,  4.60,  5.10,  2.40,  4.10,  3.60,  4.10, 
	3.10,  2.40,  3.10,  1.40,  2.60,  1.90,  2.60,  0.90,  2.60,  1.90,  2.60,  0.90,  3.10,  2.40,  3.10,  1.40,  3.10,  2.40,  3.10,  1.40,  2.60,  1.90,  2.60,  0.90, 
	2.90,  2.90,  1.70,  2.50,  2.40,  2.40,  1.20,  2.00,  2.40,  2.40,  1.20,  2.00,  2.90,  2.90,  1.70,  2.50,  2.90,  2.90,  1.70,  2.50,  2.40,  2.40,  1.20,  2.00, 
	2.90,  2.90,  3.40,  2.90,  2.40,  2.40,  2.90,  2.40,  2.40,  2.40,  2.90,  2.40,  2.90,  2.90,  3.40,  2.90,  2.90,  2.90,  3.40,  2.90,  2.40,  2.40,  2.90,  2.40, 
	0.50,  1.30,  1.70,  1.30,  0.00,  0.80,  1.20,  0.80,  0.00,  0.80,  1.20,  0.80,  0.50,  1.30,  1.70,  1.30,  0.50,  1.30,  1.70,  1.30,  0.00,  0.80,  1.20,  0.80, 
	3.60,  3.50,  3.60,  2.50,  2.70,  2.60,  2.70,  1.60,  2.10,  2.00,  2.10,  1.00,  3.50,  3.40,  3.50,  2.40,  3.60,  3.50,  3.60,  2.50,  2.70,  2.60,  2.70,  1.60, 
	2.90,  2.90,  1.70,  2.50,  2.40,  2.40,  1.20,  2.00,  2.40,  2.40,  1.20,  2.00,  2.90,  2.90,  1.70,  2.50,  2.90,  2.90,  1.70,  2.50,  2.40,  2.40,  1.20,  2.00, 
	2.90,  2.90,  3.40,  2.90,  2.40,  2.40,  2.90,  2.40,  2.40,  2.40,  2.90,  2.40,  2.90,  2.90,  3.40,  2.90,  2.90,  2.90,  3.40,  2.90,  2.40,  2.40,  2.90,  2.40, 
	2.90,  3.70,  4.10,  3.70,  2.10,  2.90,  3.30,  2.90,  1.70,  2.50,  2.90,  2.50,  2.50,  3.30,  3.70,  3.30,  2.90,  3.70,  4.10,  3.70,  2.10,  2.90,  3.30,  2.90, 
	3.00,  2.90,  3.00,  1.90,  2.50,  2.40,  2.50,  1.40,  2.50,  2.40,  2.50,  1.40,  3.00,  2.90,  3.00,  1.90,  3.00,  2.90,  3.00,  1.90,  2.50,  2.40,  2.50,  1.40, 
	1.70,  1.70,  0.50,  1.30,  1.20,  1.20,  0.00,  0.80,  1.20,  1.20,  0.00,  0.80,  1.70,  1.70,  0.50,  1.30,  1.70,  1.70,  0.50,  1.30,  1.20,  1.20,  0.00,  0.80, 
	3.50,  3.50,  4.00,  3.50,  2.40,  2.40,  2.90,  2.40,  2.00,  2.00,  2.50,  2.00,  3.40,  3.40,  3.90,  3.40,  3.50,  3.50,  4.00,  3.50,  2.40,  2.40,  2.90,  2.40, 
	1.70,  2.50,  2.90,  2.50,  1.20,  2.00,  2.40,  2.00,  1.20,  2.00,  2.40,  2.00,  1.70,  2.50,  2.90,  2.50,  1.70,  2.50,  2.90,  2.50,  1.20,  2.00,  2.40,  2.00, 
	3.60,  3.50,  3.60,  2.50,  2.50,  2.40,  2.50,  1.40,  2.10,  2.00,  2.10,  1.00,  3.50,  3.40,  3.50,  2.40,  3.60,  3.50,  3.60,  2.50,  2.50,  2.40,  2.50,  1.40, 
	4.10,  4.10,  2.90,  3.70,  3.10,  3.10,  1.90,  2.70,  2.50,  2.50,  1.30,  2.10,  3.60,  3.60,  2.40,  3.20,  4.10,  4.10,  2.90,  3.70,  3.10,  3.10,  1.90,  2.70, 
	2.90,  2.90,  3.40,  2.90,  2.40,  2.40,  2.90,  2.40,  2.40,  2.40,  2.90,  2.40,  2.90,  2.90,  3.40,  2.90,  2.90,  2.90,  3.40,  2.90,  2.40,  2.40,  2.90,  2.40, 
	2.90,  3.70,  4.10,  3.70,  1.90,  2.70,  3.10,  2.70,  1.30,  2.10,  2.50,  2.10,  2.40,  3.20,  3.60,  3.20,  2.90,  3.70,  4.10,  3.70,  1.90,  2.70,  3.10,  2.70, 
	2.00,  1.90,  2.00,  0.90,  1.50,  1.40,  1.50,  0.40,  1.50,  1.40,  1.50,  0.40,  2.00,  1.90,  2.00,  0.90,  2.00,  1.90,  2.00,  0.90,  1.50,  1.40,  1.50,  0.40, 
	2.90,  2.90,  1.70,  3.10,  2.40,  2.40,  1.20,  2.60,  2.40,  2.40,  1.20,  2.60,  2.90,  2.90,  1.70,  3.10,  2.90,  2.90,  1.70,  3.10,  2.40,  2.40,  1.20,  2.60, 
	2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40,  2.40,  2.40,  3.30,  2.40,  2.90,  2.90,  3.80,  2.90,  2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40, 
	0.50,  1.70,  1.70,  1.70,  0.00,  1.20,  1.20,  1.20,  0.00,  1.20,  1.20,  1.20,  0.50,  1.70,  1.70,  1.70,  0.50,  1.70,  1.70,  1.70,  0.00,  1.20,  1.20,  1.20, 
	4.20,  3.50,  4.20,  2.50,  3.30,  2.60,  3.30,  1.60,  2.70,  2.00,  2.70,  1.00,  4.10,  3.40,  4.10,  2.40,  4.20,  3.50,  4.20,  2.50,  3.30,  2.60,  3.30,  1.60, 
	2.90,  2.90,  1.70,  3.10,  2.40,  2.40,  1.20,  2.60,  2.40,  2.40,  1.20,  2.60,  2.90,  2.90,  1.70,  3.10,  2.90,  2.90,  1.70,  3.10,  2.40,  2.40,  1.20,  2.60, 
	2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40,  2.40,  2.40,  3.30,  2.40,  2.90,  2.90,  3.80,  2.90,  2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40, 
	2.90,  4.10,  4.10,  4.10,  2.10,  3.30,  3.30,  3.30,  1.70,  2.90,  2.90,  2.90,  2.50,  3.70,  3.70,  3.70,  2.90,  4.10,  4.10,  4.10,  2.10,  3.30,  3.30,  3.30, 
	3.60,  2.90,  3.60,  1.90,  3.10,  2.40,  3.10,  1.40,  3.10,  2.40,  3.10,  1.40,  3.60,  2.90,  3.60,  1.90,  3.60,  2.90,  3.60,  1.90,  3.10,  2.40,  3.10,  1.40, 
	1.70,  1.70,  0.50,  1.90,  1.20,  1.20,  0.00,  1.40,  1.20,  1.20,  0.00,  1.40,  1.70,  1.70,  0.50,  1.90,  1.70,  1.70,  0.50,  1.90,  1.20,  1.20,  0.00,  1.40, 
	3.50,  3.50,  4.40,  3.50,  2.40,  2.40,  3.30,  2.40,  2.00,  2.00,  2.90,  2.00,  3.40,  3.40,  4.30,  3.40,  3.50,  3.50,  4.40,  3.50,  2.40,  2.40,  3.30,  2.40, 
	1.70,  2.90,  2.90,  2.90,  1.20,  2.40,  2.40,  2.40,  1.20,  2.40,  2.40,  2.40,  1.70,  2.90,  2.90,  2.90,  1.70,  2.90,  2.90,  2.90,  1.20,  2.40,  2.40,  2.40, 
	4.20,  3.50,  4.20,  2.50,  3.10,  2.40,  3.10,  1.40,  2.70,  2.00,  2.70,  1.00,  4.10,  3.40,  4.10,  2.40,  4.20,  3.50,  4.20,  2.50,  3.10,  2.40,  3.10,  1.40, 
	4.10,  4.10,  2.90,  4.30,  3.10,  3.10,  1.90,  3.30,  2.50,  2.50,  1.30,  2.70,  3.60,  3.60,  2.40,  3.80,  4.10,  4.10,  2.90,  4.30,  3.10,  3.10,  1.90,  3.30, 
	2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40,  2.40,  2.40,  3.30,  2.40,  2.90,  2.90,  3.80,  2.90,  2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40, 
	2.90,  4.10,  4.10,  4.10,  1.90,  3.10,  3.10,  3.10,  1.30,  2.50,  2.50,  2.50,  2.40,  3.60,  3.60,  3.60,  2.90,  4.10,  4.10,  4.10,  1.90,  3.10,  3.10,  3.10, 
	2.60,  1.90,  2.60,  0.90,  2.10,  1.40,  2.10,  0.40,  2.10,  1.40,  2.10,  0.40,  2.60,  1.90,  2.60,  0.90,  2.60,  1.90,  2.60,  0.90,  2.10,  1.40,  2.10,  0.40, 
	3.40,  3.40,  2.20,  4.00,  2.90,  2.90,  1.70,  3.50,  2.90,  2.90,  1.70,  3.50,  3.40,  3.40,  2.20,  4.00,  3.40,  3.40,  2.20,  4.00,  2.90,  2.90,  1.70,  3.50, 
	3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90,  2.90,  2.90,  4.10,  2.90,  3.40,  3.40,  4.60,  3.40,  3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90, 
	1.00,  2.80,  2.20,  2.80,  0.50,  2.30,  1.70,  2.30,  0.50,  2.30,  1.70,  2.30,  1.00,  2.80,  2.20,  2.80,  1.00,  2.80,  2.20,  2.80,  0.50,  2.30,  1.70,  2.30, 
	5.20,  4.00,  5.20,  3.00,  4.30,  3.10,  4.30,  2.10,  3.70,  2.50,  3.70,  1.50,  5.10,  3.90,  5.10,  2.90,  5.20,  4.00,  5.20,  3.00,  4.30,  3.10,  4.30,  2.10, 
	3.40,  3.40,  2.20,  4.00,  2.90,  2.90,  1.70,  3.50,  2.90,  2.90,  1.70,  3.50,  3.40,  3.40,  2.20,  4.00,  3.40,  3.40,  2.20,  4.00,  2.90,  2.90,  1.70,  3.50, 
	3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90,  2.90,  2.90,  4.10,  2.90,  3.40,  3.40,  4.60,  3.40,  3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90, 
	3.40,  5.20,  4.60,  5.20,  2.60,  4.40,  3.80,  4.40,  2.20,  4.00,  3.40,  4.00,  3.00,  4.80,  4.20,  4.80,  3.40,  5.20,  4.60,  5.20,  2.60,  4.40,  3.80,  4.40, 
	4.60,  3.40,  4.60,  2.40,  4.10,  2.90,  4.10,  1.90,  4.10,  2.90,  4.10,  1.90,  4.60,  3.40,  4.60,  2.40,  4.60,  3.40,  4.60,  2.40,  4.10,  2.90,  4.10,  1.90, 
	2.20,  2.20,  1.00,  2.80,  1.70,  1.70,  0.50,  2.30,  1.70,  1.70,  0.50,  2.30,  2.20,  2.20,  1.00,  2.80,  2.20,  2.20,  1.00,  2.80,  1.70,  1.70,  0.50,  2.30, 
	4.00,  4.00,  5.20,  4.00,  2.90,  2.90,  4.10,  2.90,  2.50,  2.50,  3.70,  2.50,  3.90,  3.90,  5.10,  3.90,  4.00,  4.00,  5.20,  4.00,  2.90,  2.90,  4.10,  2.90, 
	2.20,  4.00,  3.40,  4.00,  1.70,  3.50,  2.90,  3.50,  1.70,  3.50,  2.90,  3.50,  2.20,  4.00,  3.40,  4.00,  2.20,  4.00,  3.40,  4.00,  1.70,  3.50,  2.90,  3.50, 
	5.20,  4.00,  5.20,  3.00,  4.10,  2.90,  4.10,  1.90,  3.70,  2.50,  3.70,  1.50,  5.10,  3.90,  5.10,  2.90,  5.20,  4.00,  5.20,  3.00,  4.10,  2.90,  4.10,  1.90, 
	4.60,  4.60,  3.40,  5.20,  3.60,  3.60,  2.40,  4.20,  3.00,  3.00,  1.80,  3.60,  4.10,  4.10,  2.90,  4.70,  4.60,  4.60,  3.40,  5.20,  3.60,  3.60,  2.40,  4.20, 
	3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90,  2.90,  2.90,  4.10,  2.90,  3.40,  3.40,  4.60,  3.40,  3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90, 
	3.40,  5.20,  4.60,  5.20,  2.40,  4.20,  3.60,  4.20,  1.80,  3.60,  3.00,  3.60,  2.90,  4.70,  4.10,  4.70,  3.40,  5.20,  4.60,  5.20,  2.40,  4.20,  3.60,  4.20, 
	3.60,  2.40,  3.60,  1.40,  3.10,  1.90,  3.10,  0.90,  3.10,  1.90,  3.10,  0.90,  3.60,  2.40,  3.60,  1.40,  3.60,  2.40,  3.60,  1.40,  3.10,  1.90,  3.10,  0.90, 
	2.90,  2.90,  1.70,  3.10,  2.40,  2.40,  1.20,  2.60,  2.40,  2.40,  1.20,  2.60,  2.90,  2.90,  1.70,  3.10,  2.90,  2.90,  1.70,  3.10,  2.40,  2.40,  1.20,  2.60, 
	2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40,  2.40,  2.40,  3.30,  2.40,  2.90,  2.90,  3.80,  2.90,  2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40, 
	0.50,  1.70,  1.70,  1.70,  0.00,  1.20,  1.20,  1.20,  0.00,  1.20,  1.20,  1.20,  0.50,  1.70,  1.70,  1.70,  0.50,  1.70,  1.70,  1.70,  0.00,  1.20,  1.20,  1.20, 
	4.20,  3.50,  4.20,  2.50,  3.30,  2.60,  3.30,  1.60,  2.70,  2.00,  2.70,  1.00,  4.10,  3.40,  4.10,  2.40,  4.20,  3.50,  4.20,  2.50,  3.30,  2.60,  3.30,  1.60, 
	2.90,  2.90,  1.70,  3.10,  2.40,  2.40,  1.20,  2.60,  2.40,  2.40,  1.20,  2.60,  2.90,  2.90,  1.70,  3.10,  2.90,  2.90,  1.70,  3.10,  2.40,  2.40,  1.20,  2.60, 
	2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40,  2.40,  2.40,  3.30,  2.40,  2.90,  2.90,  3.80,  2.90,  2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40, 
	2.90,  4.10,  4.10,  4.10,  2.10,  3.30,  3.30,  3.30,  1.70,  2.90,  2.90,  2.90,  2.50,  3.70,  3.70,  3.70,  2.90,  4.10,  4.10,  4.10,  2.10,  3.30,  3.30,  3.30, 
	3.60,  2.90,  3.60,  1.90,  3.10,  2.40,  3.10,  1.40,  3.10,  2.40,  3.10,  1.40,  3.60,  2.90,  3.60,  1.90,  3.60,  2.90,  3.60,  1.90,  3.10,  2.40,  3.10,  1.40, 
	1.70,  1.70,  0.50,  1.90,  1.20,  1.20,  0.00,  1.40,  1.20,  1.20,  0.00,  1.40,  1.70,  1.70,  0.50,  1.90,  1.70,  1.70,  0.50,  1.90,  1.20,  1.20,  0.00,  1.40, 
	3.50,  3.50,  4.40,  3.50,  2.40,  2.40,  3.30,  2.40,  2.00,  2.00,  2.90,  2.00,  3.40,  3.40,  4.30,  3.40,  3.50,  3.50,  4.40,  3.50,  2.40,  2.40,  3.30,  2.40, 
	1.70,  2.90,  2.90,  2.90,  1.20,  2.40,  2.40,  2.40,  1.20,  2.40,  2.40,  2.40,  1.70,  2.90,  2.90,  2.90,  1.70,  2.90,  2.90,  2.90,  1.20,  2.40,  2.40,  2.40, 
	4.20,  3.50,  4.20,  2.50,  3.10,  2.40,  3.10,  1.40,  2.70,  2.00,  2.70,  1.00,  4.10,  3.40,  4.10,  2.40,  4.20,  3.50,  4.20,  2.50,  3.10,  2.40,  3.10,  1.40, 
	4.10,  4.10,  2.90,  4.30,  3.10,  3.10,  1.90,  3.30,  2.50,  2.50,  1.30,  2.70,  3.60,  3.60,  2.40,  3.80,  4.10,  4.10,  2.90,  4.30,  3.10,  3.10,  1.90,  3.30, 
	2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40,  2.40,  2.40,  3.30,  2.40,  2.90,  2.90,  3.80,  2.90,  2.90,  2.90,  3.80,  2.90,  2.40,  2.40,  3.30,  2.40, 
	2.90,  4.10,  4.10,  4.10,  1.90,  3.10,  3.10,  3.10,  1.30,  2.50,  2.50,  2.50,  2.40,  3.60,  3.60,  3.60,  2.90,  4.10,  4.10,  4.10,  1.90,  3.10,  3.10,  3.10, 
	2.60,  1.90,  2.60,  0.90,  2.10,  1.40,  2.10,  0.40,  2.10,  1.40,  2.10,  0.40,  2.60,  1.90,  2.60,  0.90,  2.60,  1.90,  2.60,  0.90,  2.10,  1.40,  2.10,  0.40, 
	3.40,  3.40,  2.20,  4.00,  2.90,  2.90,  1.70,  3.50,  2.90,  2.90,  1.70,  3.50,  3.40,  3.40,  2.20,  4.00,  3.40,  3.40,  2.20,  4.00,  2.90,  2.90,  1.70,  3.50, 
	3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90,  2.90,  2.90,  4.10,  2.90,  3.40,  3.40,  4.60,  3.40,  3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90, 
	1.00,  2.80,  2.20,  2.80,  0.50,  2.30,  1.70,  2.30,  0.50,  2.30,  1.70,  2.30,  1.00,  2.80,  2.20,  2.80,  1.00,  2.80,  2.20,  2.80,  0.50,  2.30,  1.70,  2.30, 
	5.20,  4.00,  5.20,  3.00,  4.30,  3.10,  4.30,  2.10,  3.70,  2.50,  3.70,  1.50,  5.10,  3.90,  5.10,  2.90,  5.20,  4.00,  5.20,  3.00,  4.30,  3.10,  4.30,  2.10, 
	3.40,  3.40,  2.20,  4.00,  2.90,  2.90,  1.70,  3.50,  2.90,  2.90,  1.70,  3.50,  3.40,  3.40,  2.20,  4.00,  3.40,  3.40,  2.20,  4.00,  2.90,  2.90,  1.70,  3.50, 
	3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90,  2.90,  2.90,  4.10,  2.90,  3.40,  3.40,  4.60,  3.40,  3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90, 
	3.40,  5.20,  4.60,  5.20,  2.60,  4.40,  3.80,  4.40,  2.20,  4.00,  3.40,  4.00,  3.00,  4.80,  4.20,  4.80,  3.40,  5.20,  4.60,  5.20,  2.60,  4.40,  3.80,  4.40, 
	4.60,  3.40,  4.60,  2.40,  4.10,  2.90,  4.10,  1.90,  4.10,  2.90,  4.10,  1.90,  4.60,  3.40,  4.60,  2.40,  4.60,  3.40,  4.60,  2.40,  4.10,  2.90,  4.10,  1.90, 
	2.20,  2.20,  1.00,  2.80,  1.70,  1.70,  0.50,  2.30,  1.70,  1.70,  0.50,  2.30,  2.20,  2.20,  1.00,  2.80,  2.20,  2.20,  1.00,  2.80,  1.70,  1.70,  0.50,  2.30, 
	4.00,  4.00,  5.20,  4.00,  2.90,  2.90,  4.10,  2.90,  2.50,  2.50,  3.70,  2.50,  3.90,  3.90,  5.10,  3.90,  4.00,  4.00,  5.20,  4.00,  2.90,  2.90,  4.10,  2.90, 
	2.20,  4.00,  3.40,  4.00,  1.70,  3.50,  2.90,  3.50,  1.70,  3.50,  2.90,  3.50,  2.20,  4.00,  3.40,  4.00,  2.20,  4.00,  3.40,  4.00,  1.70,  3.50,  2.90,  3.50, 
	5.20,  4.00,  5.20,  3.00,  4.10,  2.90,  4.10,  1.90,  3.70,  2.50,  3.70,  1.50,  5.10,  3.90,  5.10,  2.90,  5.20,  4.00,  5.20,  3.00,  4.10,  2.90,  4.10,  1.90, 
	4.60,  4.60,  3.40,  5.20,  3.60,  3.60,  2.40,  4.20,  3.00,  3.00,  1.80,  3.60,  4.10,  4.10,  2.90,  4.70,  4.60,  4.60,  3.40,  5.20,  3.60,  3.60,  2.40,  4.20, 
	3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90,  2.90,  2.90,  4.10,  2.90,  3.40,  3.40,  4.60,  3.40,  3.40,  3.40,  4.60,  3.40,  2.90,  2.90,  4.10,  2.90, 
	3.40,  5.20,  4.60,  5.20,  2.40,  4.20,  3.60,  4.20,  1.80,  3.60,  3.00,  3.60,  2.90,  4.70,  4.10,  4.70,  3.40,  5.20,  4.60,  5.20,  2.40,  4.20,  3.60,  4.20, 
	3.60,  2.40,  3.60,  1.40,  3.10,  1.90,  3.10,  0.90,  3.10,  1.90,  3.10,  0.90,  3.60,  2.40,  3.60,  1.40,  3.60,  2.40,  3.60,  1.40,  3.10,  1.90,  3.10,  0.90 
};
// check the direction of X and Y 
/*Special cases for internal loop 2x2.

  Y 
  ---------------------------------------------------------------- 
  (X)   A   A   A   A   C   C   C   C   G   G   G   G   U   U   U   U 
  A   C   G   U   A   C   G   U   A   C   G   U   A   C   G   U 
  ---------------------------------------------------------------- 
  5' ------> 3' 
  A \/ \_/ A 
  U /\  |  U 
  3' <------ 5' 

  (AA)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (AC)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (AG)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (AU)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (CA)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (CC)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (CG)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (CU)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (GA)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (GC)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (GG)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (GU)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (UA)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (UC)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (UG)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
  (UU)  .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   
 */
float sint2_2[] = {
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.40,  3.40,  2.20,  4.60,  3.40,  3.40,  4.00,  3.40,  2.20,  4.60,  3.40,  4.60,  4.00,  3.40,  4.00,  2.40, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.70,  3.70,  2.50,  4.90,  3.70,  3.70,  4.30,  3.70,  2.50,  4.90,  3.70,  4.90,  4.30,  3.70,  4.30,  2.70, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.40,  3.40,  2.20,  4.60,  3.40,  3.40,  4.00,  3.40,  2.20,  4.60,  3.40,  4.60,  4.00,  3.40,  4.00,  2.40, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.40,  3.40,  2.20,  4.60,  3.40,  3.40,  4.00,  3.40,  2.20,  4.60,  3.40,  4.60,  4.00,  3.40,  4.00,  2.40, 
	3.60,  3.60,  2.40,  4.80,  3.60,  3.60,  4.20,  3.60,  2.40,  4.80,  3.60,  4.80,  4.20,  3.60,  4.20,  2.60, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.60,  3.60,  2.40,  4.80,  3.60,  3.60,  4.20,  3.60,  2.40,  4.80,  3.60,  4.80,  4.20,  3.60,  4.20,  2.60, 
	1.90,  1.90,  0.70,  3.10,  1.90,  1.90,  2.50,  1.90,  0.70,  3.10,  1.90,  3.10,  2.50,  1.90,  2.50,  0.90, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  2.90,  2.90,  1.70,  3.80,  2.90,  3.60,  3.10,  2.90,  2.90,  1.90, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.20,  3.20,  2.00,  3.90,  3.20,  3.20,  3.20,  3.20,  2.00,  4.10,  3.20,  3.90,  3.40,  3.20,  3.20,  2.20, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  2.90,  2.90,  1.70,  3.80,  2.90,  3.60,  3.10,  2.90,  2.90,  1.90, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  2.90,  2.90,  1.70,  3.80,  2.90,  3.60,  3.10,  2.90,  2.90,  1.90, 
	3.10,  3.10,  1.90,  3.80,  3.10,  3.10,  3.10,  3.10,  1.90,  4.00,  3.10,  3.80,  3.30,  3.10,  3.10,  2.10, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.10,  3.10,  1.90,  3.80,  3.10,  3.10,  3.10,  3.10,  1.90,  4.00,  3.10,  3.80,  3.30,  3.10,  3.10,  2.10, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.40,  1.40,  0.20,  2.30,  1.40,  2.10,  1.60,  1.40,  1.40,  0.40, 

	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	1.20,  1.20,  0.00,  1.30,  1.20,  1.20,  0.80,  1.20,  0.00,  1.70,  1.20,  1.30,  0.80,  1.20,  0.80,  0.20, 
	2.90,  2.90,  1.70,  3.00,  2.90,  2.90,  2.50,  2.90,  1.70,  3.40,  2.90,  3.00,  2.50,  2.90,  2.50,  1.90, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	3.20,  3.20,  2.00,  3.30,  3.20,  3.20,  2.80,  3.20,  2.00,  3.70,  3.20,  3.30,  2.80,  3.20,  2.80,  2.20, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	1.20,  1.20,  0.00,  1.30,  1.20,  1.20,  0.80,  1.20,  0.00,  1.70,  1.20,  1.30,  0.80,  1.20,  0.80,  0.20, 
	2.90,  2.90,  1.70,  3.00,  2.90,  2.90,  2.50,  2.90,  1.70,  3.40,  2.90,  3.00,  2.50,  2.90,  2.50,  1.90, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	2.90,  2.90,  1.70,  3.00,  2.90,  2.90,  2.50,  2.90,  1.70,  3.40,  2.90,  3.00,  2.50,  2.90,  2.50,  1.90, 
	3.10,  3.10,  1.90,  3.20,  3.10,  3.10,  2.70,  3.10,  1.90,  3.60,  3.10,  3.20,  2.70,  3.10,  2.70,  2.10, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	3.10,  3.10,  1.90,  3.20,  3.10,  3.10,  2.70,  3.10,  1.90,  3.60,  3.10,  3.20,  2.70,  3.10,  2.70,  2.10, 
	1.40,  1.40,  0.20,  1.50,  1.40,  1.40,  1.00,  1.40,  0.20,  1.90,  1.40,  1.50,  1.00,  1.40,  1.00,  0.40, 

	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	1.70,  1.70,  0.50,  2.40,  1.70,  1.70,  2.20,  1.70,  0.50,  2.50,  1.70,  2.40,  2.20,  1.70,  2.20,  0.70, 
	3.40,  3.40,  2.20,  4.10,  3.40,  3.40,  3.90,  3.40,  2.20,  4.20,  3.40,  4.10,  3.90,  3.40,  3.90,  2.40, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	3.70,  3.70,  2.50,  4.40,  3.70,  3.70,  4.20,  3.70,  2.50,  4.50,  3.70,  4.40,  4.20,  3.70,  4.20,  2.70, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	1.70,  1.70,  0.50,  2.40,  1.70,  1.70,  2.20,  1.70,  0.50,  2.50,  1.70,  2.40,  2.20,  1.70,  2.20,  0.70, 
	3.40,  3.40,  2.20,  4.10,  3.40,  3.40,  3.90,  3.40,  2.20,  4.20,  3.40,  4.10,  3.90,  3.40,  3.90,  2.40, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	3.40,  3.40,  2.20,  4.10,  3.40,  3.40,  3.90,  3.40,  2.20,  4.20,  3.40,  4.10,  3.90,  3.40,  3.90,  2.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  4.10,  3.60,  2.40,  4.40,  3.60,  4.30,  4.10,  3.60,  4.10,  2.60, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  4.10,  3.60,  2.40,  4.40,  3.60,  4.30,  4.10,  3.60,  4.10,  2.60, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  2.40,  1.90,  0.70,  2.70,  1.90,  2.60,  2.40,  1.90,  2.40,  0.90, 

	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.40,  3.40,  2.20,  4.60,  3.40,  3.40,  4.00,  3.40,  2.20,  4.60,  3.40,  4.60,  4.00,  3.40,  4.00,  2.40, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.70,  3.70,  2.50,  4.90,  3.70,  3.70,  4.30,  3.70,  2.50,  4.90,  3.70,  4.90,  4.30,  3.70,  4.30,  2.70, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.40,  3.40,  2.20,  4.60,  3.40,  3.40,  4.00,  3.40,  2.20,  4.60,  3.40,  4.60,  4.00,  3.40,  4.00,  2.40, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.40,  3.40,  2.20,  4.60,  3.40,  3.40,  4.00,  3.40,  2.20,  4.60,  3.40,  4.60,  4.00,  3.40,  4.00,  2.40, 
	3.60,  3.60,  2.40,  4.80,  3.60,  3.60,  4.20,  3.60,  2.40,  4.80,  3.60,  4.80,  4.20,  3.60,  4.20,  2.60, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.60,  3.60,  2.40,  4.80,  3.60,  3.60,  4.20,  3.60,  2.40,  4.80,  3.60,  4.80,  4.20,  3.60,  4.20,  2.60, 
	1.90,  1.90,  0.70,  3.10,  1.90,  1.90,  2.50,  1.90,  0.70,  3.10,  1.90,  3.10,  2.50,  1.90,  2.50,  0.90, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  2.90,  2.90,  1.70,  3.80,  2.90,  3.60,  3.10,  2.90,  2.90,  1.90, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.20,  3.20,  2.00,  3.90,  3.20,  3.20,  3.20,  3.20,  2.00,  4.10,  3.20,  3.90,  3.40,  3.20,  3.20,  2.20, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  2.90,  2.90,  1.70,  3.80,  2.90,  3.60,  3.10,  2.90,  2.90,  1.90, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  2.90,  2.90,  1.70,  3.80,  2.90,  3.60,  3.10,  2.90,  2.90,  1.90, 
	3.10,  3.10,  1.90,  3.80,  3.10,  3.10,  3.10,  3.10,  1.90,  4.00,  3.10,  3.80,  3.30,  3.10,  3.10,  2.10, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.10,  3.10,  1.90,  3.80,  3.10,  3.10,  3.10,  3.10,  1.90,  4.00,  3.10,  3.80,  3.30,  3.10,  3.10,  2.10, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.40,  1.40,  0.20,  2.30,  1.40,  2.10,  1.60,  1.40,  1.40,  0.40, 

	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.00,  2.00,  0.80,  3.20,  2.00,  2.00,  2.60,  2.00,  0.80,  3.20,  2.00,  3.20,  2.60,  2.00,  2.60,  1.00, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.00,  2.00,  0.80,  3.20,  2.00,  2.00,  2.60,  2.00,  0.80,  3.20,  2.00,  3.20,  2.60,  2.00,  2.60,  1.00, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.00,  2.00,  0.80,  3.20,  2.00,  2.00,  2.60,  2.00,  0.80,  3.20,  2.00,  3.20,  2.60,  2.00,  2.60,  1.00, 
	2.50,  2.50,  1.30,  3.70,  2.50,  2.50,  3.10,  2.50,  1.30,  3.70,  2.50,  3.70,  3.10,  2.50,  3.10,  1.50, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.50,  2.50,  1.30,  3.70,  2.50,  2.50,  3.10,  2.50,  1.30,  3.70,  2.50,  3.70,  3.10,  2.50,  3.10,  1.50, 
	1.40,  1.40,  0.20,  2.60,  1.40,  1.40,  2.00,  1.40,  0.20,  2.60,  1.40,  2.60,  2.00,  1.40,  2.00,  0.40, 

	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	1.50,  1.50,  0.30,  2.20,  1.50,  1.50,  1.50,  1.50,  0.30,  2.40,  1.50,  2.20,  1.70,  1.50,  1.50,  0.50, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	1.50,  1.50,  0.30,  2.20,  1.50,  1.50,  1.50,  1.50,  0.30,  2.40,  1.50,  2.20,  1.70,  1.50,  1.50,  0.50, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.50,  1.50,  0.30,  2.20,  1.50,  1.50,  1.50,  1.50,  0.30,  2.40,  1.50,  2.20,  1.70,  1.50,  1.50,  0.50, 
	2.00,  2.00,  0.80,  2.70,  2.00,  2.00,  2.00,  2.00,  0.80,  2.90,  2.00,  2.70,  2.20,  2.00,  2.00,  1.00, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.00,  2.00,  0.80,  2.70,  2.00,  2.00,  2.00,  2.00,  0.80,  2.90,  2.00,  2.70,  2.20,  2.00,  2.00,  1.00, 
	0.90,  0.90, -0.30,  1.60,  0.90,  0.90,  0.90,  0.90, -0.30,  1.80,  0.90,  1.60,  1.10,  0.90,  0.90, -0.10, 

	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	0.70,  0.70, -0.50,  0.80,  0.70,  0.70,  0.30,  0.70, -0.50,  1.20,  0.70,  0.80,  0.30,  0.70,  0.30, -0.30, 
	1.50,  1.50,  0.30,  1.60,  1.50,  1.50,  1.10,  1.50,  0.30,  2.00,  1.50,  1.60,  1.10,  1.50,  1.10,  0.50, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	0.70,  0.70, -0.50,  0.80,  0.70,  0.70,  0.30,  0.70, -0.50,  1.20,  0.70,  0.80,  0.30,  0.70,  0.30, -0.30, 
	1.50,  1.50,  0.30,  1.60,  1.50,  1.50,  1.10,  1.50,  0.30,  2.00,  1.50,  1.60,  1.10,  1.50,  1.10,  0.50, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.50,  1.50,  0.30,  1.60,  1.50,  1.50,  1.10,  1.50,  0.30,  2.00,  1.50,  1.60,  1.10,  1.50,  1.10,  0.50, 
	2.00,  2.00,  0.80,  2.10,  2.00,  2.00,  1.60,  2.00,  0.80,  2.50,  2.00,  2.10,  1.60,  2.00,  1.60,  1.00, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	2.00,  2.00,  0.80,  2.10,  2.00,  2.00,  1.60,  2.00,  0.80,  2.50,  2.00,  2.10,  1.60,  2.00,  1.60,  1.00, 
	0.90,  0.90, -0.30,  1.00,  0.90,  0.90,  0.50,  0.90, -0.30,  1.40,  0.90,  1.00,  0.50,  0.90,  0.50, -0.10, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.70,  1.20,  0.00,  2.00,  1.20,  1.90,  1.70,  1.20,  1.70,  0.20, 
	2.00,  2.00,  0.80,  2.70,  2.00,  2.00,  2.50,  2.00,  0.80,  2.80,  2.00,  2.70,  2.50,  2.00,  2.50,  1.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.70,  1.20,  0.00,  2.00,  1.20,  1.90,  1.70,  1.20,  1.70,  0.20, 
	2.00,  2.00,  0.80,  2.70,  2.00,  2.00,  2.50,  2.00,  0.80,  2.80,  2.00,  2.70,  2.50,  2.00,  2.50,  1.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.00,  2.00,  0.80,  2.70,  2.00,  2.00,  2.50,  2.00,  0.80,  2.80,  2.00,  2.70,  2.50,  2.00,  2.50,  1.00, 
	2.50,  2.50,  1.30,  3.20,  2.50,  2.50,  3.00,  2.50,  1.30,  3.30,  2.50,  3.20,  3.00,  2.50,  3.00,  1.50, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.50,  2.50,  1.30,  3.20,  2.50,  2.50,  3.00,  2.50,  1.30,  3.30,  2.50,  3.20,  3.00,  2.50,  3.00,  1.50, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.90,  1.40,  0.20,  2.20,  1.40,  2.10,  1.90,  1.40,  1.90,  0.40, 

	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.00,  2.00,  0.80,  3.20,  2.00,  2.00,  2.60,  2.00,  0.80,  3.20,  2.00,  3.20,  2.60,  2.00,  2.60,  1.00, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.00,  2.00,  0.80,  3.20,  2.00,  2.00,  2.60,  2.00,  0.80,  3.20,  2.00,  3.20,  2.60,  2.00,  2.60,  1.00, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.00,  2.00,  0.80,  3.20,  2.00,  2.00,  2.60,  2.00,  0.80,  3.20,  2.00,  3.20,  2.60,  2.00,  2.60,  1.00, 
	2.50,  2.50,  1.30,  3.70,  2.50,  2.50,  3.10,  2.50,  1.30,  3.70,  2.50,  3.70,  3.10,  2.50,  3.10,  1.50, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.50,  2.50,  1.30,  3.70,  2.50,  2.50,  3.10,  2.50,  1.30,  3.70,  2.50,  3.70,  3.10,  2.50,  3.10,  1.50, 
	1.40,  1.40,  0.20,  2.60,  1.40,  1.40,  2.00,  1.40,  0.20,  2.60,  1.40,  2.60,  2.00,  1.40,  2.00,  0.40, 

	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	1.50,  1.50,  0.30,  2.20,  1.50,  1.50,  1.50,  1.50,  0.30,  2.40,  1.50,  2.20,  1.70,  1.50,  1.50,  0.50, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	1.50,  1.50,  0.30,  2.20,  1.50,  1.50,  1.50,  1.50,  0.30,  2.40,  1.50,  2.20,  1.70,  1.50,  1.50,  0.50, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.50,  1.50,  0.30,  2.20,  1.50,  1.50,  1.50,  1.50,  0.30,  2.40,  1.50,  2.20,  1.70,  1.50,  1.50,  0.50, 
	2.00,  2.00,  0.80,  2.70,  2.00,  2.00,  2.00,  2.00,  0.80,  2.90,  2.00,  2.70,  2.20,  2.00,  2.00,  1.00, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.00,  2.00,  0.80,  2.70,  2.00,  2.00,  2.00,  2.00,  0.80,  2.90,  2.00,  2.70,  2.20,  2.00,  2.00,  1.00, 
	0.90,  0.90, -0.30,  1.60,  0.90,  0.90,  0.90,  0.90, -0.30,  1.80,  0.90,  1.60,  1.10,  0.90,  0.90, -0.10, 

	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.60,  2.60,  1.40,  3.80,  2.60,  2.60,  3.20,  2.60,  1.40,  3.80,  2.60,  3.80,  3.20,  2.60,  3.20,  1.60, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.30,  3.30,  2.10,  4.50,  3.30,  3.30,  3.90,  3.30,  2.10,  4.50,  3.30,  4.50,  3.90,  3.30,  3.90,  2.30, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.10,  3.10,  1.90,  4.30,  3.10,  3.10,  3.70,  3.10,  1.90,  4.30,  3.10,  4.30,  3.70,  3.10,  3.70,  2.10, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.10,  3.10,  1.90,  4.30,  3.10,  3.10,  3.70,  3.10,  1.90,  4.30,  3.10,  4.30,  3.70,  3.10,  3.70,  2.10, 
	1.40,  1.40,  0.20,  2.60,  1.40,  1.40,  2.00,  1.40,  0.20,  2.60,  1.40,  2.60,  2.00,  1.40,  2.00,  0.40, 

	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	2.10,  2.10,  0.90,  2.80,  2.10,  2.10,  2.10,  2.10,  0.90,  3.00,  2.10,  2.80,  2.30,  2.10,  2.10,  1.10, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.80,  2.80,  1.60,  3.50,  2.80,  2.80,  2.80,  2.80,  1.60,  3.70,  2.80,  3.50,  3.00,  2.80,  2.80,  1.80, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  2.60,  2.60,  1.40,  3.50,  2.60,  3.30,  2.80,  2.60,  2.60,  1.60, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  2.60,  2.60,  1.40,  3.50,  2.60,  3.30,  2.80,  2.60,  2.60,  1.60, 
	0.90,  0.90, -0.30,  1.60,  0.90,  0.90,  0.90,  0.90, -0.30,  1.80,  0.90,  1.60,  1.10,  0.90,  0.90, -0.10, 

	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	0.70,  0.70, -0.50,  0.80,  0.70,  0.70,  0.30,  0.70, -0.50,  1.20,  0.70,  0.80,  0.30,  0.70,  0.30, -0.30, 
	2.10,  2.10,  0.90,  2.20,  2.10,  2.10,  1.70,  2.10,  0.90,  2.60,  2.10,  2.20,  1.70,  2.10,  1.70,  1.10, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	2.80,  2.80,  1.60,  2.90,  2.80,  2.80,  2.40,  2.80,  1.60,  3.30,  2.80,  2.90,  2.40,  2.80,  2.40,  1.80, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	0.70,  0.70, -0.50,  0.80,  0.70,  0.70,  0.30,  0.70, -0.50,  1.20,  0.70,  0.80,  0.30,  0.70,  0.30, -0.30, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	2.60,  2.60,  1.40,  2.70,  2.60,  2.60,  2.20,  2.60,  1.40,  3.10,  2.60,  2.70,  2.20,  2.60,  2.20,  1.60, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	2.60,  2.60,  1.40,  2.70,  2.60,  2.60,  2.20,  2.60,  1.40,  3.10,  2.60,  2.70,  2.20,  2.60,  2.20,  1.60, 
	0.90,  0.90, -0.30,  1.00,  0.90,  0.90,  0.50,  0.90, -0.30,  1.40,  0.90,  1.00,  0.50,  0.90,  0.50, -0.10, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.70,  1.20,  0.00,  2.00,  1.20,  1.90,  1.70,  1.20,  1.70,  0.20, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  3.10,  2.60,  1.40,  3.40,  2.60,  3.30,  3.10,  2.60,  3.10,  1.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	3.30,  3.30,  2.10,  4.00,  3.30,  3.30,  3.80,  3.30,  2.10,  4.10,  3.30,  4.00,  3.80,  3.30,  3.80,  2.30, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.70,  1.20,  0.00,  2.00,  1.20,  1.90,  1.70,  1.20,  1.70,  0.20, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	3.10,  3.10,  1.90,  3.80,  3.10,  3.10,  3.60,  3.10,  1.90,  3.90,  3.10,  3.80,  3.60,  3.10,  3.60,  2.10, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	3.10,  3.10,  1.90,  3.80,  3.10,  3.10,  3.60,  3.10,  1.90,  3.90,  3.10,  3.80,  3.60,  3.10,  3.60,  2.10, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.90,  1.40,  0.20,  2.20,  1.40,  2.10,  1.90,  1.40,  1.90,  0.40, 

	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.60,  2.60,  1.40,  3.80,  2.60,  2.60,  3.20,  2.60,  1.40,  3.80,  2.60,  3.80,  3.20,  2.60,  3.20,  1.60, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.30,  3.30,  2.10,  4.50,  3.30,  3.30,  3.90,  3.30,  2.10,  4.50,  3.30,  4.50,  3.90,  3.30,  3.90,  2.30, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.10,  3.10,  1.90,  4.30,  3.10,  3.10,  3.70,  3.10,  1.90,  4.30,  3.10,  4.30,  3.70,  3.10,  3.70,  2.10, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.10,  3.10,  1.90,  4.30,  3.10,  3.10,  3.70,  3.10,  1.90,  4.30,  3.10,  4.30,  3.70,  3.10,  3.70,  2.10, 
	1.40,  1.40,  0.20,  2.60,  1.40,  1.40,  2.00,  1.40,  0.20,  2.60,  1.40,  2.60,  2.00,  1.40,  2.00,  0.40, 

	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	2.10,  2.10,  0.90,  2.80,  2.10,  2.10,  2.10,  2.10,  0.90,  3.00,  2.10,  2.80,  2.30,  2.10,  2.10,  1.10, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.80,  2.80,  1.60,  3.50,  2.80,  2.80,  2.80,  2.80,  1.60,  3.70,  2.80,  3.50,  3.00,  2.80,  2.80,  1.80, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  2.60,  2.60,  1.40,  3.50,  2.60,  3.30,  2.80,  2.60,  2.60,  1.60, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  2.60,  2.60,  1.40,  3.50,  2.60,  3.30,  2.80,  2.60,  2.60,  1.60, 
	0.90,  0.90, -0.30,  1.60,  0.90,  0.90,  0.90,  0.90, -0.30,  1.80,  0.90,  1.60,  1.10,  0.90,  0.90, -0.10, 

	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	1.90,  1.90,  0.70,  3.10,  1.90,  1.90,  2.50,  1.90,  0.70,  3.10,  1.90,  3.10,  2.50,  1.90,  2.50,  0.90, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.40,  1.40,  0.20,  2.30,  1.40,  2.10,  1.60,  1.40,  1.40,  0.40, 

	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	1.20,  1.20,  0.00,  1.30,  1.20,  1.20,  0.80,  1.20,  0.00,  1.70,  1.20,  1.30,  0.80,  1.20,  0.80,  0.20, 
	3.00,  3.00,  1.80,  3.10,  3.00,  3.00,  2.60,  3.00,  1.80,  3.50,  3.00,  3.10,  2.60,  3.00,  2.60,  2.00, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	3.60,  3.60,  2.40,  3.70,  3.60,  3.60,  3.20,  3.60,  2.40,  4.10,  3.60,  3.70,  3.20,  3.60,  3.20,  2.60, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	1.20,  1.20,  0.00,  1.30,  1.20,  1.20,  0.80,  1.20,  0.00,  1.70,  1.20,  1.30,  0.80,  1.20,  0.80,  0.20, 
	3.00,  3.00,  1.80,  3.10,  3.00,  3.00,  2.60,  3.00,  1.80,  3.50,  3.00,  3.10,  2.60,  3.00,  2.60,  2.00, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	3.00,  3.00,  1.80,  3.10,  3.00,  3.00,  2.60,  3.00,  1.80,  3.50,  3.00,  3.10,  2.60,  3.00,  2.60,  2.00, 
	3.60,  3.60,  2.40,  3.70,  3.60,  3.60,  3.20,  3.60,  2.40,  4.10,  3.60,  3.70,  3.20,  3.60,  3.20,  2.60, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	3.60,  3.60,  2.40,  3.70,  3.60,  3.60,  3.20,  3.60,  2.40,  4.10,  3.60,  3.70,  3.20,  3.60,  3.20,  2.60, 
	1.40,  1.40,  0.20,  1.50,  1.40,  1.40,  1.00,  1.40,  0.20,  1.90,  1.40,  1.50,  1.00,  1.40,  1.00,  0.40, 

	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	1.70,  1.70,  0.50,  2.40,  1.70,  1.70,  2.20,  1.70,  0.50,  2.50,  1.70,  2.40,  2.20,  1.70,  2.20,  0.70, 
	3.50,  3.50,  2.30,  4.20,  3.50,  3.50,  4.00,  3.50,  2.30,  4.30,  3.50,  4.20,  4.00,  3.50,  4.00,  2.50, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	4.10,  4.10,  2.90,  4.80,  4.10,  4.10,  4.60,  4.10,  2.90,  4.90,  4.10,  4.80,  4.60,  4.10,  4.60,  3.10, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	1.70,  1.70,  0.50,  2.40,  1.70,  1.70,  2.20,  1.70,  0.50,  2.50,  1.70,  2.40,  2.20,  1.70,  2.20,  0.70, 
	3.50,  3.50,  2.30,  4.20,  3.50,  3.50,  4.00,  3.50,  2.30,  4.30,  3.50,  4.20,  4.00,  3.50,  4.00,  2.50, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	3.50,  3.50,  2.30,  4.20,  3.50,  3.50,  4.00,  3.50,  2.30,  4.30,  3.50,  4.20,  4.00,  3.50,  4.00,  2.50, 
	4.10,  4.10,  2.90,  4.80,  4.10,  4.10,  4.60,  4.10,  2.90,  4.90,  4.10,  4.80,  4.60,  4.10,  4.60,  3.10, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	4.10,  4.10,  2.90,  4.80,  4.10,  4.10,  4.60,  4.10,  2.90,  4.90,  4.10,  4.80,  4.60,  4.10,  4.60,  3.10, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  2.40,  1.90,  0.70,  2.70,  1.90,  2.60,  2.40,  1.90,  2.40,  0.90, 

	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	1.90,  1.90,  0.70,  3.10,  1.90,  1.90,  2.50,  1.90,  0.70,  3.10,  1.90,  3.10,  2.50,  1.90,  2.50,  0.90, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.40,  1.40,  0.20,  2.30,  1.40,  2.10,  1.60,  1.40,  1.40,  0.40, 

	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.60,  2.60,  1.40,  3.80,  2.60,  2.60,  3.20,  2.60,  1.40,  3.80,  2.60,  3.80,  3.20,  2.60,  3.20,  1.60, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.30,  3.30,  2.10,  4.50,  3.30,  3.30,  3.90,  3.30,  2.10,  4.50,  3.30,  4.50,  3.90,  3.30,  3.90,  2.30, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.10,  3.10,  1.90,  4.30,  3.10,  3.10,  3.70,  3.10,  1.90,  4.30,  3.10,  4.30,  3.70,  3.10,  3.70,  2.10, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.10,  3.10,  1.90,  4.30,  3.10,  3.10,  3.70,  3.10,  1.90,  4.30,  3.10,  4.30,  3.70,  3.10,  3.70,  2.10, 
	1.40,  1.40,  0.20,  2.60,  1.40,  1.40,  2.00,  1.40,  0.20,  2.60,  1.40,  2.60,  2.00,  1.40,  2.00,  0.40, 

	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	2.10,  2.10,  0.90,  2.80,  2.10,  2.10,  2.10,  2.10,  0.90,  3.00,  2.10,  2.80,  2.30,  2.10,  2.10,  1.10, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.80,  2.80,  1.60,  3.50,  2.80,  2.80,  2.80,  2.80,  1.60,  3.70,  2.80,  3.50,  3.00,  2.80,  2.80,  1.80, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  2.60,  2.60,  1.40,  3.50,  2.60,  3.30,  2.80,  2.60,  2.60,  1.60, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  2.60,  2.60,  1.40,  3.50,  2.60,  3.30,  2.80,  2.60,  2.60,  1.60, 
	0.90,  0.90, -0.30,  1.60,  0.90,  0.90,  0.90,  0.90, -0.30,  1.80,  0.90,  1.60,  1.10,  0.90,  0.90, -0.10, 

	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	0.70,  0.70, -0.50,  0.80,  0.70,  0.70,  0.30,  0.70, -0.50,  1.20,  0.70,  0.80,  0.30,  0.70,  0.30, -0.30, 
	2.10,  2.10,  0.90,  2.20,  2.10,  2.10,  1.70,  2.10,  0.90,  2.60,  2.10,  2.20,  1.70,  2.10,  1.70,  1.10, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	2.80,  2.80,  1.60,  2.90,  2.80,  2.80,  2.40,  2.80,  1.60,  3.30,  2.80,  2.90,  2.40,  2.80,  2.40,  1.80, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	0.70,  0.70, -0.50,  0.80,  0.70,  0.70,  0.30,  0.70, -0.50,  1.20,  0.70,  0.80,  0.30,  0.70,  0.30, -0.30, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	2.60,  2.60,  1.40,  2.70,  2.60,  2.60,  2.20,  2.60,  1.40,  3.10,  2.60,  2.70,  2.20,  2.60,  2.20,  1.60, 
	1.90,  1.90,  0.70,  2.00,  1.90,  1.90,  1.50,  1.90,  0.70,  2.40,  1.90,  2.00,  1.50,  1.90,  1.50,  0.90, 
	2.60,  2.60,  1.40,  2.70,  2.60,  2.60,  2.20,  2.60,  1.40,  3.10,  2.60,  2.70,  2.20,  2.60,  2.20,  1.60, 
	0.90,  0.90, -0.30,  1.00,  0.90,  0.90,  0.50,  0.90, -0.30,  1.40,  0.90,  1.00,  0.50,  0.90,  0.50, -0.10, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.70,  1.20,  0.00,  2.00,  1.20,  1.90,  1.70,  1.20,  1.70,  0.20, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  3.10,  2.60,  1.40,  3.40,  2.60,  3.30,  3.10,  2.60,  3.10,  1.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	3.30,  3.30,  2.10,  4.00,  3.30,  3.30,  3.80,  3.30,  2.10,  4.10,  3.30,  4.00,  3.80,  3.30,  3.80,  2.30, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.70,  1.20,  0.00,  2.00,  1.20,  1.90,  1.70,  1.20,  1.70,  0.20, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	3.10,  3.10,  1.90,  3.80,  3.10,  3.10,  3.60,  3.10,  1.90,  3.90,  3.10,  3.80,  3.60,  3.10,  3.60,  2.10, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.90,  2.40,  1.20,  3.20,  2.40,  3.10,  2.90,  2.40,  2.90,  1.40, 
	3.10,  3.10,  1.90,  3.80,  3.10,  3.10,  3.60,  3.10,  1.90,  3.90,  3.10,  3.80,  3.60,  3.10,  3.60,  2.10, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.90,  1.40,  0.20,  2.20,  1.40,  2.10,  1.90,  1.40,  1.90,  0.40, 

	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.60,  2.60,  1.40,  3.80,  2.60,  2.60,  3.20,  2.60,  1.40,  3.80,  2.60,  3.80,  3.20,  2.60,  3.20,  1.60, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.30,  3.30,  2.10,  4.50,  3.30,  3.30,  3.90,  3.30,  2.10,  4.50,  3.30,  4.50,  3.90,  3.30,  3.90,  2.30, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	1.20,  1.20,  0.00,  2.40,  1.20,  1.20,  1.80,  1.20,  0.00,  2.40,  1.20,  2.40,  1.80,  1.20,  1.80,  0.20, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.10,  3.10,  1.90,  4.30,  3.10,  3.10,  3.70,  3.10,  1.90,  4.30,  3.10,  4.30,  3.70,  3.10,  3.70,  2.10, 
	2.40,  2.40,  1.20,  3.60,  2.40,  2.40,  3.00,  2.40,  1.20,  3.60,  2.40,  3.60,  3.00,  2.40,  3.00,  1.40, 
	3.10,  3.10,  1.90,  4.30,  3.10,  3.10,  3.70,  3.10,  1.90,  4.30,  3.10,  4.30,  3.70,  3.10,  3.70,  2.10, 
	1.40,  1.40,  0.20,  2.60,  1.40,  1.40,  2.00,  1.40,  0.20,  2.60,  1.40,  2.60,  2.00,  1.40,  2.00,  0.40, 

	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	2.10,  2.10,  0.90,  2.80,  2.10,  2.10,  2.10,  2.10,  0.90,  3.00,  2.10,  2.80,  2.30,  2.10,  2.10,  1.10, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.80,  2.80,  1.60,  3.50,  2.80,  2.80,  2.80,  2.80,  1.60,  3.70,  2.80,  3.50,  3.00,  2.80,  2.80,  1.80, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	0.70,  0.70, -0.50,  1.40,  0.70,  0.70,  0.70,  0.70, -0.50,  1.60,  0.70,  1.40,  0.90,  0.70,  0.70, -0.30, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  2.60,  2.60,  1.40,  3.50,  2.60,  3.30,  2.80,  2.60,  2.60,  1.60, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  1.90,  1.90,  0.70,  2.80,  1.90,  2.60,  2.10,  1.90,  1.90,  0.90, 
	2.60,  2.60,  1.40,  3.30,  2.60,  2.60,  2.60,  2.60,  1.40,  3.50,  2.60,  3.30,  2.80,  2.60,  2.60,  1.60, 
	0.90,  0.90, -0.30,  1.60,  0.90,  0.90,  0.90,  0.90, -0.30,  1.80,  0.90,  1.60,  1.10,  0.90,  0.90, -0.10, 

	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	1.90,  1.90,  0.70,  3.10,  1.90,  1.90,  2.50,  1.90,  0.70,  3.10,  1.90,  3.10,  2.50,  1.90,  2.50,  0.90, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.40,  1.40,  0.20,  2.30,  1.40,  2.10,  1.60,  1.40,  1.40,  0.40, 

	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	1.20,  1.20,  0.00,  1.30,  1.20,  1.20,  0.80,  1.20,  0.00,  1.70,  1.20,  1.30,  0.80,  1.20,  0.80,  0.20, 
	3.00,  3.00,  1.80,  3.10,  3.00,  3.00,  2.60,  3.00,  1.80,  3.50,  3.00,  3.10,  2.60,  3.00,  2.60,  2.00, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	3.60,  3.60,  2.40,  3.70,  3.60,  3.60,  3.20,  3.60,  2.40,  4.10,  3.60,  3.70,  3.20,  3.60,  3.20,  2.60, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	1.20,  1.20,  0.00,  1.30,  1.20,  1.20,  0.80,  1.20,  0.00,  1.70,  1.20,  1.30,  0.80,  1.20,  0.80,  0.20, 
	3.00,  3.00,  1.80,  3.10,  3.00,  3.00,  2.60,  3.00,  1.80,  3.50,  3.00,  3.10,  2.60,  3.00,  2.60,  2.00, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	3.00,  3.00,  1.80,  3.10,  3.00,  3.00,  2.60,  3.00,  1.80,  3.50,  3.00,  3.10,  2.60,  3.00,  2.60,  2.00, 
	3.60,  3.60,  2.40,  3.70,  3.60,  3.60,  3.20,  3.60,  2.40,  4.10,  3.60,  3.70,  3.20,  3.60,  3.20,  2.60, 
	2.40,  2.40,  1.20,  2.50,  2.40,  2.40,  2.00,  2.40,  1.20,  2.90,  2.40,  2.50,  2.00,  2.40,  2.00,  1.40, 
	3.60,  3.60,  2.40,  3.70,  3.60,  3.60,  3.20,  3.60,  2.40,  4.10,  3.60,  3.70,  3.20,  3.60,  3.20,  2.60, 
	1.40,  1.40,  0.20,  1.50,  1.40,  1.40,  1.00,  1.40,  0.20,  1.90,  1.40,  1.50,  1.00,  1.40,  1.00,  0.40, 

	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	1.70,  1.70,  0.50,  2.40,  1.70,  1.70,  2.20,  1.70,  0.50,  2.50,  1.70,  2.40,  2.20,  1.70,  2.20,  0.70, 
	3.50,  3.50,  2.30,  4.20,  3.50,  3.50,  4.00,  3.50,  2.30,  4.30,  3.50,  4.20,  4.00,  3.50,  4.00,  2.50, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	4.10,  4.10,  2.90,  4.80,  4.10,  4.10,  4.60,  4.10,  2.90,  4.90,  4.10,  4.80,  4.60,  4.10,  4.60,  3.10, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	1.70,  1.70,  0.50,  2.40,  1.70,  1.70,  2.20,  1.70,  0.50,  2.50,  1.70,  2.40,  2.20,  1.70,  2.20,  0.70, 
	3.50,  3.50,  2.30,  4.20,  3.50,  3.50,  4.00,  3.50,  2.30,  4.30,  3.50,  4.20,  4.00,  3.50,  4.00,  2.50, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	3.50,  3.50,  2.30,  4.20,  3.50,  3.50,  4.00,  3.50,  2.30,  4.30,  3.50,  4.20,  4.00,  3.50,  4.00,  2.50, 
	4.10,  4.10,  2.90,  4.80,  4.10,  4.10,  4.60,  4.10,  2.90,  4.90,  4.10,  4.80,  4.60,  4.10,  4.60,  3.10, 
	2.90,  2.90,  1.70,  3.60,  2.90,  2.90,  3.40,  2.90,  1.70,  3.70,  2.90,  3.60,  3.40,  2.90,  3.40,  1.90, 
	4.10,  4.10,  2.90,  4.80,  4.10,  4.10,  4.60,  4.10,  2.90,  4.90,  4.10,  4.80,  4.60,  4.10,  4.60,  3.10, 
	1.90,  1.90,  0.70,  2.60,  1.90,  1.90,  2.40,  1.90,  0.70,  2.70,  1.90,  2.60,  2.40,  1.90,  2.40,  0.90, 

	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	1.70,  1.70,  0.50,  2.90,  1.70,  1.70,  2.30,  1.70,  0.50,  2.90,  1.70,  2.90,  2.30,  1.70,  2.30,  0.70, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	3.50,  3.50,  2.30,  4.70,  3.50,  3.50,  4.10,  3.50,  2.30,  4.70,  3.50,  4.70,  4.10,  3.50,  4.10,  2.50, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	2.90,  2.90,  1.70,  4.10,  2.90,  2.90,  3.50,  2.90,  1.70,  4.10,  2.90,  4.10,  3.50,  2.90,  3.50,  1.90, 
	4.10,  4.10,  2.90,  5.30,  4.10,  4.10,  4.70,  4.10,  2.90,  5.30,  4.10,  5.30,  4.70,  4.10,  4.70,  3.10, 
	1.90,  1.90,  0.70,  3.10,  1.90,  1.90,  2.50,  1.90,  0.70,  3.10,  1.90,  3.10,  2.50,  1.90,  2.50,  0.90, 

	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	1.20,  1.20,  0.00,  1.90,  1.20,  1.20,  1.20,  1.20,  0.00,  2.10,  1.20,  1.90,  1.40,  1.20,  1.20,  0.20, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.00,  3.00,  1.80,  3.70,  3.00,  3.00,  3.00,  3.00,  1.80,  3.90,  3.00,  3.70,  3.20,  3.00,  3.00,  2.00, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	2.40,  2.40,  1.20,  3.10,  2.40,  2.40,  2.40,  2.40,  1.20,  3.30,  2.40,  3.10,  2.60,  2.40,  2.40,  1.40, 
	3.60,  3.60,  2.40,  4.30,  3.60,  3.60,  3.60,  3.60,  2.40,  4.50,  3.60,  4.30,  3.80,  3.60,  3.60,  2.60, 
	1.40,  1.40,  0.20,  2.10,  1.40,  1.40,  1.40,  1.40,  0.20,  2.30,  1.40,  2.10,  1.60,  1.40,  1.40,  0.40 
};
////////////////////////////////////////////////////

float calculate_energy(int start,int end);
float eHairpin_size(int len);
float eBuldge_size(int len);
float eInternal_size(int len);
float eHairpin(int i, int j);
float eTerminal_mismatch(int i, int j);
float eDangle(int i,int j,int d);
float eInternal(int i,int j,int i2,int j2);
float eTriloop(int i);
float eTetraloop(int i);
float eloop(int branches, int i[],int j[]);
float eMultiloop(int i,int j);
float eStack(int i, int j);

int isPairing(char A, char B) {
    char A1;
	char B1;
	A1 = toupper(A);
	B1 = toupper(B);

	if (A1 == 'A' && B1 == 'U') 
		return 0;
	else if (A1 == 'U' && B1 == 'A') 
		return 3;
	else if (A1 == 'G' && B1 == 'U') 
		return 4;
	else if (A1 == 'U' && B1 == 'G') 
		return 5;
	else if (A1 == 'G' && B1 == 'C') 
		return 2;
	else if (A1 == 'C' && B1 == 'G')
		return 1;
	else
		return -1;
}
float eTriloop(int i){ return 0; }

float eTetraloop(int i){
	int j;
	char *str; 
	str = (char *)malloc(sizeof(char)*6);
	for (j=0; j<6; j++){
		str[j] = str[i+j];
	}

	if (strcmp(str,"GGGGAC") == 0) return -3.00;
	else if (strcmp(str,"GGUGAC") == 0) return -3.00; 
	else if (strcmp(str,"CGAAAG") == 0) return -3.00; 
	else if (strcmp(str,"GGAGAC") == 0) return -3.00; 
	else if (strcmp(str,"CGCAAG") == 0) return -3.00; 
	else if (strcmp(str,"GGAAAC") == 0) return -3.00; 
	else if (strcmp(str,"CGGAAG") == 0) return -3.00; 
	else if (strcmp(str,"CUUCGG") == 0) return -3.00; 
	else if (strcmp(str,"CGUGAG") == 0) return -3.00; 
	else if (strcmp(str,"CGAAGG") == 0) return -2.50; 
	else if (strcmp(str,"CUACGG") == 0) return -2.50; 
	else if (strcmp(str,"GGCAAC") == 0) return -2.50; 
	else if (strcmp(str,"CGCGAG") == 0) return -2.50; 
	else if (strcmp(str,"UGAGAG") == 0) return -2.50; 
	else if (strcmp(str,"CGAGAG") == 0) return -2.00; 
	else if (strcmp(str,"AGAAAU") == 0) return -2.00;
	else if (strcmp(str,"CGUAAG") == 0) return -2.00; 
	else if (strcmp(str,"CUAACG") == 0) return -2.00; 
	else if (strcmp(str,"UGAAAG") == 0) return -2.00; 
	else if (strcmp(str,"GGAAGC") == 0) return -1.50; 
	else if (strcmp(str,"GGGAAC") == 0) return -1.50; 
	else if (strcmp(str,"UGAAAA") == 0) return -1.50; 
	else if (strcmp(str,"AGCAAU") == 0) return -1.50; 
	else if (strcmp(str,"AGUAAU") == 0) return -1.50; 
	else if (strcmp(str,"CGGGAG") == 0) return -1.50; 
	else if (strcmp(str,"AGUGAU") == 0) return -1.50; 
	else if (strcmp(str,"GGCGAC") == 0) return -1.50; 
	else if (strcmp(str,"GGGAGC") == 0) return -1.50; 
	else if (strcmp(str,"GUGAAC") == 0) return -1.50; 
	else if (strcmp(str,"UGGAAA") == 0) return -1.50;
	else return 0;
}

float min(float a,float b){
	if (a <b)
		return a;
	else
		return b;
}

int isBase(char A){
	A = toupper(A);

	if (A == 'A') 
		return 0;
	else if (A == 'U') 
		return 1;
	else if (A == 'G') 
		return 2;
	else if (A == 'C')
		return 3;
	else
		return -1;
}

float eHairpin_size(int len){
	if (len == 0)
		return INF;
	else if(len ==1)
		return INF;
	else if(len ==2)
		return INF;
	else if(len ==3)
		return 5.7;
	else if(len ==4)
		return 5.6;
	else if(len ==5)
		return 5.6;
	else if(len ==6)
		return 5.4;
	else if(len ==7)
		return 5.9;
	else if(len ==8)
		return 5.6;
	else if(len ==9)
		return 6.4;
	else if(len ==10)
		return 6.5;
	else if(len ==11)
		return 6.6;
	else if(len ==12)
		return 6.7;
	else if(len ==13)
		return 6.8;
	else if(len ==14)
		return 6.9;
	else if(len ==15)
		return 6.9;
	else if(len ==16)
		return 7.0;
	else if(len ==17)
		return 7.1;
	else if(len ==18)
		return 7.1;
	else if(len ==19)
		return 7.2;
	else if(len ==20)
		return 7.2;
	else if(len ==21)
		return 7.3;
	else if(len ==22)
		return 7.3;
	else if(len ==23)
		return 7.4;
	else if(len ==24)
		return 7.4;
	else if(len ==25)
		return 7.5;
	else if(len ==26)
		return 7.5;
	else if(len ==27)
		return 7.5;
	else if(len ==28)
		return 7.6;
	else if(len ==29)
		return 7.6;
	else if(len >=30)
		return 7.7/*+(1.079*log(len/30))*/;
	else
		return INF;
}

float eBuldge_size(int len){

	if (len == 0)
		return INF;
	else if(len ==1)
		return 3.8;
	else if(len ==2)
		return 2.8;
	else if(len ==3)
		return 3.2;
	else if(len ==4)
		return 3.6;
	else if(len ==5)
		return 4;
	else if(len ==6)
		return 4.4;
	else if(len ==7)
		return 4.6;
	else if(len ==8)
		return 4.7;
	else if(len ==9)
		return 4.8;
	else if(len ==10)
		return 4.9;
	else if(len ==11)
		return 5;
	else if(len ==12)
		return 5.1;
	else if(len ==13)
		return 5.2;
	else if(len ==14)
		return 5.3;
	else if(len ==15)
		return 5.4;
	else if(len ==16)
		return 5.4;
	else if(len ==17)
		return 5.5;
	else if(len ==18)
		return 5.5;
	else if(len ==19)
		return 5.6;
	else if(len ==20)
		return 5.7;
	else if(len ==21)
		return 5.7;
	else if(len ==22)
		return 5.8;
	else if(len ==23)
		return 5.8;
	else if(len ==24)
		return 5.8;
	else if(len ==25)
		return 5.9;
	else if(len ==26)
		return 5.9;
	else if(len ==27)
		return 6.0;
	else if(len ==28)
		return 6.0;
	else if(len ==29)
		return 6.0;
	else if(len >=30)
		return 6.1/*+(1.079*ln(d))*/;

	else
		return INF;
}

float eInternal_size(int len){
	if (len == 0)
		return INF;
	else if(len ==1)
		return INF;
	else if(len ==2)
		return INF;
	else if(len ==3)
		return INF;
	else if(len ==4)
		return 1.7;
	else if(len ==5)
		return 1.8;
	else if(len ==6)
		return 2.0;
	else if(len ==7)
		return 2.2;
	else if(len ==8)
		return 2.3;
	else if(len ==9)
		return 2.4;
	else if(len ==10)
		return 2.5;
	else if(len ==11)
		return 2.6;
	else if(len ==12)
		return 2.7;
	else if(len ==13)
		return 2.8;
	else if(len ==14)
		return 2.9;
	else if(len ==15)
		return 3.0;
	else if(len ==16)
		return 3.0;
	else if(len ==17)
		return 3.1;
	else if(len ==18)
		return 3.1;
	else if(len ==19)
		return 3.2;
	else if(len ==20)
		return 3.3;
	else if(len ==21)
		return 3.3;
	else if(len ==22)
		return 3.4;
	else if(len ==23)
		return 3.4;
	else if(len ==24)
		return 3.4;
	else if(len ==25)
		return 3.5;
	else if(len ==26)
		return 3.5;
	else if(len ==27)
		return 3.6;
	else if(len ==28)
		return 3.6;
	else if(len ==29)
		return 3.6;
	else if(len >=30)
		return 3.7/*+(1.079*log(len/30))*/;
	else
		return INF;
}

float eHairpin(int i, int j){
	int len;
	float bonus = 0;
	len = j-i-1;
	if (len == 4) bonus = eTetraloop(i);
	return eHairpin_size(len) + eTerminal_mismatch(i,j) + bonus;
}

float eTerminal_mismatch(int i, int j){

	// if i-j basepair; i<j
	if (i+1 > n || j-1 < 1)
		return 0.0; // not sure now

	if (  (str[i] == 'X' && str[j] == 'Y') ||(str[i] == 'X' && str[j] == 'Y')){
		return 0;
	}
	if (  (str[i] == 'A' && str[j] == 'U') ||(str[i] == 'U' && str[j] == 'A')){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return 0.7;
		else if (str[i+1] == 'U' && str[j-1] == 'C')
			return 0.7;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return 0.7;
		else if (str[i+1] == 'U' && str[j-1] == 'U')
			return 0;
		else if (str[i+1] == 'G' && str[j-1] == 'A')
			return -0.4;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return 0.7;
		else if (str[i+1] == 'G' && str[j-1] == 'G')
			return 0.7;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return 0.7;
		else if (str[i+1] == 'C' && str[j-1] == 'A')
			return 0.7;
		else if (str[i+1] == 'C' && str[j-1] == 'C')
			return 0.7;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return 0.7;
		else if (str[i+1] == 'C' && str[j-1] == 'U')
			return 0.7;
		else if (str[i+1] == 'A' && str[j-1] == 'A')
			return 0.7;
		else if (str[i+1] == 'A' && str[j-1] == 'C')
			return 0.7;
		else if (str[i+1] == 'A' && str[j-1] == 'G')
			return -0.4;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return 0.7;
	}

	if ( (str[i] == 'C' && str[j] == 'G') || (str[i] == 'G' && str[j] == 'C')){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return 0;
		else if (str[i+1] == 'U' && str[j-1] == 'C')
			return 0;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return 0;
		else if (str[i+1] == 'U' && str[j-1] == 'U')
			return -0.7;
		else if (str[i+1] == 'G' && str[j-1] == 'A')
			return -1.1;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return 0;
		else if (str[i+1] == 'G' && str[j-1] == 'G')
			return 0;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return 0;
		else if (str[i+1] == 'C' && str[j-1] == 'A')
			return 0;
		else if (str[i+1] == 'C' && str[j-1] == 'C')
			return 0;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return 0;
		else if (str[i+1] == 'C' && str[j-1] == 'U')
			return 0;
		else if (str[i+1] == 'A' && str[j-1] == 'A')
			return 0;
		else if (str[i+1] == 'A' && str[j-1] == 'C')
			return 0;
		else if (str[i+1] == 'A' && str[j-1] == 'G')
			return -1.1;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return 0;
	}

	if ( (str[i] == 'G' && str[j] == 'U') || (str[i] == 'U' && str[j] == 'G')){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return 0.7;
		else if (str[i+1] == 'U' && str[j-1] == 'C')
			return 0.7;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return 0.7;
		else if (str[i+1] == 'U' && str[j-1] == 'U')
			return 0;
		else if (str[i+1] == 'G' && str[j-1] == 'A')
			return -0.4;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return 0.7;
		else if (str[i+1] == 'G' && str[j-1] == 'G')
			return 0.7;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return 0.7;
		else if (str[i+1] == 'C' && str[j-1] == 'A')
			return 0.7;
		else if (str[i+1] == 'C' && str[j-1] == 'C')
			return 0.7;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return 0.7;
		else if (str[i+1] == 'C' && str[j-1] == 'U')
			return 0.7;
		else if (str[i+1] == 'A' && str[j-1] == 'A')
			return 0.7;
		else if (str[i+1] == 'A' && str[j-1] == 'C')
			return 0.7;
		else if (str[i+1] == 'A' && str[j-1] == 'G')
			return -0.4;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return 0.7;
	}
	else
		return INF;
}

float eDangle(int i,int j,int d){

	if ( d == i+1){

		if ( (str[i]=='X' && str[j]=='Y') || (str[i]=='Y' && str[j]=='X')){
			return 0;
		}

		if (str[i]=='A' && str[j]=='U'){
			if(str[d] == 'A')
				return -0.8;
			else if(str[d] == 'C')
				return -0.5;
			else if(str[d] == 'G')
				return -0.8;
			else if(str[d] == 'U')
				return -0.6;
		}
		else if (str[i]=='C' && str[j]=='G'){
			if(str[d] == 'A')
				return -1.7;
			else if(str[d] == 'C')
				return -0.8;
			else if(str[d] == 'G')
				return -1.7;
			else if(str[d] == 'U')
				return -1.2;
		}
		else if (str[i]=='G' && str[j]=='C'){
			if(str[d] == 'A')
				return -1.1;
			else if(str[d] == 'C')
				return -0.4;
			else if(str[d] == 'G')
				return -1.3;
			else if(str[d] == 'U')
				return -0.6;
		}
		else if (str[i]=='G' && str[j]=='U'){
			if(str[d] == 'A')
				return -0.8;
			else if(str[d] == 'C')
				return -0.5;
			else if(str[d] == 'G')
				return -0.8;
			else if(str[d] == 'U')
				return -0.6;
		}
		else if (str[i]=='U' && str[j]=='A'){
			if(str[d] == 'A')
				return -0.7;
			else if(str[d] == 'C')
				return -0.1;
			else if(str[d] == 'G')
				return -0.7;
			else if(str[d] == 'U')
				return -0.1;
		}
		else if (str[i]=='U' && str[j]=='G'){
			if(str[d] == 'A')
				return -0.7;
			else if(str[d] == 'C')
				return -0.1;
			else if(str[d] == 'G')
				return -0.7;
			else if(str[d] == 'U')
				return -0.1;
		}
	}
	else if(d=i-1){

		if ( (str[i]=='X' && str[j]=='Y') || (str[i]=='Y' && str[j]=='X')){
			return 0;
		}
		if (str[i]=='A' && str[j]=='U'){
			if(str[d] == 'A')
				return -0.3;
			else if(str[d] == 'C')
				return -0.1;
			else if(str[d] == 'G')
				return -0.2;
			else if(str[d] == 'U')
				return -0.2;
		}
		else if (str[i]=='C' && str[j]=='G'){
			if(str[d] == 'A')
				return -0.2;
			else if(str[d] == 'C')
				return -0.3;
			else if(str[d] == 'G')
				return 0;
			else if(str[d] == 'U')
				return 0;
		}
		else if (str[i]=='G' && str[j]=='C'){
			if(str[d] == 'A')
				return -0.5;
			else if(str[d] == 'C')
				return -0.3;
			else if(str[d] == 'G')
				return -0.2;
			else if(str[d] == 'U')
				return -0.1;
		}
		else if (str[i]=='G' && str[j]=='U'){
			if(str[d] == 'A')
				return -0.3;
			else if(str[d] == 'C')
				return -0.1;
			else if(str[d] == 'G')
				return -0.2;
			else if(str[d] == 'U')
				return -0.2;
		}
		else if (str[i]=='U' && str[j]=='A'){
			if(str[d] == 'A')
				return -0.3;
			else if(str[d] == 'C')
				return -0.3;
			else if(str[d] == 'G')
				return -0.4;
			else if(str[d] == 'U')
				return -0.2;
		}
		else if (str[i]=='U' && str[j]=='G'){
			if(str[d] == 'A')
				return -0.3;
			else if(str[d] == 'C')
				return -0.3;
			else if(str[d] == 'G')
				return -0.4;
			else if(str[d] == 'U')
				return -0.2;
		}
	}
	return INF;
}

float eInternal(int i,int j,int i2,int j2){
	float part1 = 0;
	float part2 = 0;
	float part3 = 0;
	float part4 = 0;

	int alen;
	int len;
	int len1;
	int len2;

	len1 = i2-i-1; 
	len2 = j-j2-1; 
	len = len1+len2;
	if (len1 > len2) 
		alen = len1 - len2;
	else
		alen = len2 - len1;

	if( len1 == 0 || len2 == 0 ) // it's a buldge
		return eBuldge_size(len);
	else if( len == 2 )
		return 0.80;
	else if ( len1 == 1 && len2 == 2)
		return sint1_2[84*4*isPairing(str[i],str[j])+84*isBase(str[j2+1])+24*isBase(str[i+1])+4*isPairing(str[i2],str[j2])+isBase(str[j2+2])];
	else if ( len1 == 2 && len2 == 1) //rahele: I dstr it by assumption. I was not sure!
		return sint1_2[84*4*isPairing(str[i2],str[j2])+84*isBase(str[i+1])+24*isBase(str[j2+1])+4*isPairing(str[i],str[j])+isBase(str[i+2])];
	else if ( len1 == 2 && len2 == 2)
		return sint2_2[6*16*16*isPairing(str[i],str[j])+16*16*isPairing(str[i2],str[j2])+
			4*16*isBase(str[i+1])+16*isBase(str[i+2])+4*isBase(str[j2+1])+isBase(str[j2+2])];
	else if (len >= 4){
		part1 = eInternal_size(len);
		part2 = eTerminal_mismatch(i,j);
		part3 = eTerminal_mismatch(i2,j2);
		/* check it*/
		part4 = min(3,alen*0.3);
	}
	return part1+part2+part3+part4;
}

float eMultiloop(int i,int j){
	int is[5];
	int js[5];
	int i1,i2,i3,i4;
	int j1,j2,j3,j4;
	float tmpe;
	float density;
	float tmpd;
	float energy;
	float tmp;
	int be = 0;

	is[0] = i;
	js[0] = j;

	for(i1=i+1;i1<j-10;i1++){
		for(j1=i1+4;j1<j-6;j1++){
			if( isPairing(str[i1],str[j1]) < 0) continue;
			is[1] = i1;
			js[1] = j1;
			for(i2=j1+1;i2<j-5;i2++){
				for(j2=i2+4;j2<j-1;j2++){
					if( isPairing(str[i2],str[j2]) < 0) continue;
					is[2] = i2;
					js[2] = j2;
					tmpe = eloop(2,is,js);
					density = tmpe+Es[i1][j1];
					density = density + Es[i2][j2];
					density = density /(j-i+1);
					density = density *((i1-i)+(i2-j1)+(j-j2));
					tmpd = (ELCs[i1][j1]-(sigma*Es[i1][j1])) + (ELCs[i2][j2]-(sigma*Es[i2][j2]));
					density = density+tmpd;
					energy = sigma*(Es[i1][j1]+Es[i2][j2]+tmpe);
					tmp = density+energy;

					if(tmp < ELCm[i][j]){
						ELCm[i][j]=tmp;
						be = tmpe;
						solM[i][j].type = MULTIS;
						solM[i][j].i = i1;
					}

					for(i3=j2+1;i3<j-5;i3++){
						for(j3=i3+4;j3<j-1;j3++){
							if( isPairing(str[i3],str[j3]) < 0) continue;
							is[3] =i3;
							js[3] =j3;
							tmpe = eloop(3,is,js);
							density = ((tmpe+Es[i1][j1]+Es[i2][j2]+Es[i3][j3])/(j-i+1))*((i1-i)+(i2-j1)+(i3-j2)+(j-j3));
							tmpd = (ELCs[i1][j1]-(sigma*Es[i1][j1])) + (ELCs[i2][j2]-(sigma*Es[i2][j2])) + (ELCs[i3][j3]-(sigma*Es[i3][j3]));
							density = density+tmpd;
							energy = sigma*(Es[i1][j1]+Es[i2][j2]+Es[i3][j3]+tmpe);
							tmp = density+energy;

							if(tmp < ELCm[i][j]){
								ELCm[i][j]=tmp;
								be = tmpe;
								//solM[i][j].i = i1;
							}

							for(i4=j3+1;i4<j-5;i4++){
								for(j4=i4+4;j4<j-1;j4++){
									if( isPairing(str[i4],str[j4]) < 0) continue;
									is[4] = i4;
									js[4] = j4;
									tmpe = eloop(4,is,js);
									density = ((tmpe+Es[i1][j1]+Es[i2][j2]+Es[i3][j3]+Es[i4][j4])/(j-i+1))*((i1-i)+(i2-j1)+(i3-j2)+(i4-j3)+(j-j4));
									tmpd = (ELCs[i1][j1]-(sigma*Es[i1][j1])) + (ELCs[i2][j2]-(sigma*Es[i2][j2])) + (ELCs[i3][j3]-(sigma*Es[i3][j3])) + (ELCs[i4][j4]-(sigma*Es[i4][j4]));
									density = density+tmpd;
									energy = sigma*(Es[i1][j1]+Es[i2][j2]+Es[i3][j3]+Es[i4][j4]+tmpe);
									tmp = density+energy;

									if(tmp < ELCm[i][j]){
										ELCm[i][j]=tmp;
										be = tmpe;
										//solM[i][j].i = i1;
									}
								}}}}}}}}
	return be;
}

float eloop(int branches, int i[], int j[]){

	int s_size = 0;
	int  d_size;
	float part1;
	float part2;
	int index;

	d_size = branches;
	s_size = (i[1]-i[0]-1)+(j[0]-j[branches-1]-1);
	for(index = 1; index < branches; index++)
		s_size = s_size +(i[index+1]-j[index]-1);

	part1 = (float)a+(float)s_size*b+(float)d_size*c;

	part2 = 0;

	for(index = 1; index < branches; index++)
		if(j[index]+1 != i[index+1])
			if(j[index]+1 == i[index+1]-1){
				part2 = part2 + min(eDangle(i[index],j[index],j[index]+1),eDangle(i[index+1],j[index+1],i[index+1]-1 ));
			}
			else{
				part2 = part2 + eDangle(i[index],j[index],j[index]+1);
				part2 = part2 + eDangle(i[index+1],j[index+1],i[index+1]-1);
			}

	if(i[0]+1 != i[1])
		if(i[0]+1 == i[1]-1){
			part2 = part2 + min(eDangle(i[0],j[0],i[0]+1),eDangle(i[1],j[1],i[1]-1));
		}
		else{
			part2 = part2 + eDangle(i[0],j[0],i[0]+1);
			part2 = part2 + eDangle(i[1],j[1],i[1]-1);
		}


	if(j[branches-1]+1!=j[0])
		if(j[branches-1]+1 == j[0]-1){
			part2 = part2 + min(eDangle(i[branches-1],j[branches-1],j[branches-1]+1),eDangle(j[0],i[0],j[0]-1));
		}
		else{
			part2 = part2 + eDangle(i[branches-1],j[branches-1],j[branches-1]+1);
			part2 = part2 + eDangle(j[0],i[0],j[0]-1);
		}

	return part1+part2;
}

float eStack(int i, int j){
	// if i-j basepair; i<j
	if (i+1 > n || j-1 < 1)
		return 0.0; // not sure now


	if ((str[i] == 'X' && str[j] == 'Y') || (str[i] == 'Y' && str[j] == 'X')){
		if ((str[i+1] == 'X' && str[j-1] == 'Y') || (str[i+1] == 'Y' && str[j-1] == 'X')){
			return  -10;
		}
		else
			return INF;
	}

	if (str[i] == 'A' && str[j] == 'U'){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return -1.10;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return -2.10;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return -2.20;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return -0.90;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return -1.40;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return -0.60;  
		else 
			return INF;
	}
	else if (str[i] == 'U' && str[j] == 'A'){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return -0.90;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return -2.10;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return -2.40;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return -1.30;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return -1.30;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return -1.00;
		else
			return INF;
	}

	else if (str[i] == 'C' && str[j] == 'G'){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return -2.10;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return -2.40;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return -3.30;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return -2.10;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return -2.10;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return -1.40;
		else 
			return INF;
	}

	else if (str[i] == 'G' && str[j] == 'C'){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return -2.20;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return -3.30;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return -3.40;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return -2.40;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return -2.50;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return -1.50;
		else
			return INF;
	}  

	else if (str[i] == 'U' && str[j] == 'G'){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return -0.60;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return -1.40;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return -1.50;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return -1.00;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return -0.50;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return -0.30;
		else
			return INF;
	}

	else if (str[i] == 'G' && str[j] == 'U'){
		if (str[i+1] == 'U' && str[j-1] == 'A')
			return -1.40;
		else if (str[i+1] == 'G' && str[j-1] == 'C')
			return -2.10;
		else if (str[i+1] == 'C' && str[j-1] == 'G')
			return -2.50;
		else if (str[i+1] == 'A' && str[j-1] == 'U')
			return -1.30;
		else if (str[i+1] == 'U' && str[j-1] == 'G')
			return -1.30;
		else if (str[i+1] == 'G' && str[j-1] == 'U')
			return -0.50;
		else 
			return INF;
	}
	else
		return INF;
}

void initializeTables(int n){
	ELC = new float[n];

	ELCs = new float*[n];
	for(int i = 0; i < n; i++)
		ELCs[i] = new float[n];

	ELCbi = new float*[n];
	for(int i = 0; i < n; i++)
		ELCbi[i] = new float[n];
	
	ELCm = new float*[n];
	for(int i = 0; i < n; i++)
		ELCm[i] = new float[n];
	
	Es = new float*[n];
	for(int i = 0; i < n; i++)
		Es[i] = new float[n];
	
	Ebi = new float*[n];
	for(int i = 0; i < n; i++)
		Ebi[i] = new float[n];
	
	Em = new float*[n];
	for(int i = 0; i < n; i++)
		Em[i] = new float[n];

	traceback_ELC = new data[n];
	
	traceback_multi = new data*[n];
	for(int i = 0; i < n; i++)
		traceback_multi[i] = new data[n];

	traceback_bulge = new data*[n];
	for(int i = 0; i < n; i++)
		traceback_bulge[i] = new data[n];

	traceback_ELCs = new data*[n];
	for(int i = 0; i < n; i++)
		traceback_ELCs[i] = new data[n];

	hairpin_number = new int*[n];
	for(int i = 0; i < n; i++)
		hairpin_number[i] = new int[n];
}

float calculate_energy(int begin, int end){
	int i, j, k, l, h;
	float min = INF;
	vector<vector<float> > Em_dash ( n+1, std::vector<float> ( n+1, 0 ) );
	vector<vector<float> > Em_hat ( n+1, std::vector<float> ( n+1, 0 ) );
	for(int k = 0; k < n+1; k++) 
		Em_dash[k][k] = b;
	
	int len = end-begin+1;
	for(i = 1; i < len; i++){
		for(k = begin; k <= n-i; k++){
			l = k+i;
			min = Es[k][l] + c + eDangle(k, l, k-1) + eDangle(l, k, l+1); // nuk e di si perdoren keto
			float tmp = INF;
			for(h = k; h < l; h++){
				if(Em_dash[k][h] + Em_dash[h+1][l] < tmp)
					tmp = Em_dash[k][h] + Em_dash[h+1][l];
			}
			if (tmp < min)
				min = tmp; 
			Em_dash[k][l] = min;
		}
	}

	min = INF;
	for(k = begin + 1; k < end; k++)
		if ( Em_dash[begin][k] + Em_dash[k+1][end] < min)
			min = Em_dash[begin][k] + Em_dash[k+1][end];
	return a + min;
}

void updateBIEnergyTables(int i, int j){
	// ELCbi and Ebi update
	ELCbi[i][j] = INF;
	Ebi[i][j] = INF;
	float tmp, density;
	if( isPairing(str[i], str[j]) >= 0){
		for(int i_prime = i+1; i_prime < j-1; i_prime++){
			for(int j_prime = i_prime+1; j_prime < j; j_prime++){
				if(isPairing(str[i_prime], str[j_prime]) < 0) continue;
				float eBI_temp = eInternal(i, j, i_prime, j_prime);
				density = 
					(( eBI_temp + Es[i_prime][j_prime])/(j-i+1)) * ((i_prime-i)+(j-j_prime));
				tmp = density + ELCs[i_prime][j_prime] + sigma*eBI_temp;
				if(tmp < ELCbi[i][j]){
					ELCbi[i][j] = tmp;
					Ebi[i][j] = eBI_temp + Es[i_prime][j_prime];
					traceback_bulge[i][j].i = i_prime; // -------- ~~ traceback_bulge ~~ ---------
					traceback_bulge[i][j].j = j_prime; 
				}
			}
		}
	}
}

int updateMEnergyTables(int i, int j){
	ELCm[i][j] = INF;
	Em[i][j] = INF;
	int flag = 0;
	int h, l, keep, k;
	int hairpin_n = -1;
	float density, energy, b_prime;

	if( j-i-1 > 4){ 
		if(isPairing(str[i], str[j]) >= 0){
			// init and find energy
			ELCm_dash = new float*[n+1];
			for(int i = 0; i < n+1; i++)
				ELCm_dash[i] = new float[n+1];

			solM = new data*[n+1];
			for(int i = 0; i < n+1; i++)
				solM[i] = new data[n+1];

			energy = calculate_energy(i, j); // this is the E-hat
			b_prime = energy / (j-i+1);

			for(int m = i; m <= j; m++){ 
				ELCm_dash[m][m]= b_prime + sigma*b;
				solM[m][m].type = MULTIS; // ~~ traceback solM ~~
			}
			// find ELCm min
			float min1, min2;
			for( int l_temp = 1; l_temp < j-i+1; l_temp++){
				for(h = i; h <= j - l_temp; h++){
					l = h + l_temp;
					flag = 2;
					// case 1
					min1 = ELCs[h][l] + sigma*(c + eDangle(h,l,h-1) + eDangle(l,h,l+1));
					// case 2
					min2 = INF;
					for(int hh = h; hh < l; hh++){
						if(ELCm_dash[h][hh] + ELCm_dash[hh+1][l] < min2){
							min2 = ELCm_dash[h][hh] + ELCm_dash[hh+1][l];
							keep = hh;
						}
					}

					if(min1 < min2){ 
						ELCm_dash[h][l] = min1; // ~~ traceback ELCm-dash ~~
						solM[h][l].type = MULTI1;
					}else{
						ELCm_dash[h][l] = min2;
						solM[h][l].type = MULTI2;
						solM[h][l].i = keep;
					}
				}
			}
			// update accordingly
			float min = INF;
			for( k = i+1; k < j; k++){
				if(ELCm_dash[i][k] + ELCm_dash[k+1][j] < min){
					flag = 1;
					min = ELCm_dash[i][k] + ELCm_dash[k+1][j];
					keep = k;
				}
			}

			ELCm[i][j] = (sigma*a) + min;
			traceback_multi[i][j].type = MULTI1;
			traceback_multi[i][j].i = keep;	
		}
	}

	if(flag == 0){
		Em[i][j] = INF;
	}
	else if(flag == 1){
		hairpin_n = 0;
		Em[i][j] = eMultiloop(i, j);
	}

	// release mem
	if(flag > 0){
		for(int i = 0; i < n+1; i++) {
			delete [] ELCm_dash[i];
		}
		delete [] ELCm_dash;

		for(int i = 0; i < n+1; i++) {
			delete [] solM[i];
		}
		delete [] solM;
	}
	return hairpin_n;
}

void updateETables(int i, int j, int hn){
	int flag = -1;
	float min = INF;
	float tmp, density;
	if(isPairing(str[i], str[j]) >= 0){
		// case 2
		tmp = (1+sigma) * eHairpin(i,j);
		if(tmp < min){
			flag = 1;
			min = tmp;
		}
		// case 3
		tmp = INF;
		if(j-i >= 4){
			density=2*((eStack(i, j) + Es[i+1][j-1]) / (j-i+1));
			tmp = density + ELCs[i+1][j-1] + sigma*eStack(i, j);
		}
		if(tmp < min){
			flag = 2;
			min = tmp;
		}
		// case 4
		if(ELCbi[i][j] < min){
			min = ELCbi[i][j];
			flag = 3;
		}
		// case 5
		if(ELCm[i][j] < min){
			min = ELCm[i][j];
			flag = 4;
		}
		
	}

	ELCs[i][j] = min;

	if(flag == -1)
		Es[i][j] = INF;
	else if(flag == 1){
		Es[i][j] = min - 5.7; // taken from densityfold old implementation
		traceback_ELCs[i][j].type = HAIRPIN; // ~~ traceback ~~
		traceback_ELCs[i][j].i = i;
		traceback_ELCs[i][j].j = j;
		hairpin_number[i][j] = 1;
	}
	else if(flag == 2){
		Es[i][j] = (eStack(i, j) + Es[i+1][j-1]);
		traceback_ELCs[i][j].type = STACK;
		traceback_ELCs[i][j].i = i+1;
		traceback_ELCs[i][j].j = j-1;
		hairpin_number[i][j] = hairpin_number[i+1][j-1];
	}
	else if(flag == 3){
		Es[i][j] = Ebi[i][j];
		traceback_ELCs[i][j].type = BULGE;
		traceback_ELCs[i][j].i = traceback_bulge[i][j].i;
		traceback_ELCs[i][j].j = traceback_bulge[i][j].j;
		hairpin_number[i][j] = hairpin_number[traceback_bulge[i][j].i][traceback_bulge[i][j].j];
	}
	else if(flag == 4){
		Es[i][j] = Em[i][j];
		traceback_ELCs[i][j].type = MULTILOOP;
		hairpin_number[i][j] = hn;
	}else{ //} if(flag == 0){
		Es[i][j] = INF;
		traceback_ELCs[i][j].type = MISC;
		traceback_ELCs[i][j].i = 0;
	}
}

void dumpToFile(float** matrix){
	ofstream myfile;
	myfile.open ("energy_file");
	for( int i = 0; i < n+1; i++){
		for( int j = 0; j < n+1; j++){
			myfile << "" << i << " " << j << " " << matrix[i][j] << "\n";
		}
	}
	// myfile << "Writing this to a file.\n";
	myfile.close();
}

void print_matrix(float** A){
	for (int i = 0; i < n+1; i++){
		for (int j = 0; j < n+1; j++){
			// if(i < j)
				cout << A[i][j] << "\t";
			// else cout << "-" << "\t"; 
		}
		cout << endl;
	}
}

/* 	
* Implementation of density fold 
* algorithm for RNA structure prediction. 
*/
int main(int argc, char **argv){
	sigma = atof(argv[1]);
	b = atof(argv[2]);
	cout << "" << sigma << " " << b << endl;
	initializeTables(n+1);

	clock_t begin = clock();

	int i, j, k;
	int hn;
	float min;

	for(k = 1; k < n; k++){
		for(i = 1;i <= n-k; i++){
			j=i+k;
			if(i==1 && j==15)
				printf("test\n");
			// Ebi and ELCbi update
			updateBIEnergyTables(i, j);
			// Em and ELCm update
			hn = updateMEnergyTables(i, j);
			// ELCs and Es update
			updateETables(i, j, hn);	
		}
	}

	clock_t filling_in = clock();
  	double filling_in_clock = double(filling_in - begin) / CLOCKS_PER_SEC;
  	cout << "" << filling_in_clock << endl;
	float tmp;
	ELC[0] = INF;
	for(j = 1; j <= n; j++){
		// case 1
		ELC[j] = ELC[j-1];
		traceback_ELC[j].i = j-1; // ~~ traceback elc ~~
		traceback_ELC[j].type = 0;
		
		// case 2
		for(i=1; i <= j-1; i++){
			tmp = ELC[i-1] + ELCs[i][j];
			if(tmp < ELC[j]){
				ELC[j]=tmp;
				traceback_ELC[j].i = i;
				traceback_ELC[j].type = 1;
			}
		}
	}

	tbELC(n, &result, sigma);
	printStructure(result);
	
	clock_t end = clock();
	double elapsed_secs = double(end - begin) / CLOCKS_PER_SEC;
	cout << "Overall: " << elapsed_secs << endl;
	// print_matrix(ELCs);
	dumpToFile(ELCs);
	return 0;
}

void tbELC(int j, std::vector<std::pair<int, int> >* result, float sigma){
	if(traceback_ELC[j].type == 0){
		tbELC(j-1, result, sigma);
	}else{
		int i = traceback_ELC[j].i;
		if(i > 1 && i <= j-1){
			tbELCs(i, j, result, sigma);
			tbELC(i-1, result, sigma);
		}
	}
}

void tbELCs(int i, int j, std::vector<std::pair<int, int> >* result, float sigma){
	float density;
	float tmp;
	float b_prime;
	float tmpe;
	float energy;
	int h;
	int l;
	data** solM;
	float** ELCml;
		
	if(traceback_ELCs[i][j].type == HAIRPIN){
		result->push_back(std::make_pair(traceback_ELCs[i][j].i, traceback_ELCs[i][j].j));
	}
	else if(traceback_ELCs[i][j].type == STACK){
		result->push_back(std::make_pair(traceback_ELCs[i][j].i, traceback_ELCs[i][j].j));
		tbELCs(i+1, j-1, result, sigma);
	}
	else if(traceback_ELCs[i][j].type == BULGE){
		result->push_back(std::make_pair(traceback_ELCs[i][j].i, traceback_ELCs[i][j].j));
		tbELCs((traceback_ELCs)[i][j].i,(traceback_ELCs)[i][j].j,result, sigma);
	}
	else if(traceback_ELCs[i][j].type == MULTILOOP){
		// result->push_back(std::make_pair(i, j));

		ELCm_dash = new float*[n+1];
		for(int i = 0; i < n+1; i++)
			ELCm_dash[i] = new float[n+1];

		solM = new data*[n+1];
		for(int i = 0; i < n+1; i++)
			solM[i] = new data[n+1];

		energy = calculate_energy(i,j);
		b_prime = energy / (j-i+1);

		for(int m = i; m <= j; m++){ 
			ELCml[m][m]= b_prime + sigma*b;
			solM[h][h].type = MULTIS;
		}

		int keep = -1;
		float min1, min2;
		int hh;
		int k;
		for( int l_temp = 1; l_temp < j-i+1; l_temp++){
			for(h = i; h <= j - l_temp; h++){
				l = h + l_temp;
				min1 = ELCs[h][l] + sigma*(c + eDangle(h,l,h-1) + eDangle(l,h,l+1));
				min2 = INF;
				for(hh = h; hh < l; hh++){
					if(ELCml[h][hh] + ELCml[hh+1][l] < min2){
						min2 = ELCml[h][hh] + ELCml[hh+1][l];
						keep = hh;
					}
				}

				if(min1 < min2){ 
					ELCml[h][l] = min1;
					solM[h][l].type = MULTI1;
				}else{
					ELCml[h][l] = min2;
					solM[h][l].type = MULTI2;
					solM[h][l].i = keep;
				}
			}
		}
		k = (traceback_multi)[i][j].i;

		reportloop(i, j, i+1, k, &solM, result, sigma);
		reportloop(i, j, k+1, j-1, &solM, result, sigma);

		
		for(int i = 0; i < n+1; i++) {
    		delete [] ELCm_dash[i];
		}
		delete [] ELCm_dash;

		for(int i = 0; i < n+1; i++) {
    		delete [] solM[i];
		}
		delete [] solM;
	
	}
}

void reportloop(int i,int j,int k,int h, struct data ***solM, 
	std::vector<std::pair<int, int> >* result, float sigma){
	int m;
	if((solM)[k][h]->type == MULTI1){
		result->push_back(std::make_pair(k, h));
		tbELCs(k+1, h-1, result, sigma);
		return;
	}
	else{
		if (k == h || k > h )
			return ;
		else{
			m = (*solM)[k][h].i;
			reportloop(i, j, k, m, solM, result, sigma);
			reportloop(i, j, m+1, h, solM, result, sigma);
		}
	}
}

void printStructure(const std::vector<std::pair<int, int> >& result){	
	cout << "Printing structure"<< endl;
	for (size_t i = 0; i < result.size(); i++){
		// we actually keep the index values in result. 
		cout << "("<< result[i].first + 1 << ", " << result[i].second + 1 << ") ";
	}
	cout << endl;

	char res[n];
	for(int index = 0; index < n; index++){
		res[index] = '.';
	}
	for(size_t in = 0; in < result.size(); in++){
		// int wee = result[in].first;
		// int s = result[in].second;
		// res[wee] = '(';
		// res[s] = ')';
		res[result[in].first+1] = '(';
		res[result[in].second+1] = ')';
	}
	string r = res;
	cout << r.substr(0, n) << endl;
}