hdu 3341 Lost's revenge

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3341

题目大意：给出一些病毒串，一个原串，问怎么调整原串的顺序使得跟最多的病毒串匹配。

题目思路：用dp[a][b][c][d][e]表示分别用了a,b,c,d个碱基，到达e结点的最值。开始加了一个总长度，后来想到既然有总长度，则有一维可以去，但是还是超时，哎，居然没有想到不用总长度直接用4种碱基的长度来表示状态，这样就又可以降低常数了，居然一下就刷到第一了，嘿嘿。。

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<string>
#include<queue>
#include<algorithm>
#include<vector>
#include<stack>
#include<list>
#include<iostream>
#include<map>
using namespace std;
#define inf 0x3f3f3f3f
#define M 110
inline int max(int a,int b)
{
return a>b?a:b;
}
int min(int a,int b)
{
return a<b?a:b;
}
int num[4];
int cnt;
int dp[2][42][42][42][550];
char s[100],str[100];
int mp[1000],q[55*10];
struct node
{
int cnt,fail;
int next[4];
void init()
{
cnt=fail=0;
memset(next,0,sizeof(next));
}
}tri[55*10];
void insert(char *s)
{
int i,p,x;
p=0;
for(i=0;s[i];i++)
{
x=mp[s[i]];
if(!tri[p].next[x])
{
tri[++cnt].init();
tri[p].next[x]=cnt;
}
p=tri[p].next[x];//不要等于cnt,只有开辟新节点才等于cnt;
}
tri[p].cnt++;
}
void bfs()
{
int i,p,head,tail,suf;
p=head=tail=0;
for(i=0;i<4;i++)
{
if(tri[0].next[i])
{
q[tail++]=tri[0].next[i];
tri[q[tail-1]].fail=0;
}
}
while(head<tail)
{
p=q[head++];suf=tri[p].fail;
tri[p].cnt+=tri[suf].cnt;
for(i=0;i<4;i++)
{
if(tri[p].next[i])
{
q[tail++]=tri[p].next[i];
tri[q[tail-1]].fail=tri[suf].next[i];
}
else
tri[p].next[i]=tri[suf].next[i];
}
}
}
void set(int a)
{
int j,b,c,d,e;
for(b=0;b<=num[1];b++)
for(c=0;c<=num[2];c++)
for(d=0;d<=num[3];d++)
for(e=0;e<=cnt;e++)
{
dp[a][b][c][d][e]=-1;
}
}
void solve(int n)
{
int tmp,i,j,a,b,c,d,e,pos,ans;
tmp=0;
set(0);
dp[0][0][0][0][0]=0;
for(a=0;a<=num[0];a++)
{
set(tmp^1);
for(b=0;b<=num[1];b++)
for(c=0;c<=num[2];c++)
for(d=0;d<=num[3];d++)
{
for(e=0;e<=cnt;e++)
{
if(dp[tmp][b][c][d][e]==-1) continue;
pos=tri[e].next[0];
dp[tmp^1][b][c][d][pos]=max(dp[tmp^1][b][c][d][pos],dp[tmp][b][c][d][e]+tri[pos].cnt);
pos=tri[e].next[1];
dp[tmp][b+1][c][d][pos]=max(dp[tmp][b+1][c][d][pos],dp[tmp][b][c][d][e]+tri[pos].cnt);
pos=tri[e].next[2];
dp[tmp][b][c+1][d][pos]=max(dp[tmp][b][c+1][d][pos],dp[tmp][b][c][d][e]+tri[pos].cnt);
pos=tri[e].next[3];
dp[tmp][b][c][d+1][pos]=max(dp[tmp][b][c][d+1][pos],dp[tmp][b][c][d][e]+tri[pos].cnt);
}
}
tmp^=1;
}
ans=0;
for(i=0;i<=cnt;i++)
{
ans=max(ans,dp[tmp^1][num[1]][num[2]][num[3]][i]);
}
printf("%d\n",ans);
}
int main()
{
int n,i,len,count=1;
mp['A']=0;mp['C']=1;mp['G']=2;mp['T']=3;
while(scanf("%d",&n),n)
{
cnt=0;
tri[0].init();
for(i=0;i<n;i++)
{
scanf("%s",str);
insert(str);
}
scanf("%s",s);
len=strlen(s);
memset(num,0,sizeof(num));
for(i=0;i<len;i++)
{
num[mp[s[i]]]++;
}
bfs();
printf("Case %d: ",count++);
solve(len);

}

}