【codeforces】2014-2015 ACM-ICPC, Central Europe Regional Contest (CERC 14) F Vocabulary【dp】

最新推荐文章于 2020-05-12 14:15:15 发布
最新推荐文章于 2020-05-12 14:15:15 发布
阅读量2k 收藏
点赞数
分类专栏： DP 文章标签： codeforces
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u013368721/article/details/42032563
版权

传送门：【codeforces】2014-2015 ACM-ICPC, Central Europe Regional Contest (CERC 14) F Vocabulary


题目分析：蒽。。。。这个dp嘛其实也是满简单易懂的，就是要考虑的地方略多一点。。。

设字典序最小的为a串，中间的为b串，最大的为c串。

首先我们预处理出'?'个数在前i个字符位置的个数num[i]（三个字符串的都算上）。

然后我们分三种情况讨论：

1.第i个位置之前三个串的字典序相同（每个位置三串的字母都相同），第i+1个位置a[i + 1] < b[i + 1] < c[i + 1]。

2.第i个位置之前a、b串的字典序相同，且b串的字典序小于c串，第i+1个位置a[i + 1] < b[i + 1]。

3.第i个位置之前b、c串的字典序相同，且a串的字典序小于b串，第i+1个位置b[i + 1] < c[i + 1]。

当第i个位置三个字符串的大小可以区分时，后面的'?'都可以随便取了，所以要乘上26^(num[n-1]-num[i])，其中n为三个字符串长度的最大值，字符串下标从0开始。

对于每一种情况，第i个位置开始区分字典序时，就是ans+=dp[x][i-1]*(当前的情况数)*(26^(num[n-1]-num[i]))。

dp[x][i - 1]*(当前的情况数)表示情况x在第i个位置符合情况x的定义的方案数。


剩下的就是无比繁琐的if了。。。。（对我来说是这样。。整个人都不好了，写了两天都要写晕了。。。各种细节问题）


代码如下：


#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;

typedef long long LL ;

#define rep( i , a , b ) for ( int i = ( a ) ; i <  ( b ) ; ++ i )
#define For( i , a , b ) for ( int i = ( a ) ; i <= ( b ) ; ++ i )
#define rev( i , a , b ) for ( int i = ( a ) ; i >= ( b ) ; -- i )
#define clr( a , x ) memset ( a , x , sizeof a )

const int MAXN = 1000005 ;
const int mod = 1e9 + 9 ;

int pow[MAXN * 3] ;
int num[MAXN] ;
char a[MAXN] ;
char b[MAXN] ;
char c[MAXN] ;
int na , nb , nc , n ;
LL dpa[MAXN] , dpb[MAXN] , dpc[MAXN] ;

int sigma ( int a , int b ) {
	return ( a + b ) * ( b - a + 1 ) / 2 ;
}

void solve () {
	int ans = 0 ;
	scanf ( "%s%s%s" , a , b , c ) ;
	na = strlen ( a ) ;
	nb = strlen ( b ) ;
	nc = strlen ( c ) ;
	n = max ( max ( na , nb ) , nc ) ;
	rep ( i , 0 , n ) {
		num[i] = 0 ;
		dpa[i] = 0 ;
		dpb[i] = 0 ;
		dpc[i] = 0 ;
	}
	dpa[0] = 1 ;
	rep ( i , 0 , n ) {
		if ( i < na ) {
			if ( a[i] == '?' ) ++ num[i] ;
			else a[i] -= 'a' - 1 ;
		} else a[i] = 0 ;
		if ( i < nb ) {
			if ( b[i] == '?' ) ++ num[i] ;
			else b[i] -= 'a' - 1 ;
		} else b[i] = 0 ;
		if ( i < nc ) {
			if ( c[i] == '?' ) ++ num[i] ;
			else c[i] -= 'a' - 1 ;
		} else c[i] = 0 ;
	}
	rep ( i , 1 , n ) num[i] += num[i - 1] ;
	int ans1 = 0 , ans2 = 0 , ans3 = 0 ;
	rep ( i , 0 , n ) {
		int Pow = pow[num[n - 1] - num[i]] ;
		int cnt = ( a[i] == '?' ) + ( b[i] == '?' ) + ( c[i] == '?' ) ;
		if ( cnt == 3 ) {
			ans1 = ( ans1 + dpa[i] * 2600 % mod * Pow ) % mod ;
			dpa[i + 1] = dpa[i] * 26 % mod ;
		}
		if ( cnt == 2 ) {
			if ( a[i] != '?' ) {
				ans1 = ( ans1 + dpa[i] * sigma ( 1 , 26 - a[i] - 1 ) % mod * Pow ) % mod ;
				if ( a[i] == 0 ) break ;
				dpa[i + 1] = dpa[i] ;
			}
			if ( b[i] != '?' ) {
				if ( b[i] == 0 ) break ;
				ans1 = ( ans1 + dpa[i] * ( b[i] - 1 ) * ( 26 - b[i] ) % mod * Pow ) % mod ;
				dpa[i + 1] = dpa[i] ;
			}
			if ( c[i] != '?' ) {
				if ( c[i] == 0 ) break ;
				ans1 = ( ans1 + dpa[i] * sigma ( 1 , c[i] - 2 ) % mod * Pow ) % mod ;
				dpa[i + 1] = dpa[i] ;
			}
		}
		if ( cnt == 1 ) {
			if ( a[i] == '?' ) {
				if ( b[i] == 0 || c[i] == 0 ) break ;
				if ( b[i] < c[i] ) ans1 = ( ans1 + dpa[i] * ( b[i] - 1 ) % mod * Pow ) % mod ;
				if ( b[i] != c[i] ) break ;
				dpa[i + 1] = dpa[i] ;
			}
			if ( b[i] == '?' ) {
				if ( c[i] == 0 ) break ;
				if ( a[i] < c[i] ) ans1 = ( ans1 + dpa[i] * ( c[i] - a[i] - 1 ) % mod * Pow ) % mod ;
				if ( a[i] == 0 ) break ;
				if ( a[i] != c[i] ) break ;
				dpa[i + 1] = dpa[i] ;
			}
			if ( c[i] == '?' ) {
				if ( b[i] == 0 ) break ;
				if ( a[i] < b[i] ) ans1 = ( ans1 + dpa[i] * ( 26 - b[i] ) % mod * Pow ) % mod ;
				if ( a[i] == 0 ) break ;
				if ( a[i] != b[i] ) break ;
				dpa[i + 1] = dpa[i] ;
			}
		}
		if ( cnt == 0 ) {
			if ( a[i] <  b[i] && b[i] <  c[i] ) ans1 = ( ans1 + dpa[i] * Pow ) % mod ;
			if ( a[i] != b[i] || b[i] != c[i] ) break ;
			dpa[i + 1] = dpa[i] ;
		}
	}
	rep ( i , 0 , n ) {
		int Pow = pow[num[n - 1] - num[i]] ;
		int cnt = ( a[i] == '?' ) + ( b[i] == '?' ) + ( c[i] == '?' ) ;
		if ( cnt == 3 ) {
			ans2 = ( ans2 + dpb[i] * sigma ( 1 , 25 ) * 26 % mod * Pow ) % mod ;
			dpb[i + 1] = ( dpb[i] * 26 * 26 + dpa[i] * sigma ( 1 , 25 ) ) % mod ;
		}
		if ( cnt == 2 ) {
			if ( a[i] != '?' ) {
				ans2 = ( ans2 + dpb[i] * ( 26 - a[i] ) * 26 % mod * Pow ) % mod ;
				if ( a[i] == 0 ) break ;
				dpb[i + 1] = ( dpb[i] * 26 + dpa[i] * ( 26 - a[i] ) ) % mod ;
			}
			if ( b[i] != '?' ) {
				if ( b[i] == 0 ) break ;
				ans2 = ( ans2 + dpb[i] * ( b[i] - 1 ) * 26 % mod * Pow ) % mod ;
				dpb[i + 1] = ( dpb[i] * 26 + dpa[i] * ( 26 - b[i] ) ) % mod ;
			}
			if ( c[i] != '?' ) {
				ans2 = ( ans2 + dpb[i] * sigma ( 1 , 25 ) % mod * Pow ) % mod ;
				if ( c[i] > 1 ) dpb[i + 1] = ( dpb[i] * 26 + dpa[i] * ( c[i] - 1 ) ) % mod ;
				else dpb[i + 1] = dpb[i] * 26 % mod ;
			}
		}
		if ( cnt == 1 ) {
			if ( a[i] == '?' ) {
				if ( b[i] == 0 ) break ;
				ans2 = ( ans2 + dpb[i] * ( b[i] - 1 ) % mod * Pow ) % mod ;
				if ( b[i] < c[i] ) dpb[i + 1] = ( dpb[i] + dpa[i] ) % mod ;
				else dpb[i + 1] = dpb[i] ;
			}
			if ( b[i] == '?' ) {
				ans2 = ( ans2 + dpb[i] * ( 26 - a[i] ) % mod * Pow ) % mod ;
				if ( a[i] == 0 ) break ;
				if ( a[i] < c[i] ) dpb[i + 1] = ( dpb[i] + dpa[i] ) % mod ;
				else dpb[i + 1] = dpb[i] ;
			}
			if ( c[i] == '?' ) {
				if ( a[i] <  b[i] ) ans2 = ( ans2 + dpb[i] * 26 % mod * Pow ) % mod ;
				if ( a[i] != b[i] ) break ;
				dpb[i + 1] = ( dpb[i] * 26 + dpa[i] * ( 26 - b[i] ) ) % mod ;
			}
		}
		if ( cnt == 0 ) {
			if ( a[i] <  b[i] ) ans2 = ( ans2 + dpb[i] * Pow ) % mod ;
			if ( a[i] != b[i] ) break ;
			if ( b[i] <  c[i] ) dpb[i + 1] = ( dpb[i] + dpa[i] ) % mod ;
			else dpb[i + 1] = dpb[i] ;
		}
	}
	rep ( i , 0 , n ) {
		int Pow = pow[num[n - 1] - num[i]] ;
		int cnt = ( a[i] == '?' ) + ( b[i] == '?' ) + ( c[i] == '?' ) ;
		if ( cnt == 3 ) {
			ans3 = ( ans3 + dpc[i] * sigma ( 1 , 25 ) * 26 % mod * Pow ) % mod ;
			dpc[i + 1] = ( dpc[i] * 26 * 26 + dpa[i] * sigma ( 1 , 25 ) ) % mod ;
		}
		if ( cnt == 2 ) {
			if ( a[i] != '?' ) {
				ans3 = ( ans3 + dpc[i] * sigma ( 1 , 25 ) % mod * Pow ) % mod ;
				dpc[i + 1] = ( dpc[i] * 26 + dpa[i] * ( 26 - a[i] ) ) % mod ;
			}
			if ( b[i] != '?' ) {
				ans3 = ( ans3 + dpc[i] * 26 * ( 26 - b[i] ) % mod * Pow ) % mod ;
				if ( b[i] == 0 ) break ;
				dpc[i + 1] = ( dpc[i] * 26 + dpa[i] * ( b[i] - 1 ) ) % mod ;
			}
			if ( c[i] != '?' ) {
				if ( c[i] == 0 ) break ;
				ans3 = ( ans3 + dpc[i] * 26 * ( c[i] - 1 ) % mod * Pow ) % mod ;
				dpc[i + 1] = ( dpc[i] * 26 + dpa[i] * ( c[i] - 1 ) ) % mod ;
			}
		}
		if ( cnt == 1 ) {
			if ( a[i] == '?' ) {
				if ( b[i] < c[i] ) ans3 = ( ans3 + dpc[i] * 26 % mod * Pow ) % mod ;
				if ( b[i] != c[i] ) break ;
				dpc[i + 1] = ( dpc[i] * 26 + dpa[i] * ( b[i] - 1 ) ) % mod ;
			}
			if ( b[i] == '?' ) {
				if ( c[i] == 0 ) break ;
				ans3 = ( ans3 + dpc[i] * ( c[i] - 1 ) % mod * Pow ) % mod ;
				if ( a[i] < c[i] ) dpc[i + 1] = ( dpc[i] + dpa[i] ) % mod ;
				else dpc[i + 1] = dpc[i] ;
			}
			if ( c[i] == '?' ) {
				ans3 = ( ans3 + dpc[i] * ( 26 - b[i] ) % mod * Pow ) % mod ;
				if ( b[i] == 0 ) break ;
				if ( a[i] < b[i] ) dpc[i + 1] = ( dpc[i] + dpa[i] ) % mod ;
				else dpc[i + 1] = dpc[i] ;
			}
		}
		if ( cnt == 0 ) {
			if ( b[i] <  c[i] ) ans3 = ( ans3 + dpc[i] * Pow ) % mod ;
			if ( b[i] != c[i] ) break ;
			if ( a[i] <  b[i] ) dpc[i + 1] = ( dpc[i] + dpa[i] ) % mod ;
			else dpc[i + 1] = dpc[i] ;
		}
	}
	/*rep ( i , 0 , n ) printf ( "%d " , dpa[i] ) ;
	printf ( "\n" ) ;
	rep ( i , 0 , n ) printf ( "%d " , dpb[i] ) ;
	printf ( "\n" ) ;
	rep ( i , 0 , n ) printf ( "%d " , dpc[i] ) ;
	printf ( "\n" ) ;
	printf ( "%d %d %d\n" , ans1 , ans2 , ans3 ) ;*/
	ans = ( ( LL ) ans1 + ans2 + ans3 ) % mod ;
	printf ( "%d\n" , ans ) ;
}

int main () {
	//freopen ( "1.in" , "r" , stdin ) ;
	//freopen ( "12.out" , "w" , stdout ) ;
	int T ;
	scanf ( "%d" , &T ) ;
	pow[0] = 1 ;
	For ( i , 1 , 3000000 ) pow[i] = pow[i - 1] * 26LL % mod ;
	while ( T -- ) solve () ;
	return 0 ;
}


poursoul
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
2014-2015 ACM-ICPC, Central Europe Regional Contest (CERC 14) [Gym-100543G]
qq965194745的博客
07-18 700
2014-2015 ACM-ICPC, Central Europe Regional Contest (CERC 14) [Gym-100543G]最后一题，训练完看了别人的代码，发现是区间dp，dp的是时间区间，先要将时间离散化，然后是时间的区间的选择问题，当时就是不知道要怎么组合区间的选择，画图往贪心去想只能推出一个最少要攻击多少次的次数，没有想到可以用区间dp考虑这种麻烦选择。#inclu
2014-2015 ACM-ICPC Central Europe Regional Contest (CERC 2014)
weixin_30481087的博客
05-04 264
原标题应该叫“SCUT五四青年节软件集训队快乐训练赛”。本来应该很快乐才对，但是没想到毒瘤BPM136居然挑了套区域赛题(虽然不是很难)，而且队友基本有事导致变成个人solo，这就非常不快乐了。 而且题风有点诡异，做完四道签到之后剩下的要么看不懂要么就是各种dp。 题目链接：http://codeforces.com/gym/100543 C： 一点意思都没有的签到题，因为没留意爆了i...
BZOJ 4043 [Cerc2014] Vocabulary
dilei7762的博客
11-26 141
题目链接 https://www.lydsy.com/JudgeOnline/problem.php?id=4043 题解 首先补全字符串。 设f[i][j]f[i][j]f[i][j]表示字符串的前iii位，jjj为000代表三个串字典序都相等，111代表前两个串字典序相等而第二个<第三个，222代表第一个<第二个=第三个，3...
[ CERC 2014 ] Vocabulary
weixin_34273479的博客
10-07 120
\(\\\) \(Description\) 给出三个长度分别为 \(lenA,lenB,lenC\) 的三个字符串 \(A,B,C\) ，其中字符集只包括所有小写字母以及 \(?\) 号。 现在将所有 \(?\) 号改为任意小写字母，问有多少种方案，使得 \(A,B,C\) 三个字符串按照字典序排列\((\)相同不算\()?\) 输出方案数对 \(10^9+9\) 取模后的结果。多组询问。 ...
cerc2014 Vocabulary
芒果街上的小屋（此博客暂时停用）
08-25 578
给你三个字符串，这些字符串有些单词模糊不可认了，用"?"来代表。 现在你可以用任意英文小写字母来代表它们。要求是使得给定的三个字符串中 所有的"?"被你认定的字母代替后，各不相同且按字典序出现。问有多少种方式。 Input 先给出一个数字N，代表数据组数。 接下来3*N行，每行给出一个字符串。长度 Output 输出结果 Mod 10^9+9
ACM-ICPC-制备：ACM-ICPC制备指南
02-02
在计算机科学领域，ACM国际大学生程序设计竞赛（ACM International Collegiate Programming Contest，简称ACM-ICPC）是一项备受瞩目的赛事，它挑战参赛者的编程技能、算法理解和团队协作能力。对于希望在ACM-ICPC中...
TrainACM-ICPC-OLP:来自Kattis和Codeforce的ACM ICPC竞赛解决方案
04-09
17/03/2021-10/11问题） 链接比赛代码 2020-2021年ICPC-墨西哥大首府-雷佩查（日期：10/03/2021-12/14问题） 链接比赛代码 ACM-ICPC在2018-2019和2019-2020赛季链接到2018年至2020年训练竞赛ACM-ICPC（文件README2....
ICPC-ACM_C++_Codeforces_
09-30
【ICPC-ACM_C++_Codeforces_] 这个标题揭示了我们即将探讨的是一个与ACM（国际大学生程序设计竞赛）相关的主题，而使用的编程语言是C++，并且这个话题与Codeforces平台上的算法比赛题目有关。Codeforces是一个全球...
Codeforces-Solution
05-03
"Codeforces-Solution"很可能是一个包含针对Codeforces上问题的解题代码库，主要使用C++语言，部分可能涉及cpp14标准。 在这个压缩包中，"Codeforces-Solution-master"可能是这个代码库的主分支或者最新版本，这...
codeforces-ACM竞赛题目-2833道.tgz
10-25
codeforces-ACM竞赛题目-2833道.tgz
【CERC2014】【BZOJ4043】Vocabulary
CreationAugust is 14 years old forever
03-26 1359
Description 给你三个字符串，这些字符串有些单词模糊不可认了，用”?”来代表。 现在你可以用任意英文小写字母来代表它们。要求是使得给定的三个字符串中 所有的”?”被你认定的字母代替后，各不相同且按字典序出现。问有多少种方式。 Input 先给出一个数字N，代表数据组数。 接下来3*N行，每行给出一个字符串。长度 Output 输出结果 Mod 10^9+9 Sa
P4761 [CERC2014]Vocabulary（DP）
tomjobs的博客
05-12 302
题意： 给你三个字符串，这些字符串有些单词模糊不可认了，用"?“来代表。 现在你可以用任意英文小写字母来代表它们。要求是使得给定的三个字符串中 所有的”?"被你认定的字母代替后，各不相同且按字典序出现。问有多少种方式。 输入：先给出一个数字N，代表数据组数。 接下来3*N行，每行给出一个字符串。长度<=1000 000 输出：输出结果 Mod 10^9+9 参考自：https://www.cnblogs.com/SGCollin/p/9751265.html 思路： 将这三个字符串一起处理，’?'特殊
【bzoj4043】【Cerc2014】【Vocabulary】【dp+预处理】
sunshinezff的专栏
03-27 691
Description 给你三个字符串，这些字符串有些单词模糊不可认了，用"?"来代表。  现在你可以用任意英文小写字母来代表它们。要求是使得给定的三个字符串中  所有的"?"被你认定的字母代替后，各不相同且按字典序出现。问有多少种方式。  Input 先给出一个数字N，代表数据组数。  接下来3*N行，每行给出一个字符串。长度 Output
bzoj3802 cerc2014 Vocabulary
u012647218的专栏
12-18 1100
3802: Vocabulary Time Limit: 10 Sec Memory Limit: 128 MB Submit: 54 Solved: 11 [Submit][Status] Description 给你三个字符串，这些字符串有些单词模糊不可认了，用"?"来代表。 现在你可以用任意英文小写字母来代表它们。要求是使得给定的三个字符串中 所有的"?"被你认定的字母
CERC 2014 B Mountainous landscape - 线段树 - 凸包 - 三分
Mys_C_K的博客
08-28 295
题目大意：给你平面上n个横坐标严格递增的点，对于每个i∈[1,n)i∈[1,n)i\in[1,n)，求一个jjj，满足PjPjP_j严格在直线PiPi+1PiPi+1P_iP_{i+1}上方。点坐标10910910^9，有10510510^5个点。 题解：答案可能不在凸包上，但如果凸包上都没有解那么就不可能有解。因此可以判断是否有解，扩展下去，可以先二分答案，然后看区间凸包是否有解；区间可以用线...
Cerc2014 Pork barre
abc473848880_的专栏
12-26 647
题目大意：给你n个点m条边，再给你q个询问，查询边权在l~r之中，且覆盖点最多的边权和最小是多少 n 这题比较神，看了题想了一个小时都还没有思路，最后还是看了一下题解才弄懂 本质就是求边权在l~r内的最小生成树 先把所有的边按照边权从大到小排序，然后每次把边权相同的插入图中维护好最小生成树，再用一个可持久化线段树维护权值就OK了 #include #include
2015-2016 ACM-ICPC, Central Europe Regional Contest (CERC 15)
banhe8035的博客
11-28 336
地址 Rank Solved A B C D E F G H I J K L 62/217 5/12 O O . O Ø O . . . Ø O . O: 当场通过 Ø: 赛后通过 .: 尚未通过 A ASCII Addition solved by chelly chelly's solution B Book Borders solved by chelly ...
Cerc2014 Outer space invaders
u012647218的专栏
12-25 1318
Cerc2014 Outer spaceinvaders      有N个外星人，第i个外星人会在ai时间出现，离你距离di，并且必须在bi时间之前被消灭。你有一把很NB的武器，攻击范围是个半径为R的圆，R可以任意调整，不过你以R的范围每攻击一次就要消耗R单位能量。外星人被攻击一次就会死掉。求需要消灭所有外星人的最小消耗能量。  范围  N 题解 我们首先将时间离散化，然后定义F[i]
2015-2016 ACM-ICPC, Central Europe Regional Contest (CERC 15)F
Soul
05-03 536
题意：给你F[1,i]和F[i,1] 接着有递推式F[i, j] = a ∗ F[i, j − 1] + b ∗ F[i − 1, j] + c.让你计算F[n,n] 题解：这题F[1,i]的贡献很显然是,F(i,1)的贡献同理 对于c的贡献我们每个点都计算很显然不好求我们可以做如下操作 并且那么 我们解得并且那么我们就可以得出以下式子 然后打个表就可以O(n)计算每个点的贡献 ...
CodeForces - 982C 树形DP
最新发布
07-25
引用\[1\]中提到了一种树形动态规划的方法来解决CodeForces - 982C问题。在这个问题中，subtree指的是子连通块，而不是子树。为了使cnt_white - cnt_black尽可能大，可以使用两次树形动态规划来求解。第一次是自底...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值