Exemplo n.º 1
0
	LCP_DAC_VAR::LCP_DAC_VAR(TextIndex *csa, char *text, size_t length){
		lcp_type = DAC_VAR;
		n = length;	
		uint *lcp;
		uint *frec_acu = new uint[n];
		uint max_lcp = 0;
		ushort *kvalues;
		uint nkvalues;
		/*create the LCP array and frec_acu array*/
		lcp = create_lcp(csa,text,n,32);
		for(uint i=0; i<n; i++)
			frec_acu[i]=0;

		frec_acu[0]++;
		for(uint i=0; i < n; i++){
			frec_acu[lcp[i]]++;
			if(max_lcp<lcp[i])
				max_lcp = lcp[i];
		}
		for(uint i=1; i<=max_lcp; i++)
			frec_acu[i] = frec_acu[i] + frec_acu[i-1];

		/*calculate the kvalues*/
		kvalues = NULL;
		calc_kvalues(max_lcp, frec_acu, &kvalues, &nkvalues);
		rep = new factorization_var(lcp, n, kvalues, nkvalues);
		next_p_var = new size_t[(uint)rep->nLevels -1];
		delete [] kvalues;
		delete [] frec_acu;
		delete [] lcp;
	}
Exemplo n.º 2
0
	LCP_Sad::LCP_Sad(TextIndex *csa, char *text, size_t n, size_t op_rs) {
		if(op_rs!=BRW32_HDR && op_rs!=DARRAY_HDR) {
			cout << "Error: op_rs must be BRW32_HDR or DARRAY_HDR\n" << endl;
			exit(1);
		}
		lcp_type = SAD_GON_OS;
		U_type = op_rs;
		uint *lcp = create_lcp(csa,text,n,32);
		uint pos =0;
		uint lcp_prev, lcp_act;
		long long nb = 1;
		nb = (nb*2*n+W-1)/W;
		if(nb > MAXINT) {
			cout << "Memory limit excess (in LCP)" << endl;
			exit(1);
		}
		uint *S = new uint[(uint)nb];
		for(uint i=0; i < (uint)nb; i++)
			S[i]=0;
		cds_utils::bitset(S,0);

		lcp_prev = lcp[csa->getISA(0)];
		pos += lcp_prev +1;
		//calculate S and create U
		for(uint i=1; i<n ;i++) {
			cds_utils::bitset(S,pos);
			lcp_act = lcp[csa->getISA(i)];
			pos+= lcp_act - lcp_prev + 2;
			lcp_prev = lcp_act;
		}
		cds_utils::bitset(S,pos);
		if(op_rs==BRW32_HDR)
								 //(bitarray, length, factor)
			U = new BitSequenceRG(S, pos+1, 4);
		else
			U = new BitSequenceDArray(S, pos+1);
		U_length = (size_t)(pos+1);
		delete [] S;
		delete [] lcp;
	}