hdu 4057

http://acm.hdu.edu.cn/showproblem.php?pid=4057


我不喜欢把它叫做ac自动机+dp ,个人觉得叫做ac自动机上的dp更好,原因就在于我们用ac自动机的一个个节点的状态来构建状态和状态之间的转移,至于非要为什么用ac自动机来构建状态,个人认为它包含了题目给出的说有字串的转移之间的信息(保证正确),而且状态数很少(保证了时间复查度小)

dp[node ][ state ] 表示当前长度到达 node节点,包含state字串(用二进制状态压缩),因为题目说了每个字串只能计算一次,否则的话 直接dp【node】 每次取最大值就可以了

#include <cstdlib>
#include <iostream>
#include <algorithm> 

using namespace std;
const int maxn=1010;
int chd[maxn][4],fail[maxn],word[maxn],sz,sw[128];
bool dp1[1010][1050],dp2[1010][1050];
int n,len;
char s[110];
int val[20];

void Ins(char *a,int val)
{
     int p=0;
    for(;*a;a++)     
    {
        int id=sw[*a];                 
        if(!chd[p][id])
        {
             memset(chd[sz],0,sizeof(chd[sz]));               
             word[sz]=0;
             chd[p][id]=sz++;
        }
        p=chd[p][id];
    }
    word[p]|=(1<<val);
}
int Que[maxn];
void ac()
{
    int *s=Que,*e=Que; 
    for(int i=0;i<4;i++)     
      if(chd[0][i])
      {
          *e++=chd[0][i];
          fail[chd[0][i]]=0;
      }
    while(s!=e)  
    {
         int p=*s++;
         for(int i=0;i<4;i++)             
           if(chd[p][i])
           {
                *e++=chd[p][i];
                fail[chd[p][i]]=chd[fail[p]][i];
                word[chd[p][i]]|=word[fail[chd[p][i]]];
           }
           else
             chd[p][i]=chd[fail[p]][i];
    }
} 
void solve()
{
    memset(dp1,0,sizeof(dp1));  
    memset(dp2,0,sizeof(dp2));  
    dp1[0][0]=1;
    int Max=(1<<n);
    for(int k=0;k<len;k++)
    {
         for(int i=0;i<sz;i++)        
         for(int j=0;j<Max;j++)
          if(dp1[i][j])
          {
               for(int r=0;r<4;r++)             
                dp2[chd[i][r]][j|word[chd[i][r]]]=1;
          }
        memcpy(dp1,dp2,sizeof(dp2));  
        memset(dp2,0,sizeof(dp2));  
    }
    long long ans=-1;
    for(int i=0;i<sz;i++)        
    for(int j=0;j<Max;j++)
      if(dp1[i][j]) 
      {
           long long tmp=0;         
           for(int r=0;r<n;r++)              
            if(j&(1<<r)) tmp+=(long long)val[r];
           if(ans<tmp) ans=tmp; 
      }
     if(ans>=0) cout<<ans<<endl; 
     else cout<<"No Rabbit after 2012!"<<endl;
}

int main(int argc, char *argv[])
{
    sw['A']=0;
    sw['T']=1;
    sw['C']=2;
    sw['G']=3;
    while(cin>>n>>len)
    {
          sz=1;
          memset(chd[0],0,sizeof(chd[0]));
          for(int i=0;i<n;i++)  
          {
             cin>>s>>val[i];
             Ins(s,i);
          }
         ac();
         solve();  
    }
   // system("PAUSE");
    return EXIT_SUCCESS;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值