题意:给出两个数列(可以看成串),问最多能在串1中找出几个不重叠的substring,使得这些substring的元素的相对大小关系与串2相同(显然,长度也要与串2相同)。
解题思路:kmp。。比赛的时候有想到这个思路,但一直不知道怎么搞定这个相对大小关系。其实很简单。。
count(Ai<An+1)=count(Bi<Bn+1),i=1..n
count(Ai=An+1)=count(Bi=Bn+1),i=1..n
count(Ai>An+1)=count(Bi>Bn+1),i=1..n
这样不就行了吗?那怎么弄这个count呢?暴力记录就好了,开两个二维数组g1[i][j]表示串1中,到j为止,有多少个i。g2[i][j]表示串2中,到j为止,有多少个i。其中i的范围1-25完全可以接受。然后就是kmp的匹配过程了,这个匹配并不是绝对的匹配,而是说,我们的相对大小关系相同。比如,我们要看串1的x位置,能否和串2的y位置匹配,那我们就看串1里面,[x-y+1,x-1]这段区间内(注意这个区间的长度是y-1,我们要匹配的是x,y,所以前y-1个长度肯定是匹配好的),小于a[x]的数的个数,以及等于a[x]的数的个数,是否相应的与串2里面,[1,y-1]这段区间内,小于b[y]的数的个数,等于b[y]的数的个数相等。这个就用刚才记录的数组暴力算好了(当然也可以用树状数组,但是字符集只有25,完全可以暴力),其他就跟普通的kmp一样了。在匹配的过程中,记录下匹配点,然后去枚举,一旦找到匹配点,就ans+1,然后往后跳m长度去找匹配点,这样贪心。
(我的代码里,把串1,串2对调了,所以g1,g2也对调了)。
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std ;
const int maxn = 111111 ;
int g1[26][maxn] , g2[26][maxn] ;
int p[maxn] , a[maxn] , b[maxn] ;
int vis[maxn] ;
int n , m , t ;
int ok ( int x , int y ) {
int c1 = 0 , i , c2 = 0 , c3 = 0 ;
if ( g1[a[x]][x-1] != g1[a[y]][y-1] - g1[a[y]][y-x] ) return 0 ;
for ( i = 1 ; i < a[x] ; i ++ ) c1 += g1[i][x-1] ;
for ( i = 1 ; i < a[y] ; i ++ ) c2 += g1[i][y-1] ;
for ( i = 1 ; i < a[y] ; i ++ ) c3 += g1[i][y-x] ;
if ( c1 != c2 - c3 ) return 0 ;
return 1 ;
}
void get_p () {
int i , j = 0 ;
p[1] = 0 ;
for ( i = 2 ; i <= n ; i ++ ) {
while ( j && !ok ( j + 1 , i ) ) j = p[j] ;
if ( ok ( j + 1 , i ) ) j ++ ;
p[i] = j ;
}
}
int judge ( int x , int y ) {
int c1 = 0 , i , c2 = 0 , c3 = 0 ;
if ( g1[a[x]][x-1] != g2[b[y]][y-1] - g2[b[y]][y-x] ) return 0 ;
for ( i = 1 ; i < a[x] ; i ++ ) c1 += g1[i][x-1] ;
for ( i = 1 ; i < b[y] ; i ++ ) c2 += g2[i][y-1] ;
for ( i = 1 ; i < b[y] ; i ++ ) c3 += g2[i][y-x] ;
if ( c1 != c2 - c3 ) return 0 ;
return 1 ;
}
void kmp () {
int i , j = 0 ;
for ( i = 1 ; i <= m ; i ++ ) {
while ( j && !judge ( j + 1 , i ) ) j = p[j] ;
if ( judge ( j + 1 , i ) ) j ++ ;
if ( j == n ) {
vis[i] = 1 ;
j = p[j] ;
}
}
// for ( i = 1 ; i <= m ; i ++ ) printf ( "%d " , vis[i] ) ; puts ( "" ) ;
}
int main () {
int i , j , k ;
while ( scanf ( "%d%d%d" , &m , &n , &t ) != EOF ) {
memset ( g1 , 0 , sizeof ( g1 ) ) ;
memset ( g2 , 0 , sizeof ( g2 ) ) ;
memset ( vis , 0 , sizeof ( vis ) ) ;
for ( i = 1 ; i <= m ; i ++ ) {
scanf ( "%d" , &b[i] ) ;
for ( j = 1 ; j <= t ; j ++ ) g2[j][i] = g2[j][i-1] ;
g2[b[i]][i] ++ ;
}
for ( i = 1 ; i <= n ; i ++ ) {
scanf ( "%d" , &a[i] ) ;
for ( j = 1 ; j <= t ; j ++ ) g1[j][i] = g1[j][i-1] ;
g1[a[i]][i] ++ ;
}
get_p () ;
kmp () ;
int ans = 0 ;
i = 1 ;
while ( i <= m ) {
// printf ( "i = %d , ans = %d\n" , i , ans ) ;
if ( vis[i] )
ans ++ , i += n ;
else i ++ ;
}
printf ( "%d\n" , ans ) ;
}
}
count(Ai=An+1)=count(Bi=Bn+1),i=1..n
count(Ai>An+1)=count(Bi>Bn+1),i=1..n