题意:给你只包含‘A’,‘G’,‘T’,‘C’四个字母的n个模板串和1个文本串,问你文本串改变多少个字符就可以使得文本串中没有一个模板串

解题思路:

我们可以知道  dp[i][j] 为文本串到 第i 个字符  AC自动机状态为j的最少的变换次数(这里为什么要用AC自动机,因为end数组可以记录哪一个状态是结束的,而且处理以后可以知道那些后缀等于前缀--也就是不能到达,因为如果能够到达的话那么状态更新就会产生错误。),这样dp即可

解题代码:

  1 // File Name: temp.cpp
  2 // Author: darkdream
  3 // Created Time: 2014年09月11日 星期四 15时18分4秒
  4
  5 #include<vector>
  6 #include<list>
  7 #include<map>
  8 #include<set>
  9 #include<deque>
 10 #include<stack>
 11 #include<bitset>
 12 #include<algorithm>
 13 #include<functional>
 14 #include<numeric>
 15 #include<utility>
 16 #include<sstream>
 17 #include<iostream>
 18 #include<iomanip>
 19 #include<cstdio>
 20 #include<cmath>
 21 #include<cstdlib>
 22 #include<cstring>
 23 #include<ctime>
 24 #include<queue>
 25 #define LL long long
 26 #define maxn 1005
 27
 28 using namespace std;
 29 int dp[maxn][maxn];
 30 struct Trie
 31 {
 32     int next[maxn][4],fail[maxn],end[maxn];
 33     int root, L;
 34     int newnode()
 35     {
 36         memset(next[L],-1,sizeof(next[L]));
 37         end[L++] = 0 ;
 38         return L-1;
 39     }
 40     void init()
 41     {
 42         L = 0 ;
 43         root = newnode();
 44     }
 45     void insert(int buf[],int len)
 46     {
 47         int now = root;
 48         for(int i = 0 ;i < len ;i ++)
 49         {
 50             if(next[now][buf[i]] ==  -1)
 51             {
 52                 next[now][buf[i]] = newnode();
 53             }
 54             now = next[now][buf[i]];
 55         }
 56         end[now] = 1;
 57     }
 58     void build()
 59     {
 60         queue<int> Q;
 61         fail[root] = root;
 62         for(int i = 0;i < 4;i ++)
 63         {
 64             if(next[root][i] == -1)
 65             {
 66                 next[root][i] = root;  //指向root 是没有问题的,我们主要是通过end 数组对个数进行计数的。
 67             }else{
 68                 fail[next[root][i]] = root;
 69                 Q.push(next[root][i]);
 70             }
 71         }
 72         while(!Q.empty())
 73         {
 74             int now = Q.front();
 75             Q.pop();
 76             for(int i = 0 ;i < 4;i ++)
 77             {
 78                 if(next[now][i] == -1)
 79                 {
 80                     next[now][i] =  next[fail[now]][i];
 81                 }
 82                 else{
 83                     fail[next[now][i]] = next[fail[now]][i];
 84                     if(end[next[fail[now]][i]] == 1)
 85                         end[next[now][i]] = 1;
 86                     Q.push(next[now][i]);
 87                 }
 88             }
 89         }
 90     }
 91     void query(int buf[],int len,int ca)
 92     {
 93         //end[3] = 1;
 94         memset(dp,2,sizeof(dp));
 95         dp[0][0] = 0;
 96         for(int i = 0 ;i < len ;i ++)
 97         {
 98             for(int j = 0 ;j < L ;j ++ )
 99             {
100              if(dp[i][j] < len && end[j] != 1)
101                for(int s = 0 ;s < 4;s ++)
102                {
103                  if(end[next[j][s]] != 1)
104                  dp[i+1][next[j][s]] = min(dp[i+1][next[j][s]],dp[i][j] + (buf[i] == s?0:1));
105                }
106             }
107         }
108        int ans = dp[len][0];
109        for(int i = 0 ;i < L ;i ++)
110        {
111            if(end[i] != 1)
112            ans = min(dp[len][i],ans);
113        }
114        if(ans <= len)
115        {
116           printf("Case %d: %d\n",ca,ans);
117        }else{
118           printf("Case %d: -1\n",ca);
119        }
120
121     }
122 };
123 char buf[3000];
124 int  temp[3000];
125 Trie ac;
126 void change(int len)
127 {
128     for(int i = 0;i < len;i ++)
129     {
130        if(buf[i] == 'A')
131        {
132          temp[i] = 0 ;
133        }else if(buf[i] == 'T'){
134          temp[i] = 1;
135        }else if(buf[i] == 'G')
136        {
137           temp[i] = 2;
138        }else if(buf[i] == 'C')
139        {
140           temp[i] = 3;
141        }
142     }
143 }
144 int main(){
145     int n;
146     int ca= 0 ;
147     //freopen("in","r",stdin);
148     while(scanf("%d",&n) != EOF,n)
149     {
150         ca ++ ;
151         ac.init();
152         int len = strlen(buf);
153         for(int i =1 ;i <= n;i ++)
154         {
155           scanf("%s",buf);
156           len = strlen(buf);
157           change(len);
158           ac.insert(temp,len);
159         }
160         ac.build();
161         scanf("%s",buf);
162         len = strlen(buf);
163         change(len);
164         ac.query(temp,len,ca);
165     }
166     return 0;
167 }

View Code

转载于:https://www.cnblogs.com/zyue/p/3973867.html

POJ 3691 DNA repair AC自动机 + DP相关推荐

  1. HDU - 2457 DNA repair(AC自动机+dp)

    题目链接:点击查看 题目大意:给出 n 个匹配串,再给出一个模式串,问最少修改模式串中多少个字母可以使得模式串中不含有任意一个匹配串 题目分析:因为又是模式串与匹配串的题目,虽然与一般意义上的匹配不太 ...

  2. HDU 2457 DNA repair (AC自动机+DP)

    题意:给N个串,一个大串,要求在最小的改变代价下,得到一个不含上述n个串的大串. 思路:dp,f[i][j]代表大串中第i位,AC自动机上第j位的最小代价. 1 #include<algorit ...

  3. POJ 2778 DNA Sequence [AC自动机 + 矩阵快速幂]

    http://poj.org/problem?id=2778 题意:给一些只由ACGT组成的模式串,问有多少种长度为n且不含有给出的模式串的DNA序列. 自动机的状态转换可以看成一个有向图(有重边的) ...

  4. POJ 2778 DNA Sequence (自动机DP+矩阵快速幂)

    题意:给出m个致病DNA片段,求长为n且不含致病片段的DNA序列共有多少种. 数据范围:0 <= m <= 10,1 <= n <=2000000000 这题初看起来与上一题差 ...

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

    题目链接:点击查看 题目大意:给出 n 个长度不大于 10 的字符串表示病毒串,再给出一个长度 len ,问长度为 len 的字符串中,有多少个字符串不含有病毒串作为子串 题目分析:因为病毒串的长度和 ...

  6. [POJ 1625] Censored! (AC自动机+DP+高精度)

    链接 POJ 1625 题意 给出P个模式串,问长度为M且不含有P中任何一个为子串的字符串有多少种. 给出了大小为N的一个字符集,属于ASCII但不一定为英文字母. 最终答案不进行取模,所以可能非常大 ...

  7. uvalive4842(AC自动机+DP)

    题意: 给出猴子打字时打某个字母的概率,猴子最多可以敲键盘m次,问得到的长度是m的单词包含模式串的概率. 思路: AC自动机+dp. 首先,我们用模式串构造一个AC自动机,用dp[i]][j]表示当前 ...

  8. HDU 2296 Ring AC自动机 + DP

    题意:给你n个模式串,每个模式串有一个得分,让你构造出一个长度为N之内且分数最高的文本串;输出字典序列最小的. 解题思路:  AC自动机 + DP , 不过要输出字典序列最小,多开一个 一个三维字符串 ...

  9. bzoj 1030: [JSOI2007]文本生成器(AC自动机+DP)

    1030: [JSOI2007]文本生成器 Time Limit: 1 Sec  Memory Limit: 162 MB Submit: 5187  Solved: 2136 [Submit][St ...

最新文章

  1. 15个目标检测开源数据集汇总
  2. python2.0教程_django2.0入门教程第一节
  3. 5月22日 格式与布局
  4. 根据文件路径检测文件大小并检测是否正在被占用
  5. js如何往div里写div
  6. 在program A里访问program B local class的两种方式
  7. 【渝粤教育】电大中专建设工程法规1作业 题库
  8. 人工智能AI实战100讲(一)-机器人语义建图(上)
  9. springMVC环境搭建
  10. redis的list类型
  11. x-scan mysql_X-Scan使用教程
  12. 大赛评委共论技术趋势,大咖观点干货满满!
  13. 使用Nginx实现负载均衡
  14. 获取datagrid中编辑列combobox的value值与text值
  15. Webstorm克隆Git后npm install报错
  16. adpcb 添加差分对_在AD中PCB设计常用规则——差分规则设置?
  17. 默认关闭IDEA的注释doc的rendered view模式
  18. 学习小甲鱼Python入门(二)习题笔记-列表
  19. 高并发如何处理,解决方案
  20. 计算机网络·计算子网地址,子网掩码,广播地址,根据子网掩码求网络地址详解

热门文章

  1. Vue报错:Error while running task D:\vue\forum:build with message‘spawn vue-cli-service ENOENT‘
  2. java jetty eclipse_用Eclipse+Maven+Jetty构建Java Web开发环境(详细笔记)
  3. python中的datatype啥意思_案例中使用的是dataType,但是用在联系上面dataType不可用,必须改写成type:..._慕课问答...
  4. python 写入excel_使用python将数据写入excel
  5. html实现上下三角,css3实现三角形(上下左右)
  6. bootstrap五星评分_星型打分插件 bootstrap-rating-input
  7. python列表的小东西_Python---列表相关操作
  8. 速算C语言程序设计,C语言程序设计课程设计-速算24.doc
  9. apk可以解压再复制到手机吗_不行了,这个打通手机和电脑的神器,必须得安利给你们...
  10. 工具类的方法怎么引用_Hutool中那些常用的工具类和方法