1. 程式人生 > >LIS-nefuoj-1209&nefuoj-1424-nlongn演算法&LIS的一些小變形

LIS-nefuoj-1209&nefuoj-1424-nlongn演算法&LIS的一些小變形

求最長公共子序列的長度,用普通的方法會超時。
用二分的方法找他插入的位置,如果在最後,那麼就插入到最後,否則就在第一個大於它的地方插入(嚴格遞增子序列這樣求,如果在第一個大於它的地方就是不下降子序列了,對應兩個方法)
挑戰程式設計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;
}