数组最长递增子序列

来源:互联网 发布:动漫壁纸软件 编辑:程序博客网 时间:2024/05/22 15:05

这个问题和矩阵相乘比较相似,但因为它只需要计算A[1...j]的代价,因此比矩阵相乘少了一层循环,矩阵相称的复杂度是n的3次方,这个问题的复杂度是n平方的。要把数据打印出来,破费了一点周折,需要好好理解,其中递归退出的条件是,数组记录的位置和自身的位置相同,但这个位置的数据也要输出,因此在else语句中输出最前面的哪个数据,同时最后一个数据需要保存在临时变量,压入栈中,待递归返回时输出数据,整个编写思路模拟矩阵相乘的思路。程序如下:

// LIS[1,j] = max{LIS[1,k] + 1, 1}#include <stdio.h>#include <vector>#include <boost/scoped_ptr.hpp>void Print(int value) {  printf("%d ", value);} template<typename T, typename VisitFun>void PrintSubarray(const std::vector<T>& array, const size_t* array_offset, size_t position, VisitFun visit_fun) {  if (position != array_offset[position]) {    size_t to_print = position;    PrintSubarray(array, array_offset, array_offset[position], visit_fun);    visit_fun(array[to_print]);  } else {    visit_fun(array[position]);  }}template<typename T>void LongestIncreaseSubarray(const std::vector<T>& array) {  size_t max = 0;  size_t max_position = 0;  boost::scoped_ptr<size_t> max_length_buffer(new size_t[array.size()]);  size_t* max_length = max_length_buffer.get();  boost::scoped_ptr<size_t> array_offset_buffer(new size_t[array.size()]);  size_t* array_offset = array_offset_buffer.get();  for (size_t length = 1; length <= array.size(); ++length) {    max_length[length - 1] = 1;    array_offset[length -1] = length - 1;    for (size_t k = 0; k < length; ++k) {      if (array[k] < array[length - 1] && max_length[k] + 1 > max_length[length - 1]) {          max_length[length - 1] = max_length[k] + 1;          array_offset[length - 1] = k;          if (max_length[length - 1] > max) {            max = max_length[length - 1];            max_position = length -1;            printf("%zd\n", max_position);          }      }    }    printf("--%zd\n", max_length[length - 1]);        }  for (int i = 0; i < array.size(); ++i) {    printf("%zd ", array_offset[i]);  }  printf("\n");  PrintSubarray(array, array_offset, max_position, Print);}int main(int argc, char** argv) {  int array_buffer[] = {1, 0, 1,4,5,2,6,1,7,2,8};  std::vector<int> array(array_buffer, array_buffer + sizeof(array_buffer) / sizeof(int));  LongestIncreaseSubarray(array);}


原创粉丝点击