西安十五日游 day12 练习 【AC 自动机】

最新推荐文章于 2024-02-22 00:23:59 发布
最新推荐文章于 2024-02-22 00:23:59 发布
阅读量371 收藏 1
点赞数
分类专栏： AC自动机 # 系列 比赛题解 矩阵DP TSP 文章标签： AC自动机
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/nike0good/article/details/52203413
版权
比赛题解 同时被 3 个专栏收录
115 篇文章 0 订阅
订阅专栏
系列
38 篇文章 0 订阅
订阅专栏
AC自动机
6 篇文章 0 订阅
订阅专栏

版子题

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<sstream>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 


#define MAXNode (1000000)
#define Sigma_size (26)
class Aho_Corasick_Automata
{
public:
    int ch[MAXNode][Sigma_size];
    int v[MAXNode],siz;
    // AC自动机 
    int f[MAXNode],last[MAXNode];
    Aho_Corasick_Automata(int _siz=0):siz(_siz){MEM(ch) MEM(v) MEM(f) MEM(last)}
    void mem(int _siz=0){siz=_siz; MEM(ch) MEM(v) MEM(f) MEM(last)  }
    int idx(char c){return c-'a';}
    void insert(char *s) //val不为0 表示str末尾 
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            int c=idx(s[i]);
            if (!ch[u][c])
            {
                ++siz;
                MEM(ch[siz]);
                ch[u][c]=siz;
            }
            u=ch[u][c];
        }
        v[u]++;
    }
    void getFail()
    {
        queue<int> q;
        Rep(c,Sigma_size)
        {
            int u=ch[0][c];
            if (u) q.push(u),last[u]=0;
        }
        while (!q.empty())
        {
            int r=q.front();q.pop();  //r--c-->u
            Rep(c,Sigma_size)
            {
                int u=ch[r][c];
                if (!u) {ch[r][c]=ch[f[r]][c]; continue;} 
                q.push(u);
                f[u]=ch[f[r]][c];
                last[u]=v[f[u]]?f[u]:last[f[u]];
            }
        }
    } 
    int ans;
    void print(int j) //打印全串中所有以j为末尾的str 
    {
        if (j)
        {
//          printf("%d %d\n",j,v[j]);
            ans+=v[j]; v[j]=0;
            print(last[j]);
        }
    } 
    void find(char *s)
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            int c=idx(s[i]);
            u=ch[u][c];
            if (v[u]) print(u);
            else if (last[u]) print(u);
        }
    }

}T;

char s[1000000+10];
int main()
{
//  freopen("A.in","r",stdin);
//  freopen(".out","w",stdout);
    int t=read();
    while(t--) {
        int n=read();
        T.mem();
        For(i,n) {
            scanf("%s",s);
            T.insert(s);
        }
        T.getFail();
        scanf("%s",s);
        T.ans=0;
        T.find(s);
        cout<<T.ans<<endl;
    }


    return 0;
}

病毒侵袭

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<sstream>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define MAXN (500+10)
bool b[MAXN];

#define MAXNode (70000)
#define Sigma_size (100)
class Aho_Corasick_Automata
{
public:
    int ch[MAXNode][Sigma_size];
    int v[MAXNode],siz;
    // AC自动机 
    int f[MAXNode],last[MAXNode];
    Aho_Corasick_Automata(int _siz=0):siz(_siz){MEM(ch) MEM(v) MEM(f) MEM(last)}
    void mem(int _siz=0){siz=_siz; MEM(ch) MEM(v) MEM(f) MEM(last)  }
    int idx(char c){return c-31;}
    void insert(char *s,int val) //val不为0 表示str末尾 
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            int c=idx(s[i]);
            if (!ch[u][c])
            {
                ++siz;
                MEM(ch[siz]);
                ch[u][c]=siz;
            }
            u=ch[u][c];
        }
        v[u]=val;
    }
    void getFail()
    {
        queue<int> q;
        Rep(c,Sigma_size)
        {
            int u=ch[0][c];
            if (u) q.push(u),last[u]=0;
        }
        while (!q.empty())
        {
            int r=q.front();q.pop();  //r--c-->u
            Rep(c,Sigma_size)
            {
                int u=ch[r][c];
                if (!u) {ch[r][c]=ch[f[r]][c]; continue;} 
                q.push(u);
                f[u]=ch[f[r]][c];
                last[u]=v[f[u]]?f[u]:last[f[u]];
            }
        }
    } 
    void print(int j) //打印全串中所有以j为末尾的str 
    {
        if (j)
        {
            if (v[j]) {
                b[v[j]]=1; b[0]=1;
            }
            j=last[j];
        }
    } 
    void find(char *s)
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            int c=idx(s[i]);
            u=ch[u][c];
            if (v[u]) print(u);
            else if (last[u]) print(u);
        }
    }

}T;

char s[10000+10];
int n,m;
int main()
{
//  freopen("B.in","r",stdin);
//  freopen(".out","w",stdout);
//  cout<<sizeof(T)/1024<<endl;
    while(scanf("%d\n",&n)!=EOF) {
        T.mem();
        For(i,n) {
            gets(s);
            T.insert(s,i);
        }
        T.getFail();
        scanf("%d\n",&m);
        int tot=0;
        For(i,m) {
            gets(s);
            MEM(b)
            T.find(s);
            if(b[0]) {
                ++tot;
                printf("web %d:",i);
                For(i,n) if(b[i]) printf(" %d",i);
                puts("");
            }
        }
        printf("total: %d\n",tot);
    }


    return 0;
}

病毒侵袭持续中

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<sstream>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define MAXN (1000+10)
int b[MAXN];

#define MAXNode (50010)
#define Sigma_size (26)
class Aho_Corasick_Automata
{
public:
    int ch[MAXNode][Sigma_size];
    int v[MAXNode],siz;
    // AC自动机 
    int f[MAXNode],last[MAXNode];
    Aho_Corasick_Automata(int _siz=0):siz(_siz){MEM(ch) MEM(v) MEM(f) MEM(last)}
    void mem(int _siz=0){siz=_siz; MEM(ch) MEM(v) MEM(f) MEM(last)  }
    int idx(char c){return c-'A';}
    void insert(char *s,int val) //val不为0 表示str末尾 
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            int c=idx(s[i]);
            if (!ch[u][c])
            {
                ++siz;
                MEM(ch[siz]);
                ch[u][c]=siz;
            }
            u=ch[u][c];
        }
        v[u]=val;
    }
    void getFail()
    {
        queue<int> q;
        Rep(c,Sigma_size)
        {
            int u=ch[0][c];
            if (u) q.push(u),last[u]=0;
        }
        while (!q.empty())
        {
            int r=q.front();q.pop();  //r--c-->u
            Rep(c,Sigma_size)
            {
                int u=ch[r][c];
                if (!u) {ch[r][c]=ch[f[r]][c]; continue;} 
                q.push(u);
                f[u]=ch[f[r]][c];
                last[u]=v[f[u]]?f[u]:last[f[u]];
            }
        }
    } 
    void print(int j) //打印全串中所有以j为末尾的str 
    {
        while (j)
        {
            if (v[j]) {
                b[v[j]]++;
            }
            j=last[j];
        }
    } 
    void find(char *s)
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            if (!isupper(s[i])) {
                u=0;continue;
            }
            int c=idx(s[i]);
            u=ch[u][c];
            if (v[u]) print(u);
            else if (last[u]) print(u);
        }
    }

}T;
char buf[MAXN][51];
char s[2000000+10];
int n,m;
int main()
{
//  freopen("C.in","r",stdin);
//  freopen(".out","w",stdout);
    while(scanf("%d\n",&n)!=EOF) {
        T.mem();
        For(i,n) {
            gets(buf[i]);
            T.insert(buf[i],i);
        }
        T.getFail();
        gets(s);
        MEM(b)
        T.find(s);
        For(i,n) if (b[i]) {
            printf("%s: %d\n",buf[i],b[i]);
        }
    }


    return 0;
}

Word Puzzles

把查找串作为模式串，记得反过来

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<sstream>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define MAXN (1000+10)
int dir[8][2]={{-1,0},{-1,1},{0,1},{1,1},{1,0},{1,-1},{0,-1},{-1,-1}};
int n,m,q;
bool inside(int i,int j){return 0<=i&&i<n&&0<=j&&j<m;}
int ans[MAXN][3];
char buf[MAXN][MAXN];
char s[MAXN+10];
#define MAXNode (5000000)
#define Sigma_size (26)
class Aho_Corasick_Automata
{
public:
    int ch[MAXNode][Sigma_size];
    int v[MAXNode],siz;
    // AC自动机 
    int f[MAXNode],last[MAXNode];
    Aho_Corasick_Automata(int _siz=0):siz(_siz){MEM(ch) MEM(v) MEM(f) MEM(last)}
    void mem(int _siz=0){siz=_siz; MEM(ch) MEM(v) MEM(f) MEM(last)  }
    int idx(char c){return c-'A';}
    void insert(char *s,int val) //val不为0 表示str末尾 
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            int c=idx(s[i]);
            if (!ch[u][c])
            {
                ++siz;
                MEM(ch[siz]);
                ch[u][c]=siz;
            }
            u=ch[u][c];
        }
        v[u]=val;
    }
    void getFail()
    {
        queue<int> q;
        Rep(c,Sigma_size)
        {
            int u=ch[0][c];
            if (u) q.push(u),last[u]=0;
        }
        while (!q.empty())
        {
            int r=q.front();q.pop();  //r--c-->u
            Rep(c,Sigma_size)
            {
                int u=ch[r][c];
                if (!u) {ch[r][c]=ch[f[r]][c]; continue;} 
                q.push(u);
                f[u]=ch[f[r]][c];
                last[u]=v[f[u]]?f[u]:last[f[u]];
            }
        }
    } 
    void print(int j,int x,int y,int k) //打印全串中所有以j为末尾的str 
    {
        while (j)
        {
            if (v[j]) {
                ans[v[j]][0]=x;
                ans[v[j]][1]=y;
                ans[v[j]][2]=k;
            }
            j=last[j];
        }
    } 
    void find(int x,int y,int k)
    {
        int u=0;
        while(inside(x,y))
        {
            char now=buf[x][y];
            int c=idx(now);
            u=ch[u][c];
            if (v[u]) print(u,x,y,k);
            else if (last[u]) print(u,x,y,k);
            x+=dir[k][0],y+=dir[k][1];
        }
    }

}T;
int main()
{
//  freopen("D.in","r",stdin);
//  freopen(".out","w",stdout);
    int t=read();
    For(kcase,t) {
        cin>>n>>m>>q;
        T.mem();
        Rep(i,n) {
            scanf("%s",buf[i]);
        }
        For(i,q) {
            scanf("%s",s);
            reverse(s,s+strlen(s));
            T.insert(s,i);
        }
        T.getFail();
        Rep(i,m) {
            T.find(n-1,i,0);
            T.find(n-1,i,1);
            T.find(n-1,i,7);
            T.find(0,i,3);
            T.find(0,i,4);
            T.find(0,i,5);
        }
        Rep(i,n) {
            T.find(i,0,1);
            T.find(i,0,2);
            T.find(i,0,3);
            T.find(i,m-1,5);
            T.find(i,m-1,6);
            T.find(i,m-1,7);
        }

        For(i,q) {
            if (ans[i][2]<4) ans[i][2]+=4; else ans[i][2]-=4;
            printf("%d %d %c\n",ans[i][0],ans[i][1],'A'+ans[i][2]);
        }
        if(kcase<t) puts("");
    }


    return 0;
}

考研路茫茫——单词情结

长度不超过L，只由小写字母组成的，至少包含一个词根的单词，一共可能有多少个呢？
不包含的可以用AC自动机建立转移矩阵A，答案=A+A^2+A^3+..+A^L的第一行的和

#include<cstdio>
#include<iostream>
#include<functional>
#include<vector>
#include<cstring>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--))  
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define MAXN (100+10)
#define MAXNode (100+10)
#define Sigma_size (26)
#define pb push_back 
typedef unsigned long long ll;

ll pow2(ll a,ll b)
{
    if (b==1) return a;
    if (b==0) return 1;
    ll p=pow2(a,b/2);
    p=p*p;
    if (b&1) p=p*a;
    return p;
}

struct M  
{  
    int n,m;  
    ll a[MAXN][MAXN];  
    M(int _n=0){n=m=_n;MEM(a);} 
    M(int _n,int _m){n=_n,m=_m;MEM(a);}
    void mem (int _n=0){n=m=_n;MEM(a);}
    void mem (int _n,int _m){n=_n,m=_m;MEM(a);}

    friend M operator*(M a,M b)  
    {  
        M c(a.n,b.m);  
        For(k,a.m)
            For(i,a.n)  
                For(j,b.m)  
                    c.a[i][j]=(c.a[i][j]+a.a[i][k]*b.a[k][j]);  
        return c;     
    }  
    friend M operator+(M a,M b)  
    {  
        For(i,a.n)  
            For(j,a.m)  
                a.a[i][j]=(a.a[i][j]+b.a[i][j]);  
        return a;
    }  
    void make_I(int _n)  
    {  
        n=m=_n; MEM(a)
        For(i,n) a[i][i]=1;  
    }  
}A;

bool a3[1000000];  
M pow222(M a,ll b)  
{  
    M c;c.make_I(a.n);    
    int n=0;while (b) a3[++n]=b&1,b>>=1;
    c=a; b=1;
    M d=c;
    ForD(i,n-1)    
    {    
        b=b*2+a3[i];
        c=c*d+c;
        d=d*d;
        if (a3[i]) c=c*a+a,d=d*a;
    }    
    return c;    
}
ll pow222(ll a,ll b)  
{  
    ll c=1;    
    int n=0;while (b) a3[++n]=b&1,b>>=1;
    c=a; b=1;
    ll d=c;
    ForD(i,n-1)    
    {    
        b=b*2+a3[i];
        c=c*d+c;
        d=d*d;
        if (a3[i]) c=c*a+a,d=d*a;
    }    
    return c;    
}

int n;
ll m;
class Trie
{
public:
    int ch[MAXNode][Sigma_size],fa[MAXNode];
    int fail[MAXNode],go[MAXNode][Sigma_size];
    int siz;
    bool tag[MAXNode];
    void mem(){siz=0; MEM(tag) MEM(ch) MEM(fa) MEM(go) MEM(fail) }
    int idx(char c){return c-'a';}
    void insert(char *s)
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            int c=idx(s[i]);
            if (!ch[u][c])
            {
                ++siz;
                MEM(ch[siz]);
                ch[u][c]=siz;
                fa[siz]=u;
            }
            u=ch[u][c];
        }
        tag[u]=1;
    }
    int q[MAXNode+10];
    void make_fail()
    {
        int Head=1,Tail=1;
        q[Head]=0;
        For(i,siz) fail[i]=0;
        while(Head<=Tail) {
            int u=q[Head];
            Rep(c,Sigma_size)
            {
                if (ch[u][c]) {
                    q[++Tail]=ch[u][c];
                    if (u==0) fail[ch[u][c]]=0;
                    else fail[ch[u][c]]=go[ fail[u] ][c];
                    go[u][c]=ch[u][c];
                }
                else {
                    go[u][c]=go[fail[u]][c];
                }
            } 
            Head++; 
        }
        For(i,Tail) {
            int now=q[i];
            if (tag[fail[now]]) tag[now]=1;
        }

        int t=0;
        vector<int> Q;
        int h[MAXN]={0};
        Rep(i,siz+1) {
            if (!tag[i]) Q.pb(i),h[i]=++t;
        }
        A.mem(t);
        For(i,t)
        {
            int now=Q[i-1]; 
            Rep(j,Sigma_size) {
                int v=go[now][j];
                if (h[v]) A.a[i][h[v]]++; 
            } 
        } 
    }
}T;

char s[20];
ll L; 
int main()
{
//  freopen("P.in","r",stdin);
//  freopen(".out","w",stdout);
    while(cin>>n>>L) {
        T.mem();
        For(i,n) {
            scanf("%s",s);
            T.insert(s);
        }
        T.make_fail();
        A=pow222(A,L);
        ll ans=0;
        ll tot=0;
        tot=pow222(Sigma_size,L);
        For(i,A.n) ans+=A.a[1][i];
        cout<<tot-ans<<endl;
    }
    return 0;
}

Resource Archiver

题意：design a 01 string that contains all your resources (their occurrences can overlap), but none of the virus codes. To make your software smaller in size, the string should be as short as possible.
对每个单机在AC自动机上bfs一遍，求它到每个节点的最短距离，然后跑TSP。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
#include<vector>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<sstream>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (1000000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
#define Pr(kcase,ans) printf("Case #%d: %I64d\n",kcase,ans);
#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;
#define PRi2D(a,n,m) For(i,n) { \
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define MAXN (1000+10)
int b[MAXN];
int n,m;
int pos[MAXN];
int g[MAXN][MAXN];
#define MAXNode (4000000)
#define Sigma_size (2)
class Aho_Corasick_Automata
{
public:
    int ch[MAXNode][Sigma_size];
    int v[MAXNode],siz;
    // AC自动机 
    int f[MAXNode],last[MAXNode];
    Aho_Corasick_Automata(int _siz=0):siz(_siz){MEM(ch) MEM(v) MEM(f) MEM(last)}
    void mem(int _siz=0){siz=_siz; MEM(ch) MEM(v) MEM(f) MEM(last)  }
    int idx(char c){return c-'0';}
    void insert(char *s,int val) //val不为0 表示str末尾 
    {
        int u=0,n=strlen(s);
        Rep(i,n)
        {
            int c=idx(s[i]);
            if (!ch[u][c])
            {
                ++siz;
                MEM(ch[siz]);
                ch[u][c]=siz;
            }
            u=ch[u][c];
        }
        v[u]=val;
        if (val!=-1) pos[val]=u; 
    }
    void getFail()
    {
        queue<int> q;
        Rep(c,Sigma_size)
        {
            int u=ch[0][c];
            if (u) q.push(u),last[u]=0;
        }
        while (!q.empty())
        {
            int r=q.front();q.pop();  //r--c-->u
            Rep(c,Sigma_size)
            {
                int u=ch[r][c];
                if (!u) {ch[r][c]=ch[f[r]][c]; continue;} 
                q.push(u);
                f[u]=ch[f[r]][c];
                last[u]=v[f[u]]?f[u]:last[f[u]];
            }
        }
    } 

    int dis[MAXNode];
    void bfs(int u) {
        queue<int> q;
        MEMI(dis) 
        q.push(pos[u]); dis[pos[u]]=0;
        while(!q.empty()) {
            int u=q.front(); q.pop();
            Rep(c,2) {
                if (dis[ch[u][c]]==INF&&v[ch[u][c]]!=-1) {
                    dis[ch[u][c]]=dis[u]+1;
                    q.push(ch[u][c]);
                }
            }
        }
        For(i,n) g[u][i]=dis[pos[i]];
    }

}T;
char s[100000+10];
int dp[1<<11][20],p2[20],len[20];
void work() {
    Rep(i,n) p2[i+1]=1<<i;
    int S=1<<n;
    MEMI(dp) 
    For(i,n) dp[p2[i]][i]=len[i];
    Rep(st,S) {
        For(i,n) {
            if (dp[st][i]!=INF) {
                For(j,n) if (!(st&p2[j])) {
                    int ne=st|p2[j];
                    if (g[i][j]==INF) continue; 
                    if (dp[ne][j]>dp[st][i]+g[i][j]) dp[ne][j]=dp[st][i]+g[i][j];
                }
            }
        }
    }
    int ans=INF;
    For(i,n) ans = min(ans, dp[S-1][i] );
    cout<<ans<<endl;
}
int main()
{
//  freopen("Q.in","r",stdin);
//  freopen(".out","w",stdout);
    while(scanf("%d%d\n",&n,&m)!=EOF&&n&&m) {
        T.mem();
        For(i,n) {
            scanf("%s",s); 
            T.insert(s,i);
            len[i]=strlen(s);
        }
        For(i,m) {
            scanf("%s",s); 
            T.insert(s,-1);
        }
        T.getFail();
        For(i,n) T.bfs(i);
        work();
    }


    return 0;
}
nike0good
关注
专栏目录
AC自动机 算法详解(图解)及模板
bestsort
10-05 11万+
ac自动机,就是在tire树的基础上,增加一个fail指针,如果当前点匹配失败,则将指针**转移**到fail指针指向的地方,这样就不用回溯,而可以路匹配下去了.(当前模式串后缀和fail指针指向的模式串部分前缀相同,如`abce`和`bcd`,我们找到`c`发现下一个要找的不是`e`,就跳到`bcd`中的`c`处,看看此处的下一个字符(`d`)是不是应该找的那一个)
python实现AC自动机
uhauha2929的专栏
02-10 1252
ac自动机可以看成带指针的字典树，每个节点的指针指向了当前节点的最大后缀的位置。在建立字典树后，可以层次遍历字典树来构建fail指针，根节点的直接孩子（第一层节点）的fail指针肯定是指向根节点的，之后的节点需要看其父节点的fail指针指向的节点，如果父节点的fail指针指向的节点的孩子中存在当前节点，则将当前节点的fail指针指向该节点，如果不存在，需要沿着fail指针一直向上寻找直至根节点，如果根节点的孩子中也找不到，说明树中不存在该字符，将当前节点的fail指针指向根节点就行了。最后，在查找时也类似，
AC自动机
qq_45832461的博客
08-09 987
AC自动机 AC自动机的作用：给定一段包含mmm个字符的文章，再给定nnn个单词，找出有多少个单词在文章里面出现过。 AC自动机实际上是先将KMP算法和Trie字典树结合，用多个模式串构建一棵字典树，然后在这棵字典树上构建失配指针，失配指针相当于KMP算法中的next[]数组（匹配失败时的回退位置），最后将主串在Trie树上进行模式匹配。 AC自动机算法分为3步： 构建一棵字典树 构建AC自动机 进行模式匹配 这里主要讲解如何构建AC自动机和如何进行模式匹配。 构建AC自动机 在构建AC自动机时，其
字符串处理【AC自动机】 - 原理 AC自动机详解
谢谢你们的关注
12-08 2073
字符串处理【AC自动机】 - 原理 AC自动机详解
AC自动机算法详解以及Java代码实现
终生学习践行者
12-03 2805
详细介绍了AC自动机算法详解以及Java代码实现。
AC自动机算法
懒猫的脚印
08-15 1330
ac自动机，双数组
关于 AC 自动机
weixin_43501634的博客
02-22 3463
AC 自动机，全称 Aho-Corasick 自动机，是一种用于字符串搜索的算法，由 Alfred V. Aho 和 Margaret J. Corasick 在 1975 年提出。这个算法是为了解决在一个主文本字符串中查找多个模式字符串（或称为“关键词”）的问题，它可以在线性时间内完成搜索，非常高效。在计算机科学中，Aho-Corasick 算法是 Alfred V. Aho 和 Margaret J. Corasick 于 1975 年发明的一种字符串搜索算法。
Python实现多模匹配——AC自动机
GeekZW的博客
02-10 6144
Python实现多模匹配——AC自动机 目标：学习AC自动机，多模匹配。 要求：尽可能用纯Python实现，提升代码的扩展性。 一、什么是AC自动机？ AC自动机，Aho-Corasick automaton，该算法在1975年产生于贝尔实验室，是著名的多...
AC自动机原理 关键点 fail指针
stevewong的专栏
02-02 1524
前提 已经了解基础的Trie树，进一步看下AC自动机比Trie树牛皮在哪里，其实就是牛皮在fail指针可以减少匹配次数 参考 b站大法好 例子 好的例子是成功的一半 这个图很好地说明了fail指针的原理，还是看视频理解的好 ...
【算法无用系列】AC自动机敏感词过滤
叁滴水 的博客
10-23 2114
文章目录前言一、AC自动机1、原理总结 前言 鄙人不才，通过敏感词过滤这个功能，逐步引出了BM算法、BF算法、RK算法，最后发现这些算法只适合单个主串和模式串的匹配，不适合单主串和多模式串的匹配。无奈，查阅各种资料，终于命中AC自动机算法，专门用于单主串和多模式串匹配的算法，很适合敏感词过滤这个场景。 一、AC自动机 1、原理 由于例举敏感词文章审核会有问题，暂且用别的词替换。 AC自动机是一种特殊的树形结构，每一个根节点到叶节点都有一个唯一的路径，这个路径下的各个节点为一个词。假设有敏感词{} 假如
AC自动机_AC自动机模板_
09-30
AC自动机，全称为Aho-Corasick自动机，是一种字符串搜索算法，广泛应用于信息学奥林匹克竞赛、数据结构与算法竞赛以及文本处理等领域。它的主要功能是高效地在一个字符串集合中进行多模式匹配，即同时查找多个模式串...
ac自动机.pptx
07-14
要学AC自动机需要自备两个前置技能:KMP和trie树(其实个人感觉不会kmp也行,失配指针的概念并不难) 其中,KMP是用于一对一的字符串匹配,而trie虽然能用于多模式匹配,但是每次匹配失败都需要进行回溯,如果模式串很长的话...
高效Java敏感词过滤系统AC自动机算法源码，支持独立部署与集成注册中心
最新发布
03-25
本项目是一款高效的Java敏感词过滤系统，基于AC自动机算法实现。系统支持独立部署，同时可便捷集成至注册中心，为各类项目提供敏感词过滤服务。包含文件共117个，其中主要构成如下： - Java源文件：49个 - Class...
36丨AC自动机：如何用多模式串匹配实现敏感词过滤功能？1
08-03
AC自动机是一种高效的多模式串匹配算法，主要应用于文本处理中的敏感词过滤、关键词搜索等功能。它的全称为Aho-Corasick算法，由Aho和Corasick在1975年提出。AC自动机在Trie树的基础上进行了扩展，增加了失败指针，...
Codeforces Round #350 (Div. 2) 题解
nike0good |Oier&ACMer | 熟能生巧
05-07 3192
Holidays Game of Robots Cinema Magic Powder - 12 Correct Bracket Sequence Editor Restore a Number Holidays#include<bits/stdc++.h> using namespace std; #define For(i,n) for(int i=1;i<=n;i++) #define For
几何算法专题
nike0good |Oier&ACMer | 熟能生巧
07-02 2625
UVA 10652 Board Wrapping UVA 11168 Airport UVA 1396 Most Distant Point from the Sea UVA 1298 Triathlon UVA 1475 Jungle OutpostUVA 10652 Board Wrapping#include<bits/stdc++.h> using namespace std; #defin
hihoCoder挑战赛20 题解
nike0good |Oier&ACMer | 熟能生巧
05-06 2524
打折机票静态区间最大值SB题#include<bits/stdc++.h> using namespace std; #define For(i,n) for(int i=1;i<=n;i++) #define Fork(i,k,n) for(int i=k;i<=n;i++) #define Rep(i,n) for(int i=0;i<n;i++) #define ForD(i,n) for(i
几何基础专题
nike0good |Oier&ACMer | 熟能生巧
06-29 2172
UVA 11437 Triangle Fun UVA 11800 Determine the Shape 四边形判定 UVA 11646 Athletics Track UVA 11817 Tunnelling the Earth 球面距离 UVA 1473 Dome of Circus UVA 11524 InCircle UVA 11731 Ex-circles 旁切圆 UVA 12300 Sm
AC自动机算法详解：高效字符串匹配技术
资源摘要信息:"AC自动机（Aho-Corasick算法）是一种用于高效字符串匹配的算法，由Aho和Corasick提出。它可以在一个文本串中查找一组模式串的出现位置。AC自动机通过构建一个有限状态自动机（DFA）来实现，其中每个...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值