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

[Usaco2008 Feb]Eating Together麻烦的聚餐

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改卡片编号的时候，都不会挪位置。

Input

第1行: 1个整数：N 第2..N+1行: 第i+1行是1个整数，为第i头奶牛的用餐批次D_i

Output

第1行: 输出1个整数，为FJ最少要改几头奶牛卡片上的编号，才能让编号变成他设想中的样子

Sample Input

5

1

3

2

1

1

输入说明:

队列中共有5头奶牛，第1头以及最后2头奶牛被设定为第一批用餐，第2头奶牛的预设是第三批用餐，第3头则为第二批用餐。

Sample Output

1

输出说明:

如果FJ想把当前队列改成一个不下降序列，他至少要改2头奶牛的编号，一种可行的方案是：把队伍中2头编号不是1的奶牛的编号都改成1。不过，如果FJ选择把第1头奶牛的编号改成3就能把奶牛们的队伍改造成一个合法的不上升序列了。

BZOJ 擦边AC COGS T了

var
 y,x:array[0..30000]of longint;
 dp:array[0..30000,1..3,1..2]of longint;
 i,j,k:longint;
 n,ans:longint;
begin
 readln(n);
 for i:=1 to n do
  readln(x[i]);
 for i:=1 to n do
  begin
   if x[i]=1
   then dp[i,1,1]:=dp[i-1,1,1]
   else dp[i,1,1]:=dp[i-1,1,1]+1;
   if x[i]=3
   then dp[i,3,1]:=dp[i-1,3,1]
   else dp[i,3,1]:=dp[i-1,3,1]+1;
  end;
 for i:=n downto 1 do
  begin
   if x[i]=1
   then dp[i,1,2]:=dp[i+1,1,2]
   else dp[i,1,2]:=dp[i+1,1,2]+1;
   if x[i]=3
   then dp[i,3,2]:=dp[i+1,3,2]
   else dp[i,3,2]:=dp[i+1,3,2]+1;
  end;
 for i:=1 to n do
  if x[i]<>2
  then y[i]:=y[i-1]+1
  else y[i]:=y[i-1];
 ans:=maxlongint;
 for i:=0 to n do
  for j:=i+1 to n+1 do
  begin
   if dp[i,1,1]+dp[j,3,2]+y[j-1]-y[i]<ans
   then ans:=dp[i,1,1]+dp[j,3,2]+y[j-1]-y[i];
   if dp[i,3,1]+dp[j,1,2]+y[j-1]-y[i]<ans
   then ans:=dp[i,3,1]+dp[j,1,2]+y[j-1]-y[i];
  end;
 writeln(ans);
end.

咳咳，这个是O（N）的

• dp[i,j]:以j结尾1到i的最小更改
• dp[i,j]:=min(dp[i,j],dp[i-1,k]) (j=x[i]) (1<=k<=j)
• dp[i,j]:=min(dp[i,j],dp[i-1,k]+1) (j<>x[i]) (1<=k<=j)

var
 x:array[0..300001]of longint;
 dp:array[0..300011,1..3]of longint;
 i,j,k:longint;
 n,ans:longint;
function min(a,b:longint):longint;
begin
 if a>b
 then exit(b)
 else exit(a);
end;

procedure swap(a,b:longint);
var c,d:longint;
begin
 c:=x[a]; d:=x[b];
 x[a]:=d; x[b]:=c;
end;

begin
 readln(n);
 for i:=1 to n do
  readln(x[i]);
 for i:=1 to n do
  for j:=1 to 3 do
   dp[i,j]:=30000000;
 dp[0,1]:=0; dp[0,2]:=0; dp[0,3]:=0;
 for i:=1 to n do
  for j:=1 to 3 do
   for k:=1 to j do
   begin
    if x[i]=j
    then dp[i,j]:=min(dp[i,j],dp[i-1,k])
    else dp[i,j]:=min(dp[i,j],dp[i-1,k]+1);
   end;
  ans:=min(min(dp[n,1],dp[n,2]),dp[n,3]);
 for i:=1 to (n div 2) do
  swap(i,n-i+1);
 for i:=1 to n do
  for j:=1 to 3 do
   dp[i,j]:=30000000;
 dp[0,1]:=0; dp[0,2]:=0; dp[0,3]:=0;
 for i:=1 to n do
  for j:=1 to 3 do
   for k:=1 to j do
    begin
     if x[i]=j
     then dp[i,j]:=min(dp[i,j],dp[i-1,k])
     else dp[i,j]:=min(dp[i,j],dp[i-1,k]+1);
   end;
  ans:=min(min(min(dp[n,1],dp[n,2]),dp[n,3]),ans);
 writeln(ans);
end.