bzoj1609[Usaco2008 Feb]Eating Together麻烦的聚餐*

bzoj1609[Usaco2008 Feb]Eating Together麻烦的聚餐

题意：

一个序列只由1﹑2﹑3三种数组成。求最少要改变多少个数使它变成不下降序列或不上升序列。序列大小≤30000

题解：

DP。设f[i][j]表示正在考虑第i个数，上一个数是j。求不下降序列最少改变个数方程：

f[i][j]=min(f[i+1][k]+1,k∈[j,3]),a[i]<j
          min(f[i+1][a[i]],f[i+1][k]+1,k∈[j,3]且k!=a[i])a[i]>=j

求不上升同理。

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 50000
#define inc(i,j,k) for(int i=j;i<=k;i++)
#define dec(i,j,k) for(int i=j;i>=k;i--)
#define INF 0x3fffffff
using namespace std;

int n,a[maxn],x[5],y[5],ans;
int main(){
    scanf("%d",&n); inc(i,1,n)scanf("%d",&a[i]); memset(x,0,sizeof(x));
    dec(i,n,1){
        memset(y,0,sizeof(y));
        inc(j,1,3)if(a[i]<j){
            y[j]=INF; inc(k,j,3)y[j]=min(y[j],x[k]+1);
        }else{
            y[j]=x[a[i]]; inc(k,j,3)if(k!=a[i])y[j]=min(y[j],x[k]+1);
        }
        swap(x,y);
    }
    ans=INF; inc(i,1,3)ans=min(ans,x[i]); memset(x,0,sizeof(x));
    dec(i,n,1){
        memset(y,0,sizeof(y));
        inc(j,1,3)if(a[i]>j){
            y[j]=INF; dec(k,j,1)y[j]=min(y[j],x[k]+1);
        }else{
            y[j]=x[a[i]]; dec(k,j,1)if(k!=a[i])y[j]=min(y[j],x[k]+1);
        }
        swap(x,y);
    }
    inc(i,1,3)ans=min(ans,x[i]); printf("%d",ans); return 0;
}

20160729