当前位置:编程学习 > C/C++ >>

poj 2185 Milking Grid

求最小重复矩阵。

        首先得证明一个结论,最小重复矩阵一定是靠着左上角的。这个可以利用黏贴技术证明。如果有一个不靠左上角的重复矩阵,必然可以通过黏贴得到一个大小相同的靠左上角的矩阵,这个自己画画就明白了。

        第二个就是s=len-next[len]等于一个字符串的最小循环节。至于它是循环节利用画图可以看出来。

        至于是最小的可以这么考虑,如果存在更小的循环节,必然能得到更大的next值,与next数组矛盾。

        还有一点不是s的整数倍的必然不是循环节。

        但是此题麻烦的地方在于,可以黏贴在空白的地方。这个需要对多出来的部分再次求循环节。


[cpp]
#include <iostream>  
#include <cstdio>  
#include <cstring>  
using namespace std; 
const int maxn=1e4+9,maxm=80; 
char a[maxn][maxm],b[maxm][maxn]; 
int count[maxn],next[maxn]; 
int n,m; 
 
void kmp(char a[],int n) 

    next[0]=0; 
    next[1]=0; 
    for(int i=2;i<=n;i++) 
    { 
        int t=i-1; 
        while(t&&a[i]!=a[next[t]+1]) 
        t=next[t]; 
        next[i]=next[t]+(a[i]==a[next[t]+1]); 
    } 

 
void work(char a[],int n) 

 
    kmp(a,n); 
    int t=n-next[n]; 
    int ret=0,len=n; 
    while(ret+t!=n) 
    { 
        for(int i=t;i+ret<=n;i+=t) 
        count[i+ret]++; 
        ret+=len/t*t; 
        len=len-len/t*t; 
        t=len-next[len]; 
    } 
    count[n]++; 

 
void solve() 

    int c,r; 
    memset(count,0,sizeof(count)); 
    for(int i=1;i<=n;i++) 
    work(a[i],m); 
    for(int i=1;i<=m;i++) 
    if(count[i]>=n) 
    { 
        c=i; 
        break; 
    } 
    memset(count,0,sizeof(count)); 
    for(int i=1;i<=m;i++) 
    work(b[i],n); 
    for(int i=1;i<=n;i++) 
    if(count[i]>=m) 
    { 
        r=i; 
        break; 
    } 
 
    printf("%d\n",c*r); 

int main() 

//    freopen("in.txt","r",stdin);  
    while(scanf("%d %d",&n,&m)!=EOF) 
    { 
        for(int i=1;i<=n;i++) 
        scanf("%s",a[i]+1); 
 
        for(int i=1;i<=n;i++) 
        for(int j=1;j<=m;j++) 
        b[j][i]=a[i][j]; 
 
        solve(); 
    } 
    return 0; 

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=1e4+9,maxm=80;
char a[maxn][maxm],b[maxm][maxn];
int count[maxn],next[maxn];
int n,m;

void kmp(char a[],int n)
{
    next[0]=0;
    next[1]=0;
    for(int i=2;i<=n;i++)
    {
        int t=i-1;
        while(t&&a[i]!=a[next[t]+1])
        t=next[t];
        next[i]=next[t]+(a[i]==a[next[t]+1]);
    }
}

void work(char a[],int n)
{

    kmp(a,n);
    int t=n-next[n];
    int ret=0,len=n;
    while(ret+t!=n)
    {
        for(int i=t;i+ret<=n;i+=t)
        count[i+ret]++;
        ret+=len/t*t;
        len=len-len/t*t;
        t=len-next[len];
    }
    count[n]++;
}

void solve()
{
    int c,r;
    memset(count,0,sizeof(count));
    for(int i=1;i<=n;i++)
    work(a[i],m);
    for(int i=1;i<=m;i++)
    if(count[i]>=n)
    {
        c=i;
        break;
    }
    memset(count,0,sizeof(count));
    for(int i=1;i<=m;i++)
    work(b[i],n);
    for(int i=1;i<=n;i++)
    if(count[i]>=m)
    {
        r=i;
        break;
    }

    printf("%d\n",c*r);
}
int main()
{
//    freopen("in.txt","r",stdin);
    while(scanf("%d %d",&n,&m)!=EOF)
    {
        for(int i=1;i<=n;i++)
        scanf("%s",a[i]+1);

        for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        b[j][i]=a[i][j];

        solve();
    }
    return 0;
}


 

补充:软件开发 , C++ ,
CopyRight © 2022 站长资源库 编程知识问答 zzzyk.com All Rights Reserved
部分文章来自网络,