hrbust 1598 求最长公共子序列

当数组长度p,q范围高达10000所以O(pq)的算法会超时。注意到序列中所有元素各不相同，因此可以把A串重新编号为1~p，例如：

A={1 7 5 4 8 3 2},B={1 4 3 5 6 2 8 7}  把A重新编号后:

A={1 2 3 4 5 6 7},B={1 4 6 3 0 7 5 2} 其中0表示在A中没有出现过，因为没有出现过所以就不用考虑了，肯定不在最长公共子序列中。对于新的A,B序列只需求B的最长递增子序列即可。因为n的范围较大，所以需要用到O(nlogn)的方法求最长递增子序列(二分)。

 

#include<stdio.h>
#include<string.h>

const int MAX = 10001;

int main(void)
{
    int t;
    int p;
    int q;
    int A[MAX];
    int B[MAX];
    int res;
    scanf("%d", &t);

    while (t--)
    {
        int x;
        int provi[MAX];
        scanf("%d%d", &p, &q);

        memset(provi, 0, sizeof(provi));

        for (int i=1; i<=p; i++)
        {
            scanf("%d", &x);
            provi[x] = i;
        }
        int k = 0;
        for (int i=1; i<=q; i++)//为什么这样写自己的好好想想写写看
        {
            scanf("%d", &x);
            if (provi[x] != 0)
            {
                A[k++] = provi[x];
            }
        }
        B[0] = A[0];
        res = 1;
        for (int i=1; i<k; i++)
        {
            int left = 0;
            int right = res;
            while (left < right)
            {
                int mid = (left + right) / 2;
                if (A[i] > B[mid])
                {
                    left = mid + 1;
                }
                else
                {
                    right = mid;
                }
            }
            B[left] = A[i];
            if (left >= res)
            {
                res++;
            }
        }
        printf("%d\n", res);
    }
    return 0;
}