JZOJ5897. 【NOIP2018模拟10.6】密钥

题解

字母要配对，而且配对的方案很多，不可能枚举。
因为有用的信息就是相同字母的相对位置。
对于每一个位置，二分一个长度，
维护前缀的哈希值，然后通过加加减减可以得出每一段的。
按照字母出现的先后顺序枚举这一段的字母，暴力跟T中的字母建立映射，
判断是否合法。

code

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string.h>
#include <cmath>
#define N 100003
#define ll long long
#define P putchar
#define G getchar
using namespace std;
const int mo=233333339,w=3;
char ch,s[N],t[N];
int n,m,M,ans,p[30],fi[30][N],id;
ll sum[30][N],vt[30][N],z[N],ny[N],V;
int l,r,mid;
bool pd;
void write(int x){if(x>9) write(x/10);P(x%10+'0');}

ll ksm(ll x,int y)
{
	ll s=1;
	for(;y;y>>=1,x=x*x%mo)
		if(y&1)s=s*x%mo;
	return s;
}

ll calc(int ch,int x,int len)
{
	return (sum[ch][x+len-1]-sum[ch][x-1]*z[len]%mo+mo)%mo;
}

int main()
{
	freopen("key.in","r",stdin);
	freopen("key.out","w",stdout);
	
	scanf("%d%d",&n,&m);
	M=max(n,m),z[0]=1;
	for(int i=1;i<=M;i++)z[i]=z[i-1]*w%mo;
	for(ch=G();ch<'a' || ch>'z';ch=G());
	for(int i=1;i<=n;i++)s[i]=ch-'a'+1,ch=G();
	for(;ch<'a' || ch>'z';ch=G());
	for(int i=1;i<=m;i++)t[i]=ch-'a'+1,ch=G();
	
	for(int j=1;j<27;j++)fi[j][n+1]=j;
	for(int i=n;i;i--)
	{
		fi[id=1][i]=s[i];
		for(int j=1;j<27;j++)
			if(fi[j][i+1]!=s[i])fi[++id][i]=fi[j][i+1];
	}
	
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<27;j++)
			sum[j][i]=sum[j][i-1]*w%mo;
		sum[s[i]][i]++;
	}
	
	for(int i=1;i<=m;i++)
	{
		for(int j=1;j<27;j++)
			vt[j][i]=vt[j][i-1]*w%mo;
		vt[t[i]][i]++;
	}
	
	for(int i=1;i<=n;i++)
	{
		for(ans=l=1,r=min(m,n-i+1);l<=r;)
		{
			mid=(l+r)>>1;pd=1;
			memset(p,0,sizeof(p));
			for(int j=1;j<27 && pd;j++)
			{
				V=calc(fi[j][i],i,mid);
				if(p[fi[j][i]]==0)
				{
					for(int k=1;k<27;k++)
						if(vt[k][mid]==V && p[k]==0)
						{
							p[k]=fi[j][i],p[fi[j][i]]=k;
							break;
						}
				}
				if(p[fi[j][i]]==0 || vt[p[fi[j][i]]][mid]!=V)pd=0;
			}
			if(pd)ans=mid,l=mid+1;else r=mid-1;
		}
		write(ans),P('\n');
	}
	
	return 0;
}