arc119f

自动机写法。开始在做的时候题解没讲每个节点代表什么状态,自己推了一遍,记录一下。

思路

计数,求有多少种替换方式使得 存在一条长度小于等于 的路径。

可以做 的 dp。设 表示前 个位置,最近的 分别在 。官方题解是优化 ,发现有意义状态下 的差不超过一定数值。

观察发现,在 处,肯定选择从 直接跳过去。具体的,当有连续的 后,会选择直接从两端的 越过。所以有效的状态并不多。

考虑建一个自动机来求出任意状态对应的下一个状态。由于 可以取反得到,先只考虑当前在 ,在 同理。以下几种有用状态。

  • 状态 。当前在 ,上一个位置为 ,记为

  • 状态 即状态 反过来,记为

  • 状态 。当前在 ,上一个位置为 ,下一个位置为 ,记为

  • 状态 。当前在 ,上一个位置为 ,下两个位置为 ,记为

  • 状态 。当前在 ,上一个位置为 ,下三个位置为 ,记为 。此时如果往后有 ,直接先后退越过;否则走向自己,不会有 的状态。注意到此时无论如何都会后退,所以稍微调整一下。定义状态 为当前在 ,接下来有连续 且一定会选择跳过。

  • 状态 。当前在 ,上一个位置不重要,下一个位置为 ,记为 。此时可以跳到下一个 ,也可以爬到 ,但不知道走一步后要去 还是

  • 状态 。当前既可以当作在 ,也可以当作在 ,记为 ,称为状态 。状态 表示当前在 ,下一个为 ,记为

综上,有 种有用状态。

  • 上面 种状态取反得到当前在 的另外 种。

  • ,,

考虑转移。记 表示当前状态为 ,加入 ,走向哪一个状态, 为转移的代价。

  • 状态 。如果加入 ,走向状态 ,点没有移动,代价为 。如果加入 ,走向状态 ,点没有移动,代价为

  • 状态 。如果加入 ,走向状态 ,点没有移动,代价为 。如果加入 ,走向状态 ,点向后移动一格,代价为

  • 状态 。如果加入 ,走向状态 ,点后退一步到上一个 ,代价为 。如果加入 ,可以先向后再向前越过 到下一个 ,也可以爬两步 ,即走向状态 ,代价为

  • 状态 。注意这里状态不同于以前,当前点为 。如果加入 ,走向状态自己 ,点不动,代价为 。如果加入 ,向前越过 到下一个 ,即走向状态 ,代价为

  • 状态 。如果加入 ,可以向后走到 ,也可以跳到下一个 ,即走向状态 ,代价为 。如果加入 ,虽然没有连续 ,但一定会从 跳过这段 ,即状态 ,代价为

  • 状态 。如果加入 ,虽然没有连续 ,但一定会把 当作 跳过这段 ,即状态 ,代价为 。如果加入 ,可以向后走到 ,也可以把 当作 跳到下一个 ,即走向状态 ,代价为

  • 状态 。如果加入 ,走向状态 ,代价为 。如果加入 ,走向状态 ,代价为

分析完自动机状态的转移,dp 即可。设 表示加入前 位,走了 步,当前状态为

注意到加入一位并不是走到一位。状态设定的当前走到的点不是状态的最后一位。记录 表示状态的最后一位是 ,状态设定的当前走到的点离终点的距离,稍微处理即可。

复杂度

code

int n,m,ans;
char c[maxn];
int tree[15][2],val[15][2];
int dis[7];
int dp[maxn][maxn][15];
signed main(){
	n=read();m=read();
	scanf("%s",c+1);
	tree[0][0]=2;val[0][0]=0;tree[0][1]=8;val[0][1]=0;
	tree[2][0]=4;val[2][0]=0;tree[2][1]=8;val[2][1]=1;
	tree[4][0]=6;val[4][0]=1;tree[4][1]=12;val[4][1]=2;
	tree[6][0]=6;val[6][0]=0;tree[6][1]=1;val[6][1]=1;
	tree[8][0]=12;val[8][0]=1;tree[8][1]=7;val[8][1]=0;
	tree[10][0]=6;val[10][0]=0;tree[10][1]=12;val[10][1]=1;
	tree[12][0]=10;val[12][0]=0;tree[12][1]=11;val[12][1]=0;
	for(int i=0;i<=5;i++){//状态取反。
		for(int k=0;k<=1;k++){
			if(tree[i*2][k^1]==12)tree[i*2+1][k]=tree[i*2][k^1];
			else tree[i*2+1][k]=tree[i*2][k^1]^1;
			val[i*2+1][k]=val[i*2][k^1];
		}
	}
//	for(int i=0;i<=12;i++)cout<<tree[i][0]<<" "<<val[i][0]<<" "<<tree[i][1]<<" "<<val[i][1]<<"\n";
	dis[0]=dis[3]=dis[4]=dis[5]=dis[6]=1;
	dis[1]=dis[2]=2;
	dp[0][0][12]=1;
	for(int i=0;i<=n-2;i++){
		for(int j=0;j<=m;j++){
			for(int k=0;k<=12;k++){
				if(dp[i][j][k]){
					if(c[i+1]=='A'){
						dp[i+1][j+val[k][0]][tree[k][0]]+=dp[i][j][k];
						dp[i+1][j+val[k][0]][tree[k][0]]%=mod;
					}
					else if(c[i+1]=='B'){
						dp[i+1][j+val[k][1]][tree[k][1]]+=dp[i][j][k];
						dp[i+1][j+val[k][1]][tree[k][1]]%=mod;
					}
					else{
						dp[i+1][j+val[k][0]][tree[k][0]]+=dp[i][j][k];
						dp[i+1][j+val[k][0]][tree[k][0]]%=mod;
						dp[i+1][j+val[k][1]][tree[k][1]]+=dp[i][j][k];
						dp[i+1][j+val[k][1]][tree[k][1]]%=mod;
					}
				}
			}
		}
	}
	for(int j=0;j<=m;j++){
		for(int k=0;k<=12;k++)if(j+dis[k>>1]<=m){
			ans+=dp[n-1][j][k];
			ans%=mod;
		}
	}
	printf("%lld\n",ans);
}