【ZOJ 3228】Searching the String 【AC自动机】

题意

  给出n个模式串和一个文本串,输出各个模式串在文本串中出现的次数。模式串有两种类型,0类型代表可以有重叠,1类型代表不能有重叠。模式串可能出现重复。

分析

  算是AC自动机的模板题?

  因为模式串可以重复,所以如果直接插入并且用val数组来保存模式串的编号的话,后面出现的会把前面出现的给覆盖。所以我这里用了一个map来保存每个模式串在trie中的编号。

  如何处理1类型不能有重叠的情况?对于1类型的每个模式串,记录一下它的长度和上次匹配到的位置。当再次匹配到这个模式串的时候,看一下这次的位置和上次位置的差有没有大于它的长度,如果大于,则说明这个可以选择不会重叠。

  我们在插入模式串的时候不区分是哪种类型,在进行find的时候也不进行区分,找到一个模式串以后,同时更新两种类型。只在最后输出的时候区分一下。

  

  1 #include <cstdio>
  2 #include <cstring>
  3 #include <algorithm>
  4 #include <iostream>
  5 #include <queue>
  6 #include <map>
  7 
  8 using namespace std;
  9 const int maxnode=700000;
 10 const int sigma_size=26;
 11 const int maxs=100000+100;
 12 char T[maxs],P[maxs][20],type[maxs];
 13 map<string,int>ms;
 14 int kase;
 15 struct AC_Automata{
 16     int ch[maxnode][sigma_size],val[maxnode],cnt[maxnode][2],tim[maxnode];
 17     int f[maxnode],last[maxnode],len[maxnode];
 18     int sz;
 19     void init(){
 20         sz=1;
 21         memset(ch[0],0,sizeof(ch[0]));
 22         memset(cnt,0,sizeof(cnt));
 23         memset(tim,0,sizeof(tim));
 24         val[0]=0;
 25         ms.clear();
 26     }
 27     void insert(char *s){
 28         int n=strlen(s),u=0;
 29         for(int i=0;i<n;i++){
 30             int c=s[i]-'a';
 31             if(!ch[u][c]){
 32                 ch[u][c]=sz;
 33                 memset(ch[sz],0,sizeof(ch[sz]));
 34                 val[sz++]=0;
 35             }
 36             u=ch[u][c];
 37         }
 38         val[u]=1;
 39         len[u]=n;
 40         string S=(string)s;
 41         ms[S]=u;
 42     }
 43     void getFail(){
 44         queue<int>q;
 45         f[0]=last[0]=0;
 46         for(int i=0;i<sigma_size;i++){
 47             int u=ch[0][i];
 48             if(u){
 49                 q.push(u);
 50                 f[u]=last[u]=0;
 51             }
 52         }
 53         while(!q.empty()){
 54             int r=q.front();q.pop();
 55             for(int i=0;i<sigma_size;i++){
 56                 int u=ch[r][i];
 57                 if(!u)continue;
 58                 q.push(u);
 59                 int v=f[r];
 60                 while(v&&!ch[v][i])v=f[v];
 61                 f[u]=ch[v][i];
 62                 last[u]=val[f[u]]?f[u]:last[f[u]];
 63             }
 64         }
 65     }
 66     void print(int i,int pos){
 67         if(val[i]){
 68             cnt[i][0]++;
 69             if(tim[i]+len[i]<=pos){
 70                 cnt[i][1]++;
 71                 tim[i]=pos;
 72             }
 73             print(last[i],pos);
 74         }
 75     }
 76     void find(char *s){
 77         int n=strlen(s),j=0;
 78         for(int i=0;i<n;i++){
 79             int c=s[i]-'a';
 80             while(j&&!ch[j][c])j=f[j];
 81             j=ch[j][c];
 82             if(val[j])
 83                 print(j,i+1);
 84             else if(last[j])
 85                 print(last[j],i+1);
 86         }
 87     }
 88 }ac;
 89 int n;
 90 int main(){
 91     kase=0;
 92     while(scanf("%s",T)!=EOF){
 93         ++kase;
 94         scanf("%d",&n);
 95         ac.init();
 96         for(int i=1;i<=n;i++){
 97             scanf("%d %s",&type[i],P[i]);
 98             ac.insert(P[i]);
 99         }
100         ac.getFail();
101         ac.find(T);
102         printf("Case %d\n",kase);
103         for(int i=1;i<=n;i++){
104             string S=(string)P[i];
105             int u=ms[S];
106             printf("%d\n",ac.cnt[u][type[i]]);
107           // printf("%d\n",u);
108         }
109         printf("\n");
110     }
111 return 0;
112 }
View Code

 

转载于:https://www.cnblogs.com/LQLlulu/p/9330112.html

以下是ZOJ1626的C++ AC代码,使用了旋转卡壳算法: ```c++ #include <iostream> #include <cstdio> #include <algorithm> #include <cmath> #include <cstring> #define MAXN 100010 #define eps 1e-8 #define INF 1e20 using namespace std; struct point { double x,y; friend point operator -(point a,point b) { point res; res.x=a.x-b.x; res.y=a.y-b.y; return res; } friend bool operator <(point a,point b) { if(fabs(a.x-b.x)<eps) return a.y<b.y; return a.x<b.x; } friend double operator *(point a,point b) { return a.x*b.y-a.y*b.x; } friend double dis(point a,point b) { return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); } }a[MAXN],b[MAXN],st[MAXN]; int n; double ans=INF; int cmp(point a,point b) { double tmp=(a-b)*(a[1]-b); if(fabs(tmp)<eps) return dis(a,a[1])-dis(b,a[1])<0; return tmp>0; } int main() { while(~scanf("%d",&n) && n) { for(int i=1;i<=n;i++) scanf("%lf%lf",&a[i].x,&a[i].y); sort(a+1,a+n+1); int tot=0; for(int i=1;i<=n;i++) { while(tot>=2 && (st[tot]-st[tot-1])*(a[i]-st[tot])<0) tot--; st[++tot]=a[i]; } int k=tot; for(int i=n-1;i>=1;i--) { while(tot>k && (st[tot]-st[tot-1])*(a[i]-st[tot])<0) tot--; st[++tot]=a[i]; } tot--; for(int i=1;i<=tot;i++) b[i]=st[i]; int tmp=1; for(int i=2;i<=tot;i++) if(b[i].y<b[tmp].y) tmp=i; swap(b[1],b[tmp]); sort(b+2,b+tot+1,cmp); st[1]=b[1]; st[2]=b[2]; k=2; for(int i=3;i<=tot;i++) { while(k>1 && (st[k]-st[k-1])*(b[i]-st[k])<=0) k--; st[++k]=b[i]; } double ans=0; if(k==2) ans=dis(st[1],st[2]); else { st[k+1]=st[1]; for(int i=1;i<=k;i++) for(int j=1;j<=k;j++) ans=max(ans,dis(st[i],st[j])); } printf("%.2lf\n",ans/2); } return 0; } ``` 其中,结构体 `point` 表示二维平面上的一个点,包含了点的坐标和一些基本操作。函数 `cmp` 是旋转卡壳算法中的比较函数,按照点到起点的极角从小到大排序。在主函数中,先使用 Graham 扫描法求出点集的凸包,然后按照旋转卡壳的步骤,求出凸包上的最远点对距离作为最小直径。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值