Good Bye 2016 E. New Year and Old Subsequence

题意

自己看英文

题解

左老师在KOI讲的题，没看题就听题解了QAQ
还没想过呢。。
线段树维护矩阵，好像很难啊。。
然后今天做了一下，发现也不怎么难
就线段树每个节点维护一个DP数组，表示从i状态到j状态所需的代价
0:什么都没有
1:2
2:20
3:201
4:2017
CODE:

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int MAX=1<<28;
const int N=200005;
char ss[N];
struct matrix
{
    int s[5][5];
    matrix operator + (const matrix x) const
    {
        matrix a;
        for (int u=0;u<5;u++)
            for (int i=0;i<5;i++)
            {
                a.s[u][i]=MAX;
                for (int j=0;j<5;j++)
                {
                //  printf("u:%d i:%d j:%d s[u][j]:%d s[j][i]:%d\n",u,i,j,s[u][j]);
                    a.s[u][i]=min(a.s[u][i],s[u][j]+x.s[j][i]);
                }
            }
        return a;
    }
};
struct qq
{
    int l,r;
    int s1,s2;
    matrix c;
}tr[N<<1];int num;
void bt (int l,int r)
{
    int a=++num;
    tr[a].l=l;tr[a].r=r;
    if (l==r)
    {
        for (int u=0;u<5;u++)
            for (int i=0;i<5;i++)
                tr[a].c.s[u][i]=(u==i)?0:MAX;
        if (ss[l]=='2') {tr[a].c.s[0][0]=1;tr[a].c.s[0][1]=0;}
        if (ss[l]=='0') {tr[a].c.s[1][1]=1;tr[a].c.s[1][2]=0;}
        if (ss[l]=='1') {tr[a].c.s[2][2]=1;tr[a].c.s[2][3]=0;}
        if (ss[l]=='7') {tr[a].c.s[3][3]=1;tr[a].c.s[3][4]=0;}
        if (ss[l]=='6') {tr[a].c.s[3][3]=1;tr[a].c.s[4][4]=1;}
    /*  for (int u=0;u<4;u++)
    {
        for (int i=0;i<4;i++)
            printf("%d ",tr[a].c.s[u][i]);
        printf("\n");
    }*/
        return ;
    }
    int mid=(l+r)>>1;
    tr[a].s1=num+1;bt(l,mid);
    tr[a].s2=num+1;bt(mid+1,r);
    tr[a].c=tr[tr[a].s1].c+tr[tr[a].s2].c;
/*  printf("%d %d\n",l,r);
    for (int u=0;u<4;u++)
    {
        for (int i=0;i<4;i++)
            printf("%d ",tr[a].c.s[u][i]);
        printf("\n");
    }
    system("pause");*/
}
matrix find (int now,int l,int r)
{
/*  printf("%d %d %d %d %d\n",now,tr[now].l,tr[now].r,l,r);
    system("pause");*/
    if (tr[now].l==l&&tr[now].r==r) return tr[now].c;
    int s1=tr[now].s1,s2=tr[now].s2;
    int mid=(tr[now].l+tr[now].r)>>1;
    if (r<=mid) return find(s1,l,r);
    else if (l>mid) return find(s2,l,r);
    else return find(s1,l,mid)+find(s2,mid+1,r);
}
int main()
{
    int n,q;
    scanf("%d%d",&n,&q);
    scanf("%s",ss+1);
    bt(1,n);
    for (int u=1;u<=q;u++)
    {
        int l,r;
        scanf("%d%d",&l,&r);
        matrix ans=find(1,l,r);
        if (ans.s[0][4]==MAX) printf("-1\n");
        else printf("%d\n",ans.s[0][4]);
    }
    return 0;
}