【CF696D】Legen...(AC自动机)(矩阵快速幂)

题目描述
Barney was hanging out with Nora for a while and now he thinks he may have feelings for her. Barney wants to send her a cheesy text message and wants to make her as happy as possible.

Initially, happiness level of Nora is 0 0 . Nora loves some pickup lines like "I‘m falling for you" and stuff. Totally, she knows n n pickup lines, each consisting only of lowercase English letters, also some of them may be equal (in writing, but different in pronouncing or meaning though). Every time Nora sees i i -th pickup line as a consecutive subsequence of Barney‘s text message her happiness level increases by a_{i} a
i
? . These substrings may overlap, for example, Nora will see the pickup line aa twice and the pickup line ab once in text message aaab.

Due to texting app limits, Barney‘s text may have up to l l characters.

Barney asked you to help him make Nora as much happy as possible, it‘s gonna be legen...

Barney与Nora一起玩的时候,觉得自己喜欢上Nora了。他想给她发一段信息让她开心。

Nora的初始快乐值是0 。Nora喜欢“我深深地爱上你了”这样的情话。她一共知道n句情话,每句仅包含小写英文字母,其中的一些可能相同(写法相同,但读音或意思是不同的)。每次Nora在Barney的信息中看到第i句情话,她的快乐值就会增加a_ia
i
? 。这些情话在信息中可能重叠。例如,在aaab 中,Nora会看到aa 两次,看到ab 一次。

因为短信的长度限制,Barney的短信最长含有lll 个字符。Barney想让你帮他让Nora尽可能开心。

输入输出样例
输入 #1

3 6
3 2 1
heart
earth
art

输出 #1
6
输入 #2

3 6
3 2 8
heart
earth
art

输出 #2
16

这道题可以先把fail树建出来,应为权值是可以重叠的,所以在建fail树的同时也要更新val。

然后可以考虑DP。

\(dp[i][j]表示从fail树上点i到点j所可以带来的贡献\)

显然可以求得\(dp[i][j]\),但是N过大,所以我们考虑用矩阵快速幂加速DP。

然后求从根节点到每个节点的贡献,取最大值即可。

#include<bits/stdc++.h>
#define M 210
using namespace std;
struct data
{
    long long ans[M][M];
}ans1;
int tree[M][26],cnt,val[M],sum[M*26],fail[M*26],m;
long long n,maxn;
char s[M];
queue<int> q;
void work()
{
    for(int i=0;i<=cnt;i++)
    {
        for(int j=0;j<=25;j++)
        {
            ans1.ans[i][tree[i][j]]=sum[tree[i][j]];
        }
    }
}
data operator *(data a,data b)
{
    data c;
    memset(c.ans,-1,sizeof(c.ans));
    for(int k=0;k<=cnt;k++)
    {
        for(int i=0;i<=cnt;i++)
        {
            if(a.ans[i][k]!=-1)
            {
                for(int j=0;j<=cnt;j++)
                {
                    if(b.ans[k][j]!=-1)
                    {
                        c.ans[i][j]=max(c.ans[i][j],a.ans[i][k]+b.ans[k][j]);
                    }
                }
            }
        }
    }
    return c;
}
void insert(int x)
{
    int rt=0,len=strlen(s);
    for(int i=0;i<len;i++)
    {
        if(!tree[rt][s[i]-'a'])
        {
            tree[rt][s[i]-'a']=++cnt;
        }
        rt=tree[rt][s[i]-'a'];
    }
    sum[rt]+=val[x];
}
void getfail()
{
    for(int i=0;i<=25;i++)
    {
        if(tree[0][i])
        {
            q.push(tree[0][i]);
        }
    }
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=0;i<=25;i++)
        {
            int &v=tree[u][i];
            if(!v)
            {
                v=tree[fail[u]][i];
                continue;
            }
            fail[v]=tree[fail[u]][i];
            sum[v]+=sum[fail[v]];
            q.push(v);
        }
    }
}
data fastpow(data a,int b)
{
    data sum=a;
    while(b)
    {
        if(b&1)
        {
            sum=sum*a;
        }
        b>>=1;
        a=a*a;
    }
    return sum;
}
signed main()
{
    memset(ans1.ans,-1,sizeof(ans1.ans));
    scanf("%lld%lld",&m,&n);
    for(int i=1;i<=m;i++)
    {
        scanf("%lld",&val[i]);
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%s",s);
        insert(i);
    }
    getfail();
    work();
    ans1=fastpow(ans1,n-1);
    for(int i=0;i<=cnt;i++)
    {
        maxn=max(maxn,ans1.ans[0][i]);
    }
    printf("%lld\n",maxn);
    return 0;
}

原文地址:https://www.cnblogs.com/2017gdgzoi44/p/11780500.html

时间: 2024-10-17 17:38:00

【CF696D】Legen...(AC自动机)(矩阵快速幂)的相关文章

poj 2778 AC自动机 + 矩阵快速幂

// poj 2778 AC自动机 + 矩阵快速幂 // // 题目链接: // // http://poj.org/problem?id=2778 // // 解题思路: // // 建立AC自动机,确定状态之间的关系,构造出,走一步 // 能到达的状态矩阵,然后进行n次乘法,就可以得到状态间 // 走n步的方法数. // 精髓: // 1):这个ac自动机有一些特别,根节点是为空串,然而 // 每走一步的时候,如果没法走了,这时候,不一定是回到根 // 节点,因为有可能单个的字符时病毒,这样

POJ POJ 2778 DNA Sequence AC自动机 + 矩阵快速幂

首先建立Trie和失败指针,然后你会发现对于每个节点 i 匹配AGCT时只有以下几种情况: i 节点有关于当前字符的儿子节点 j 且安全,则i 到 j找到一条长度为 1的路. i 节点有关于当前字符的儿子节点 j 且 不安全,则i 到 j没有路. i 节点没有关于当前字符的儿子节点 但是能通过失败指针找到一个安全的节点j,那么 i 到 j 找到一条长度为1的路. 关于节点安全的定义: 当前节点不是末节点且当前节点由失败指针指回跟节点的路径上不存在不安全节点,那么这个节点就是安全节点. 然后问题就

hdu 2243 AC自动机 + 矩阵快速幂

// hdu 2243 AC自动机 + 矩阵快速幂 // // 题目大意: // // 给你一些短串,问在长度不超过k的任意串,包含至少一个这些短串的其中 // 一个.问这样的串有多少个. // // 解题思路: // // 首先, 包含和不包含是一种互斥关系,包含+不包含 = 全集u.全集的答案就是 // 26 ^ 1 + 26 ^ 2 + .... + 26 ^ k.不包含的比较好求.构建一个自动机,得到 // 一个转移矩阵A.表示状态i能到状态j的方法数.而这些状态中都是不包含所给的 //

HDU 2243 考研路茫茫――单词情结 (AC自动机 + 矩阵快速幂)

题目链接:考研路茫茫――单词情结 做本题前,个人建议先做一下POJ 2778 http://blog.csdn.net/u013446688/article/details/47378255 POJ2778 是求长度为n,不包含模式串的字符串个数. 而本题是求长度为n,包含模式串的字符串个数.直接用字符串总数减去不包含模式串的字符串个数即为所求. 同样是AC自动机 + 矩阵快速幂.但是还是有所不同的. 因为对2^64取模,所以定义数据类型为unsigned long long就可以了,这样就实现

HDU 2243 考研路茫茫――单词情结 ——(AC自动机+矩阵快速幂)

和前几天做的AC自动机类似. 思路简单但是代码200余行.. 假设solve_sub(i)表示长度为i的不含危险单词的总数. 最终答案为用总数(26^1+26^2+...+26^n)减去(solve_sub(1)+solve(2)+...+solve_sub(n)).前者构造f[i]=f[i-1]*26+26然后矩阵快速幂即可(当然也可以分治的方法).后者即构造出dp矩阵p,然后计算(p^1+p^2+...+p^n),对其分治即可. 代码如下: 1 #include <stdio.h> 2 #

poj2778DNA Sequence (AC自动机+矩阵快速幂)

转载请注明出处: http://www.cnblogs.com/fraud/          ——by fraud DNA Sequence Time Limit: 1000MS   Memory Limit: 65536K Description It's well known that DNA Sequence is a sequence only contains A, C, T and G, and it's very useful to analyze a segment of DN

poj 2778 DNA Sequence(AC自动机+矩阵快速幂)

题目链接:poj 2778 DNA Sequence 题目大意:给定一些含有疾病的DNA序列,现在给定DNA长度,问有多少种不同的DNA序列是健康的. 解题思路:对DNA片段建立AC自动机,因为最多10个串,每个串最长为10,所以最多可能有100个节点,在长度为n时 以每个节点终止的健康字符串个数形成一个状态集,通过AC自动机形成的边可以推导出n+1的状态集,走到单词节点是 非法的,所以同样的我们可以先走到单词节点,但是从单词节点不向后转移.这样可以构造一个矩阵,剩下的就是矩阵 快速幂.注意的一

poj2778 ac自动机+矩阵快速幂

给m个子串,求长度为n的不包含子串的母串数,最直接的应该是暴搜,肯定tle,考虑用ac自动机 将子串建成字典树,通过next表来构造矩阵,然后用矩阵快速幂求长度为n的数量 邻接矩阵https://wenku.baidu.com/view/d7b9787f1711cc7931b716b0.html 对于a(i,j)^k  是指从i到j经过k个点的所有情况数 注意对于End数组,如果某个节点如果fail指针End数组为1,那么End[该节点]也是1 string要开全局变量,不然不能运行= = #i

POJ 2778 DNA Sequence (AC自动机 + 矩阵快速幂)

题目链接:DNA Sequence 解析:AC自动机 + 矩阵加速(快速幂). 这个时候AC自动机 的一种状态转移图的思路就很透彻了,AC自动机就是可以确定状态的转移. AC代码: #include <iostream> #include <cstdio> #include <queue> #include <cstring> using namespace std; const int MOD = 100000; struct Matrix{ int ma

HDU 2243考研路茫茫——单词情结 (AC自动机+矩阵快速幂)

背单词,始终是复习英语的重要环节.在荒废了3年大学生涯后,Lele也终于要开始背单词了. 一天,Lele在某本单词书上看到了一个根据词根来背单词的方法.比如"ab",放在单词前一般表示"相反,变坏,离去"等. 于是Lele想,如果背了N个词根,那这些词根到底会不会在单词里出现呢.更确切的描述是:长度不超过L,只由小写字母组成的,至少包含一个词根的单词,一共可能有多少个呢?这里就不考虑单词是否有实际意义. 比如一共有2个词根 aa 和 ab ,则可能存在104个长度不