poj-1625 Censored![ac自动机+dp+高精度]

题目地址
先把病毒串丢进ac自动机里面。
dp[i][j]表示长度为i的从trie图的根节点到j满足条件的串的数量。
因为答案很大，要用到高精度。
注意的是建图的时候，ed[u] |= ed[f[u]]，如果失配指针的位置是病毒的话那么u也不能匹配。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn = 1e5+5;
struct BigInter {private:int num[50], len;int base;
public:     BigInter() {memset(num,0,sizeof(num));len = 1;base = 10000;}BigInter(int v) {memset(num,0,sizeof(num));base = 10000;len = 0;do {num[++len] = v%base;v /= base;}while(v);}BigInter operator=(int x) {return *this = BigInter(x);}BigInter operator=(long long x) {return *this = BigInter(x);}BigInter operator+(BigInter &b) {BigInter res;res.len = max(this->len, b.len);for(int i = 1; i <= res.len; i++) {res.num[i] += this->num[i] + b.num[i];res.num[i+1] += res.num[i]/base;res.num[i] = res.num[i]%base;}if(res.num[res.len+1]) res.len++;return res;}BigInter operator+(long long x) {BigInter res;BigInter b = BigInter(x);res.len = max(this->len, b.len);for(int i = 1; i <= res.len; i++) {res.num[i] += this->num[i] + b.num[i];res.num[i+1] += res.num[i]/base;res.num[i] = res.num[i]%base;}if(res.num[res.len+1]) res.len++;return res;}BigInter operator*(BigInter &b) {BigInter res;res.len = this->len+b.len;for(int i = 1; i <= this->len; i++) {for(int j = 1; j <= b.len; j++) {res.num[j+i-1] += this->num[i] * b.num[j];res.num[j+i] += res.num[j+i-1]/base;res.num[j+i-1] %= base;}}while(!res.num[res.len]) res.len--;return res;}BigInter operator*(long long x) {BigInter res;BigInter b(x);res.len = this->len+b.len;for(int i = 1; i <= this->len; i++) {for(int j = 1; j <= b.len; j++) {res.num[j+i-1] += this->num[i] * b.num[j];res.num[j+i] += res.num[j+i-1]/base;res.num[j+i-1] %= base;}}while(!res.num[res.len]) res.len--;return res;}void operator +=(BigInter &b) {*this = *this + b;}void operator *=(BigInter &b) {*this = *this * b;}void operator +=(long long x) {BigInter b = BigInter(x);*this = *this + b;}void operator *=(long long x) {BigInter b = BigInter(x);*this = *this * b;}BigInter operator++() {*this = *this+1;return *this;}BigInter operator ++(int) {BigInter old = *this;++(*this);return old;}void print() {printf("%d", num[len]);for(int i = len-1; i >= 1; i--) {printf("%04d", this->num[i]);}puts("");}
}dp[105][305];
int indx[5005];
struct node{int nex[maxn][305], tot, root;int f[maxn], ed[maxn];int newnode() {for(int i = 0; i < 300; i++) {nex[tot][i] = -1;}ed[tot] = 0;return tot++;}void init() {tot = 0;root = newnode();}void insert(char *s) {int len = strlen(s), u = root;for(int i = 0; i < len; i++) {int ch = indx[s[i]+100];if(nex[u][ch] == -1) nex[u][ch] = newnode();u = nex[u][ch];}ed[u] = 1;}void getfail() {queue<int>Q;f[root] = root;for(int i = 0; i < 300; i++) {if(nex[root][i] == -1) nex[root][i] = root;else {f[nex[root][i]] = root;Q.push(nex[root][i]);}}while(!Q.empty()) {int u = Q.front();Q.pop();ed[u] |= ed[f[u]];for(int i = 0; i < 300; i++) {if(nex[u][i] == -1) nex[u][i] = nex[f[u]][i];else {f[nex[u][i]] = nex[f[u]][i];Q.push(nex[u][i]);}}}}void query(int m, int n) {dp[0][0] = 1;for(int i = 1; i <= m; i++) {for(int j = 0; j < tot; j++) {if(ed[j] == 1) continue;for(int k = 0; k < n; k++) {if(ed[nex[j][k]] != 1) {dp[i][nex[j][k]] += dp[i-1][j];}}}}BigInter res;for(int i = 0; i < tot; i++) res = res + dp[m][i];res.print();}
}ac;
char ss[maxn];
int main() {int n, m, p;scanf("%d%d%d", &n, &m, &p);scanf("%s", ss);ac.init();int len = strlen(ss);for(int i = 0; i < len; i++) indx[ss[i]+100] = i;for(int i = 1; i <= p; i++) {scanf("%s", ss);ac.insert(ss);}ac.getfail();ac.query(m, n);
}

poj-1625 Censored![ac自动机+dp+高精度]相关推荐

POJ - 1625 Censored!(AC自动机+dp+高精度运算)
题目链接:点击查看题目大意:给出一个含有 n 个不同字符的字符集,接着规定所有单词的长度为 m ,再给出 k 个病毒串,问有多少个字符串中不含有病毒串题目分析:这个题目和之前做过的DNA的那个题有 ...
Censored! POJ - 1625（AC自动机 + dp +高精度模板）
题目链接题目大意:给你一个字母表,给定一些敏感字符串,问长度为m且不含任意敏感字符串的串有多少个.(字符全部来自字母表) 思路:首先第一个坑点是输入的字符是unsigned char,可能出现负的A ...
POJ 1625 Censored ( Trie图 DP 高精度 )
题意 : 给出 n 个单词组成的字符集以及 p 个非法串,问你用字符集里面的单词构造长度为 m 的单词的方案数有多少种? 分析 : 与 POJ 2778 非常相似的一道题目,如果没有做过就尝试去了解 ...
[POJ 1625] Censored! （AC自动机+DP+高精度）
链接 POJ 1625 题意给出P个模式串,问长度为M且不含有P中任何一个为子串的字符串有多少种. 给出了大小为N的一个字符集,属于ASCII但不一定为英文字母. 最终答案不进行取模,所以可能非常大 ...
Censored! :ac自动机 + DP
分析如果字符串的长度大一点可以用矩阵快速幂来做,因为数据范围比较小,可以用DP来代替我们去枚举每一位在a自动机中的位置,再去枚举下一位的位置,判断是否打过标记,然后进行状态转移需要注意的是,如果 ...
uvalive4842(AC自动机+DP)
题意: 给出猴子打字时打某个字母的概率,猴子最多可以敲键盘m次,问得到的长度是m的单词包含模式串的概率. 思路: AC自动机+dp. 首先,我们用模式串构造一个AC自动机,用dp[i]][j]表示当前 ...
HDU 2296 Ring AC自动机 + DP
题意:给你n个模式串,每个模式串有一个得分,让你构造出一个长度为N之内且分数最高的文本串;输出字典序列最小的. 解题思路: AC自动机 + DP , 不过要输出字典序列最小,多开一个一个三维字符串 ...
bzoj 1030: [JSOI2007]文本生成器（AC自动机+DP）
1030: [JSOI2007]文本生成器 Time Limit: 1 Sec Memory Limit: 162 MB Submit: 5187 Solved: 2136 [Submit][St ...
【BZOJ】4861: [Beijing2017]魔法咒语 AC自动机+DP+矩阵快速幂
[题意]给定n个原串和m个禁忌串,要求用原串集合能拼出的不含禁忌串且长度为L的串的数量.(60%)n,m<=50,L<=100.(40%)原串长度为1或2,L<=10^18. [算法 ...

poj-1625 Censored![ac自动机+dp+高精度]

poj-1625 Censored![ac自动机+dp+高精度]相关推荐

最新文章

热门文章