[BZOJ4043/CERC2014]Vocabulary

Description
给你三个字符串，这些字符串有些单词模糊不可认了，用"?"来代表。
现在你可以用任意英文小写字母来代表它们。要求是使得给定的三个字符串中
所有的"?"被你认定的字母代替后，各不相同且按字典序出现。问有多少种方式。

Input
先给出一个数字N，代表数据组数。
接下来3*N行，每行给出一个字符串。长度<=1000 000

Output
输出结果 Mod 10^9+9

Sample Input
3
?heoret?cal
c?mputer
?cience
jagiellonia
?niversity
kra?ow
?
b
c

Sample Output
42562
52
1

---

我们令三个字符串为A,B,C，设\(f[i][0/1/2/3]\)表示当前扫到第i个字符串，0表示A=B=C，1表示A=B<C，2表示A<B=C，3表示A<B<C的方案数，但是这样转移需要枚举当前填的字符和上一个填的字符，复杂度太高

所以我们可以预处理出\(g[x][y][i][j][k]\)，表示当前字符为i,j,k，状态由x转移到y的方案数，然后就可以大力转移f了

/*problem from Wolfycz*/
#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline char gc(){
    static char buf[1000000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
inline int frd(){
    int x=0,f=1; char ch=gc();
    for (;ch<'0'||ch>'9';ch=gc())   if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=gc()) x=(x<<3)+(x<<1)+ch-'0';
    return x*f;
}
inline int read(){
    int x=0,f=1; char ch=getchar();
    for (;ch<'0'||ch>'9';ch=getchar())  if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=getchar())    x=(x<<3)+(x<<1)+ch-'0';
    return x*f;
}
inline void print(int x){
    if (x<0)    putchar('-');
    if (x>9)    print(x/10);
    putchar(x%10+'0');
}
const int N=1e6,p=1e9+9;
int g[5][5][30][30][30];
void prepare(){
    for (int i=0;i<=27;i++){
        for (int j=0;j<=27;j++){
            for (int k=0;k<=27;k++){
                int li=(i==27)?1:i,ri=(i==27)?26:i;
                for (int x=li;x<=ri;x++){
                    int lj=(j==27)?1:j,rj=(j==27)?26:j;
                    for (int y=lj;y<=rj;y++){
                        int lk=(k==27)?1:k,rk=(k==27)?26:k;
                        for (int z=lk;z<=rk;z++){
                            if (x==y&&y==z) g[0][0][i][j][k]++;
                            if (x==y&&y< z) g[0][1][i][j][k]++;
                            if (x< y&&y==z) g[0][2][i][j][k]++;
                            if (x< y&&y< z) g[0][3][i][j][k]++;
                            if (x==y)   g[1][1][i][j][k]++;
                            if (y==z)   g[2][2][i][j][k]++;
                            if (x< y)   g[1][3][i][j][k]++;
                            if (y< z)   g[2][3][i][j][k]++;
                            g[3][3][i][j][k]++;
                        }
                    }
                }
            }
        }
    }
}
char s[N+10];
int v[5][N+10];
int f[5][N+10];
void init(int &len,int *a){
    scanf("%s",s+1);
    len=strlen(s+1);
    for (int i=1;i<=len;i++)    a[i]=(s[i]=='?'?27:s[i]-'a'+1);
}
int main(){
    prepare();
    for (int Data=read();Data;Data--){
        int n1,n2,n3,len;
        init(n1,v[1]),init(n2,v[2]),init(n3,v[3]);
        len=max(max(n1,n2),n3);
        for (int i=n1+1;i<=len;i++) v[1][i]=0;
        for (int i=n2+1;i<=len;i++) v[2][i]=0;
        for (int i=n3+1;i<=len;i++) v[3][i]=0;
        f[0][0]=1;
        for (int i=1;i<=len;i++){
            for (int j=0;j<4;j++)   f[j][i]=0;
            for (int j=0;j<4;j++)
                for (int k=0;k<4;k++)
                    f[j][i]=(f[j][i]+1ll*f[k][i-1]*g[k][j][v[1][i]][v[2][i]][v[3][i]]%p)%p;
        }
        printf("%d\n",f[3][len]);
    }
    return 0;
}

转载于:https://www.cnblogs.com/Wolfycz/p/10030234.html