SPOJ
来源:互联网 发布:两位字母域名 编辑:程序博客网 时间:2024/06/03 08:30
题意:让你求这个字符串能组合出来多少个不重复的子串
思路:利用后缀数组height【i】的性质就可以求出每个位置开始能组成的字符串个数为n-sa[i]-height[i] , n为字符串的长度
#include<iostream>#include<cstring>#include<algorithm>#include<cstdio>using namespace std;const int mx = 50005;int wa[mx],wb[mx],wv[mx],c[mx];int sa[mx],rk[mx],r[mx],height[mx];char s[mx];bool cmp(int *r,int a,int b,int i){ return r[a] == r[b]&&r[a+i] == r[b+i];}void da(int *r,int *sa,int n,int m){ int *x = wa,*y = wb; for(int i = 0; i < m; i++) c[i] = 0; for(int i = 0; i < n; i++) c[x[i]=r[i]]++; for(int i = 1; i < m; i++) c[i] += c[i-1]; for(int i = n-1; i >= 0; i--) sa[--c[x[i]]] = i; for(int j = 1,p = 1; p < n;j<<=1,m = p){ p = 0; for(int i = n-j; i < n; i++) y[p++] = i; for(int i = 0; i < n; i++) if(sa[i]>=j) y[p++] = sa[i]-j; for(int i = 0; i < n; i++) wv[i] = x[y[i]]; for(int i = 0; i < m; i++) c[i] = 0; for(int i = 0; i < n; i++) c[wv[i]]++; for(int i = 1; i < m; i++) c[i] += c[i-1]; for(int i = n-1; i >= 0; i--) sa[--c[wv[i]]] = y[i]; swap(x,y),p = 1,x[sa[0]] = 0; for(int i = 1; i < n; i++) x[sa[i]] = cmp(y,sa[i],sa[i-1],j)?p-1:p++; }}void getheight(int *r,int n){ int k = 0,j; for(int i = 1; i <= n; i++) rk[sa[i]] = i; for(int i = 0; i < n; height[rk[i++]] = k) for(j = sa[rk[i]-1],k?k--:0;r[i+k] == r[j+k]; k++);}int main(){ int t; scanf("%d",&t); while(t--){ scanf("%s",s); int n = strlen(s); for(int i = 0; i < n; i++) r[i] = s[i]; r[n] = 0; da(r,sa,n+1,256); getheight(r,n); int ans = 0; for(int i = 1; i <= n; i++){ ans += n-sa[i]-height[i]; //printf("%d %d %d\n",n,sa[i],height[i]); } cout<<ans<<endl; } return 0;}
阅读全文
0 0
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- SPOJ
- Java中的多线程你只要看这一篇就够了
- spot 状压dp
- Android 在本地使用SharedPreferences,保存大量的数据
- Unicode编码的byte 数组转为String
- 交换排序——快速排序
- SPOJ
- PHP时间戳转日期
- 解决jenkins运行磁盘满的问题
- 二.HighGUI图形用户界面初步(1.Read_Show_Write)
- HDU3486-RMQ+二分
- windows技术----HOOK技术学习
- navicat与mysql的连接
- 高仿微信APP实战(一)-Actionbar的制作与应用
- 微信小程序 模拟加载请求数据(六)