题意:
给定一个长为 n 的只含数字的串 s,给有 q 次询问,每次询问 s[l…r] 子串中至少删去多少数字才能使得不含 8102 子序列但含有 9102 子序列。(n, q <= 2e5)
链接:
https://nanti.jisuanke.com/t/41350
题解:
https://vjudge.net/problem/CodeForces-750E
原题为求不含 2016 但含 2017,这题只是把字符串反转了,其他的都一样。
X / Y 表示接收 X 输入后跳转的代价为 Y,对于串 s,令 Ms[ i ][ j ] 为接收 s 后从状态 i 跳转到状态 j 的最小代价,则对s = sp + sq,Ms = Mp * Mq(重载乘法,类似 floyed 的转移)。对于每个询问,只要得到区间的转移矩阵 M,则答案为 M[S][S2019],可用线段树区间合并来维护。
参考代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
#define pb push_back
#define sz(a) ((int)a.size())
#define mem(a, b) memset(a, b, sizeof a)
#define lson (rt << 1)
#define rson (rt << 1 | 1)
#define gmid (l + r >> 1)
const int maxn = 2e5 + 5;
const int maxm = 2e6 + 5;
const int mod = 1e9 + 7;
const int inf = 0x3f3f3f3f;
struct Ma{
int a[5][5];
Ma(int v = 0) {
mem(a, 0x3f);
for(int i = 0; i < 5; ++i) a[i][i] = v;
}
Ma operator * (const Ma &o) const{
Ma ret(inf);
for(int i = 0; i < 5; ++i)
for(int j = 0; j < 5; ++j)
for(int k = 0; k < 5; ++k)
ret.a[i][j] = min(ret.a[i][j], a[i][k] + o.a[k][j]);
return ret;
}
};
char s[maxn];
Ma sum[maxn << 2];
int n, q;
void add(int rt, int l){
if(s[l] == '2') sum[rt].a[0][0] = 1, sum[rt].a[0][1] = 0;
else if(s[l] == '0') sum[rt].a[1][1] = 1, sum[rt].a[1][2] = 0;
else if(s[l] == '1') sum[rt].a[2][2] = 1, sum[rt].a[2][3] = 0;
else if(s[l] == '9') sum[rt].a[3][3] = 1, sum[rt].a[3][4] = 0;
else if(s[l] == '8') sum[rt].a[3][3] = 1, sum[rt].a[4][4] = 1;
}
void pushUp(int rt){
sum[rt] = sum[lson] * sum[rson];
}
void build(int l, int r, int rt){
if(l == r) { add(rt, l); return; }
int mid = gmid;
build(l, mid, lson);
build(mid + 1, r, rson);
pushUp(rt);
}
void query(int l, int r, int rt, int L, int R, Ma &ret){
if(l >= L && r <= R){
ret = ret * sum[rt];
return;
}
int mid = gmid;
if(L <= mid) query(l, mid, lson, L, R, ret);
if(R > mid) query(mid + 1, r, rson, L, R, ret);
}
int main(){
scanf("%d%d%s", &n, &q, s + 1);
reverse(s + 1, s + 1 + n);
build(1, n, 1);
while(q--){
int l, r; scanf("%d%d", &l, &r);
l = n - l + 1, r = n - r + 1; swap(l, r);
Ma ret;
query(1, n, 1, l, r, ret);
printf("%d\n", ret.a[0][4] != inf ? ret.a[0][4] : -1);
}
}