poj 1147 BWT压缩方法

题目描述

问题源自BWT压缩方法。详见http://acm.pku.edu.cn/JudgeOnline/problem?id=1147。

该问题就是BWT压缩方法的逆过程。关于BWT压缩方法见http://en.wikipedia.org/wiki/Burrows%E2%80%93Wheeler_transform。

 

解题过程

有的时候感觉智商不够用。如果没有答案估计我是想不出答案了。

我解题过程写了3套代码。

第一套代码

按照wiki上说的建立一个表，问题的规模是n=3000，表的规模是n*n，所以空间满足需求。这里面有一个排序的过程，我用STL的sort，就是快速排序，问题的规模是n=3000，每次排序平均需比较1000左右的字符，快速排序的复杂度为n*lgn，所以总的时间为O(n*n*lgn*n/2)。所以可想而知，CPU越要计算10^10，当然是超时了。

第二套代码

发现排序的时候有个规律，在第i步之前，原来的表是从小到大排序好的，利用这个规律改进排序方法，排序用O(n)就可以了，所以总共的时间复杂度是O(n*n), 最后空间是9000k，时间是900ms.

第三套代码

继而发现表都不用建立，那么填充表的过程也可以省略，最终空间是180k，时间是100ms.

改进所带来的乐趣，只有努力AC的人才能体会！没有最快，只有更快!

 

代码

#include <cstdio>#include <cmath>#include <algorithm>#include <float.h>using namespace std;#define M 3004int num;char g_in[M];int g_mapIdx[M];void CalcMap(){int idx=0;for (int i=0; i<num; i++){if (g_in[i]==0) {g_mapIdx[idx++] = i;}}for (int i=0; i<num; i++){if (g_in[i]==1) {g_mapIdx[idx++] = i;}}}int main(){ //freopen("ReadMe.txt", "r", stdin);int *pIdxSrc = new int[M];int *pIdxDst = new int[M];int nTmp;scanf("%d", &num);for (int i=0; i<num; i++){scanf("%d", &nTmp);g_in[i] = nTmp;pIdxSrc[i] = i;}CalcMap();for (int i=num-1; i>=0; i--){// sortfor (int i=0; i<num; i++){pIdxDst[i] = pIdxSrc[g_mapIdx[i] ];}printf("%d ", g_in[pIdxDst[0] ]);swap(pIdxSrc, pIdxDst);}printf("/n");delete []pIdxSrc;delete []pIdxDst; return 0;}