首页 > 代码库 > hdu 4758 Walk Through Squares(AC自动机+状态压缩DP)
hdu 4758 Walk Through Squares(AC自动机+状态压缩DP)
题目链接:hdu 4758 Walk Through Squares
题意:
给你一个n*m的网格,现在你要从(1,1)走到(n,m),每次只能向右走或者向下走,走完后会形成一个包含R,D的序列,这个序列必须要包含题目给出的两个字符串,问有多少种方案。
题解:
由于要从(1,1)走到(n,m),所以这个形成的字符串包含R和D的数量是一定的。
现在考虑dp[u][i][j][k],表示包含两种串的组合状态,走了i个D,走了j个R,走到了k这个AC自动机的节点的方案数,然后dp一下就行了。复杂度为O(3*n*m*len(s1)*len(s2))
1 #include<bits/stdc++.h> 2 #define mst(a,b) memset(a,b,sizeof(a)) 3 #define F(i,a,b) for(int i=(a);i<=(b);++i) 4 using namespace std; 5 6 const int AC_N=300+7,tyn=2; 7 int t,n,m,dp[4][111][111][212],P=1e9+7,ans; 8 char a[101],b[101]; 9 inline void up(int &a,int b){a=(a+b)%P;} 10 11 struct AC_automation{ 12 int tr[AC_N][tyn],cnt[AC_N],Q[AC_N],fail[AC_N],tot; 13 inline int getid(char x){ 14 if(x==‘D‘)return 0; 15 return 1; 16 } 17 void nw(){cnt[++tot]=0,fail[tot]=0;mst(tr[tot],0);} 18 void init(){tot=-1,fail[0]=-1,nw();} 19 void insert(char *s,int idx,int x=0){ 20 for(int len=strlen(s),i=0,w;i<len;x=tr[x][w],i++) 21 if(!tr[x][w=getid(s[i])])nw(),tr[x][w]=tot; 22 cnt[x]=1<<idx; 23 } 24 void build(int head=1,int tail=0){ 25 for(int i=0;i<tyn;i++)if(tr[0][i])Q[++tail]=tr[0][i]; 26 while(head<=tail)for(int x=Q[head++],i=0;i<tyn;i++) 27 if(tr[x][i]) 28 { 29 fail[tr[x][i]]=tr[fail[x]][i],Q[++tail]=tr[x][i]; 30 cnt[tr[x][i]]|=cnt[tr[fail[x]][i]]; 31 } 32 else tr[x][i]=tr[fail[x]][i]; 33 } 34 void solve() 35 { 36 mst(dp,0),dp[0][0][0][0]=1; 37 F(u,0,3)F(i,0,n)F(j,0,m)F(k,0,tot) 38 { 39 if(!dp[u][i][j][k])continue; 40 up(dp[u|cnt[tr[k][0]]][i+1][j][tr[k][0]],dp[u][i][j][k]); 41 up(dp[u|cnt[tr[k][1]]][i][j+1][tr[k][1]],dp[u][i][j][k]); 42 } 43 ans=0; 44 F(k,0,tot)up(ans,dp[3][n][m][k]); 45 printf("%d\n",ans); 46 } 47 }AC; 48 49 50 int main() 51 { 52 scanf("%d",&t); 53 while(t--) 54 { 55 scanf("%d%d%s%s",&m,&n,a,b); 56 AC.init(),AC.insert(a,0),AC.insert(b,1); 57 AC.build(),AC.solve(); 58 } 59 return 0; 60 }
hdu 4758 Walk Through Squares(AC自动机+状态压缩DP)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。