国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發設計 > 正文

bzoj 3172: [Tjoi2013]單詞 后綴數組+rmq+二分

2019-11-08 01:44:19
字體:
來源:轉載
供稿:網友

題意

某人讀論文,一篇論文是由許多單詞組成。但他發現一個單詞會在論文中出現很多次,現在想知道每個單詞分別在論文中出現多少次。每個單詞由小寫字母組成,N<=200,文章長度不超過10^6

分析

AC自動機模板,但由于很久沒打sa了想要復習一下,于是就用了sa+rmq+二分的方法來做。

代碼

#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#include<cmath>#define N 1010005using namespace std;int sa[N],rank[N],height[N],b[N],c[N],d[N],m,pos[205],rmq[N][30],s[N],len[205];char ch[N];void getsa(int n,int m){ for (int i=0;i<=m;i++) b[i]=0; for (int i=1;i<=n;i++) b[s[i]]++; for (int i=1;i<=m;i++) b[i]+=b[i-1]; for (int i=n;i>=1;i--) c[b[s[i]]--]=i; int t=0; for (int i=1;i<=n;i++) { if (s[c[i]]!=s[c[i-1]]) t++; rank[c[i]]=t; } int j=1; while (j<=n) { for (int i=0;i<=n;i++) b[i]=0; for (int i=1;i<=n;i++) b[rank[i+j]]++; for (int i=1;i<=n;i++) b[i]+=b[i-1]; for (int i=n;i>=1;i--) c[b[rank[i+j]]--]=i; for (int i=0;i<=n;i++) b[i]=0; for (int i=1;i<=n;i++) b[rank[i]]++; for (int i=1;i<=n;i++) b[i]+=b[i-1]; for (int i=n;i>=1;i--) d[b[rank[c[i]]]--]=c[i]; int t=0; for (int i=1;i<=n;i++) { if (rank[d[i]]!=rank[d[i-1]]||rank[d[i]]==rank[d[i-1]]&&rank[d[i]+j]!=rank[d[i-1]+j]) t++; c[d[i]]=t; } for (int i=1;i<=n;i++) rank[i]=c[i]; if (t==n) break; j*=2; } for (int i=1;i<=n;i++) sa[rank[i]]=i;}void getheight(int n){ int k=0; for (int i=1;i<=n;i++) { if (k) k--; int j=sa[rank[i]-1]; while (i+k<=n&&j+k<=n&&s[i+k]==s[j+k]) k++; height[rank[i]]=k; }}void get_rmq(int n){ int lg=log(n)/log(2); for (int i=1;i<=n;i++) rmq[i][0]=height[i]; for (int j=1;j<=lg;j++) for (int i=1;i<=n-(1<<j)+1;i++) rmq[i][j]=min(rmq[i][j-1],rmq[i+(1<<(j-1))][j-1]);}int get_min(int l,int r){ int lg=log(r-l+1)/log(2); return min(rmq[l][lg],rmq[r-(1<<lg)+1][lg]);}int efl(int pos,int len,int n){ int l=1,r=pos-1; while (l<=r) { int mid=(l+r)/2; if (get_min(mid+1,pos)>=len) r=mid-1; else l=mid+1; } return r+1;}int efr(int pos,int len,int n){ int l=pos+1,r=n; while (l<=r) { int mid=(l+r)/2; if (get_min(pos+1,mid)>=len) l=mid+1; else r=mid-1; } return l-1;}int main(){ scanf("%d",&m); int n=0,tot=123; for (int i=1;i<=m;i++) { scanf("%s",ch); len[i]=strlen(ch);pos[i]=n+1; for (int j=0;j<len[i];j++) s[++n]=ch[j]; s[++n]=tot++; } getsa(n,tot); getheight(n); get_rmq(n); for (int i=1;i<=m;i++) { int l=efl(rank[pos[i]],len[i],n),r=efr(rank[pos[i]],len[i],n);;
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 昭苏县| 满城县| 洱源县| 柘城县| 类乌齐县| 乡城县| 古田县| 施甸县| 家居| 张北县| 交口县| 鄂伦春自治旗| 长葛市| 铜梁县| 剑阁县| 长丰县| 乳山市| 临沭县| 长子县| 乌什县| 安岳县| 辽源市| 拉萨市| 东乌珠穆沁旗| 思茅市| 富平县| 神池县| 晋宁县| 绥滨县| 神农架林区| 曲沃县| 依安县| 买车| 新宾| 太保市| 平利县| 进贤县| 武宣县| 天津市| 云梦县| 北碚区|