自动机写法。开始在做的时候题解没讲每个节点代表什么状态,自己推了一遍,记录一下。
思路
计数,求有多少种替换方式使得 到 存在一条长度小于等于 的路径。
可以做 的 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);
}