首页 > 代码库 > 【HDU3341】AC自动机状态压缩DP,或者说hash枚举DP,-------出题人卡常数都是狗!!!!!
【HDU3341】AC自动机状态压缩DP,或者说hash枚举DP,-------出题人卡常数都是狗!!!!!
题意:给若干种个串,再给个主串,然后把主串打乱顺序,使得包含子串尽量多(一种可以有多个,两个之间可以部分重叠)。如第一组数据,ACGT,包含AC、CG、GT,三个,输出3。第二组数据A1A2A3,包含A1A2和A2A3两个“AA”,答案为2。
其实我并没有AC。我被卡常数TLE了。。。实在不想写这种没意义的东西了。
贴代码,待填坑。
#include <queue> #include <cstdio> #include <cstring> #include <iostream> #define max(a,b) (a>b?a:b) #define N 41 #define P 505 #define M 15000// 14641 #define T 4 #define inf 0x3f3f3f3f using namespace std; char ttt[50]; int f[P][M],g[P][M],crs[200]; int cn[5],mul[5],ed; struct ACAUTO { int next[P][T],num[P],fail[P],cnt,stack[P]; queue<int>q; void init() { while(!q.empty())q.pop(); memset(num,0,sizeof(num)); memset(fail,0,sizeof(fail)); } void insert() { static int i,x,alp; x = 0; scanf("%s",ttt); for(i=0;ttt[i];i++) { alp=crs[ttt[i]]; if(!next[x][alp])next[x][alp]=++cnt; x=next[x][alp]; } ++num[x]; } bool build() { int n,i,u,v,temp; int tcnt=0; scanf("%d",&n); if(!n)return 0; init(); for(i=1;i<=n;++i)insert(); q.push(0); while(!q.empty())/*維護fail指針*/ { // static int _ = 0; // cout << ++_ << endl; u=q.front(),q.pop(); stack[++tcnt]=u; for(i=0;i<T;++i)if(v=next[u][i]) { if(u) { temp=fail[u]; while(temp&&!next[temp][i])temp=fail[temp]; fail[v]=next[temp][i]; num[v]+=num[fail[v]]; } q.push(v); } } return 1; } void dp() { int i,j,u,v; int zt; for(u=0;u<=cnt;u++) { for(i=0;i<ed;i++) { for(j=0;j<T;j++)if(next[u][j]&&((j?i%mul[j-1]:i)/mul[j]+1)<cn[j]) { v=next[u][j]; zt=i+mul[j]; f[v][zt]=max(f[v][zt],g[u][i]+num[v]); // f[v][zt]=max(f[v][zt],f[u][i]+num[v]); } } } for(i=cnt+1;i;i--) { u=stack[i]; v=fail[u]; // for(j=0;j<ed;j++)f[v][j]=max(f[v][j],f[u][j]); for(j=0;j<ed;j++)g[u][j]=max(g[u][j],f[u][j]),g[v][j]=max(g[v][j],g[u][j]); } } }acauto; void work(int cas) { int i,j,k; scanf("%s",ttt); memset(cn,0,sizeof(cn)); memset(f,0xc0,sizeof(f)); memset(g,0xc0,sizeof(g)); f[0][0]=g[0][0]=0; for(i=0;ttt[i];i++)cn[crs[ttt[i]]]++; for(i=0;i<T;i++)cn[i]++; for(mul[T-1]=1,i=T-2;i>=0;i--)mul[i]=cn[i+1]*mul[i+1]; ed=cn[0]*cn[1]*cn[2]*cn[3]; int ans=0,llen=strlen(ttt); for(i=0;i<llen;i++) acauto.dp(); for(i=0;i<=acauto.cnt;i++)for(j=0;j<ed;j++)ans=max(ans,g[i][j]); printf("Case %d: %d\n",cas,ans); } int main() { // freopen("test.in","r",stdin); int cas=0; crs['A']=0,crs['C']=1,crs['G']=2,crs['T']=3; while(acauto.build())work(++cas); return 0; } /*2896 3065 2243 2825 3341 */
【HDU3341】AC自动机状态压缩DP,或者说hash枚举DP,-------出题人卡常数都是狗!!!!!
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。