2018 湖南多校（2）----CSU 2037 Mars (后缀自动机 + DFS）

题目链接：http://acm.csu.edu.cn/csuoj/problemset/problem?pid=2037

由于博主水平有限，无法解出这道题！查阅了 国科大大佬alpc_qleonardo的博文，感觉写的非常好，顿时思路变清晰，感觉对后缀自动机了解了一点，所以这篇文章只是对大佬文章的带有个人的理解的译文！

题意： 给你一个长度不超过10000的01字符串，不超过1000次查询区间[l,r]的字符串更改最少个字符后的新字符串不为原本串的子串.如果不存在，输出"Impossible".

1、 按照模板将原01字符串建立自动机数据结构，仅需要O(N)的存储空间；

2、 自动机的每条路径都是一个后缀，路径的每个前缀（未到终点），都是原
字符串的子串；另外用数组将查询区间储存，开始对其进行从左至右DFS遍历，
因为遍历到每一位字符，都有两种可能：更改和不改，对应影响是下一边是
1或0，若当前状态存在此出边，即遍历下一位，若不存在，即表明此前路径表示
的新字符串的前缀已经不是本串的子串了，剩下的字符不改变（答案求改变最少）
结束。

3、 DFS到了更改次数超过了此前遍历的current更改次数，退出，剪枝！

时间复杂度并不是很明显，但个人感觉应该有个乐观的结果！

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <vector>
#include <queue>
#include <stack>
using namespace std;

const int Size=20010;
struct Suffix_Automation{
    struct node_t{
        int pre;//指上一个可接收的Right集合的state节点的idx；
        int step;//建自动机时需用
        int son[2];//表明后面是否添i元素到达下一状态；
    }Node[Size];

    int cur,toUsed;

    void init(){
        cur=toUsed=1;
        memset(Node,0,sizeof(Node));
    }

    void add(char a,int idx){
        int sn=a-'0';
        int last=cur; cur=++toUsed;//此时cur为当前添加进来的节点的idx，last指向原来的cur；
        Node[cur].step=idx;

        //遍历所有可以接纳sn的状态；
        while(last&&Node[last].son[sn]==0){
            Node[last].son[sn]=cur;
            last=Node[last].pre;
        }

        //第一种情况，所有last后面都没有sn
        if(last==0) {Node[cur].pre=1;return ;}
        //第二种情况，last后面有sn
        int q=Node[last].son[sn];
        //直接跳到cur，节省空间
        if(Node[last].step==Node[q].step-1){
            Node[cur].pre=q;
            return ;
        }
        //不直连，last、q之间的u、v可能也通过此节点跳下去；所以我们分出一个q做nq来接
        int nq=++toUsed;
        memcpy(Node+nq,Node+q,sizeof(node_t));
        Node[cur].pre=Node[q].pre=nq;
        Node[nq].step=Node[last].step+1;
        while(last&&Node[last].son[sn]==q){
            Node[last].son[sn]=nq;
            last=Node[last].pre;
        }
    }
}SAM;

int ans,len;
char Q[10010],s[10010];
void dfs(int x,int pos,int val){
    if(x>=len) return;
    if(val>=ans)return;

    int sn=Q[x]-'0';
    //不更改
    int nxt=SAM.Node[pos].son[sn];
    if(nxt==0){//结束，此时的dfs的询问串所取得前缀不为文本串的前缀
        ans=min(ans,val);
    }else dfs(x+1,nxt,val);

    //更改
    nxt=SAM.Node[pos].son[sn^1];
    if(nxt==0)
        ans=min(ans,val+1);
    else dfs(x+1,nxt,val+1);
}
int main(){
    int n,q;
    while(~scanf("%d%d",&n,&q)&&n&&q){
        scanf("%s",s);
        SAM.init();
        for(int i=0;i<n;++i)
            SAM.add(s[i],i+1);
        while(q--){
            int l,r;
            scanf("%d%d",&l,&r);
            len=r-l+1;
            l--;
            for(int i=0;i<len;++i)
                Q[i]=s[l+i];
            ans=0x3f3f3f3f;
            dfs(0,1,0);
            if(ans==0x3f3f3f3f)
                puts("Impossible");
            else printf("%d\n",ans);
        }
    }

    return 0;
}