首页 > 其他 > 详细

HNU 13108-Just Another Knapsack Problem (ac自动机上的dp)

时间:2015-09-03 23:10:32      阅读:228      评论:0      收藏:0      [点我收藏+]

题意:

给你一个母串,多个模式串及其价值,求用模式串拼接成母串(不重叠不遗漏),能获得的最大价值。

分析:

ac自动机中,在字典树上查找时,用dp,dp[i]拼成母串以i为结尾的子串,获得的最大价值,dp[i]=max(dp[i],dp[i-len]+val[tmp])。,len是模式串的长度,val[tmp]为其价值。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
#include <algorithm>
#include<iostream>
using namespace std;
const int maxnode = 300010;
const int sigma_size = 26;
struct Trie{
    int dp[100010];
    int ch[maxnode][sigma_size];
    int val[maxnode];     //该单词在模式串中出现的次数
    int last[maxnode];
    int f[maxnode];         //失配数组
    int num[maxnode];     //该单词出现在文本串的次数
    int pre[maxnode];    //该单词的前驱
    int len[maxnode];    //以该单词结尾的单词长度
    int Char[maxnode];     //该单词对应的字母
    int road[maxnode];   //路径压缩优化 针对计算模式串出现的种数
    int sz;
    int Newnode()
    {
        val[sz] = f[sz] = last[sz] = len[sz] = num[sz] = 0;
        memset(ch[sz], 0, sizeof ch[sz]);
        return sz++;
    }
    void init(){
        sz=0;
        Newnode();
    }
    int idx(char c){ return c-a; }
    void build(char *s,int v){
        int u = 0;
        for(int i = 0, c; s[i] ;i++){
            c = idx(s[i]);
            if(!ch[u][c])
                ch[u][c] = Newnode();
            pre[ch[u][c]] = u;
            Char[ch[u][c]] = s[i];
            len[ch[u][c]] = len[u]+1;
            road[ch[u][c]] = 1;
            u = ch[u][c];
        }
        val[u] = max(val[u],v);
        //num[u] = 0;
        //return u;
    }
   void getFail(){
        queue<int> q;
        for(int i = 0; i<sigma_size; i++)
            if(ch[0][i]) q.push(ch[0][i]);
        int r, c, u, v;
        while(!q.empty()){
            r = q.front(); q.pop();
            for(c = 0; c<sigma_size; c++){
                u = ch[r][c];
                if(!u)continue;
                q.push(u);
                v = f[r];
                while(v && ch[v][c] == 0) v = f[v]; //沿失配边走上去 如果失配后有节点 且 其子节点c存在则结束循环
                f[u] = ch[v][c];
            }
        }
    }
    void find(char *T){
    //计算模式串出现的个数:(每种多次出现算多次)
        int j = 0,l=strlen(T);
        for(int i = 0, c, temp; i<l ; i++){
                dp[i]=0;
            c = idx(T[i]);
            while(j && ch[j][c]==0) j = f[j];
            j = ch[j][c];

            temp = j;
            while(temp){
                if(val[temp]){
                    if(i-len[temp]<0)
                        dp[i]=max(dp[i],val[temp]);
                    else if(dp[i-len[temp]])
                        dp[i]=max(dp[i],dp[i-len[temp]]+val[temp]);
                }
                temp = f[temp];
            }
        }
        printf("%d\n",dp[l-1]);
    }
    int find_kind(char *T){
    //计算种数, 重复出现的不再计算(若多个询问则要在此处加for(i=0->sz)lu[i]=1;
        int j = 0, i, c, temp,ans=0;
        for(i = 0; T[i]; i++){
            c = idx(T[i]);
            while(j && ch[j][c] == 0) j = f[j];
            j = ch[j][c];
            temp = j;
            while(temp && road[temp]){
                if(val[temp])
                {
                    ++ans;
                    val[temp] = 0;
                }
                road[temp] = 0;
                temp = f[temp];
            }
        }
        return ans;
    }
}ac;
int n,v;
char T[100100],P[310];
int main() {
    while(~scanf("%s",T)){
            ac.init();
        scanf("%d",&n);
        while(n--){
            scanf("%s%d",P,&v);
            ac.build(P,v);
        }
        ac.getFail();
        ac.find(T);
    }
}

 

HNU 13108-Just Another Knapsack Problem (ac自动机上的dp)

原文:http://www.cnblogs.com/zsf123/p/4780904.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!