BZOJ_1609_[Usaco2008_Feb]_Eating_Together_麻烦的聚餐_(LIS)

Posted 2020-07-12 晴歌。

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了BZOJ_1609_[Usaco2008_Feb]_Eating_Together_麻烦的聚餐_(LIS)相关的知识，希望对你有一定的参考价值。

描述

http://www.lydsy.com/JudgeOnline/problem.php?id=1609

给出一串由1,2,3组成的数,求最少需要改动多少个数,使其成为不降或不升序列.

分析

法1:改动一些数字后变为不升(不降)序列,那么除了需要改动的数字以外,其他的数字本身满足不升(不降),所以求最长不升(不降)子序列即可.O(nlogn)

法2:考虑搜索的思路,枚举当前位置的值,如果和原来的值相等,那么不许改动,否则改动数+1,然后搜索下一个位置,值要大于等于当前位置的改动数.这样会有重叠子问题,所以可以记忆化,倒过来就是dp.

dp[i][j]表示前i个,结尾为j的改动数,则有:

如果a[i]==j,则dp[i][j]=dp[i-1][k](k<=j)

如果a[i]!=j,则dp[i][j]=dp[i-1][k]+1(k<=j) O(n).

LIS:

O(n^2)的算法:

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 
 4 const int maxn=30000+5;
 5 int n,ans;
 6 int a[maxn],dp[maxn];
 7 
 8 int main(){
 9     scanf("%d",&n);
10     for(int i=1;i<=n;i++) scanf("%d",&a[i]);
11     for(int i=1;i<=n;i++){
12         dp[i]=1;
13         for(int j=1;j<i;j++)
14             if(a[j]<=a[i]) dp[i]=max(dp[j]+1,dp[i]);
15         ans=max(ans,dp[i]);
16     }
17     for(int i=1;i<=n;i++){
18         dp[i]=1;
19         for(int j=1;j<i;j++)
20             if(a[j]>=a[i]) dp[i]=max(dp[j]+1,dp[i]);
21         ans=max(ans,dp[i]);
22     }
23     printf("%d\\n",n-ans);
24     return 0;
25 }

View Code

O(nlogn)的算法:

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 
 4 const int maxn=30000+5,INF=0x7fffffff;
 5 int n,ans;
 6 int a[maxn],b[maxn],dp[maxn];
 7 
 8 int main(){
 9     scanf("%d",&n);
10     for(int i=1;i<=n;i++) scanf("%d",&a[i]), b[n+1-i]=a[i];
11     fill(dp+1,dp+n+1,INF);
12     for(int i=1;i<=n;i++)
13         *upper_bound(dp+1,dp+n+1,a[i])=a[i];
14     ans=lower_bound(dp+1,dp+n+1,INF)-(dp+1);
15     fill(dp+1,dp+n+1,INF);
16     for(int i=1;i<=n;i++)
17         *upper_bound(dp+1,dp+n+1,b[i])=b[i];
18     ans=max(ans,lower_bound(dp+1,dp+n+1,INF)-(dp+1));
19     printf("%d\\n",n-ans);
20     return 0;
21 }

View Code

快得不是一点点...

O(n)的算法:

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 
 4 const int maxn=30000+5,INF=0x3fffffff;
 5 int n,ans;
 6 int a[maxn];
 7 int dp[maxn][5];
 8 int main(){
 9     scanf("%d",&n);
10     for(int i=1;i<=n;i++) scanf("%d",&a[i]);
11     ans=INF;
12     for(int i=1;i<=n;i++)for(int j=1;j<=3;j++) dp[i][j]=INF;
13     for(int i=1;i<=n;i++)for(int j=1;j<=3;j++)for(int k=1;k<=j;k++)
14         if(a[i]==j) dp[i][j]=min(dp[i][j],dp[i-1][k]);
15         else dp[i][j]=min(dp[i][j],dp[i-1][k]+1);
16     for(int i=1;i<=3;i++) ans=min(ans,dp[n][i]);
17     for(int i=1;i<=n;i++)for(int j=1;j<=3;j++) dp[i][j]=INF;
18     for(int i=n;i>=1;i--)for(int j=1;j<=3;j++)for(int k=1;k<=j;k++)
19         if(a[i]==j) dp[i][j]=min(dp[i][j],dp[i+1][k]);
20         else dp[i][j]=min(dp[i][j],dp[i+1][k]+1);
21     for(int i=1;i<=3;i++) ans=min(ans,dp[1][i]);
22     printf("%d\\n",ans);
23 }

View Code

这里用到了STL里的二分查找.对于严格上升子序列,找到dp数组中a[i]的下界(满足dp[k]>=a[i]的最小的k(上界也可以,因为没有重复所以不会取等))即可,这样的话如果和以前的相同就覆盖了.对于这个问题,求的是不降子序列,所以相同的不能覆盖,所以我们找到dp数组中a[i]的上界(满足dp[k]>a[i]的最小的k),这样如果和以前相同的话就继续往后排,把原来上升序列中更大的结尾更新小.

对于求严格下降子序列,就需要找满足dp[k]<=a[i]的最小的k,不严格的就是dp[k]<a[i]的最小的k.

1609: [Usaco2008 Feb]Eating Together麻烦的聚餐

Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 1326 Solved: 801
[Submit][Status][Discuss]

Description

为了避免餐厅过分拥挤，FJ要求奶牛们分3批就餐。每天晚饭前，奶牛们都会在餐厅前排队入内，按FJ的设想所有第3批就餐的奶牛排在队尾，队伍的前端由设定为第1批就餐的奶牛占据，中间的位置就归第2批就餐的奶牛了。由于奶牛们不理解FJ的安排，晚饭前的排队成了一个大麻烦。第i头奶牛有一张标明她用餐批次D_i(1 <= D_i <= 3)的卡片。虽然所有N(1 <= N <= 30,000)头奶牛排成了很整齐的队伍但谁都看得出来，卡片上的号码是完全杂乱无章的。在若干次混乱的重新排队后，FJ找到了一种简单些的方法：奶牛们不动，他沿着队伍从头到尾走一遍把那些他认为排错队的奶牛卡片上的编号改掉，最终得到一个他想要的每个组中的奶牛都站在一起的队列，例如111222333或者333222111。哦，你也发现了，FJ不反对一条前后颠倒的队列，那样他可以让所有奶牛向后转，然后按正常顺序进入餐厅。你也晓得，FJ是个很懒的人。他想知道，如果他想达到目的，那么他最少得改多少头奶牛卡片上的编号。所有奶牛在FJ改卡片编号的时候，都不会挪位置。