【HDU5469】Antonidas（点分治，字符串哈希）
题面
题解
啊哈？什么垃圾一眼点分治+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;
}