【模板时间】◆模板·III◆ 单调子序列
◆模板·III◆ 单调子序列
以前只知道DP用 O(n2) 的做法,现在才发现求单调子序列方法好多……
◇ 模板简述
单调子序列包括 升序/降序/非升序/非降序 子序列。主要题型如下:
①在原串中找到一个最长的单调子序列;
②将原串分解为若干个单调子序列;
③通过修改元素使原串变为单调序列。
Tab: 子序列在原串中可以断开,也就是说若原串为A{a[1]~a[n]},则其子序列可以是 A'{a[b[1]]~a[b[n]]}满足 b[1]<b[2]<...<b[n]
方法的确很多,包括STL(lower_bound),DP……接着讲吧……
◇ 求原串中最长的单调子序列
【OpenJudge 1759】 +传送门+
一道非常经典的模板题——最长上升子序列。由于n最大1000,O(n2) 的算法是可以通过的。
解决这类问题的基础算法是DP。对于第i个数,我们可以找到一条以 i 结尾的上升子序列,记以 i 结尾的最长上升子序列的长度为 dp[i],也就是我们的DP状态。
如何转移?
对于第 i 个数,若存在第j个数(j<i)小于第i个数,则第i个数可以继续连接在以j结尾的最长上升子序列上,因此转移可以写成下列:
没有多大难度……其实只是数据规模不算大!
先附上代码(一个在 n≤3000 的情况下不会超时的版):
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 /*Lucky_Glass*/ 2 #include<cstdio> 3 #include<algorithm> 4 using namespace std; 5 int A[1005],dp[1005],n,ans=1; 6 int DP(int x) 7 { 8 if(dp[x]) return dp[x]; 9 dp[x]=1; //只有 第i个元素 本身 10 for(int i=0;i<x;i++) 11 if(A[i]<A[x]) 12 dp[x]=max(dp[x],DP(i)+1); 13 return dp[x]; 14 } 15 int main() 16 { 17 scanf("%d",&n); 18 for(int i=0;i<n;i++) scanf("%d",&A[i]); 19 dp[0]=1; 20 for(int i=n-1;i>=1;i--) //枚举对于以每一个元素为结尾的序列 21 ans=max(DP(i),ans); 22 printf("%d\n",ans); 23 return 0; 24 }
难度↑【POJ 1631】Bridging signals +传送门+
只是数据规模上升到40000了……O(n2)算法会炸……我们需要一种比DP更优秀的算法(DP有时并不是最优算法)。
不难发现,我们有时会沿用上一个最长上升子序列的前部分(也就是小于等于当前元素的最大元素的后面)。所以我们可以把以当前元素结尾的最长上升子序列给存下来(seq[])。由于seq[]存储的是一个单调序列,所以我们可以在seq里二分查找(lower_bound)!
具体怎么做?
先初始化seq的所有值为正无穷(memset(seq,0x3f,sizeof seq)),这样才能在第一次查找时找到seq[1]。
若当前要求以num[i]结尾的最长上升子序列。先通过lower_bound找到seq[1~n-1]中小于等于num[i]的元素 seq[pos-1] (如果找不到此元素,则应找到seq[1]),然后将 seq[pos] 替换为num[i],此时seq[1~pos]则是以 num[i] 结尾的最长上升子序列,且序列长度为 pos,此时就可以更新答案了。
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 /*Lucky_Glass*/ 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 using namespace std; 6 const int MAXN=40000; 7 int num[MAXN+5],seq[MAXN+5]; 8 int main() 9 { 10 int T;scanf("%d",&T); 11 while(T--) 12 { 13 int n;scanf("%d",&n); 14 for(int i=0;i<n;i++) scanf("%d",&num[i]); 15 memset(seq,0x3f,sizeof seq); //初始化为正无穷 16 int ans=0; 17 for(int i=0;i<n;i++) 18 { 19 int pos=lower_bound(seq+1,seq+n,num[i])-seq; //pos-1是小于等于num[i]的最大元素 20 seq[pos]=num[i]; //更新序列 21 ans=max(ans,pos); //更新答案 22 } 23 printf("%d\n",ans); 24 } 25 return 0; 26 }
◇ 通过修改原串元素形成单调序列
【Codeforces 13C】Sequence +传送门+
将元素i从a改成b的花费为 |a-b| ,现给出一个序列,修改一些元素,使得序列形成单调序列(不下降/不上升),求最小花费。
一个简单结论:对于每一个数,修改其值为先前的数或不修改,一定可以使原串形成单调序列。
于是我们定义dp[i][j]为第i个数值为j时的最小花费。
从小到大枚举i先前的每一个数值(包括i,此时相当于不修改)j,表示第i个数将修改为j,则花费为|j-num[i]|。由于我们是从小到大枚举的值,我们可以储存Min,表示dp[i-1][1~j]中最小的值,以此更新dp[i][j],也就是说:
但是我们发现j太大了(1e9),存不下,于是可以用离散化来优化一下……最坏情况下,原序列每一个元素都不相等,第二维为5000,仍然存不下。回过头看我们的转移式,我们其实只需要用到dp[i-1],所以我们只需要滚动数组就可以了。
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 /*Lucky_Glass*/ 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #include<vector> 6 #include<cmath> 7 using namespace std; 8 const int MAXN=5000; 9 typedef long long ll; 10 int n; 11 int hgt[MAXN+5],fhgt[MAXN+5]; 12 vector<int> uni; 13 ll dp[2][MAXN+5]; 14 int main() 15 { 16 scanf("%d",&n); 17 for(int i=1;i<=n;i++) 18 scanf("%d",&hgt[i]),uni.push_back(hgt[i]); 19 sort(uni.begin(),uni.end()); 20 uni.erase(unique(uni.begin(),uni.end()),uni.end()); //离散化 21 ll ans=ll(1e16); 22 for(int i=1;i<=n;i++) 23 { 24 ll Min=ll(1e16); 25 for(int j=0;j<uni.size();j++) 26 { 27 Min=min((ll)Min,dp[(i-1)%2][j]); //取1~j的最小值 28 dp[i%2][j]=Min+fabs(1.0*hgt[i]-uni[j]); //%2 滚动数组 29 if(i==n) ans=min(ans,dp[i%2][j]); //更新答案 30 } 31 } 32 printf("%lld\n",ans); 33 return 0; 34 }
(Tab:以上每一篇代码都是模板,可以直接套上用的~)
The End
Thanks for reading!
- Lucky_Glass