题目链接:[JSOI2009]密码

我们先看第一问:输出方案数

我们把所有给出来的串丢到AC自动机里面去,然后在建出来的\(trie\)图上跑dp

由于\(n\leq 10\)我们很自然的就想到了状压

记\(dp[i][j][sta]\)表示原串匹配到了第\(i\)位,在AC自动机里走到了第\(j\)个节点,已经出现了\(sta\)个单词(压缩状态)的方案数

注意到如果有两个串\(s_i\)和\(s_j\)满足\(s_i\)是\(s_j\)的子串,那么我们所构建的串并不一定必须要有\(s_i\),因为如果有了\(s_j\)那么一定满足了必须有\(s_i\)

那么我们可以确定最终统计答案时的\(sta\),暴力累加即可

接下来就是第二问:输出方案

我么先要考虑一下在什么情况需要做第二问,题目中给的条件是\(ans\leq 42\)

如果此时我们构建的字符串一定有一位“空闲”,即一定有一位可以放置任意字符,而剩下的则必须构成给定的字符串,那么方案数就是\(2*26=52>42\)(该字符可以放在开头和末尾)

因此此时一定不需要输出答案

所以当要进行第二问的时候,一定是只能用所给定的串进行组合

并且不可能出现的重复的串,因为重复的串是可以变成空闲节点的,这样就和上面一样了

预处理出对于任意两个字符串\(s_i,s_j\),\(s_i\)的后\(k\)位和\(s_j\)的前\(k\)位相等的所有\(k\)值

然后暴力dfs排序顺序即可

#include<iostream>
#include<string.h>
#include<string>
#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<vector>
#include<queue>
#include<map>
using namespace std;
const int maxd=1000000007,N=100000;
const double pi=acos(-1.0);
typedef long long ll; namespace AC{
struct node{
int ch[27];
int fail,end;
node() {memset(ch,0,sizeof(ch));fail=0;end=0;}
}point[100100];
int cnt; void init()
{
cnt=0;
point[0].fail=0;
} void insert(char *s,int id)
{
int len=strlen(s),i,now=0;
for (i=0;i<len;i++)
{
if (!point[now].ch[s[i]-'a'])
{
point[now].ch[s[i]-'a']=(++cnt);
point[cnt]=node();
}
now=point[now].ch[s[i]-'a'];
}
point[now].end=(1<<(id-1));
} void get_fail()
{
queue<int> q;
int i;
for (i=0;i<26;i++)
{
if (point[0].ch[i])
{
point[point[0].ch[i]].fail=0;
q.push(point[0].ch[i]);
}
}
while (!q.empty())
{
int u=q.front();q.pop();
for (i=0;i<26;i++)
{
if (point[u].ch[i])
{
point[point[u].ch[i]].fail=point[point[u].fail].ch[i];
q.push(point[u].ch[i]);
}
else point[u].ch[i]=point[point[u].fail].ch[i];
}
}
}
}
using namespace AC;
int n,m,len[20],restm,LINK[50][50],ansord[50];
ll dp[2][110][1030];
char s[20][20];
bool del[20]; int read()
{
int x=0,f=1;char ch=getchar();
while ((ch<'0') || (ch>'9')) {if (ch=='-') f=-1;ch=getchar();}
while ((ch>='0') && (ch<='9')) {x=x*10+(ch-'0');ch=getchar();}
return x*f;
} int getl(int u,int v)
{
int maxl,i;
for (maxl=min(len[u],len[v]);maxl>0;maxl--)
{
bool ok=1;
for (i=0;i<maxl;i++)
{
if (s[v][i]!=s[u][len[u]-maxl+i]) {ok=0;break;}
}
if (ok) return maxl;
}
return 0;
} int ord[50],tot=0;
bool vis[50];
char ansch[110][110]; void dfs(int dep)
{
if (dep>restm)
{
tot++;
int i,j,now=0;
for (i=1;i<=restm;i++)
{
for (j=LINK[ord[i-1]][ord[i]];j<len[ord[i]];j++)
{
ansch[tot][now++]=s[ord[i]][j];
if (now>n) {tot--;return;}
}
}
if (now<n) tot--;
return;
} int i;
for (i=1;i<=m;i++)
{
if ((del[i]) || (vis[i])) continue;
ord[dep]=i;vis[i]=1;
dfs(dep+1);
vis[i]=0;
}
} bool cmp(int p,int q)
{
int i;
for (i=0;i<n;i++)
if (ansch[p][i]!=ansch[q][i]) return ansch[p][i]<ansch[q][i];
} int main()
{
n=read();m=read();restm=m;
int i,j,k;
for (i=1;i<=m;i++) {scanf("%s",s[i]);len[i]=strlen(s[i]);}
for (i=1;i<=m;i++)
{
for (j=1;j<=m;j++)
{
if ((del[i]) || (del[j]) || (i==j)) continue;
if (len[i]<len[j]) continue;
if (strstr(s[i],s[j])) {del[j]=1;restm--;break;}
}
}
init();
for (i=1;i<=m;i++)
if (!del[i]) insert(s[i],i);
get_fail();
dp[0][0][0]=1;
int now=1,pre=0;
for (i=1;i<=n;i++)
{
memset(dp[now],0,sizeof(dp[now]));
for (j=0;j<=cnt;j++)
{
for (k=0;k<(1<<m);k++)
{
if (dp[pre][j][k])
{
int p;
for (p=0;p<26;p++)
{
int v=point[j].ch[p];
dp[now][v][(k|point[v].end)]+=dp[pre][j][k];
}
}
}
}
now^=1;pre^=1;
}
int all=0;
for (i=1;i<=m;i++)
if (!del[i]) all|=(1<<(i-1));
ll ans=0;
for (i=0;i<=cnt;i++) ans+=dp[pre][i][all];
printf("%lld\n",ans);
if (ans<=42)
{
for (i=1;i<=m;i++)
{
for (j=1;j<=m;j++)
{
if ((del[i]) || (del[j])) continue;
LINK[i][j]=getl(i,j);
}
}
memset(vis,0,sizeof(vis));
dfs(1);
for (i=1;i<=ans;i++) ansord[i]=i;
sort(ansord+1,ansord+1+tot,cmp);
for (i=1;i<=tot;i++)
{
for (j=0;j<n;j++) putchar(ansch[ansord[i]][j]);
puts("");
}
}
return 0;
}

bzoj

bzoj1559 [JSOI2009]密码的更多相关文章

  1. BZOJ1559[JSOI2009]密码——AC自动机+DP+搜索

    题目描述 输入 输出 样例输入 10 2 hello world 样例输出 2 helloworld worldhello 提示 这题算是一个套路题了,多个串求都包含它们的长为L的串的方案数. 显然是 ...

  2. [BZOJ1559][JSOI2009]密码(AC自动机)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1559 2009年的省选题虽然比起现在简单了不少,但对我来说还是很有挑战性的. 首先对于这种多串匹配问 ...

  3. BZOJ1559 [JSOI2009]密码 【AC自动机 + 状压dp】

    题目链接 BZOJ1559 题解 考虑到这是一个包含子串的问题,而且子串非常少,我们考虑\(AC\)自动机上的状压\(dp\) 设\(f[i][j][s]\)表示长度为\(i\)的串,匹配到了\(AC ...

  4. 【BZOJ1559】[JSOI2009]密码(AC自动机,动态规划,搜索)

    [BZOJ1559][JSOI2009]密码(AC自动机,动态规划,搜索) 题面 BZOJ 洛谷 题解 首先求方案数显然是构建\(AC\)自动机之后再状压\(dp\),似乎没有什么好讲的. 现在考虑答 ...

  5. BZOJ 1559: [JSOI2009]密码( AC自动机 + 状压dp )

    建AC自动机后, dp(x, y, s)表示当前长度为x, 在结点y, 包括的串的状态为s的方案数, 转移就在自动机上走就行了. 对于输出方案, 必定是由给出的串组成(因为<=42), 所以直接 ...

  6. P4045 [JSOI2009]密码

    题目 P4045 [JSOI2009]密码 做法 AC自动机+状压+爆搜 建AC自动机是显然的,顺便预处理\(lst_i\)表示\(i\)结点以哪些串结束(二进制) 然后跑状压\(dp[i][j][k ...

  7. [JSOI2009]密码 [AC自动机]

    题面 bzoj luogu 首先看到这题就知道随便暴枚 只要是多项式算法都能过 先常规建AC自动机 注意被别的单词包含的单词没有存在的价值 剩余单词状压 大力dp f[长度][节点编号][状态] \( ...

  8. JSOI2009 密码 和 JSOI2007 文本生成器 和 ZOJ3545 Rescue the Rabbit

    密码 众所周知,密码在信息领域起到了不可估量的作用.对于普通的登陆口令,唯一的破解 方法就是暴力破解一逐个尝试所有可能的字母组合,但这是一项很耗时又容易被发现的工 作.所以,为了获取对方的登陆口令,在 ...

  9. [BZOJ 1559] [JSOI2009] 密码 【AC自动机DP】

    题目链接:BZOJ - 1559 题目分析 将给定的串建成AC自动机,然后在AC自动机上状压DP. 转移边就是Father -> Son 或 Now -> Fail. f[i][j][k] ...

随机推荐

  1. 423. Reconstruct Original Digits from English(Medium)

    Given a non-empty string containing an out-of-order English representation of digits 0-9, output the ...

  2. Oja’s rule

    目录 Oja's rule 背景 Hebbian learning 主要的一些理论 论文里面一些主要的假设 引理1 引理2 引理3 定理1 LEMMA 3(ALL) 引理 4 定理 2 定理 3(关于 ...

  3. pycharm异常问题之Unable to save settings: Failed to save settings. Please restart PyCharm

    pycharm异常之Unable to save settings: Failed to save settings. Please restart PyCharm 今天一不小心将电脑关了,但是关机之 ...

  4. 字符串和ASCII之间的转换

    public class CharToAscii { public static void main(String[] args) { CharToAscii.AscToString(); CharT ...

  5. mysql实现首字母从A-Z排序

    1.常规排序ASC DESC ASC 正序 DESC倒叙 -- 此处不用多讲 2.自定义排序 自定义排序是根据自己想要的特定字符串(数字)顺序进行排序.主要是使用函数 FIELD(str,str1,s ...

  6. AngularJS:directive自定义的指令

    除了 AngularJS 内置的指令外,我们还可以创建自定义指令. 你可以使用 .directive 函数来添加自定义的指令. 要调用自定义指令,HTML 元素上需要添加自定义指令名. 使用驼峰法来命 ...

  7. 你不知道的JavaScript——第一章:作用域是什么?

    编译原理 JavaScript事实上是一门编译语言,但与传统的编译语言不同,它不是提前编译的,编译结果也不能在分布式系统中进行移植. 任何JavaScript代码片段在执行前都要进行编译(通常就在执行 ...

  8. AssemblyScript的测试

    详细文档介绍 export function f(x: i32): i32 { if (x === 1 || x === 2) { return 1; } return f(x - 1) + f(x ...

  9. scrapy几种反反爬策略

    一.浏览器代理 1.直接处理: 1.1在setting中配置浏览器的各类代理: user_agent_list=[ "Mozilla/5.0 (Windows NT 10.0; Win64; ...

  10. Arrays类

    package cn.zhou.com; import java.util.Arrays; /* * Arrays 类 ? * * 数组工具类 操作数组! * * Arrays的方法? * * */ ...