3167

/*

解决方法好多:
1. KMP
2. Hash


这个题目不能直接的套用KMP模板来解决
计算实时排名非常关键
首先求next数组时,不能在像往常的来判断两个字符相等了,而是要判断其动态排名是相等的。
求动态的排名是这道题目的关键。
不过求解起来确实是有点的繁琐的.
*/


// include file
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cctype>
#include <ctime>

#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <bitset>
#include <strstream>

#include <algorithm>
#include <string>
#include <vector>
#include <queue>
#include <set>
#include <list>
#include <functional>

using namespace std;

// typedef
typedef long long LL;
typedef unsigned long long ULL;

// 
#define read freopen("in.txt","r",stdin)
#define write freopen("out.txt","w",stdout)
#define FORi(a,b,c) for(int i=(a);i<(b);i+=c)
#define FORj(a,b,c) for(int j=(a);j<(b);j+=c)
#define FORk(a,b,c) for(int k=(a);k<(b);k+=c)
#define FORp(a,b,c) for(int p=(a);p<(b);p+=c)

#define FF(i,a)    for(int i=0;i<(a);i+++)
#define FFD(i,a)   for(int i=(a)-1;i>=0;i--)
#define Z(a) (a<<1)
#define Y(a) (a>>1)

const double eps = 1e-6;
const double INFf = 1e10;
const int INFi = 1000000000;
const double Pi = acos(-1.0);

template<class T> inline T sqr(T a){return a*a;}
template<class T> inline T TMAX(T x,T y)
{
	if(x>y) return x;
	return y;
}
template<class T> inline T TMIN(T x,T y)
{
	if(x<y) return x;
	return y;
}
template<class T> inline void SWAP(T &x,T &y)
{
	T t = x;
	x = y;
	y = t;
}
template<class T> inline T MMAX(T x,T y,T z)
{
	return TMAX(TMAX(x,y),z);
}


// code begin
#define MAXT 100010
#define MAXP 25110
int N,S,K;

int T[MAXT];
int P[MAXP];

int cT[MAXT][26];
int cP[MAXP][26];
int next[MAXP];
int rk[26];
int ans[MAXT],top;

void GetRankP(int l,int pos,int &rank1,int &rank2)
{
	//区间是从l到pos这个区间里pos的排名
	FORi(1,S+1,1)
	{
		rk[i] = cP[pos][i];
	}
	if(l!=0)
	{
		FORi(1,S+1,1)
		{
			rk[i] -= cP[l-1][i];
		}
	}

	// 判断pos位置的值 P[pos]所处的位置
	rank1 = 0;
	FORi(1,P[pos],1)
		rank1+=rk[i];
	rank2 = rk[P[pos]];
}

void GetRankT(int l,int pos,int &rank1,int &rank2)
{
	//区间是从l到pos这个区间里pos的排名
	FORi(1,S+1,1)
	{
		rk[i] = cT[pos][i];
	}
	if(l!=0)
	{
		FORi(1,S+1,1)
		{
			rk[i] -= cT[l-1][i];
		}
	}

	// 判断pos位置的值 P[pos]所处的位置
	rank1 = 0;
	FORi(1,T[pos],1)
		rank1+=rk[i];
	rank2 = rk[T[pos]];
}

void CalNext()
{
	next[0] = 0;
	next[1] = 0;
	int j;
	int ranka1,ranka2,rankb1,rankb2;
	FORi(2,K+1,1)
	{
		j = next[i-1]; // 此处j是指i-1时能否返回的位置,

		while(true)
		{
			// 计算 i-1的排名和j的排名
			GetRankP(0,j,ranka1,ranka2); 
			GetRankP(i-1-j,i-1,rankb1,rankb2);
			//if(P[i-1]==P[j]) //如果i-1时返回的位置和i-1的字符相等的话,那么就应该+1. 但是对于本题,这里不在是指字符相等了,而是指相应的排名相等
			if( ranka1==rankb1 && ranka2==rankb2 )
			{
				next[i] = j+1;
				break;
			}
			if(j==0)
			{
				next[i] = 0;
				break;
			}
			j = next[j];
		}
	}
}

void KMP()
{
	int i = 0; //文本
	int j = 0; //模式 
	int ranka1,ranka2,rankb1,rankb2;
	while(true)
	{
		if(i==N) break;
		GetRankT(i-j,i,ranka1,ranka2);
		GetRankP(0,j,rankb1,rankb2);
		//if( T[i]==P[j]) // 这里不再判断字符是否相等了,而是判断排名
		if(ranka1==rankb1 && ranka2==rankb2)
		{
			i++;
			j++;
			if(j==K)
			{
				ans[top++] = i-K;
				j = next[j];
			}
		}
		else if(j>0)
		{
			j = next[j];
		}
		else 
			i++;
	}
}

int main()
{
	read;
	write;
	while(scanf("%d %d %d",&N,&K,&S)!=-1)
	{
		FORi(0,N,1)
		{
			scanf("%d",T+i);
		}
		FORi(0,K,1)
		{
			scanf("%d",P+i);
		}
		if(K==1)
		{
			printf("%d\n",N);
			FORi(0,N,1)
			{
				printf("%d\n",T[i]);
			}
		}

		memset(cT,0,sizeof(cT));
		memset(cP,0,sizeof(cP));
		cT[0][T[0]] = 1;
		FORi(1,N,1)
		{
			FORj(1,S+1,1)
			{
				cT[i][j] = cT[i-1][j];
			}
			cT[i][T[i]]++;
		}
		cP[0][P[0]] = 1;
		FORi(1,K,1)
		{
			FORj(1,S+1,1)
			{
				cP[i][j] = cP[i-1][j];
			}
			cP[i][P[i]]++;
		}

		// 计算失败函数
		CalNext();
		top = 0;
		KMP();

		printf("%d\n",top);
		FORi(0,top,1)
		{
			printf("%d\n",ans[i]+1);
		}
	}
	return 0;
}

posted @ 2011-03-25 23:06  AC2012  阅读(131)  评论(0编辑  收藏  举报