链接:https://nanti.jisuanke.com/t/41350
题意:
给一个只包含数字的字符串s,长度为n。
定义好的字符串:子串中(不一定连续)包含9012,不包含8012的字符串。
有Q次查询,每次查询给出l,r,求在这个区间中至少删除几个数字可以把这一段的字符串变成好的字符串。如果不能成为好字符串,则输出-1。
思路:
线段树+区间dp。可以先把字符串翻转,那么包含2019,不包含2018的字符串。我们用状态压缩来处理。0表示空串,1表示“2”,2表示“20”,3表示“201”,4表示“2019”。那么dp[i][j]就表示从i这个状态到j状态的最小花费。所以用线段树维护dp[5][5]的矩阵。初始化全为inf。区间查询和普通线段树有所不同,但是大体相似。详细见代码。
AC代码:
#include<bits/stdc++.h>
#include <stdlib.h>
#include<algorithm>
#include<complex>
#include<iostream>
#include<iomanip>
#include<ostream>
#include<cstring>
#include<string.h>
#include<string>
#include<cstdio>
#include<cctype>
#include<vector>
#include<queue>
#include<set>
#include<stack>
#include<map>
#include<cstdlib>
#include<time.h>
#include<ctime>
#include<bitset>
// #include<ext/pb_ds/assoc_container.hpp>
// #include<ext/pb_ds/hash_policy.hpp>
using namespace std;
// using namespace __gnu_pbds;
#define pb push_back
#define _fileout freopen("C:\\Users\\zsk18\\Desktop\\out.txt","w",stdout)
#define _filein freopen("C:\\Users\\zsk18\\Desktop\\in.txt","r",stdin)
#define ok(i) printf("ok%d\n",i)
#define gcd(a,b) __gcd(a,b) ;
typedef double db;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int>PII;
const double PI = acos(-1.0);
const ll MOD=1e9+7;
const ll NEG=1e9+6;
const int MAXN=2e5+10;
const int INF=0x3f3f3f3f;
const ll ll_INF=1e15;
const double eps=1e-9;
ll qm(ll a,ll b){ll ret=1;while(b){if(b&1)ret=ret*a%MOD;a=a*a%MOD;b>>=1;}return ret;}
ll Inv(ll x){return qm(x,MOD-2);}
struct st
{
int l,r;
int dp[5][5];
st(){memset(dp,INF,sizeof(dp));}
}t[MAXN*4];
char s[MAXN];
int n,m;
st get_new(char c,int l,int r)
{
st ret;
ret.l=l,ret.r=r;
for(int i=0;i<5;i++)
ret.dp[i][i]=0;
if(c=='2')ret.dp[0][1]=0,ret.dp[0][0]=1;
if(c=='0')ret.dp[1][2]=0,ret.dp[1][1]=1;
if(c=='1')ret.dp[2][3]=0,ret.dp[2][2]=1;
if(c=='9')ret.dp[3][4]=0,ret.dp[3][3]=1;
if(c=='8')ret.dp[3][3]=1,ret.dp[4][4]=1;
return ret;
}
st merge(st a,st b)
{
st ret;
ret.l=a.l,ret.r=b.r;
for(int i=0;i<5;i++)
for(int j=i;j<5;j++)
for(int k=i;k<=j;k++)
ret.dp[i][j]=min(ret.dp[i][j],a.dp[i][k]+b.dp[k][j]);
return ret;
}
void build(int l,int r,int p)
{
if(l==r)
{
t[p]=get_new(s[l],l,r);
return;
}
int mid=(l+r)>>1;
build(l,mid,p*2);
build(mid+1,r,p*2+1);
t[p]=merge(t[p*2],t[p*2+1]);
// printf("t[%d]")
}
st query(int l,int r,int p)
{
// printf("l=%d r=%d\n",l,r);
// printf("t[%d].l=%d t[%d].r=%d\n",p,t[p].l,p,t[p].r);
if(l<=t[p].l&&r>=t[p].r)return t[p];
int mid=(t[p].l+t[p].r)>>1;
// printf("mid=%d\n",mid);
if(r<=mid)return query(l,r,p*2);
if(l>mid)return query(l,r,p*2+1);
return merge(query(l,r,p*2),query(l,r,p*2+1));
}
int main()
{
scanf("%d%d",&n,&m);
scanf("%s",s+1);
reverse(s+1,s+1+n);
build(1,n,1);
while(m--)
{
int l,r;
scanf("%d%d",&l,&r);
st ans=query(n-r+1,n-l+1,1);
printf("%d\n",ans.dp[0][4]==INF?-1:ans.dp[0][4]);
}
return 0;
}