TYVJ 1071 LCIS （蓝书刷题记录）

问题描述：熊大妈的奶牛在小沐沐的熏陶下开始研究信息题目。小沐沐先让奶牛研究了最长上升子序列，再让他们研究了最长公共子序列，现在又让他们要研究最长公共上升子序列了。小沐沐说，对于两个串A，B，如果它们都包含一段位置不一定连续的数字，且数字是严格递增的，那么称这一段数字是两个串的公共上升子串，而所有的公共上升子串中最长的就是最长公共上升子串了。奶牛半懂不懂，小沐沐要你来告诉奶牛什么是最长公共上升子串。不过，只要告诉奶牛它的长度就可以了。
题解：此题为LIS和LCS的结合，所以容易列出状态转移方程：f[i][j]表示a中的前i 个 和 b 中以 b[j] 结尾的最大LCIS。
f[i][j]=f[i-1][j], a[i]!=b[j]
f[i][j]=max(f[i-1][k])+1, a[i]=b[j],0<=k<j,b[j]>b[k]
所以容易得到一个O(n^3) 的算法
code：

 for(int i=1;i<=n;i++)
   for(int j=1;j<=m;j++)
    if(a[i]==b[j])
    {
      for(int k=0;k<j;k++)
       if(a[i]>b[k]) f[i][j]=max(f[i][j],f[i-1][k]+1);
	}
    else f[i][j]=f[i-1][j];
  for(int j=1;j<=m;j++) ans=max(f[n][j],ans);

考虑每次的最优的f[i-1][k]是确定的，对于每次的i,j 是可以提前处理出来k，不必要每次循环找k。
即定义一个val 表示最优k值，在j循环中更新当前最优的k值。
code:

//val 记录值，直接用于决策值的更新
  for(int i=1;i<=n;i++)
  {
  	int val=0;
  	for(int j=1;j<=m;j++)
  	{
  	 if(a[i]==b[j]) f[i][j]=val+1;	
      else f[i][j]=f[i-1][j];
  	 if(a[i]>b[j]) val=max(val,f[i-1][j]);
	}
  }
 //val 记录位置，用于状态中
   for(int i=1;i<=n;i++)
  {
  	int val=0;
  	for(int j=1;j<=m;j++)
  	{
  	 if(a[i]==b[j]) f[i][j]=f[i-1][val]+1;	
      else f[i][j]=f[i-1][j];
  	 if(a[i]>b[j] && f[i-1][j]>f[i-1][val]) val=j;
	}
  }

本题还有一个拓展，即输出方案。每次记录下当前状态是由哪里推过来的。最后递归输出方案即可。
code：

  for(int i=1;i<=n;i++)
  {
  	int val=0;
  	for(int j=1;j<=m;j++)
  	{
  	 if(a[i]==b[j]) f[i][j]=f[i-1][val]+1,g[i][j]=t;	
      else f[i][j]=f[i-1][j];
  	 if(a[i]>b[j] && f[i-1][j]>f[i-1][val]) val=j;
	}
  }

这样就可以优化到O(n^2)

另外，此题中的每次推导，只根据之前的状态，即f[i][j]由f[i-1][k] 而k是小于j的，所以可以优化掉一维。
code：

 for(int i=1;i<=n;i++)
  {
   int t=0;
   for(int j=1;j<=m;j++)
   {
    if(a[i]==b[j]) f[j]=f[t]+1,g[j]=t; 	
    if(a[i]>b[j] && f[t]<f[j]) t=j;
   }
  }

这道题转移部分的优化告诉我们，在实现状态转移方程时，要注意观察决策集合的范围随着状态的变化情况。对于“决策集合中的元素只增多不减少”的情形，就可以像本题一样维护一个变量来记录决策信息，避免重复扫描(lyd)