USACO22OPEN Pair Programming G

P8273 [USACO22OPEN] Pair Programming G

题目大意

一个程序由一系列指令组成,每条指令的类型如下:

  • × d \times d ×d,其中 d d d是一个 [ 0 , 9 ] [0,9] [0,9]范围内的整数
  • + s +s +s,其中 s s s是一个表示变量名称的字符串,不同操作的变量名称互不相同

指令字符串中一个数字 d d d表示操作一,一个加号 + + +表示操作二。

B B B和小 E E E各有一个有 n n n条指令的程序,交错这些程序可以得到一个有 2 n 2n 2n个指令的新程序。计算执行这些交错程序可能得到的不同表达式的数量,输出答案对 1 0 9 + 7 10^9+7 109+7取模后的值。

T T T组数据。

1 ≤ T ≤ 10 , ∑ n ≤ 2000 1\leq T\leq 10,\sum n\leq 2000 1T10,n2000


题解

首先,要对小 B B B和小 E E E的字符串做一些处理。如果有 × 0 \times 0 ×0的操作,则之前的部分都没用了,但注意这个 × 0 \times 0 ×0还有用。如果有 × 1 \times 1 ×1操作,则对答案没有影响,要将其舍去。

f i , j , 0 / 1 f_{i,j,0/1} fi,j,0/1表示小 B B B的串匹配到第 i i i位,小 E E E的串匹配到第 j j j位,最后一个是小 B B B的指令还是小 E E E的指令的答案。

我们知道,加法和乘法都满足交换律,即 a + b = b + a , a × b = b × a a+b=b+a,a\times b=b\times a a+b=b+a,a×b=b×a。也就是说,如果有两个类型相同且位置相邻的操作,则两个操作的先后顺序对最后得到的表达式没有影响。

那么,如果两个字符串当前的操作类型相同,则先执行小 E E E的操作。

由此可推得转移式
f i + 1 , j , 0 = f i , j , 0 + f i , j , 1 f_{i+1,j,0}=f_{i,j,0}+f_{i,j,1} fi+1,j,0=fi,j,0+fi,j,1
f i , j + 1 , 1 = f i , j , 0 × [ S i ≠ T j + 1 ] + f i , j , 1 f_{i,j+1,1}=f_{i,j,0}\times [S_i\neq T_{j+1}]+f_{i,j,1} fi,j+1,1=fi,j,0×[Si=Tj+1]+fi,j,1

这样的话,如果有一段相同类型的指令(其中包含小 B B B的指令和小 E E E的指令),则计算完一段小 E E E的指令再执行小 B B B的指令后,因为有 × [ S i ≠ T j + 1 ] \times [S_i\neq T_{j+1}] ×[Si=Tj+1]的限制,所以不会再倒回来执行小 E E E的指令,那么小 B B B的指令后的第一个小 E E E的指令一定是不同类型的。所以,这个转移式只会计算先执行的小 E E E的指令,再执行小 B B B的指令的贡献,这样就能保证没有表达式会被重复计算。

注意一开始 f 0 , 0 , 1 = 1 f_{0,0,1}=1 f0,0,1=1,如果将 1 1 1赋值到 f 0 , 0 , 0 f_{0,0,0} f0,0,0则会在下一次转移中被计算两次,所以要将 1 1 1赋值到 f 0 , 0 , 1 f_{0,0,1} f0,0,1才能使开始的空字符串只被计算一次。

时间复杂度为 O ( ∑ n 2 ) O(\sum n^2) O(n2)

code

#include<bits/stdc++.h>
using namespace std;
const long long mod=1e9+7;
int T,n,s1,t1;
long long f[2005][2005][2];
char s[2005],t[2005];
void dd(int &a1,char a[]){
	a1=0;
	for(int i=1;i<=n;i++){
		if(a[i]=='0') a1=0;
		else if(a[i]=='1') continue;
		if(a[i]!='+') a[i]='*';
		a[++a1]=a[i];
	}
}
int main()
{
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n);
		scanf("%s%s",s+1,t+1);
		dd(s1,s);dd(t1,t);
		f[0][0][1]=1;
		for(int i=0;i<=s1;i++){
			for(int j=0;j<=t1;j++){
				if(i<s1) f[i+1][j][0]=(f[i][j][0]+f[i][j][1])%mod;
				if(j<t1){
					f[i][j+1][1]=f[i][j][1];
					if(i&&s[i]!=t[j+1])
					f[i][j+1][1]=(f[i][j+1][1]+f[i][j][0])%mod;
				}
			}
		}
		printf("%lld\n",(f[s1][t1][0]+f[s1][t1][1])%mod);
	}
	return 0;
}
  • 3
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值