最长上升子序列的回溯 ZOJ 2432

题目大意：

找一组最长上升公共子序列，并把任意一组满足的情况输出出来

最长公共上升子序列不清楚可以先看这篇文章

http://www.cnblogs.com/CSU3901130321/p/4182618.html

然后在这基础上加回溯，我自己一开始利用两个一维数组写回溯，测了很多数据都没问题

但一直给segment fault，网上也看到有人跟我一样说不知道为什么，一维数组的代码主要函数先放在这里留待以后看能否解决，或者有大神帮忙解决

int dp[N] , a[N] , b[N] , rec[N] , fa[N] , src[N] , maxn , cnt;

void LCIS(int m , int n)
{
    memset(dp , 0 , sizeof(dp));
    memset(src , 0 , sizeof(src));
    memset(fa , 0 , sizeof(fa));
    for(int i = 1 ; i<=m ; i++){
        int  k = 0;
        for(int j = 1 ; j<=n ; j++){
            if(a[i] == b[j]){
                if(dp[j] < dp[k] + 1){
                    dp[j] = dp[k] + 1;
                    src[j] = i;
                    fa[i] = src[k];
                }
            }
            if(a[i] > b[j] && dp[k] < dp[j]) k = j;
        }
    }

    maxn = 0 , cnt = 0;
    int s;
    for(int i = 1 ; i <= n ; i++)
    {
        if(maxn < dp[i])
            maxn = dp[i] , s = src[i];
    }
    rec[cnt++] = s;
    while(fa[s]){
        rec[cnt++] = fa[s];
        s = fa[s];
    }
}

后来自己改成了二维数组来回溯

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int N = 1005;
#define max(a,b) a>b?a:b
int dp[N] , a[N] , b[N] , rec[N] , maxn , cnt;
int s[N][N]; //用来回溯,记录前一次出现最大的j的位置，因为那个位置一定是会出现b[pos] = 某个a[i]的

void TraceBack(int i , int j)
{
    if(i < 1 || j < 1) return ;
   // cout<<"here: "<<i<<" "<<s[i][j]<<endl;
    if(s[i][j] >= 0){
        rec[cnt++] = i;

        TraceBack(i-1 , s[i][j]);
    }else TraceBack(i-1 , j);
}

void LCIS(int m , int n)
{
    memset(dp , 0 , sizeof(dp));
    memset(s , -1 , sizeof(s));
    for(int i = 1 ; i<=m ; i++){
        int  k = 0;
        for(int j = 1 ; j<=n ; j++){
            if(a[i] == b[j]){
                if(dp[j] < dp[k] + 1){
                    dp[j] = dp[k] + 1;
                    s[i][j] = k;//记录上一次出现在最长子序列中能够进行匹配的j的位置
                }
            }
            if(a[i] > b[j] && dp[k] < dp[j]) k = j;
        }
    }

    maxn = 0 , cnt = 0;
    int pos ;
    //我自己写的函数原因，所以必须找到第一个出现最大值的位置pos，保证在这个位置会出现某个a[i]与其匹配
    /*这里从后往前找和从前往后找效果一样，但是输出的序列可能不同，
    但是题目要求只输出一种情况所以也没问题,方向找，输出的正好是样例的结果
    for(int i = 1 ; i<=n ; i++) 也确实AC了没问题
    */
    for(int i = n ; i >= 1 ; i--)
    {
        if(maxn < dp[i])
            maxn = dp[i] , pos = i;
    }
    TraceBack(m , pos);
}

int main()
{
    int m , n , T;
    scanf("%d" , &T);
    while(T--){
        scanf("%d" , &m);
        for(int i = 1 ; i<=m ; i++)
            scanf("%d" , a+i);

        scanf("%d" , &n);
        for(int i= 1 ; i<=n ; i++)
            scanf("%d" , b+i);

        LCIS(m , n);

        printf("%d
" , maxn);
        for(int i = cnt - 1 ; i>=0 ; i--)
            printf("%d " , a[rec[i]]);
        printf("
");
        if(T>0) puts("");
    }
    return 0;
}