Bzoj 3409: [Usaco2009 Oct]Barn Echoes 牛棚回声

SAM练手题

反正就是构造两个SAM然后匹配就好了。。。。

#include<cstdio  >居然会错..

#include<cstdio>才行。。。

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
struct Node
{
	Node *ch[200];
	Node *last;
	int len;
	bool end;
	Node(){for(len=0;len<=199;len++)ch[len]=NULL;len=0;end=false;last=NULL;}
};
Node *root=new Node;
Node *last=root;

inline void insert(int data)
{
	Node *now=new Node,*tp;
	now->len=last->len+1;
	for(tp=last;tp&&!tp->ch[data];tp=tp->last)tp->ch[data]=now;
	if(!tp)last=now,now->last=root;
	  else if(tp->len==tp->ch[data]->len-1)now->last=tp->ch[data];
	  else
	   {
	   	  Node *a=tp->ch[data],*b=new Node;
	   	  *b=*a;
	   	  b->len=a->len+1;
	   	  a->last=now->last=b;
	   	  for(tp;tp&&tp->ch[data]==a;tp=tp->last)tp->ch[data]=b;
	   }
   last=now;
}
inline void up(Node *a)
{
	for(a;a;a=a->last)a->end=true;
}
char a[10001],b[100001];
int main()
{
	int  i;
	scanf("%s",a);
	scanf("%s",b);
	int len1=strlen(a),len2=strlen(b);
	for(i=1;i<=len1;i++) 
	  insert(a[i-1]);
	up(last);
	int ans=0,data;
	Node *ll,*now=root,*tp;
	int nowlen=0;
	for(i=1;i<=len2;i++)
	  {
	  	data=b[i-1];
          for(tp=now;tp&&!tp->ch[data];tp=tp->last);
	      if(tp==NULL){break;} 
		  if(tp!=now)
		  nowlen=tp->len+1;
		  else nowlen++;
	      now=tp->ch[data];
	      if(now->end)
	         ans=max(ans,nowlen);
	  }
	  
	  
	  
	root=last=new Node;
	for(i=1;i<=len2;i++) 
	  insert(b[i-1]);
	up(last);
	now=root;
	nowlen=0;
	for(i=1;i<=len1;i++)
	  {
	  	data=a[i-1];
          for(tp=now;tp&&!tp->ch[data];tp=tp->last);
	      if(tp==NULL){break;} 
		  if(tp!=now)
		  nowlen=tp->len+1;
		  else nowlen++;
	      now=tp->ch[data];
	      if(now->end)
	         ans=max(ans,nowlen);
	  }
	printf("%d\n",ans);
	return 0;
}


你可能感兴趣的:(后缀自动机,sam,bzoj)