我在SPOJ上尝试了一个问题,我们必须简单地找到给定数组 A的最长递增子序列的长度。
我已经使用动态编程 O(n^2) 算法解决了这个问题,并且该解决方案被接受了。这是被接受的代码:
void LIS(int *A,int A_Length)
{
int Seq[MAX];
for(int i=0;i<A_Length;++i)
{
int maxima=0;
for(int j=0;j<i;++j)
{
if(A[i]>A[j])
{
maxima=max(Seq[j],maxima);
}
}
Seq[i]=maxima+1;
//cout<<Seq[i]<<endl;
}
cout<<*max_element(Seq,Seq+A_Length)<<endl;
}
但是当我尝试使用第二种方法(LINK)解决它时,它是::
A simple way of finding the longest increasing subsequence is
to use the Longest Common Subsequence (Dynamic Programming) algorithm.
[1]Make a sorted copy of the sequence A, denoted as B. O(nlog(n)) time.
[2]Use Longest Common Subsequence on with A and B. O(n2) time.
,我得到了错误的答案。
这是我的 C++ 代码
//Global Variable
int A[100],B[100];
int DP[100][100];
//This function Finds the Longest common subsequce of Array A[1,2,3...,N] and B[1,2,3...,N]
void LIS(int N)
{
sort((B+1),(B+1)+N);//STL SORT sort from index 1 to N of Array B.
int i,j;
//Base Cases
for(i=0;i<=N;++i)
DP[i][0]=0;
for(j=0;j<=N;++j)
DP[0][j]=0;
for(i=1;i<=N;++i)
{
for(j=1;j<=N;++j)
{
if(A[i]==B[j])
DP[i][j]=DP[i-1][j-1]+1;
else
DP[i][j]=max(DP[i-1][j],DP[i][j-1]);
}
}
printf("%d\n",DP[N][N]);
}
int main()
{
int N,i;
scanf("%d",&N);
for(i=1;i<=N;++i)
{
scanf("%d",&A[i]);
B[i]=A[i];
}
LIS(N);
return 0;
}
我不知道为什么我得到错误的答案。你能帮我找到错误吗?或者站点中给出的 LIS by LCS Algorithm不正确?