程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
 程式師世界 >> 編程語言 >> C語言 >> C++ >> C++入門知識 >> SPOJ 7258 SUBLEX (SAM)

SPOJ 7258 SUBLEX (SAM)

編輯:C++入門知識

由於SAM能遍歷所有的子串,只要預處理出某個結點的後繼中有多少個不同的子串就可以了。

首先以每個結點為終態算一個子串,所以初始化計數為1。

然後按照拓撲序,用後繼更新pre。類似數DP那種。

SPOJ很卡時,需要各種優化,首先是經典的拓撲。

一次預處理,把後繼全部存好,以為相應的字母。

有個地方需要注意,在查詢的時候,有個k--。因為以當前字母結束也是一個子串,所以需要減掉


[cpp]
#include<iostream>  
#include<cstdio>  
#include<map>  
#include<cstring>  
#include<cmath>  
#include<vector>  
#include<algorithm>  
#include<set>  
#include<string>  
#include<queue>  
#define inf 100000005  
#define M 40  
#define N 200015  
#define maxn 300005  
#define eps 1e-10  
#define zero(a) fabs(a)<eps  
#define Min(a,b) ((a)<(b)?(a):(b))  
#define Max(a,b) ((a)>(b)?(a):(b))  
#define pb(a) push_back(a)  
#define mp(a,b) make_pair(a,b)  
#define mem(a,b) memset(a,b,sizeof(a))  
#define LL unsigned long long  
#define MOD 1000000007  
#define lson step<<1  
#define rson step<<1|1  
#define sqr(a) ((double)(a)*(a))  
#define Key_value ch[ch[root][1]][0]  
#define test puts("OK");  
#pragma comment(linker, "/STACK:1024000000,1024000000")  
using namespace std; 
struct SAM 

    SAM *pre,*son[26]; 
    int len,cnt; 
}*root,*tail,que[N],*b[N]; 
char str[N/2]; 
int son[N][26],cnt[N]; 
int tot=0; 
int c[N]; 
char ch[N][26]; 
void add(int c,int l) 

    SAM *p=tail,*np=&que[tot++]; 
    np->len=l; 
    while(p&&p->son[c]==NULL) p->son[c]=np,p=p->pre; 
    if(p==NULL) np->pre=root; 
    else 
    { 
        SAM *q=p->son[c]; 
        if(p->len+1==q->len) np->pre=q; 
        else 
        { 
            SAM *nq=&que[tot++]; 
            *nq=*q; 
            nq->len=p->len+1; 
            np->pre=q->pre=nq; 
            while(p&&p->son[c]==q) p->son[c]=nq,p=p->pre; 
        } 
    } 
    tail=np; 

void slove(int k) 

    int l=0; 
    int now=0; 
    while(k) 
    { 
        for(int i=0;i<cnt[now];i++) 
        { 
            if(k>que[son[now][i]].cnt) 
            { 
                k-=que[son[now][i]].cnt; 
            } 
            else 
            { 
                str[l++]=ch[now][i]; 
                now=son[now][i]; 
                k--; 
                break; 
            } 
        } 
    } 
    str[l]='\0'; 
    puts(str); 

int main() 

    root=tail=&que[tot++]; 
    scanf("%s",str); 
    int l=strlen(str); 
    for(int i=0;str[i];i++) add(str[i]-'a',i+1); 
    for(int i=0;i<tot;i++) c[que[i].len]++; 
    for(int i=1;i<tot;i++) c[i]+=c[i-1]; 
    for(int i=0;i<tot;i++) b[--c[que[i].len]]=&que[i]; 
    for(int i=0;i<tot;i++) que[i].cnt=1; 
    for(int i=tot-1;i>=0;i--) 
    { 
        SAM *p=b[i]; 
        for(int j=0;j<26;j++) 
        { 
            if(p->son[j]) 
            { 
                int u=p-que,v=p->son[j]-que; 
                son[u][cnt[u]]=v; 
                ch[u][cnt[u]++]=j+'a'; 
                p->cnt+=p->son[j]->cnt; 
            } 
        } 
    } 
    int q,k; 
    scanf("%d",&q); 
    while(q--) 
    { 
        scanf("%d",&k); 
        slove(k); 
    } 
    return 0; 

#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
#include<set>
#include<string>
#include<queue>
#define inf 100000005
#define M 40
#define N 200015
#define maxn 300005
#define eps 1e-10
#define zero(a) fabs(a)<eps
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define pb(a) push_back(a)
#define mp(a,b) make_pair(a,b)
#define mem(a,b) memset(a,b,sizeof(a))
#define LL unsigned long long
#define MOD 1000000007
#define lson step<<1
#define rson step<<1|1
#define sqr(a) ((double)(a)*(a))
#define Key_value ch[ch[root][1]][0]
#define test puts("OK");
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
struct SAM
{
    SAM *pre,*son[26];
    int len,cnt;
}*root,*tail,que[N],*b[N];
char str[N/2];
int son[N][26],cnt[N];
int tot=0;
int c[N];
char ch[N][26];
void add(int c,int l)
{
    SAM *p=tail,*np=&que[tot++];
    np->len=l;
    while(p&&p->son[c]==NULL) p->son[c]=np,p=p->pre;
    if(p==NULL) np->pre=root;
    else
    {
        SAM *q=p->son[c];
        if(p->len+1==q->len) np->pre=q;
        else
        {
            SAM *nq=&que[tot++];
            *nq=*q;
            nq->len=p->len+1;
            np->pre=q->pre=nq;
            while(p&&p->son[c]==q) p->son[c]=nq,p=p->pre;
        }
    }
    tail=np;
}
void slove(int k)
{
    int l=0;
    int now=0;
    while(k)
    {
        for(int i=0;i<cnt[now];i++)
        {
            if(k>que[son[now][i]].cnt)
            {
                k-=que[son[now][i]].cnt;
            }
            else
            {
                str[l++]=ch[now][i];
                now=son[now][i];
                k--;
                break;
            }
        }
    }
    str[l]='\0';
    puts(str);
}
int main()
{
    root=tail=&que[tot++];
    scanf("%s",str);
    int l=strlen(str);
    for(int i=0;str[i];i++) add(str[i]-'a',i+1);
    for(int i=0;i<tot;i++) c[que[i].len]++;
    for(int i=1;i<tot;i++) c[i]+=c[i-1];
    for(int i=0;i<tot;i++) b[--c[que[i].len]]=&que[i];
    for(int i=0;i<tot;i++) que[i].cnt=1;
    for(int i=tot-1;i>=0;i--)
    {
        SAM *p=b[i];
        for(int j=0;j<26;j++)
        {
            if(p->son[j])
            {
                int u=p-que,v=p->son[j]-que;
                son[u][cnt[u]]=v;
                ch[u][cnt[u]++]=j+'a';
                p->cnt+=p->son[j]->cnt;
            }
        }
    }
    int q,k;
    scanf("%d",&q);
    while(q--)
    {
        scanf("%d",&k);
        slove(k);
    }
    return 0;
}

 

  1. 上一頁:
  2. 下一頁:
Copyright © 程式師世界 All Rights Reserved