CDOJ1300 LCIS

这道题当初想当然的认为LCIS只要求出LCS记录路径然后求LCS的LIS就可以，发现根本不是这回事儿。。

这题很明显的dp，但是转移方程写的各种戳。。。。- -

 

最裸的dp方程：

dp[i][j]=max{dp[ki][kj]}+1; ki<i&&kj<j&&a[i]==a[ki]b[j]==b[kj]

dp[i][j]表示以a[i]和b[j]结尾的LCIS

这样转移正确性就不做证明了，很明显，但是转移费用太高了，是O(n^2)的转移，这样总的时间复杂度就变成了O(n^4)。对于1000的数据不可接受

改进版本：

dp[i][j]=max{dp[ki][j-1]}+1; ki<i&&a[i]==b[j]

这时候dp[i][j]表示以a[i]结尾，b[j]或其之前元素结尾的LCIS

这样转移就可以写成O(n)的转移，加上O(n^2)的状态，算法复杂度是O(n^3)

注意，这时候还有一个很显著的优化，记录一个下标p使得a[p]<=b[j]并且dp[p][j-1]是dp[0][j-1]到dp[i-1][j]中最大的，这样转移就是O(1)的

dp[i][j]=dp[p][j-1]+1;

至此，我们的O(n^2)的dp方程就已经完成，总的复杂度就是O(n^2)，对于1000的数据量正好适合

 

 

我的代码：

#include <cstdio> #include <cstring> #include <cstdlib> #include <algorithm> using namespace std; const int MAX=1100; int dp[MAX][MAX],a[MAX],b[MAX]; int get() { char ch; int flag = 0, tmp = 0; for (ch = getchar(); ch < 48 || ch > 57; ch = getchar()) if (ch == int('-')) break; if (ch == int('-')) flag = 1; else tmp = int(ch) - 48; for (ch = getchar(); 48 <= ch && ch <= 57; ch = getchar()) tmp = tmp * 10 + int(ch) - 48; return (flag) ? -tmp : tmp; } int main(){ int t; int lena,lenb; t=get(); while(t--){ lena=get(); lenb=get(); for(int i=1;i<=lena;i++)a[i]=get(); for(int i=1;i<=lenb;i++)b[i]=get(); for(int i=0;i<=lena;i++)dp[i][0]=0; for(int j=1;j<=lenb;j++){ dp[0][j]=0; int p=0; for(int i=1;i<=lena;i++){ dp[i][j]=dp[i][j-1]; if(a[i]==b[j]&&dp[i][j]<dp[p][j-1]+1){ dp[i][j]=dp[p][j-1]+1; }else if(a[i]<b[j]&&dp[i][j]>dp[p][j-1]){ p=i; } } } int ret=0; for(int i=0;i<=lena;i++){ ret=max(ret,dp[i][lenb]); } printf("%d/n",ret); } return 0; }