hud2243 考研路茫茫——单词情结
考研路茫茫——单词情结
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Problem Description
背单词,始终是复习英语的重要环节。在荒废了3年大学生涯后,Lele也终于要开始背单词了。
一天,Lele在某本单词书上看到了一个根据词根来背单词的方法。比如"ab",放在单词前一般表示"相反,变坏,离去"等。
于是Lele想,如果背了N个词根,那这些词根到底会不会在单词里出现呢。更确切的描述是:长度不超过L,只由小写字母组成的,至少包含一个词根的单词,一共可能有多少个呢?这里就不考虑单词是否有实际意义。
比如一共有2个词根 aa 和 ab ,则可能存在104个长度不超过3的单词,分别为
(2个) aa,ab,
(26个)aaa,aab,aac...aaz,
(26个)aba,abb,abc...abz,
(25个)baa,caa,daa...zaa,
(25个)bab,cab,dab...zab。
这个只是很小的情况。而对于其他复杂点的情况,Lele实在是数不出来了,现在就请你帮帮他。
Input
本题目包含多组数据,请处理到文件结束。
每组数据占两行。
第一行有两个正整数N和L。(0<N<6,0<L<2^31)
第二行有N个词根,每个词根仅由小写字母组成,长度不超过5。两个词根中间用一个空格分隔开。
Output
对于每组数据,请在一行里输出一共可能的单词数目。
由于结果可能非常巨大,你只需要输出单词总数模2^64的值。
Sample Input
2 3
aa ab
1 2
a
Sample Output
104
52
ac自动机上dp。。。我们才用反向思考,求出没有出现的个数。
dp[i][j] 表示节点 i ,走了 j 步,没有出现的个数。
用矩阵优化 j 这一维就好啦~
但是我们要求了是j 1~L 的前缀和。
所以还有一个矩阵
a 1\\
0 1 \\
\end{matrix}\right]
\]
这个就是求前缀和。。。如果是矩阵前缀和就把a换成矩阵,1换成等大小的单位矩阵就好了。
#include<bits/stdc++.h>
#define LL unsigned long long
using namespace std;
const int maxn = 120;
struct lpl{
int n, m;
LL a[maxn][maxn];
}op, bas, ini;
struct ld{int son[30]; int fail; bool f;}node[maxn];
bool vis[30];
int N, L, cnt;
LL all, ans;
char s[maxn];
inline lpl operator * (lpl &A, lpl &B){
lpl ret; ret.n = A.n; ret.m = B.m;
for(int i = 1; i <= ret.n; ++i)
for(int j = 1; j <= ret.m; ++j){
ret.a[i][j] = 0;
for(int k = 1; k <= ret.n; ++k) ret.a[i][j] += A.a[i][k] * B.a[k][j];
}
return ret;
}
inline lpl Fpow(lpl A, int t){
lpl ret; lpl tmp = A; ret.n = ret.m = A.n;
memset(ret.a, 0, sizeof(ret.a)); memset(ini.a, 0, sizeof(ini.a));
for(int i = 1; i <= ret.n; ++i) ret.a[i][i] = 1;
while(t){
if(t & 1) ret = ret * tmp;
tmp = tmp * tmp; t >>= 1;
}
return ret;
}
inline void putit(){
int len, qwe, now; cnt = 1; memset(node, 0, sizeof(node));
for(int i = 1; i <= N; ++i){
scanf("%s", s + 1); len = strlen(s + 1); now = 1;
for(int j = 1; j <= len; ++j){
qwe = s[j] - 'a' + 1;
if(!node[now].son[qwe]) node[now].son[qwe] = ++cnt;
now = node[now].son[qwe]; if(j == len) node[now].f = true;
}
}
lpl lin; lin.n = lin.m = 2; memset(lin.a, 0, sizeof(lin.a));
lin.a[1][1] = 26; lin.a[1][2] = lin.a[2][2] = 1;
lin = Fpow(lin, L); all = lin.a[1][1] + lin.a[1][2] - 1;
}
inline void AC(){
queue<int> q; queue<int> qqq;
q.push(1); qqq.push(1);
while(!q.empty()){
int now = q.front(); q.pop();
for(int i = 1; i <= 26; ++i){
if(!node[now].son[i]) continue;
q.push(node[now].son[i]); qqq.push(node[now].son[i]);
if(now == 1){node[node[now].son[i]].fail = 1; continue;}
int t = node[now].fail;
while(!node[t].son[i] && t) t = node[t].fail;
if(!t) node[node[now].son[i]].fail = 1;
else node[node[now].son[i]].fail = node[t].son[i];
}
}
while(!qqq.empty()){
int now = qqq.front(); qqq.pop();
int qwe = now;
while(qwe){
if(node[qwe].f){node[now].f = true; break;}
qwe = node[qwe].fail;
}
}
}
inline void workk(){
op.n = cnt, op.m = 1; memset(op.a, 0, sizeof(op.a)); op.a[1][1] = 1;
memset(bas.a, 0, sizeof(bas.a)); bas.n = bas.m = cnt;
for(int i = 1; i <= cnt; ++i){
int now = i, lft = 26; memset(vis, false, sizeof(vis));
while(now){
for(int j = 1; j <= 26; ++j){
if(!node[now].son[j]) continue;
if(!vis[j]){
vis[j] = true; lft--;
bas.a[node[now].son[j]][i] = 1;
}
}
now = node[now].fail;
}
bas.a[1][i] = lft;
}
for(int i = 1; i <= cnt; ++i)
if(node[i].f)
for(int j = 1; j <= cnt; ++j) bas.a[i][j] = 0;
ini.n = ini.m = cnt * 2;
for(int i = 1; i <= cnt; ++i){
ini.a[i][i + cnt] = ini.a[i + cnt][i + cnt] = 1;
for(int j = 1; j <= cnt; ++j) ini.a[i][j] = bas.a[i][j];
}
ini = Fpow(ini, L);
for(int i = 1; i <= cnt; ++i)
for(int j = 1; j <= cnt; ++j){
bas.a[i][j] = ini.a[i][j] + ini.a[i][j + cnt];
if(i == j) bas.a[i][j]--;
}
op = bas * op; ans = 0;
for(int i = 1; i <= cnt; ++i) ans += op.a[i][1];
ans = all - ans;
cout << ans << endl;
}
int main()
{
while(~scanf("%d%d", &N, &L)){
putit();
AC();
workk();
}
return 0;
}
hud2243 考研路茫茫——单词情结的更多相关文章
- HDU 2243 考研路茫茫——单词情结(AC自动机+矩阵)
考研路茫茫——单词情结 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- hdu_2243_考研路茫茫——单词情结(AC自动机+矩阵)
题目链接:hdu_2243_考研路茫茫——单词情结 题意: 让你求包含这些模式串并且长度不小于L的单词种类 题解: 这题是poj2788的升级版,没做过的强烈建议先做那题. 我们用poj2778的方法 ...
- hdu 2243 考研路茫茫——单词情结(AC自动+矩阵)
考研路茫茫——单词情结 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- hdu2243 考研路茫茫——单词情结【AC自动机】【矩阵快速幂】
考研路茫茫——单词情结 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- 考研路茫茫――单词情结 HDU - 2243(ac自动机 + 矩阵快速幂)
考研路茫茫——单词情结 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- HDU2243 考研路茫茫——单词情结 ——AC自动机、矩阵优化
题目链接:https://vjudge.net/problem/HDU-2243 考研路茫茫——单词情结 Time Limit: 2000/1000 MS (Java/Others) Memor ...
- HDU 2243 考研路茫茫——单词情结
考研路茫茫——单词情结 Time Limit: 1000ms Memory Limit: 32768KB This problem will be judged on HDU. Original ID ...
- HDU-2243 考研路茫茫——单词情结(AC自动机)
题目大意:给n个单词,长度不超过L的单词有多少个包含n个单词中的至少一个单词. 题目分析:用长度不超过L的单词书目减去长度在L之内所有不包含任何一个单词的书目. 代码如下: # include< ...
- HDU2243 考研路茫茫——单词情结(AC自动机+矩阵快速幂)
与POJ2778一样.这题是求长度不超过n且包含至少一个词根的单词总数. 长度不超过n的单词总数记为Sn,长度不超过n不包含词根的单词总数记为Tn. 答案就是,Sn-Tn. Sn=26+262+263 ...
随机推荐
- 1145. Hashing - Average Search Time (25)
The task of this problem is simple: insert a sequence of distinct positive integers into a hash tabl ...
- mysql数据精度丢失问题深入探讨
不要盲目的说float和double精度可能发生丢失,而是说在存取时因为精度不一致会发生丢失,当然这里的丢失指的是扩展或者截断了,丢失了原有的精度.decimal是好,但不是说不会发生任何精度丢失.如 ...
- Flsak中的socket是基于werkzeug实现的。
from werkzeug.serving import run_simple from werkzeug.wrappers import Request,Respinse @Request.appl ...
- Git关联JIRA的issue
指导文章 http://www.51testing.com/html/30/n-3724930.html http://{$host_url}/help/user/project/integratio ...
- vue对特殊特性的研究
key 预期:number | string key 的特殊属性主要用在 Vue 的虚拟 DOM 算法,在新旧 nodes 对比时辨识 VNodes.如果不使用 key,Vue 会使用一种最大限度减少 ...
- Telegraf根据配置文件启动(Influxdb的数据收集)
1.创建一个telegraf.config文件 telegraf -sample-config -input-filter cpu:disk:diskio:net:system:mem -output ...
- 浅谈MySQL优化
本文整理了一些MySQL的通用优化方法,做个简单的总结分享,旨在帮助那些没有专职MySQL DBA的企业做好基本的优化工作,至于具体的SQL优化,大部分通过加适当的索引即可达到效果,更复杂的就需要具体 ...
- JS自定义随机键盘
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- SpringBoot怎么访问html文件
pom.xml <dependency> <groupId>org.springframework.boot</groupId> <artifactId> ...
- Python_010(迭代器)
一.函数名的运用 1.函数名的内存地址 def func(): print("英雄联盟") print(func) #输出结果: <function func at 0x00 ...