LIS-nefuoj-1209&nefuoj-1424-nlongn演算法&LIS的一些小變形
阿新 • • 發佈:2018-12-24
求最長公共子序列的長度,用普通的方法會超時。
用二分的方法找他插入的位置,如果在最後,那麼就插入到最後,否則就在第一個大於它的地方插入(嚴格遞增子序列這樣求,如果在第一個大於它的地方就是不下降子序列了,對應兩個方法)
挑戰程式設計2上的程式碼實現更好。哈哈
順便講一下LIS的變形。
1 正常求LIS,這誰不會啊,直接闆闆美滋滋,具體思路是,用暴力的方法求LIS,之所以不可以就是因為諸如 1 7 2 3 4 之類的,選7直接矇蔽,構成子序列 1 7 顯然不是最大的,最好的方法如果7大於2,這時在子序列裡二分第一個大於他的(很猛先這是不嚴格的上升),就替換。
2 求下降子序列,就把序列反轉一下就行,當然你求的dp陣列也要反轉哦。
3 求以i為開頭的 LIS,先反轉,在變成負數。
4 求以i為開頭的下降, 就先反轉,在反轉,在變成負數(反轉個毛線。。)
具體看圖。圖1時正常LIS。圖3時以i為底的LIS
給個例題
https://blog.csdn.net/qq_35781950/article/details/71273679
#include <iostream>
#include <cstdio>
using namespace std;
/* 開心,看了兩天終於看懂了。
之前 不怎麼確定的地方有下面這些方面,
這道題其實特別簡單,求一下嚴格的最長公共子序列,
然後在求一下嚴格的下降子序列(咋求下降子序列啊哈哈)
*/
#define INF 0x3f3f3f3f
int main()
{ int m;
int a[50006];
int dp[50006];
while (cin>>m)
{for(int i=0;i<m;i++)
cin>>a[i];
fill(dp,dp+m,INF);//全部初始化為無窮大,不用比較就找到了端點。
for(int i=0;i<m;i++)
*lower_bound(dp,dp+m,a[i])=a[i];
printf("%d\n",lower_bound(dp,dp+m,INF)-dp);
}
return 0;
}
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
using namespace std;
int main()
{ int t;
int m;
int a[50006];
int dp[50006];
memset(dp,0,sizeof(dp));
while(cin>>m)
{
for(int i=0;i<m;i++)
{ cin>>a[i];
}
vector<int>q;
for(int i=0;i<m;i++)
{ int u=upper_bound(q.begin(),q.end(),a[i])-q.begin();
if(u==q.size())
{ q.push_back(a[i]);
dp[i]=q.size();
}
else
{ q[u]=a[i];
dp[i]=u+1;
}
}
int big=-1;
for(int i=0;i<m;i++)
big=max(big,dp[i]);
printf("%d\n",big);
//for(int i=0;i<q.size();i++)
//printf("%d ",q[i]);
}
return 0;
}
nefu1424
這個,,我錯了好幾次,粗心。。
你們看程式碼就懂了,我渾身難受。
#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxn=1e5+5;
int dp[maxn];
int dp2[maxn];
int m;
void solve(vector<int>q){
vector<int>a;
for(int i=0;i<q.size();i++){
int u=lower_bound(a.begin(),a.end(),q[i])-a.begin();
if(u==a.size()){
a.push_back(q[i]);
dp[i]=max(dp[i],(int)a.size());
}
else{
a[u]=q[i];
dp[i]=max(u+1,dp[i]);//維護max無意義,直接等就好
}
}
}
void solve2(vector<int>q){
vector<int>a;
for(int i=0;i<q.size();i++){
int u=lower_bound(a.begin(),a.end(),q[i])-a.begin();
if(u==a.size()){
a.push_back(q[i]);
dp2[i]=max((int)a.size(),dp2[i]);
}
else{
a[u]=q[i];
dp2[i]=max(u+1,dp2[i]);
}
}
reverse(dp2,dp2+m);
}
int main()
{ int a;
vector<int>q;
while(cin>>m){
//if(!m)break;
q.clear();
for(int i=0;i<m;i++)
{ dp[i]=dp2[i]=0;
}
for(int i=0;i<m;i++){
cin>>a;
q.push_back(a);
}
solve(q);
reverse(q.begin(),q.end());
solve2(q);
int ans=-1;
for(int i=0;i<m;i++){
//cout<<dp[i]<<" "<<dp2[i]<<endl;
//if(dp[i]==dp2[i])
ans=max(ans,2*min(dp[i],dp2[i])-1);
}
cout<<ans<<endl;
}
return 0;
}