zoj 3430(ac自动机)

题意:给出n个模式串,是以base64编码的,m个目标串,也是base64编码,问每个目标串里有多少种模式串。base64是把普通ASCII字符二进制表示的每6位截断对应一个最大编号到63的字符表可以对应一个字符,也就是每四个字符对应三个个普通字符,不够就末尾补零。补两个0是”=”,补四个0是”==”。
题解:先把所有给出的串都先解码为普通的字符串,然后是自动机模板解法。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <map>
using namespace std;
const int N = 33000;
const int SIGMA_SIZE = 256;
const char cb64[] = "ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz0123456789+/";
unsigned char str[3000];
char str2[3000];
int Next[N][SIGMA_SIZE], val[N], fail[N], sz, n, m, vis[550];
map<char, unsigned char> mp;

int decode(char *s) {
    int cnt = 0, len = strlen(s);
    while (s[len - 1] == '=')
        len--;
    for (int i = 0; i < len; i += 4) {
        str[cnt++] = (mp[s[i]] << 2) | (mp[s[i + 1]] >> 4);
        if (i + 2 < len)
            str[cnt++] = (mp[s[i + 1]] << 4) | (mp[s[i + 2]] >> 2);
        if (i + 3 < len)
            str[cnt++] = (mp[s[i + 2]] << 6) | mp[s[i + 3]];
    }
    return cnt;
}

void init() {
    sz = 1;
    memset(Next[0], 0, sizeof(Next[0]));
    val[0] = 0;
}

void insert(unsigned char *s, int len, int v) {
    int u = 0;
    for (int i = 0; i < len; i++) {
        if (!Next[u][s[i]]) {
            memset(Next[sz], 0, sizeof(Next[sz]));
            val[sz] = 0;
            Next[u][s[i]] = sz++;
        }
        u = Next[u][s[i]];
    }
    val[u] = v;
}

void getFail() {
    queue<int> Q;
    fail[0] = 0;
    for (int i = 0; i < SIGMA_SIZE; i++) {
        if (Next[0][i]) {
            fail[Next[0][i]] = 0;
            Q.push(Next[0][i]);
        }
    }
    while (!Q.empty()) {
        int u = Q.front();
        Q.pop();
        for (int i = 0; i < SIGMA_SIZE; i++) {
            if (!Next[u][i])
                Next[u][i] = Next[fail[u]][i];
            else {
                fail[Next[u][i]] = Next[fail[u]][i];
                Q.push(Next[u][i]);
            }
        }
    }
}

int query(unsigned char *s, int len) {
    int u = 0, cnt = 0;
    memset(vis, 0, sizeof(vis));    
    for (int i = 0; i < len; i++) {
        u = Next[u][s[i]];
        int temp = u;
        while (temp != 0) {
            if (val[temp])
                vis[val[temp]] = 1;
            temp = fail[temp];
        }
    }
    for (int i = 1; i <= n; i++)
        if (vis[i])
            cnt++;
    return cnt;
}

int main() {
    for (int i = 0; i < 64; i++)
        mp[cb64[i]] = i;
    while (scanf("%d", &n) == 1) {
        init();
        for (int i = 1; i <= n; i++) {
            scanf("%s", str2);
            int len = decode(str2);
            insert(str, len, i);
        }
        getFail();
        scanf("%d", &m);
        while (m--) {
            scanf("%s", str2);
            int len = decode(str2);
            printf("%d\n", query(str, len));
        }
        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、付费专栏及课程。

余额充值