poj 2185

本文介绍了一种利用KMP算法高效求解最小覆盖重复子矩阵的方法。通过计算每行每列的最小重复串长度,利用这些长度的最小公倍数来确定最小重复子矩阵的宽度和高度。代码实现了核心算法,包括获取最小重复长度及计算最小公倍数。
思路:KMP很好的一道题。首先易证:最小覆盖子矩阵一定靠左上角。那么,我们考虑求出每一行的最小重复串长度,所有行的最小重复串的长度的lcm就是最小重复子矩阵的宽。然后我们对列也做相同的操作。于是我们就可以求得最小重复子矩阵的大小了。(这里要注意一点:当所得的宽大于原来的宽时,就让等于原来的宽,长也如此)。算法实现:算法的核心在于高效的求出每一行和每一列的最小重复串,这个可以最原串做一次KMP中的get_next()。(注:以上部分为转载)最小覆盖就等i-next[i];
代码:
#include<stdio.h>
#define N 10005
int r,l,next[N];
char s[N][N];
int getr(int n)
{
	int i,j;
	i=0;j=-1;
	next[0]=-1;
	while(i<l)
	{
		if(j==-1||s[n][i]==s[n][j])
		{
			i++;j++;
			next[i]=j;
		}
		else
			j=next[j];
	}
	return i-next[i];
}
int getl(int n)
{
	int i,j;
	i=0;j=-1;
	next[0]=-1;
	while(i<r)
	{
		if(j==-1||s[i][n]==s[j][n])
		{
			i++;j++;
			next[i]=j;
		}
		else
			j=next[j];
	}
	return i-next[i];
}
int gcd(int a,int b)
{
	return b==0 ? a : gcd(b,a%b);
}
int lcm(int a,int b)
{
	return a/gcd(a,b)*b;
}
int main()
{
	int i,j,ans1,ans2;
	while(scanf("%d%d",&r,&l)!=EOF)
	{
		for(i=0;i<r;i++)
			scanf("%s",s[i]);
		ans1=1;
		for(i=0;i<r;i++)
			ans1=lcm(ans1,getr(i));
		if(ans1>l)//注意判断
			ans1=l;
		ans2=1;
		for(j=0;j<l;j++)
			ans2=lcm(ans2,getl(j));
		if(ans2>r)
			ans2=r;
		printf("%d\n",ans1*ans2);
	}
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值