DP最长子序列子串全总结--GYM 101291 A和M

一.子序列和子串

   子序列中间不一定要连续，子串一定要连续

 

二.最长上升子序列(LIS)

①算法：（今天VJ的A题，只不过数字换成了字母，问法也有改变）

for (int i = 0; i < a.length(); ++i)
    {
        dp[i] = 1;
        for (int j = 0; j < i; ++j)
        {
            if (a[j] < a[i] && dp[j] + 1 > dp[i])
                dp[i] = dp[j] + 1;
        }
    }

    int ans = *max_element(dp,dp+a.size());//求最大值简便方法，省一个for循环

②nlogn算法：转自

我们再开一个数组low，low始终是有序的

（贪心）对于一个上升子序列，显然其结尾元素越小，越有利于在后面接其他的元素，也就越可能变得更长。

（二分）对替换（不是插入！）位置进行二分，确保替换后仍然有序

比如2  3 4，下一个要放入1则替换2，low变成1 3  4

const int maxn =300003,INF=0x7f7f7f7f;
int low[maxn],a[maxn];
int n,ans;
int binary_search(int *a,int r,int x)
//二分查找，返回a数组中第一个>=x的位置 
{
    int l=1,mid;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(a[mid]<=x)
            l=mid+1;
        else 
            r=mid-1;
    }
    return l;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) 
    {
        scanf("%d",&a[i]); 
        low[i]=INF;//由于low中存的是最小值，所以low初始化为INF 
    }
    low[1]=a[1]; 
    ans=1;//初始时LIS长度为1 
    for(int i=2;i<=n;i++)
    {
        if(a[i]>=low[ans])//若a[i]>=low数组最大的那个（注意先ans+再接），直接把a[i]接到后面 
            low[++ans]=a[i];
        else //否则找到low中第一个>=a[i]的位置，替换之 ，如 23来个1，换成13而不接到前面换成123
            low[binary_search(low,ans,a[i])]=a[i];
    }
    printf("%d\n",ans);//输出答案 
    return 0;
}

 

三.最长公共子序列（LCS）

复杂度，转自，子串只用对满足a[i]==b[j]进行操作。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
char a[1005], b[1005];
int dp[1005][1005];

void lcs(char a[], char b[]){
	int len1 = strlen(a);
	int len2 = strlen(b);
	memset(dp, 0, sizeof(dp));
	for(int i=1; i<=len1; i++)
		for(int j=1; j<=len2; j++){
			if(a[i-1]==b[j-1]) 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[len1][len2]);
}

void print(char a[], char b[]){
	stack <char> s;
	int len1 = strlen(a);
	int len2 = strlen(b);
	while(dp[len1][len2]){
		if(dp[len1][len2] == dp[len1-1][len2]) len1--;
		else if(dp[len1][len2] == dp[len1][len2-1]) len2--;
		else {
			len1--, len2--;
			s.push(a[len1]);
		}
	}
	while(!s.empty()){
		printf("%c", s.top());
		s.pop();
	}
	puts("");
}

int main(){
	while(~scanf("%s %s", a, b)){
		lcs(a, b);
		print(a, b);
	}
}

 

四.最大连续子序列和

HDU1003变式题（也是数字改成字母）

当前和为负数的话，舍去之前的数，在当前位置开始重新计算。（可以定义两个下标，方便记录）

for(int i=1;i<n;i++)
	{
		if(dp[i-1]<0)
		dp[i]=a[i];
		else
		dp[i]=dp[i-1]+a[i];
	}
	int ans=0;
	for(int i=0;i<n;i++)
	{
		ans=max(ans,dp[i]);
	}

 

五.DP无后效性与最长交替子序列

今天碰到一个最长交替子序列的题，比如21342  最长交替子序列为2132  长度为4，其实只用for一遍求拐点（大于前一个以及后一个的值），两个端点不算，有一个拐点就++即可。

以为一定可以DP搞出来，但是忽略了DP成立的条件，必须无后效性，即当前求得的结果，不会影响下面的结果，换句话说，每个将来的状态，只受当前状态影响，而与当前的前面那些状态无关。

比如这道题如果想用dp做，因为前面递增，后面一定要递减，反之亦然，那如果用dp[i]记录当前求得的结果，后面没办法得知应该增还是减