【动态规划/递推】BZOJ1806[IOI2007]- Miners

IOI历史上的著名水题，我这种蒟蒻都能写的东西。

【思路】

用1、2、3分别代替三种食物，0表示当前矿井没有食物。f[i][a][b][c][d]当前第i个食物，矿1的食物顺序由上至下为a,b；矿2的食物顺序由上至下为c,d。 

判断产物数量的方法很巧妙，由下至上a,b,c。初始时默认投入一个食物至少生产一单位，如果a为有食物且与bc不同，则加一单位；如果b为有食物且与c不同，再加一个单位。最后加一个滚动数组就可以了。

【错因】

1.因为a,b,c,d大小范围是0..3，但是我把下标范围写3！一定要写4！我就是一不小心写错了，居然改了两个晚上，完全没有考虑到错因在这里....

2.还有初始化的时候只有f[0][0][0][0][0]是0，不要把所有第一位是0的都写成0了..

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int MAXN=100000+500;
int n;
int food[MAXN];
int f[2][4][4][4][4];//f[i][a][b][c][d]当前第i个食物，矿1的食物顺序由上至下为a,b；矿2的食物顺序由上至下为c,d。 

int compare(int x,int y,int z)
{
    int t=1;
    if (x!=0 && x!=y && x!=z) t++;
    if (y!=0 && y!=z) t++;
    return t;
}

int main()
{
    scanf("%d",&n);
    getchar();
    for (int i=0;i<n;i++)
    {
        char c;
        scanf("%c",&c);
        if (c=='M') food[i]=1;
        if (c=='F') food[i]=2;
        if (c=='B') food[i]=3;
    }
    
    memset(f,-1,sizeof(f));
    f[0][0][0][0][0]=0;33     int ans=-1;
    for (int i=0;i<n;i++)
        for (int a=0;a<=3;a++)
            for (int b=0;b<=3;b++)
                for (int c=0;c<=3;c++)
                    for (int d=0;d<=3;d++)
                    {
                        if (f[i%2][a][b][c][d]==-1) continue;
                        int add1=compare(b,a,food[i]),add2=compare(d,c,food[i]);
                        
                        int x=food[i];
                        f[(i+1)%2][x][a][c][d]=max(f[(i+1)%2][x][a][c][d],f[i%2][a][b][c][d]+add1);
                        f[(i+1)%2][a][b][x][c]=max(f[(i+1)%2][a][b][x][c],f[i%2][a][b][c][d]+add2);
                        if ((i==n-1) && (f[(i+1)%2][x][a][c][d]>ans)) ans=f[(i+1)%2][x][a][c][d];
                        if ((i==n-1) && (f[(i+1)%2][a][b][x][c]>ans)) ans=f[(i+1)%2][a][b][x][c];
                    }
    cout<<ans<<endl;
    return 0;
}