Tyvj P1180(情况少分支多的过程-Dp）

这题看着就是不可解的。

但是仔细观察会发现它的情况较少，分支巨多（2^N）？？

于是我们用滚动数组+Dp显然可以用F[i][j][k][l][m]表示子结构

j,k,l,m为2个矿工最近2次的伙食（0表示没有）

任何坑爹的题目先想想能不能Dp……

记忆化搜索不能滚动还会爆栈，于是Dp的存在性证毕。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<cctype>
#include<iostream>
#include<functional>
#include<algorithm>
using namespace std;
#define MAXN (100000+10)
#define NDEBUG
int n,a[MAXN];
int f[2][4][4][4][4]; //1->M 2->F 3->B
int cost[4][4][4]={0};
int _cost(int i,int j,int k)
{
  if (k==0) return 0;
  if (i==0) i=k;
  if (j==0) j=k;  
  if (i!=j&&i!=k&&j!=k) return 3;
  if (i!=j||i!=k||j!=k) return 2;
  return 1;   
}
int main()
{
  #ifndef NDEBUG
  freopen("Vijos1386.in","r",stdin);
  #endif
  scanf("%d",&n);getchar();
  for (int i=0;i<n;i++) 
  {
    switch (getchar())
    {
    case 'M':a[i]=1;break;
    case 'F':a[i]=2;break;
    case 'B':a[i]=3;break;
    }
  } 
  for (int i=0;i<=3;i++)
    for (int j=0;j<=3;j++)
      for (int k=0;k<=3;k++)
      {
        cost[i][j][k]=_cost(i,j,k);
        #ifndef NDEBUG
        cout<<i<<' '<<j<<' '<<k<<':'<<cost[i][j][k]<<endl;
        #endif
      }
  
  
  memset(f,128,sizeof(f));
  f[0][0][0][0][0]=0;
  for (int ii=0;ii<n;ii++)
  {
    int i=ii%2;
    for (int j=0;j<=3;j++)
      for (int k=0;k<=3;k++)
        for (int l=0;l<=3;l++)
          for (int m=0;m<=3;m++)
          {
            if (f[i][j][k][l][m]>=0)
            {
              f[i^1][k][a[ii]][l][m]=max(f[i^1][k][a[ii]][l][m],f[i][j][k][l][m]+cost[j][k][a[ii]]);
              f[i^1][j][k][m][a[ii]]=max(f[i^1][j][k][m][a[ii]],f[i][j][k][l][m]+cost[l][m][a[ii]]);              
            }           
          }     
  }
  int ii=n%2,ans=0;
  for (int j=0;j<=3;j++)
    for (int k=0;k<=3;k++)
      for (int l=0;l<=3;l++)
        for (int m=0;m<=3;m++)
          ans=max(ans,f[ii][j][k][l][m]);
  cout<<ans<<endl;  
  
  
  return 0;
}