HDU 2243

原创 2012年03月30日 18:47:29

这题WA了.. 实在找不出错,CP了很多AC代码,大数据都是一样的。

哪位大神行行好帮我看看哪错了吧.....

#include<iostream>
#include<cstdio>
#include<string.h>
#define kind 26
#define MAX 100
using namespace std;

struct node{
 	   int fail,flag;
 	   int next[kind];
 	   void init(){
		   memset(next,0,sizeof(next));
		   fail=flag=0;
	   }
}Tire[MAX];

typedef unsigned __int64 ull; 
int root,allocp;
int N,L;
ull matrix[MAX][MAX];

void init(){
 	 root=allocp=0;
 	 Tire[root].init();
 	 memset( matrix,0,sizeof(matrix) );
}

void Insert( char *s )
{
 	 int i=0,k;
 	 int p=root;
 	 while( s[i] ){
	  		k=s[i++]-'a';
	  		if( !Tire[p].next[k] ){
			 	Tire[++allocp].init();
			 	Tire[p].next[k]=allocp;
		 	}else if( Tire[Tire[p].next[k]].flag )
		 		  break;
		 	p=Tire[p].next[k];
	 }
	 Tire[p].flag++;
}

void build_ac_automation()
{
 	 int queue[MAX],head,tail;
 	 head=tail=0;
 	 queue[tail++]=root;
 	 while( head!=tail ){
	  		int cur=queue[head++];
	  		for( int i=0;i<kind;i++ ){
			 	 if( Tire[cur].next[i] ){
				  	 int son=Tire[cur].next[i];
				  	 int p=Tire[cur].fail;
				  	 if( cur==root ) Tire[son].fail=root;
				  	 else Tire[son].fail=Tire[p].next[i];
				  	 if( Tire[Tire[son].fail].flag )
				  	 	 Tire[son].flag++;
				  	 queue[tail++]=son;
			   	 }
			   	 else{
				  	 int p=Tire[cur].fail;
				  	 if( cur==root ) Tire[cur].next[i]=root;
				  	 else Tire[cur].next[i]=Tire[p].next[i];
			   	 }
 	 	 	}
	 }
}

void set_matrix()
{
 	 allocp++;
 	 for( int i=0;i<allocp;i++ ){
	  	  if( Tire[i].flag )
	  	  	  continue;
	  	  for( int j=0;j<kind;j++ ){
		   	   int son=Tire[i].next[j];
		   	   if( Tire[son].flag==0 )
		   	   	   matrix[i][son]++;
  	   	  }
	 }
	 for( int i=0;i<allocp;i++ )
	 	  matrix[i][allocp+i]=matrix[allocp+i][allocp+i]=1;
}

void matriXmult( ull b[][MAX],ull c[][MAX],int r )
{
 	 ull a[MAX][MAX];
 	 memset( a,0,sizeof(a) );
 	 for( int i=0;i<r;i++ )
 	 for( int j=0;j<r;j++ )
 	 for( int k=0;k<r;k++ )
   	 	  a[i][j]+=b[i][k]*c[k][j];
 	 for( int i=0;i<r;i++ )
 	 for( int j=0;j<r;j++ )
 	 	  b[i][j]=a[i][j];
}

ull getPower(){
 	ull res[MAX][MAX],temp[MAX][MAX];
 	memset( res,0,sizeof(res) );
 	memset( temp,0,sizeof(temp) );
 	res[0][0]=res[1][1]=1;
	temp[0][0]=kind;
	temp[0][1]=temp[1][1]=1;
	temp[1][0]=0;
 	for( int i=0;i<31;i++ ){
		 if( L&(1<<i) )
		  	 matriXmult(res,temp,2);
	 	 matriXmult(temp,temp,2);
	}
	return res[0][0]+res[0][1]-1;
}

ull matrix_Power( ull matrix[][MAX] )
{
 	 ull res[MAX][MAX],temp[MAX][MAX];
 	 memset( res,0,sizeof(res) );
 	 memset( temp,0,sizeof(temp) );
 	 for( int i=0;i<2*allocp;i++ )
 	 {
 	 	  res[i][i]=1;
 	 	  for( int j=0;j<2*allocp;j++ )
 	 	  	   temp[i][j]=matrix[i][j];
     }
     for( int i=0;i<31;i++ ){
	  	  if( L&(1<<i) )
	  	  	  matriXmult( res,temp,2*allocp );
	  	  matriXmult( temp,temp,2*allocp );
	 }
	 ull ret=0;
	 for( int i=0;i<2*allocp;i++ )
	 	  ret+=res[0][i];
	 return ret-1;
}

int main()
{
 	//freopen( "my.out","w",stdout );
 	//freopen( "test.in","r",stdin );
 	char str[11];
 	while( scanf("%d%d",&N,&L)!=EOF )
 	{
	 	   init();
	 	   while( N-- ){
		   		  scanf( "%s",&str );
		   		  Insert(str);
  		   }
  		   build_ac_automation();
  		   set_matrix();
  		   ull total=getPower();
  		   ull sum=matrix_Power(matrix);
  		   printf( "%I64u\n",total-sum );
  	}
 	return 0;
}


版权声明:本文为博主原创文章,转载请写明出处。

相关文章推荐

hdu 2243 考研路茫茫——单词情结

ac自动机 + 快速矩阵连成求和#include #include #include #include #include #include using namespace st...

HDU 2243 考研路茫茫――单词情结 (AC自动机 + dp)

HDU 2243 考研路茫茫――单词情结 题意:给定一些词根,如果一个单词包含有词根,则认为是有效的。现在问长度不超过L的单词里面,有多少有效的单词? 思路:这道题和POJ 2778 是同...
  • SIOFive
  • SIOFive
  • 2014年08月20日 16:10
  • 959

AC自动机专题——P - 考研路茫茫――单词情结 HDU - 2243 矩阵快速幂+AC自动机

背单词,始终是复习英语的重要环节。在荒废了3年大学生涯后,Lele也终于要开始背单词了。  一天,Lele在某本单词书上看到了一个根据词根来背单词的方法。比如"ab",放在单词前一般表示"相反,变坏...

hdu 2243 长度为1~N的串中包含了模式串的串总共有几个

Problem Description背单词,始终是复习英语的重要环节。在荒废了3年大学生涯后,Lele也终于要开始背单词了。一天,Lele在某本单词书上看到了一个根据词根来背单词的方法。比如"ab"...

HDU 2243 考研路茫茫――单词情结 AC自动机 + 矩阵快速幂

题目大意: 就是现在给出一些单词的ci

poj2778 & hdu 2243 AC自动机+矩阵 &&BZOJ2938

AC自动机。 具体可以详看AC自动机 做这几道题需要注意一点、、如果某节点A的fail指针指向标志节点,那么这个节点A也是标志节点。如图 节点4和节点6肯定是病毒,由于节点3的fa...

HDU2243 考研路茫茫——单词情结 AC自动机DP矩阵优化

传送门:点击打开链接 考研路茫茫——单词情结 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java...

zoj 1091(pku2243,hdu1372) Knight Moves(广搜)

啥叫蛋疼……为啥蛋疼……我WA了一晚上,就因为……其实到现在我也不知道因为啥……之前一直错,在main开始run的时候,坐标录入我用的是scanf("%1s%d%1s%d",&ch1,&sj,&ch2...

【MZ】hdu 2243 AC自动机 + 矩阵加速

problem: 给你 n 个病毒,问长度为 [1, m] 的,包括病毒的有多少。全都是小写字母。 think: 和 poj 2778 有点像。不一样的是: 1)  这个是求包括的,我...

hdu2243之AC自动机+矩阵乘法

经典AC自动机+矩阵乘法
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDU 2243
举报原因:
原因补充:

(最多只允许输入30个字)