洛谷P3966 单词 [TJOI2013] AC自动机

正解:AC自动机

解题报告:

传送门!

先来提供一个40pts错解QAQ

首先看到这题就会想到AC自动机板子题2鸭!然后就照着那题的套路打一下,随便改一点儿,简单来说就是每次经过一个节点都要++,然后每搜到一个节点就暴跳所有fail,tr[fail].as+=tr[nw].as

然后你就能得到40pts的好成绩,,,另外60pts还是WA的,,,就很难过,,,

#include<bits/stdc++.h>
using namespace std;
#define ll int
#define rp(i,x,y) for(register ll i=x;i<=y;++i)

const ll N=+,M=+;
ll n,cnt,as,ed[M];
struct tre{ll as,nxt[],fail;tre(){as=;memset(nxt,,sizeof(nxt));fail=;}}tr[N];
bool gdgs=;

inline ll read()
{
    register char ch=getchar();register ll x=;register bool y=;
    while(ch!='-' && (ch>'' || ch<''))ch=getchar();
    if(ch=='-')ch=getchar(),y=;
    while(ch>='' && ch<='')x=(x<<)+(x<<)+(ch^''),ch=getchar();
    return y?x:-x;
}
inline void bd(ll nm)
{
    string str;cin>>str;ll lth=str.length()-,nw=;
    rp(i,,lth)
    {
        if(tr[nw].nxt[str[i]-'a'+]==)tr[nw].nxt[str[i]-'a'+]=++cnt;
        nw=tr[nw].nxt[str[i]-'a'+];++tr[nw].as;
    }
    ed[nm]=nw;
}
inline void fl()
{
    queue<ll>Q;
    rp(i,,)if(tr[].nxt[i])Q.push(tr[].nxt[i]);
    while(!Q.empty())
    {
        ll nw=Q.front();Q.pop();
        rp(i,,)
        {
            if(tr[nw].nxt[i])tr[tr[nw].nxt[i]].fail=tr[tr[nw].fail].nxt[i],Q.push(tr[nw].nxt[i]);
            else tr[nw].nxt[i]=tr[tr[nw].fail].nxt[i];
        }
    }
}
inline void zdj()
{
    rp(i,,cnt)
    {
    //    printf("i=%d as=%d\n",i,tr[i].as);
        for(register ll j=i;j;j=tr[j].fail)tr[tr[j].fail].as+=tr[i].as/*,printf("j=%d fl=%d\n",j,tr[j].fail)*/;
    }
}

int main()
{
    n=read();rp(i,,n)bd(i);fl();tr[].fail=;zdj();//printf("QAQ?\n");
    rp(i,,n)printf("%d\n",tr[ed[i]].as);
    return ;
}

放下sd灵巧的sd代码呜呜呜

研究辣半天回来repo一下我错哪儿辣,,,我好像是打挂辣呜呜呜

不过还是有个点有问题

就是,暴跳会T

可以在配fail的时候记录下bfs序,然后直接按bfs倒叙序每次就只有这个点的fail要跳辣而不是所有的fail都要跳一遍

然后就over辣!

话说为什么这题居然有紫,,,?病毒都比这题难一些啊我jio得,,,?

等下放代码QAQ

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define rg register
#define rp(i,x,y) for(rg ll i=x;i<=y;++i)
#define my(i,x,y) for(rg ll i=x;i>=y;--i)

const ll N=+,M=+;
ll n,cnt,as,ed[M],que[N],tmp;
struct tre{ll as,nxt[],fail;tre(){as=;memset(nxt,,sizeof(nxt));fail=;}}tr[N];

inline ll read()
{
    rg char ch=getchar();rg ll x=;rg bool y=;
    while(ch!='-' && (ch>'' || ch<''))ch=getchar();
    if(ch=='-')ch=getchar(),y=;
    while(ch>='' && ch<='')x=(x<<)+(x<<)+(ch^''),ch=getchar();
    return y?x:-x;
}
inline void bd(ll nm)
{
    string str;cin>>str;ll lth=str.length()-,nw=;
    rp(i,,lth)
    {
        if(!tr[nw].nxt[str[i]-'a'+])tr[nw].nxt[str[i]-'a'+]=++cnt;
        nw=tr[nw].nxt[str[i]-'a'+];++tr[nw].as;
    }
    ed[nm]=nw;
}
inline void fl()
{
    queue<ll>Q;
    rp(i,,)if(tr[].nxt[i])Q.push(tr[].nxt[i]),que[++tmp]=tr[].nxt[i];
    while(!Q.empty())
    {
        ll nw=Q.front();Q.pop();
        rp(i,,)
        {
            if(tr[nw].nxt[i])tr[tr[nw].nxt[i]].fail=tr[tr[nw].fail].nxt[i],Q.push(tr[nw].nxt[i]),que[++tmp]=tr[nw].nxt[i];
            else tr[nw].nxt[i]=tr[tr[nw].fail].nxt[i];
        }
    }
}

int main()
{
//    freopen("cjkcjk.in","r",stdin);
//    freopen("cjkcjk.out","w",stdout);
    n=read();rp(i,,n)bd(i);fl();tr[].fail=;my(i,tmp,)tr[tr[que[i]].fail].as+=tr[que[i]].as;
    rp(i,,n)printf("%lld\n",tr[ed[i]].as);
    return ;
}

这儿是代码吼!