poj 2185 kmp

最新推荐文章于 2021-02-09 19:24:31 发布

Eirlys_North

最新推荐文章于 2021-02-09 19:24:31 发布

阅读量256

点赞数

分类专栏： kmp&extend-kmp 文章标签： kmp

本文链接：https://blog.csdn.net/Eirlys_North/article/details/56011135

版权

kmp&extend-kmp 专栏收录该内容

5 篇文章 0 订阅

订阅专栏

题意：给你一个二维的字符矩阵，找到一个最小的子矩阵使得多次复制以后可以得到原矩阵（不要求整除）

一开始我naive的觉得最小子矩阵的宽x就是每一行的最小循环节取个最大值，然后再把每行前x个字符组成的字符串看做一个整体，跑kmp出最小的矩阵长y

然后随便来一个串就会知道，如果前l位是最小循环节，但是前（l+x）位不一定是一个循环节，所以每行最小循环节取最大值是naive的错误

后来感觉取最大不行，那就取每行最小循环节的最小公倍数，后来被讨论区的大神卡掉了 ...又是一个naive的想法，大神题解 orz

然后就去想，这个最小子矩阵必须满足的条件是什么

它的宽x，必须同时满足每一行的字符串的前x位是该字符串的循环节之一且最小

它的长y，必须同时满足每一列的字符串的前y位是该字符串的循环节之一且最小

那么对于每一行，我们可以暴力找出改行字符串的所有循环节情况，然后就找到每一行的公共最小循环节就是x

那么对于前x列，我们同样可以像对行处理那样，但是并没有这个必要，因为x是已经确定了的且所有的要求是同时满足，所以完全可以整体处理

那么我们就可以把每一行的前x位看做一个整体，跑kmp出最小循环节即可

var
        n,m,x,y,p       :longint;
        i,j,k           :longint;
        f               :boolean;
        s,map           :array[0..10010] of string;
        flag            :array[0..80] of longint;
        next            :array[0..10010] of longint;


begin
   readln(n,m);
   for i:=1 to n do
   begin
      readln(map[i]);
      for j:=1 to m do
      begin
         p:=1;
         k:=j+1;
         f:=true;
         while (k<=m) do
         begin
            if map[i][k]<>map[i][p] then
            begin
                f:=false;break;
            end;
            p:=p mod j+1;
            inc(k);
         end;
         if f then inc(flag[j]);
      end;
   end;
   //
   for i:=1 to m do
     if flag[i]=n then break;
   x:=i;
   p:=0;
   for i:=1 to n do s[i]:=copy(map[i],1,x);
   for i:=2 to n do
   begin
      while (p<>0) and (s[i]<>s[p+1]) do p:=next[p];
      if (s[i]=s[p+1]) then
      begin
         inc(p);
         next[i]:=p;
      end else next[i]:=0;
   end;
   y:=n-next[n];
   writeln(x*y);
end.

——by Eirlys