[gmoj] [GDOI 2016 Day1] 最长公共子串

[gmoj] [GDOI 2016 Day1] 最长公共子串

题目

---

解题思路

题意
给出k个区间，区间内可任意换位置
最大的公共子串最长为多少

首先不可能枚举出每种情况再找，太慢了
可以想到记录每个区间内每个字母含有多少个
匹配时用一个减去一个
区间内可以随意换，不需要给出具体换成什么样

给出了k个区间
如果区间之间有交点那么就可以合并
因为你可以将它先换到公共区域，再将它放到另一个区间
再将字符串分为一段一段
做一遍初始匹配
l[i][j]表示第i段从匹配串第j位开始，左边有多少位能匹配
r[i][j]表示第i段从匹配串第j位开始，右边有多少位能匹配
如果往右可以拿完整个区间，可以继续往右拿
用记搜，避免重复计算

---

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
struct lzf{
	int l,r;
}a[100010],b[10100],nw[2010];
int lz[2200][2010],rz[2200][2010],f[2200][2010],p[2020][31];
int lenx,leny,t=1,tot,n,ans;
char x[2020],y[2020];
bool cmp(lzf f,lzf t)
{
	 if (f.l==t.l) return f.r<t.r;
	    else return f.l<t.l;
}
int work(int i,int j)
{
	if (i>tot||j>lenx) return 0;
	if (f[i][j]) return f[i][j];
	if (rz[i][j]!=nw[i].r-nw[i].l+1) 
        return f[i][j]=rz[i][j];
	   else return f[i][j]=rz[i][j]+work(i+1,j+rz[i][j]);
}
int main()
{
	freopen("lcs.in","r",stdin);
	freopen("lcs.out","w",stdout);
	scanf("%s",x+1);
	lenx=strlen(x+1);
	scanf("%s",y+1);
	leny=strlen(y+1);  //从第一位开始存储
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
	{
	    scanf("%d%d",&a[i].l,&a[i].r);
	    a[i].l++;
	    a[i].r++;  
	}
	sort(a+1,a+n+1,cmp);  //按左边界排序
	b[t].l=a[1].l,b[t].r=a[1].r;
	for (int i=2;i<=n;i++)  //合并区间
	{
	    if (a[i].l<=b[t].r)
	       b[t].r=max(b[t].r,a[i].r);
	       else {
	       	    b[++t].l=a[i].l;
	       	    b[t].r=a[i].r;
		   }
	}
	if (b[1].l>1)  //将字符串分为一段一段，并记录所含有字母的个数
	   for (int i=1;i<b[1].l;i++) 
	       p[++tot][y[i]-'a']++,nw[tot]=(lzf){i,i};
	b[t+1].l=leny+1; 
	for (int i=1;i<=t;i++)
	{
		tot++;
		nw[tot]=(lzf){b[i].l,b[i].r};
		for (int j=b[i].l;j<=b[i].r;j++)
		    p[tot][y[j]-'a']++;
		for (int j=b[i].r+1;j<=b[i+1].l-1;j++)
		    p[++tot][y[j]-'a']++,nw[tot]=(lzf){j,j};
	}  
	for (int i=1;i<=tot;i++)  //初始匹配，一段一段的和匹配串匹配
	    for (int j=1;j<=lenx;j++)
	    {
	    	int w=0;
	    	while (j+w<=lenx&&p[i][x[j+w]-'a'])
	    	      p[i][x[j+w]-'a']--,w++;
	    	rz[i][j]=w; 
	    	for (int c=0;c<w;c++)
	    	    p[i][x[j+c]-'a']++;
	    	w=0;
	    	while (j-w>0&&p[i][x[j-w]-'a'])
	    	      p[i][x[j-w]-'a']--,w++;
	    	lz[i][j]=w;
	    	for (int c=0;c<w;c++)
	    	    p[i][x[j-c]-'a']++;
		 } 
    for (int i=1;i<=tot+1;i++)
        for (int j=1;j<=lenx+1;j++)
            ans=max(ans,lz[i-1][j-1]+work(i,j));  //左边也可来一点，尝试向右边扩展
    printf("%d",ans);
    fclose(stdin);
    fclose(stdout);
    return 0;
}