洛谷_4762:[CERC2014]Virus synthesis

题目描述:

  • 初始有一个空串,利用下面的操作构造给定串\(S\)。\(len(S)\leq10^5\)

    • 1: 串开头或末尾加一个字符。
    • 2: 串开头或末尾添加一个该串的逆串。

输入描述:

  • 给出一个\(T\)表示要处理的字符串的数量。
  • 接下来第\(2\)~\(T+1\)行每行给出一个字符串。

输出描述:

  • 对于每一个字符串,输出一个正整数表示答案。

思路:

  • 回文自动机。
  • 先建立一个回文自动机,然后记\(f(i)\)表示转移到\(i\)节点结尾代表的回文串的最少需要的次数。
  • 操作\(2\)肯定越多越好,经过操作\(2\)得到的肯定是一个回文串,那么最后的答案肯定是回文串\(+\)暴力
    • 可以知道最后答案为\(ans=min(ans,f(i)+strlen-len(i))\)
  • 对于一个串\(i\),如果在他前面和后面加上一个相同的字母可以形成回文串\(j\),则\(f(j)=f(i)+1\)。
  • 为什么是对的呢?就相当于形成\(i\)只有回文串的前一半,先在前一半的前面加上那个字母后再进行操作\(2\)就可以变成字符串\(j\)了。
  • 对应在自动机里就是如果一个节点有子节点,那么\(f(son)=f(fa)+1\)。
  • 对于遍历到的位置\(x\),他的\(trans\)指针指向\(y\),那么有\(f(x)=min(f(x),f(y)+\:(len(x)/2-len(y))\:+1)\)
  • 最后对\(trie\)图进行\(bfs\)更新答案。

代码:

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e5 +10, INF = 0x3f3f3f3f;
int T, ans; char s[maxn];
int cnt, trie[maxn][6], trans[maxn], w[110];
int len[maxn], fail[maxn], len_str, last; int f[maxn]; void init()
{
w['A'] = 0, w['T'] = 1, w['C'] = 2, w['G'] = 3;
scanf("%s", s + 1);
len_str = ans = strlen(s + 1);
cnt = 1, last = 0;
len[0] = 0, len[1] = -1;
fail[0] = 1, fail[1] = 0;
memset(trie[0], 0, sizeof(trie[0]));
memset(trie[1], 0, sizeof(trie[1]));
} inline int get_fail(int las, int i)
{
while(s[i - len[las] - 1] != s[i])
las = fail[las];
return las;
} void build_PAM()
{
for(int i = 1; i <= len_str; i++)
{
int num = w[int(s[i])];
int p = get_fail(last, i);
if(!trie[p][num])
{
len[++cnt] = len[p] + 2;
memset(trie[cnt], 0, sizeof(trie[cnt]));
fail[cnt] = trie[get_fail(fail[p], i)][num];
trie[p][num] = cnt; //--------求trans指针
if(len[cnt] <= 2) trans[cnt] = fail[cnt];
else
{
int tmp = trans[p];
while((s[i - len[tmp] - 1] != s[i]) || ( ((len[tmp] + 2) * 2) > len[cnt]))
tmp = fail[tmp];
trans[cnt] = trie[tmp][num];
}
//--------
}
last = trie[p][num];
}
} int main()
{
scanf("%d", &T);
while(T--)
{
init(); build_PAM();
for(int i = 2; i <= cnt; i++)
f[i] = len[i]; f[0] = 1;
queue<int> q; q.push(0);
while(q.size())
{
int t = q.front(); q.pop();
for(int i = 0, x, y; i <= 3; i++)
{
x = trie[t][i]; if(!x) continue;
f[x] = f[t] + 1; y = trans[x];
f[x] = min(f[x], f[y]+1+len[x]/2-len[y]);
ans = min(ans, f[x]+len_str-len[x]);
q.push(x);
}
}
printf("%d\n", ans);
}
return 0;
}
/*
4
AAAA
AGCTTGCA
AAGGGGAAGGGGAA
AAACAGTCCTGACAAAAAAAAAAAAC
*/

luogu_4762: [CERC2014]Virus synthesis的更多相关文章

  1. [CERC2014]Virus synthesis【回文自动机+DP】

    [CERC2014]Virus synthesis 初始有一个空串,利用下面的操作构造给定串 SS . 1.串开头或末尾加一个字符 2.串开头或末尾加一个该串的逆串 求最小化操作数, \(|S| \l ...

  2. bzoj4044/luoguP4762 [Cerc2014]Virus synthesis(回文自动机+dp)

    bzoj4044/luoguP4762 [Cerc2014]Virus synthesis(回文自动机+dp) bzoj Luogu 你要用ATGC四个字母用两种操作拼出给定的串: 1.将其中一个字符 ...

  3. bzoj4044 [Cerc2014] Virus synthesis

    回文自动机上dp f[x]表示形成x代表的回文串所需的最小步数, 若len[x]为奇数,f[x]=len[x],因为即使有更优的,也是直接添加,没有复制操作,那样就不用从x转移了. 若len[x]为偶 ...

  4. [CERC2014] Virus synthesis

    设f[i]为形成极长回文串i的最小操作数.答案为min f[i]+n-len[i]. 在不形成偶回文的情况下形成奇回文的最小操作数为该串长度.可以不考虑(但ans赋为len). 正确性基于: 1)奇. ...

  5. 洛谷P4762 [CERC2014]Virus synthesis(回文自动机+dp)

    传送门 回文自动机的好题啊 先建一个回文自动机,然后记$dp[i]$表示转移到$i$节点代表的回文串的最少的需要次数 首先肯定2操作越多越好,经过2操作之后的串必定是一个回文串,所以最后的答案肯定是由 ...

  6. bzoj 4044: [Cerc2014] Virus synthesis【回文自动机+dp】

    建回文自动机,注意到一个回文串是可以通过一个长度小于等于这个串长度的一半的回文串添上一些字符然后复制得到的,也就是在自动机上向fa走,相当于treedp 每次都走显然会T,记录一个up,指向祖先中最下 ...

  7. BZOJ 4044 Luogu P4762 [CERC2014]Virus Synthesis (回文自动机、DP)

    好难啊..根本不会做..基本上是抄Claris... 题目链接: (bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id=4044 (luogu) ...

  8. BZOJ4044: [Cerc2014] Virus synthesis(回文树+DP)

    Description Viruses are usually bad for your health. How about fighting them with... other viruses? ...

  9. P4762 [CERC2014]Virus synthesis

    题意 真是道回文自动机好题. 首先考虑答案必定是一个回文串+剩余部分的形式,因此可以建出回文自动机,之后考虑每个长度为偶数的回文串. 对于一个长度为偶数的回文串,设它在回文自动机上对应的节点为\(x\ ...

随机推荐

  1. SpringBootSecurity学习(25)前后端分离版之OAuth2.0 令牌中继

    增加Eureka 前面介绍的项目都是授权服务和资源服务单独两个,这样在资源服务中的 check_token 地址都是写死的地址 : 下面我们把eureka加上,这样就可以直接用服务名了.eureka服 ...

  2. 一张图看懂SharpCamera

    通过下面的图片,可以瞬间看懂整个类库的脉络.

  3. ASP.NET Core 配置文件

    在ASP.NET Core 中,应用程序配置数据可以使用JSON, XML 和 INI格式 和内置环境变量,命令行参数或内存中的集合. 1.如何获取和设置配置 ASP.NET Core配置系统针对以前 ...

  4. pands模块的妙用爬取网页中的表格

    拿我这篇为例https://www.cnblogs.com/pythonywy/p/11574340.html import pandas as pd df = pd.read_html('https ...

  5. webpack 里的 import, exports 实现原理

    在使用 webpack 对脚本进行打包, 在开发中, 每个文件中都会使用 import 语句来导入一些功能,又会使用 export 语句导出一些功能,为了研究 import 和 export 原理,研 ...

  6. 【MySQL】mysql中的锁机制

    一.分类 MySQL的锁机制不同的存储引擎支持不同的锁机制,分为表级锁.行级锁.页面锁.MyISAM和MEMORY存储引擎采用的是表级锁(table-level locking):BDB存储引擎采用的 ...

  7. django子应用

    在Web应用中,通常有一些业务功能模块是在不同的项目中都可以复用的,故在开发中通常将工程项目拆分为不同的子功能模块,各功能模块间可以保持相对的独立,在其他工程项目中需要用到某个特定功能模块时,可以将该 ...

  8. 学习笔记之自然语言处理(Natural Language Processing)

    自然语言处理 - 维基百科,自由的百科全书 https://zh.wikipedia.org/wiki/%E8%87%AA%E7%84%B6%E8%AF%AD%E8%A8%80%E5%A4%84%E7 ...

  9. vue路由切换时内容组件的滚动条回到顶部

    在使用vue的时候会出现切换路由的时候滚动条保持在原来的位置,要切换路由的时候滚动条回到顶部才有更好的用户体验 1.当页面整体都要滚动到顶部的情况 router.afterEach(() => ...

  10. spring boot 的request.getServletContext().getRealPath路径获取问题

    默认情况下springboot中request.getServletContext().getRealPath 返回的是一个临时文件夹的地址 通过查看源代码 位置在 org.springframewo ...