【ZOJ 3430 Detect the Virus】 AC自动机

ZOJ3430
题意就是给你一个加密的单词,加密方式就是先把字符按照ASCII表转为2进制,再截取6个为一段转为10进制,再去表中改为对应字符。
我们只需要把所有关键字和文本串都解码之后,就是一个最朴素的AC自动机了。
密文转换方法是向kuangbin大神学习的,我们通过思考可以发现,解码的过程就是把四个字符转换为3个字符的过程,所以我们只要对加密串进行一些二进制操作,就可以直接将24位二进制码分为3个八位二进制码。具体过程请参考代码。
注意此题字符可能有256种,所以要用unsigned char
而且 ZOJ RE和MLE均返回 Segmentation Fault ,此题卡内存比较严重,需要好好计算一下
ZOJ3430代码

#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<string.h>
#include<queue>
#include<set>
using namespace std;
const int maxn =5e4+5;
unsigned char aa[maxn];
#define dbg(x) cout<<#x<<" "<<x<<endl;
struct ACTrie
{
  int tree[maxn][256], fail[maxn],end_[maxn];
  int root,cnt,num;
  set<int> s;
  int newnode()
  {
      for(int i=0;i<256;i++)
        tree[cnt][i]=-1;
      end_[cnt]=0;
      return cnt++;
  }
  void init()
  {
      cnt=0;
      num=0;
      root=newnode();
  }
  void insert_(unsigned char str[],int len)
  {
      int pos=root;
      for(int i=0;i<len;i++)
      {
          int id=str[i];
          if(tree[pos][id]==-1) tree[pos][id]=newnode();
          pos=tree[pos][id];
      }
      end_[pos]=++num;
  }
  void build()
  {
      queue<int> que;
      fail[root]=root;
      for(int i=0;i<256;i++)
      {
          if(tree[root][i]==-1) tree[root][i]=root;
          else
          {
              fail[tree[root][i]]=root;
              que.push(tree[root][i]);
          }
      }
      while(!que.empty())
      {
          int now=que.front();
          que.pop();
          for(int i=0;i<256;i++)
          {
              if(tree[now][i]==-1) tree[now][i]=tree[fail[now]][i];
              else
              {
                  fail[tree[now][i]]=tree[fail[now]][i];
                  que.push(tree[now][i]);
              }
          }
      }
  }
  int query(unsigned char *str,int len)
  {
       s.clear();
        int now=root;
        int res=0;
        for(int i=0;i<len;i++)
        {
            int id=(int)str[i];
            now=tree[now][id];
            int temp=now;
            while(temp!=root)
            {
                if(end_[temp]!=0) s.insert(end_[temp]);
                temp=fail[temp];
            }
        }
        return (int)s.size();
  }
};
unsigned char cal(char c)
{
    if(c>='A'&&c<='Z')
    {
        return c-'A';
    }
    if(c>='a'&&c<='z')
    {
        return c-'a'+26;
    }
    if(c>='0'&&c<='9')
    {
        return c-'0'+52;
    }
    if(c=='+')  return 62;
    else return 63;
}
int change(unsigned char str[],int len)
{
    int t=0;
    for(int i=0;i<len;i+=4)//每四个字符计算一次,转换为三个字符
    {
        aa[t++]=((str[i]<<2)|(str[i+1]>>4));//截取第一个字符的全部+第二个字符前两位
        if(i+2 < len)
            aa[t++]=( (str[i+1]<<4)|(str[i+2]>>2) );//截取第二个字符后四位+第三个字符前四位
        if(i+3 < len)
            aa[t++]= ( (str[i+2]<<6)|str[i+3] );//截取第三个字符后两位+第四个字符全部
    }
    return t;
}
ACTrie ac;
char str[maxn];
unsigned char ss[maxn];
int main()
{
    int n,m;
    while(scanf("%d",&n)!=EOF)
    {
        ac.init();
        while(n--)
        {
            scanf("%s",str);
            int len=strlen(str);
            while(str[len-1]=='=') len--;
            for(int i=0;i<len;i++) ss[i]=cal(str[i]);
            int len2 = change(ss,len);
            ac.insert_(aa,len2);//转换之后插入AC自动机
        }
        ac.build();
        scanf("%d",&m);
        while(m--)
        {
            //gets(str);
            scanf("%s",str);
            int len=strlen(str);
            while(str[len-1]=='=') len--;
            for(int i=0;i<len;i++) ss[i]=cal(str[i]);
            int len2 = change(ss,len);
            printf("%d\n",ac.query(aa,len2));
        }
        printf("\n");
    }
    return 0;
}
以下是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、付费专栏及课程。

余额充值