hdu 4436 str2int ( 后缀数组 )

题意 : 给你n个数字字符串 ， 把其中的子串转化成整数，求所有不同的整数的和。

比如 101 和 123 ， 那么 结果就是 1 + 10 + 101 + 12 + 123 + 2 + 23 + 3 = 275

思路 : 首先我们先要求出这n个字符串中所有不同的子串 ， 并且这些子串不能是‘0’ 开头的

将这些字符串用没有出现过的不同的字符分隔开，然后求height数组。

从头开始遍历height数组，每遇到一个后缀，我们需要求出这个后缀能提供多少前缀。因为这个是n个字符串拼接成的，所以我们需要知道当前后缀的长度，这个在开始预处理一下就可以了。但是中间肯定会有重复，那么哪些是重复的，就是跟前面后缀的LCP，但是因为要求的子串不能以'0'开头，所以我这里前面的后缀应该是前面第一个不是以'0'开头的后缀。

然后就是如何求和的问题。假设这个后缀的长度为L，跟前面的前缀的LCP为len，那么我们要求的后面L-len个后缀的和(前面的是重复的)。

假设我们求的字符串 "12345"中的34，我们要怎么快速的求出。

我们定义了两个数组，一个是sum[]表示前缀和，一个是num[]表示数字，比如上面的例子

sum[1] = 1                                                          num[1] = 1

sum[2] = 1 + 12                                                  num[2] = 12

sum[3] = 1 + 12 + 123                                        num[3] = 123

...                                                                        num[4] = 1234

sum[5] = 1 + 12 + 123 + 1234 + 12345              num[5] = 12345

我们首先将 sum[4] - sum[2] 得到 123 + 1234 ， 而我们想要得到是 3 + 34

所以我们要减去 120 和 1200 即 num[2] * 110

用这样的方法就可以快速的求和

#include <stdio.h>  
#include <string.h>
#include <string>
#include <vector>
#include <iostream>
#include <algorithm>  
using namespace std;  

#define INF 0x3f3f3f3f
#define maxn 120005
int wa[maxn],wb[maxn],wv[maxn],wt[maxn];  

typedef long long LL ;

int cmp(int *r,int a,int b,int l)  
{return r[a]==r[b]&&r[a+l]==r[b+l];}  

void da(int *r,int *sa,int n,int m){  
	int i,j,p,*x=wa,*y=wb,*t;  
	for(i=0;i<m;i++) wt[i]=0;  
	for(i=0;i<n;i++) wt[x[i]=r[i]]++;  
	for(i=1;i<m;i++) wt[i]+=wt[i-1];  
	for(i=n-1;i>=0;i--) sa[--wt[x[i]]]=i;  
	for(j=1,p=1;p<n;j*=2,m=p){  
		for(p=0,i=n-j;i<n;i++) y[p++]=i;  
		for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;  
		for(i=0;i<n;i++) wv[i]=x[y[i]];  
		for(i=0;i<m;i++) wt[i]=0;  
		for(i=0;i<n;i++) wt[wv[i]]++;  
		for(i=1;i<m;i++) wt[i]+=wt[i-1];  
		for(i=n-1;i>=0;i--) sa[--wt[wv[i]]]=y[i];  
		for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)  
			x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;  
	}  
}  

int Rank[maxn],height[maxn];  

void calheight(int *r,int *sa,int n){  
	int i , j , k = 0 ;  
	for( i=1 ; i<=n ; i++ ) Rank[sa[i]]=i;  
	for(i=0;i<n;i++) {  
		if(k)k--;  
		int j = sa[Rank[i]-1];  
		while(r[i+k]==r[j+k]) k++ ;  
		height[Rank[i]] = k ;  
	}     
	return;  
}  

int r[maxn] , sa[maxn] ;  
char str[maxn] ;

int n ;
int length[maxn] ;
int sum[maxn] ;
int num[maxn] ;
int ten[maxn] ;

int main(){

	ten[0] = 0 ;
	for( int i = 1 ; i <= 120000 ; i ++ ) 
		ten[i] = ( ( ten[i-1] + 1 ) * 10 ) % 2012 ;


	while( scanf( "%d" , &n ) != EOF ) {
		int len = 1 ;
		memset( length , 0 , sizeof(length) ) ;
		
		sum[0] = num[0] = 0 ;

		r[0] = '9' + 1 ;
		for( int i = 1 ; i <= n; i ++ ) {
			scanf( "%s" , str ) ;
			int l = strlen( str ) ;
			for( int j = 0 ; j < l ; j ++ ) {
				r[len] = str[j] ;
				num[len] = ( num[len-1] * 10 + str[j] - '0' ) % 2012 ;
				sum[len] = ( sum[len-1] + num[len] ) % 2012 ;
				length[len++] = l - j ;
			}
			sum[len] = num[len] = 0 ;
			r[len++] = '9' + i + 1 ;
		}
		sum[len] = num[len] = 0 ;
		r[len] = 0 ;

		da( r , sa , len + 1 , 10200 ) ;
		calheight( r , sa , len ) ;
		int ans = 0 ;
		int Min = 0 ;
		if( r[sa[1]] >= '1' && r[sa[1]] <= '9' ) {
			int s = sa[1] - 1 , t1 = s + length[sa[1]] , t2 = s + Min ;
			ans += sum[t1] - sum[s] - num[s] * ten[t1-s] - ( sum[t2] - sum[s] - num[s] * ten[t2-s] ) ; 
			Min = INF ;
			ans %= 2012 ;
			ans = ( ans + 2012 ) % 2012 ;
		} 
		for( int i = 2 ; i <= len ; i ++ ) {
			Min = min( Min , height[i] ) ;
			if( r[sa[i]] >= '1' && r[sa[i]] <= '9' ) {
				int s = sa[i] - 1 , t1 = s + length[sa[i]] , t2 = s + Min ;
				ans += sum[t1] - sum[s] - num[s] * ten[t1-s] - ( sum[t2] - sum[s] - num[s] * ten[t2-s] ) ; 
				Min = INF ;
				ans %= 2012 ;
				ans = ( ans + 2012 ) % 2012 ;
			}
		}
		printf( "%d\n" , ans ) ;
	}
	return 0 ;
}