【GDOI 2016 Day1】第二题 最长公共子串

Description

抽象题意：给出两字符串S,T，求经过一些奇怪的操作后的最长公共子串

Solution

很显然一些有交集的区间是可以合并的，
一个区间内我们只需要记录每个字母的个数，
我们就可以枚举两个字符串的相对位置，然后在重叠的部分用指针维护，
设区间i,j是合法的，
我们现在把j往后移一位，j所在的区间的字母$T_j$个数就-1，
然后把再移动指针i去使其合法，i所在的区间的$T_i$个数就+1，并且i++；直到合法或移爆为止。

Code

#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int N=2500;
int read(int &n)
{
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
bool ZM(char ch){return('a'<=ch && ch<='z');}
int n,m,m1,ans;
int a[N],b[N];
int ag[N];
struct qqw{int s,t;}c[N*5];
int f[N][28],F[N][28],f0;
bool PX(qqw a,qqw b){return(a.s<b.s || a.s==b.s&&a.t<b.t);}
void solue(int q,int w)
{
    int i=q,j,I=w-q;
    for(j=q;j<=n && j+I<=m;j++)
    {
        while(!f[ag[j]][b[j+I]] && i<=j)
        {
            if(i<j)f[ag[i]][b[i+I]]++;i++;
        }
        if(i<=j)f[ag[j]][b[j+I]]--;
        ans=max(ans,j-i+1);
    }
    while(i<j)f[ag[i]][b[i+I]]++,i++;
}
int main()
{
    freopen("lcs.in","r",stdin);
    freopen("lcs.out","w",stdout);
    int q,w;char ch;
    while(!ZM(ch=getchar()));
    for(;ZM(ch);ch=getchar())b[++m]=ch-96;b[m+1]=27;
    while(!ZM(ch=getchar()));
    for(;ZM(ch);ch=getchar())a[++n]=ch-96;a[n+1]=0;
    read(m1);
    fo(i,1,m1)read(c[i].s),read(c[i].t);
    c[++m1].s=c[m1].t=n-1;
    sort(c+1,c+1+m1,PX);
    q=w=0;
    fo(i,1,m1)
    {
        c[i].s++,c[i].t++;
        if(q<c[i].s)while(q+1<=c[i].s)ag[++q]=++w;
        fo(j,q+1,c[i].t)ag[j]=w;
        q=max(q,c[i].t);
    }
    f0=w;ans=0;
    fo(i,1,n)F[ag[i]][a[i]]++;
    fo(i,1,f0)fo(j,1,26)f[i][j]=F[i][j];
    fo(i,2,n)solue(i,1);
    fo(i,1,m)solue(1,i);
    printf("%d\n",ans);
    return 0;
}