hdoj 5918 Sequence I（kmp， 暴力？）

题目地址：点击打开链接

题意：

有两个序列，a1~an与b1~bm。 给出p，问在a1~an中有几个q能满足  b1==aq，  b2==bq+p，b3==bq+2p . . . bm==aq+(m-1)p

思路：

用了最朴素的暴力。。。从后往前找，找到和bm相同的开始匹配，看看是不是匹配上。。。这样应该是n^2复杂度吧。我想应该会TLE，结果过了。。

别人写的题解是用kmp ，表示没怎么看懂。不过用他们代码交一下发现居然都比暴力跑的慢。。。。

就是在原始的kmp匹配上待匹配的串做点小变比，本来某个字符匹配成功下标+1，现在改成下标+p即可。

kmp方法博客：点击打开链接

题解：典型的模式串匹配问题，和kmp模板不同的地方就是每隔p个距离查找一次， 而

然要查找出所有满足条件的模式串， 就要进行p次匹配。  例如第一次匹配为  1 ，

1+p， 1+2*p.....    第二次匹配为  2， 2+p， 2+2*p....  第p次匹配就是 p， 2*p， 3*p .....

暴力代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn = 1e6+5;
int a[maxn], b[maxn], p;
int main(void)
{
    int t, n, m, ca = 1;
    cin >> t;
    while(t--)
    {
        int ans = 0;
        scanf("%d%d%d", &n, &m, &p);
        for(int i = 1; i <= n; i++) scanf("%d", &a[i]);
        for(int i = 1; i <= m; i++) scanf("%d", &b[i]);
        for(int i = n; i >= 1; i--)
        {
            if(a[i] == b[m])
            {
                int j;
                for(j = 1; i-j*p >= 1 && m-j >= 1; j++)
                    if(b[m-j] != a[i-j*p]) break;
                if(m-j < 1) ans++;
            }
        }
        printf("Case #%d: %d\n", ca++, ans);
    }
    return 0;
}

KMP代码：

#include<cstdio>
#include<cstring>
using namespace std;
const int N=1e6+5;
int n,m,p,next[N],a[N],b[N];
void getnext()
{
    int i=0,j=-1;
    next[i]=j;
    while(i<m)
    {
        if(j==-1||b[i]==b[j])
        {
            i++;
            j++;
            next[i]=j;
        }
        else
            j=next[j];
    }
}
int kmp()
{
    int res=0;
    getnext();
    for(int k=0;k<p;++k){//进行p次匹配 
        int i=k,j=0;
        while(i<n)
        {
            if(j==-1||a[i]==b[j])
            {
                i+=p;//在母串上每隔p个位置匹配一次 
                j++;
                if(j==m)
                	res++;
            }
            else
                j=next[j];
        }
    }
    return res;
}
int main()
{
    int T;
    scanf("%d",&T);
    for(int kase=1;kase<=T;kase++)
    {
        scanf("%d%d%d",&n,&m,&p);
        memset(next,0,sizeof(next));
        memset(b,0,sizeof(b));//注意要清空b数组，这里WA了两发，迷之耽误一个小时 
        memset(a,0,sizeof(a));
        for(int i=0;i<n;i++)scanf("%d",&a[i]);
        for(int i=0;i<m;i++)scanf("%d",&b[i]);
        printf("Case #%d: %d\n",kase,kmp());
    }
    return 0;
}

Sequence I

Time Limit: 3000/1500 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 1428    Accepted Submission(s): 534

Problem Description

Mr. Frog has two sequences  a1,a2,⋯,an  and  b1,b2,⋯,bm  and a number p. He wants to know the number of positions q such that sequence  b1,b2,⋯,bm  is exactly the sequence  aq,aq+p,aq+2p,⋯,aq+(m−1)p  where  q+(m−1)p≤n  and  q≥1 .

 

Input

The first line contains only one integer  T≤100 , which indicates the number of test cases.

Each test case contains three lines.

The first line contains three space-separated integers  1≤n≤106,1≤m≤106  and  1≤p≤106 .

The second line contains n integers  a1,a2,⋯,an(1≤ai≤109) .

the third line contains m integers  b1,b2,⋯,bm(1≤bi≤109) .

 

Output

For each test case, output one line “Case #x: y”, where x is the case number (starting from 1) and y is the number of valid q’s.

 

Sample Input

  
  

   
   2
6 3 1
1 2 3 1 2 3
1 2 3
6 3 2
1 3 2 2 3 1
1 2 3
  
  

 

Sample Output

  
  

   
   Case #1: 2
Case #2: 1
  
  

 

Source

2016中国大学生程序设计竞赛（长春）-重现赛

 

Recommend

wange2014