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) ;
    }
}

 

转载于:https://www.cnblogs.com/rayrayrainrain/p/7635471.html