codeforces750E New Year and Old Subsequence 矩阵dp + 线段树

题目传送门

思路：

　　先看一个大牛的题解

　　题解里面对矩阵的构造已经写的很清楚了，其实就是因为在每个字符串都有固定的很多中状态，刚好可以用矩阵来表达，所以$(i,j)$这种状态可以通过两个相邻的矩阵的$min(i,k)+(k,j)$得到，取最小值即可，由于这是一个区间问题，所以用线段树来维护区间的矩阵运算，这个运算就是取min的过程。

　　虽然这道原题被出在2019icpc南昌网络赛中了，但这个做法以前确实没有遇见过，开阔了思路。

　　代码和博客里的其实几乎一样。

#pragma GCC optimize (2)
#pragma G++ optimize (2)
#pragma comment(linker, "/STACK:102400000,102400000")
#include<bits/stdc++.h>
#include<unordered_map>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define dep(i,b,a) for(int i=b;i>=a;--i)
#define clr(a,b) memset(a,b,sizeof(a))
#define pb push_back
#define pii pair<int,int >
using namespace std;
typedef long long ll;
ll rd()
{
    ll x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
const int maxn = 2e5+;
int n,m;
char s[maxn];
struct matrix{
    int a[][];
    matrix(){clr(a,0x3f);}
    matrix operator *(matrix x){
        matrix res;
        rep(i,,){
            rep(j,,){
                rep(k,,){
                    res.a[i][j]=min(res.a[i][j],a[i][k]+x.a[k][j]);
                }
            }
        }
        return res;
    }
};
struct Tree{
    int l,r;
    matrix M;
}t[maxn<<];
void pushup(int x){
    t[x].M=t[x<<].M*t[x<<|].M;
}
void build(int x,int l,int r){
    t[x].l=l,t[x].r=r;int mid=(l+r)>>;
    if(l == r){
        for(int i = ; i < ; i++) t[x].M.a[i][i] = ;
        if(s[l] == '') t[x].M.a[][] = , t[x].M.a[][] = ;
        if(s[l] == '') t[x].M.a[][] = , t[x].M.a[][] = ;
        if(s[l] == '') t[x].M.a[][] = , t[x].M.a[][] = ;
        if(s[l] == '') t[x].M.a[][] = , t[x].M.a[][] = ;
        if(s[l] == '') t[x].M.a[][] = , t[x].M.a[][] = ;
        return;
    }
    build(x<<,l,mid),build(x<<|,mid+,r);
    pushup(x);
}
matrix query(int x,int l,int r){
    if(l<=t[x].l&&t[x].r<=r)return t[x].M;
    int mid=(t[x].l+t[x].r)>>;
    if(r<=mid)return query(x<<,l,r);
    if(l>mid)return query(x<<|,l,r);
    return query(x<<,l,r)*query(x<<|,l,r);
}
int main(){
    int x,y;
    cin>>n>>m>>(s+);
    build(,,n);
    rep(i,,m){
        scanf("%d%d",&x,&y);
        matrix ans=query(,x,y);
        printf("%d\n",ans.a[][]>n?-:ans.a[][]);
    }
}