uvalive 6932
Posted 天翎月
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了uvalive 6932相关的知识,希望对你有一定的参考价值。
三个串必须要一起dp
之前刚学了dfs的记忆化搜索的dp方式 觉得很舒服 现学现卖然后两个小时都没有做出来
优化1:之前在dfs中 对每一个pos都会枚举所有可能的组合 结合当前状态来产生新的状态 来决定接下来是直接算还是继续dfs
枚举的结果 只与当前状态和for的区间有关 所以提前预处理出每种情况的结果
优化完跑了个长2W的全?序列 在跑到pos=1W3左右的时候崩溃 手动加栈依然在这里崩溃 猜想是栈崩了。。
优化2:进行完优化1后代码变得好看多了 发现在每一步的结果只与下一步有关 还很严格 于是可以搞一个dp[zt][pos]出来 根据上一个优化算的状态数组直接从pos+1推出来
感想:dfs写dp是挺舒服 但是得是特定的 像数位dp就会很舒服 普通dp就不一定 大概是我接受了for写dp的设定
dfs写dp得注意深度。。之前写的dp才300层。。1e6的再敢写就剁手
for循环计数有时候可以用很舒服的预处理给做了
#include<stdio.h> #include<math.h> #include<string.h> #include<vector> #include<queue> #include<map> #include<string> #include<iostream> #include<algorithm> #include<stack> #include<set> using namespace std; #define L long long #define pb push_back #define lala printf("--------\n"); #define ph push #define rep(i, a, b) for (L i=a;i<=b;++i) #define dow(i, b, a) for (L i=b;i>=a;--i) #define fmt(i,n) if(i==n)printf("\n");else printf(" ") ; #define rnode(i,u) for(L i = head[u] ; i != -1 ; i = b[i].nex) #define fi first #define se second template<class T> inline void flc(T &A, L x){memset(A, x, sizeof(A));} L read(){L x=0,f=1;char ch=getchar();while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}return x*f;} const L mod = 1e9+9 ; char a[1000050] ; char b[1000050] ; char c[1000050] ; L la,lb,lc; L dp[3][1000050] ; L G[28][28][28][5][5] ; void init() { flc(G,0) ; rep(aa,0,27) rep(bb,0,27) rep(cc,0,27) { L La=1,Ra=26; L Lb=1,Rb=26; L Lc=1,Rc=26; if(aa!=27)La=Ra=aa; if(bb!=27)Lb=Rb=bb; if(cc!=27)Lc=Rc=cc; L (&g)[5][5] = G[aa][bb][cc] ; rep(a,La,Ra) rep(b,Lb,Rb) rep(c,Lc,Rc){ g[1][2] += a<b&&b<c ; g[1][1] += a==b&&b==c ; g[1][3] += a<b&&b==c ; g[1][4] += a==b&&b<c ; g[3][2] += b<c ; g[3][3] += b==c ; g[4][2] += a<b ; g[4][4] += a==b ; g[2][2] ++ ; } } } L jfhsb[1000050][5] ; L zhrsb() { flc(jfhsb,0) ; dow(pos,la,1) { rep(zt,1,4){ L aa,bb,cc; if(a[pos]==‘?‘)aa=27;else aa=a[pos]-‘a‘+1; if(b[pos]==‘?‘)bb=27;else bb=b[pos]-‘a‘+1; if(c[pos]==‘?‘)cc=27;else cc=c[pos]-‘a‘+1; L ztt2 = (((((dp[0][pos] * dp[1][pos])) % mod ) * dp[2][pos])%mod) ; L ztt1 = jfhsb[pos+1][1]%mod ; L ztt3 = jfhsb[pos+1][3]%mod ; L ztt4 = jfhsb[pos+1][4]%mod ; L zt1 = (ztt1 * G[aa][bb][cc][zt][1])%mod ; L zt2 = (ztt2 * G[aa][bb][cc][zt][2])%mod ; L zt3 = (ztt3 * G[aa][bb][cc][zt][3])%mod ; L zt4 = (ztt4 * G[aa][bb][cc][zt][4])%mod ; jfhsb[pos][zt] = (((((zt1+zt2)%mod)+zt3)%mod)+zt4) % mod; } } return jfhsb[1][1] ; } int main () { init() ; L t ; scanf("%lld" , &t) ; while(t--){ flc(dp,0) ; scanf("%s",a+1); scanf("%s",b+1); scanf("%s",c+1); la = strlen(a+1) ; lb = strlen(b+1) ; lc = strlen(c+1) ; L maxx = max(la,max(lb,lc)) ; rep(i,la+1,maxx) a[i] = ‘a‘ - 1 ; rep(i,lb+1,maxx) b[i] = ‘a‘ - 1 ; rep(i,lc+1,maxx) c[i] = ‘a‘ - 1 ; a[maxx+1] = ‘\0‘ ; b[maxx+1] = ‘\0‘ ; c[maxx+1] = ‘\0‘ ; la = lb = lc = maxx ; L wh = 1 ; for(L i=la;i>=1;i--){ if(a[i]==‘?‘){ wh %= mod ; dp[0][i] = wh ; wh *= 26 ; wh %= mod ; } else { dp[0][i] = wh ; } } wh = 1 ; for(L i=lb;i>=1;i--){ if(b[i]==‘?‘){ wh %= mod ; dp[1][i] = wh ; wh *= 26 ; wh %= mod ; } else { dp[1][i] = wh ; } } wh = 1 ; for(L i=lc;i>=1;i--){ if(c[i]==‘?‘){ wh %= mod ; dp[2][i] = wh ; wh *= 26 ; wh %= mod ; } else { dp[2][i] = wh ; } } L ans = zhrsb() ; printf("%lld\n" , ans) ; } }
以上是关于uvalive 6932的主要内容,如果未能解决你的问题,请参考以下文章
UVALive 4726 Average ——(斜率优化DP)
UVALive 4670 Dominating Patterns