bzoj3802 cerc2014 Vocabulary

3802: Vocabulary

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 54 Solved: 11
[ Submit][ Status]

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

HINT

Source

题解

首先数据范围是1000000，且多组数据
那么显然必须是O(n)的算法
先把所有的串补成一样长，就是加'a'-1就好了
定义f[i][0...3]
f[i][0]表示做到第i位3个字符串前面i位相同的情况数
f[i][1] 一号和二号 字符串相同且字典序比三号小的情况数
f[i][2]二号和三号字符串相同且比一号字典序小的情况数
f[i][3]已经按照字典序排列的情况数
我们首先将所有的字符串长度变成一样长
转移脑补一下就好了...是不是发现脑补出来32种转移...


所有要先做一遍预处理d[0..3][0...3][0...27][0...27][0...27]
d[i][j][a][b][c]
上一个情况为i，这一位a，b，c三个字符 //0表示补上的，27则是'?'
转移成j号情况的方案数，直接暴力6个for循环即可,不会T的，因为枚举'?'的很少的说
然后就可以转移了
上代码，自己看看吧...

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn = 1001000;
long long f[maxn][4]; int n;
int a[maxn], b[maxn], c[maxn];
char ch1[maxn], ch2[maxn], ch3[maxn]; 
int get_ch(int *a, char *ch, int n){
	int len = strlen(ch);
	for (int i = 0; i < len; i ++)
		if (ch[i] == '?') a[i + 1] = 27;
		else a[i + 1] = ch[i] - 'a' + 1;
	for (int i = len + 1; i <= n; i ++) a[i] = 0;
}
void init(){
	n = 0; scanf("%s", ch1); 
	scanf("%s", ch2);
	scanf("%s", ch3);
	n = max(strlen(ch1), max(strlen(ch2), strlen(ch3)));
	get_ch(a, ch1, n);
	get_ch(b, ch2, n);
	get_ch(c, ch3, n);
}
int d[4][4][28][28][28];
void prapare(){
	for (int i1 = 0; i1 <= 27; i1 ++)
		for (int j1 = 0; j1 <= 27; j1 ++)
			for (int k1 = 0; k1 <= 27; k1 ++){
				int st1 = i1, ed1 = i1; 
				if (i1 == 27) st1 = 1, ed1 = 26;
				for (int i2 = st1; i2 <= ed1; i2 ++){
					int st2 = j1, ed2 = j1; 
					if (j1 == 27) st2 = 1, ed2 = 26;
					for (int j2 = st2; j2 <= ed2; j2 ++){
						int st3 = k1, ed3 = k1;
						if (k1 == 27) st3 = 1, ed3 = 26;
						for (int k2 = st3; k2 <= ed3; k2 ++){
							if (i2 == j2 && j2 == k2)
								d[0][0][i1][j1][k1] ++;
							if (i2 == j2 && j2 < k2)
								d[0][1][i1][j1][k1] ++;
							if (i2 < j2 && j2 == k2)
								d[0][2][i1][j1][k1] ++;
							if (i2 < j2 && j2 < k2)
								d[0][3][i1][j1][k1] ++; 
							if (i2 == j2)
								d[1][1][i1][j1][k1] ++;
							if (i2 < j2)
								d[1][3][i1][j1][k1] ++;
							if (j2 == k2)
								d[2][2][i1][j1][k1] ++;
							if (j2 < k2)
								d[2][3][i1][j1][k1] ++;
							d[3][3][i1][j1][k1] ++;
						}
					}
				
				}
			}
}
const int mod = 1e9 + 9;
void updata(long long &g, long long s){ g = (g + s) % mod; }
void work(){
	f[0][0] = 1;
	for (int i = 1; i <= n; i ++){
		f[i][0] = f[i][1] = f[i][2] = f[i][3] = 0;
		for (int j = 0; j < 4; j ++)
			for (int k = 0; k < 4; k ++)if (f[i - 1][k])
				updata(f[i][j], f[i - 1][k] * d[k][j][a[i]][b[i]][c[i]]);
		}
	printf("%d\n", (int)f[n][3]);
}
int main(){
	prapare(); int T = 0;
	scanf("%d", &T);
	while (T--){
		init();work();
	}
	return 0;
}