【HDU5469】Antonidas（点分治，字符串哈希）

【HDU5469】Antonidas（点分治，字符串哈希）

题面

HDU
Vjudge

题解

啊哈？什么垃圾一眼点分治+Hash判断，哈哈哈哈哈，让我来码码码。
诶，怎么WA了。改改改改改。
诶，怎么很对啊，去网上蒯一个标程来拍拍拍。
诶，怎么拍不WA啊，让我来人工检测一波啊哈哈哈哈。
感觉没有问题啊？诶，我换个方式来拍把，把字符集设小点就好了。
诶，WA不了，诶诶诶诶，WA了。
WOC，数据1w让我怎么手玩。
不管了，先把剪枝删掉再试试，诶诶诶，怎么答案对了？
我之前交不加剪枝的不是TLE了吗？
诶诶诶，WOC，怎么过了，哈哈哈哈哈，果然是傻逼点分治+Hash题。
诶，怎么，，，，怎么我的晚自习就没了呢？我怎么哭了？诶诶诶诶？

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<set>
using namespace std;
#define MAX 10010
#define ull unsigned long long
const int base=19260817;
inline int read()
{
	int x=0;bool t=false;char ch=getchar();
	while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
	if(ch=='-')t=true,ch=getchar();
	while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
	return t?-x:x;
}
struct Line{int v,next;}e[MAX<<1];
int h[MAX],cnt=1;
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
int n;
int size[MAX],Size,mx,root;
bool vis[MAX];
void getroot(int u,int ff)
{
	size[u]=1;int ret=0;
	for(int i=h[u];i;i=e[i].next)
	{
		int v=e[i].v;if(v==ff||vis[v])continue;
		getroot(v,u);size[u]+=size[v];
		ret=max(ret,size[v]);
	}
	ret=max(ret,Size-size[u]);
	if(ret<mx)mx=ret,root=u;
}
int len;
char ch[MAX],s[MAX];
ull hs[MAX],pw[MAX];
set<int> S1,S2;
bool ans;
ull get(int l,int r){return hs[r]-hs[l-1]*pw[r-l+1];}
void calc(int u,int ff,int dep,ull h1,ull h2)
{
	if(ans)return;if(dep>=len)return;
	h1=h1*base+ch[u];h2=h2+pw[dep]*ch[u];
	if(dep+1==len&&h1==get(1,len)){ans=true;return;}
	if(h1==get(len-dep,len)&&S2.count(len-dep-1)){ans=true;return;}
	if(h2==get(1,dep+1)&&S1.count(dep+2)){ans=true;return;}
	for(int i=h[u];i;i=e[i].next)
		if(!vis[e[i].v]&&e[i].v!=ff)
			calc(e[i].v,u,dep+1,h1,h2);
}
void insert(int u,int ff,int dep,ull h1,ull h2)
{
	if(ans)return;if(dep>=len)return;
	h1=h1*base+ch[u];h2=h2+pw[dep]*ch[u];
	if(h1==get(len-dep,len))S1.insert(len-dep);
	if(h2==get(1,dep+1))S2.insert(dep+1);
	for(int i=h[u];i;i=e[i].next)
		if(!vis[e[i].v]&&e[i].v!=ff)
			insert(e[i].v,u,dep+1,h1,h2);
}
void Divide(int u)
{
	if(ans)return;vis[u]=true;
	S1.clear();S2.clear();
	if(ch[u]==s[1])S2.insert(1);
	if(ch[u]==s[len])S1.insert(len);
	if(len==1&&ch[u]==s[1]){ans=true;return;}
	for(int i=h[u];i;i=e[i].next)
	{
		int v=e[i].v;if(vis[v])continue;
		calc(v,u,0,0,0);insert(v,u,1,ch[u],ch[u]);
	}
	for(int i=h[u];i;i=e[i].next)
	{
		int v=e[i].v;if(vis[v])continue;
		Size=mx=size[v];getroot(v,u);
		Divide(root);
	}
}
int main()
{
	pw[0]=1;for(int i=1;i<MAX;++i)pw[i]=pw[i-1]*base;
	int T=read();
	for(int TT=1;TT<=T;++TT)
	{
		n=read();cnt=1;memset(h,0,sizeof(h));
		memset(vis,0,sizeof(vis));
		for(int i=1;i<n;++i)
		{
			int u=read(),v=read();
			Add(u,v);Add(v,u);
		}
		scanf("%s",ch+1);scanf("%s",s+1);len=strlen(s+1);
		memset(hs,0,sizeof(hs));
		for(int i=1;i<=len;++i)hs[i]=hs[i-1]*base+s[i];
		ans=false;Size=mx=n;getroot(1,0);
		Divide(root);
		printf("Case #%d: ",TT);
		ans?puts("Find"):puts("Impossible");
	}
	return 0;
}