洛谷 P1091 合唱队形

题目链接：https://www.luogu.org/problemnew/solution/P1091

思路：先求一遍最长递增子序列(dp1[i]表示以i位置结尾的LIS），再求一遍最长递减子序列(dp2[i]表示以i位置开始的LDS（最长递减子序列）,PS：相当于从后往前的最长递增子序列。答案是就是len1+len2-1。

先是我自己的超笨代码（复杂度o(n^2logn）)
写出下面这个代码是我对LIS的理解不够，我这个是对每个位置i,求一遍[0,i]的LIS，再求一遍[i,n-1]的LDS,然后再求len1+len2-1

#include <iostream>
#include <cstring>

using namespace std;

const int maxn=200;

int a[maxn],dp1[maxn],dp2[maxn],b[maxn];

int mfind(int left,int right,int x,int ta) {
    int mid;
    if((ta&&b[right]<x)||(!ta&&b[right]>x))
       return right+1;
    if(ta)
    {
       while(left<right) {  //求最长递增子序列 
         mid=(left+right)/2;
        if(b[mid]<x)left=mid+1;
        else if(b[mid]>=x)right=mid;
       }	
    }
    else  {
       while(left<right) {  //求最长递减子序列 
          mid=(left+right)/2;
       	 if(b[mid]>x)left=mid+1;
       	 else if(b[mid]<=x)right=mid;
       	 //cout <<left <<" "<<right<<endl;
       }	
    }
    return left;	
}

int main() {
    int n;
    cin >> n;
    for(int i=0;i<n;i++) cin >> a[i];
    memset(dp1,0,sizeof(dp1));
    memset(dp2,0,sizeof(dp2));
    int ans=0;
    for(int i=1;i<n-1;i++) {
          int cnt=1;
          memset(b,0,sizeof(b));
          b[0]=a[0];
          int k=0;
        for(int j=0;j<=i;j++) 
         {
         	 int p=mfind(0,cnt-1,a[j],1);
         	 //cout <<"?" <<endl;
         	 if(p<cnt) b[p]=min(b[p],a[j]);
         	 else  {
         	 	k=j;
         	 	b[cnt++]=a[j];
              }
         }
         //cout <<k <<" "<<i<<endl;
         if(cnt<=1)continue;
         memset(b,0,sizeof(b));
         b[0]=a[k];
         int cnt1=1;
         for(int j=k+1;j<n;j++) {
         	int p=mfind(0,cnt1-1,a[j],0);
         	//cout <<b[0]<<" "<<a[j]<<" "<<p <<endl;
         	if(p<cnt1) b[p]=max(b[p],a[j]);
         	 else  {
         	 	b[cnt1++]=a[j];
            }
         }
         if(cnt1<=1)continue;
        //cout <<i<<" "<<"k"<<k<<cnt <<" "<<cnt1<<endl;
            cnt=cnt+cnt1-1;
            ans=max(ans,cnt);
    }
    if(ans==0)cout <<0<<endl;
    else cout << n-ans <<endl;
    return 0;
}

改正以后的n^2复杂度代码：(因为上面已经写了nlogn的LIS(虽然并没有正确使用)了，所以就不再写了)

#include <iostream>
#include <cstring>

using namespace std;

const int maxn=200;

int a[maxn],dp1[maxn],dp2[maxn],b[maxn];



int main() {
	int n;
	cin >> n;
	for(int i=0;i<n;i++) cin >> a[i];
    for(int i=0;i<n;i++)dp1[i]=1;
    for(int i=0;i<n;i++)dp2[i]=1;
	for(int i=0;i<n-1;i++)
      for(int j=i+1;j<n;j++)
        if(a[j]>a[i])dp1[j]=max(dp1[j],dp1[i]+1);
    for(int i=n-1;i>0;i--)
      for(int j=i-1;j>=0;j--)
      	  if(a[j]>a[i])dp2[j]=max(dp2[j],dp2[i]+1);
    int ans=0;
    for(int i=0;i<n;i++) 
       ans=max(dp1[i]+dp2[i]-1,ans);
	cout <<n-ans<<endl;
	return 0;
}