luogu_4762: [CERC2014]Virus synthesis
洛谷_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的更多相关文章
- [CERC2014]Virus synthesis【回文自动机+DP】
[CERC2014]Virus synthesis 初始有一个空串,利用下面的操作构造给定串 SS . 1.串开头或末尾加一个字符 2.串开头或末尾加一个该串的逆串 求最小化操作数, \(|S| \l ...
- bzoj4044/luoguP4762 [Cerc2014]Virus synthesis(回文自动机+dp)
bzoj4044/luoguP4762 [Cerc2014]Virus synthesis(回文自动机+dp) bzoj Luogu 你要用ATGC四个字母用两种操作拼出给定的串: 1.将其中一个字符 ...
- bzoj4044 [Cerc2014] Virus synthesis
回文自动机上dp f[x]表示形成x代表的回文串所需的最小步数, 若len[x]为奇数,f[x]=len[x],因为即使有更优的,也是直接添加,没有复制操作,那样就不用从x转移了. 若len[x]为偶 ...
- [CERC2014] Virus synthesis
设f[i]为形成极长回文串i的最小操作数.答案为min f[i]+n-len[i]. 在不形成偶回文的情况下形成奇回文的最小操作数为该串长度.可以不考虑(但ans赋为len). 正确性基于: 1)奇. ...
- 洛谷P4762 [CERC2014]Virus synthesis(回文自动机+dp)
传送门 回文自动机的好题啊 先建一个回文自动机,然后记$dp[i]$表示转移到$i$节点代表的回文串的最少的需要次数 首先肯定2操作越多越好,经过2操作之后的串必定是一个回文串,所以最后的答案肯定是由 ...
- bzoj 4044: [Cerc2014] Virus synthesis【回文自动机+dp】
建回文自动机,注意到一个回文串是可以通过一个长度小于等于这个串长度的一半的回文串添上一些字符然后复制得到的,也就是在自动机上向fa走,相当于treedp 每次都走显然会T,记录一个up,指向祖先中最下 ...
- BZOJ 4044 Luogu P4762 [CERC2014]Virus Synthesis (回文自动机、DP)
好难啊..根本不会做..基本上是抄Claris... 题目链接: (bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id=4044 (luogu) ...
- BZOJ4044: [Cerc2014] Virus synthesis(回文树+DP)
Description Viruses are usually bad for your health. How about fighting them with... other viruses? ...
- P4762 [CERC2014]Virus synthesis
题意 真是道回文自动机好题. 首先考虑答案必定是一个回文串+剩余部分的形式,因此可以建出回文自动机,之后考虑每个长度为偶数的回文串. 对于一个长度为偶数的回文串,设它在回文自动机上对应的节点为\(x\ ...
随机推荐
- Python之路【第二十五篇】:数据库之pymysql模块
数据库进阶 一.pymysql模块 pymysql是Python中操作Mysql的模块,其使用的方法和py2的MySQLdb几乎相同. 二.pymysql模块安装 pip install pymysq ...
- Python 3 + Selenium 3 实现汉堡王客户调查提交
用Python 3 + Selenium 3实现汉堡王客户调查的自动填写,可以用来作为 python selenium的入门学习实现脚本,列举了几个比较不太好弄的知识点. 上代码: from sele ...
- AGC038
Contest Page 开题开错翻车场.jpg A sol $A > \frac{W}{2}$或者$B > \frac{H}{2}$的时候无解,否则构造方法长下面这样 #include& ...
- Replication:distribution 中一直在运行 waitfor delay @strdelaytime 语句
Replication 自动创建来一个 Job:Replication monitoring refresher for distribution,这个Agent执行一个sp: dbo.sp_repl ...
- spring Boot 学习(八、Spring Boot与与监控管理)
一.监控管理通过引入spring-boot-starter-actuator,可以使用Spring Boot为我们提供的准 生产环境下的应用监控和管理功能.我们可以通过HTTP,JMX,SSH协议来进 ...
- Geany——Python配置
Geany是一个很不错的编辑器,操作很简单,这里记录一下Geany的入手设置(在下是一个Python程序猿,就以Python为例): 1:新建:选择 下拉菜单中的 main.py ,然后就能生成P ...
- vue-cli项目中使用vw——相比flexible更原生的移动端解决方案
安装命令行输入: yarn add postcss-px-to-viewport 或 npm i postcss-px-to-viewport -save -dev 配置package.json中,在 ...
- 使用EF Core+CodeFirst建立ASP.NET Core MVC项目
本篇随笔介绍如何使用.NET Core+EF Core创建Web应用程序 首先借用官网的话简单介绍一下ASP.NET Core ASP.NET Core 是一个跨平台的高性能开源框架,用于生成基于云且 ...
- vue-cli2和3中的config
最近在网上找了个vue搭建的后台管理的框架,在使用的时候发现没有了config和build文件夹,所以当时就蒙圈了,以为是作者自己改了什么东西,所以感觉自己不知道从何下手了,不过通过查资料发现原来是v ...
- js学习之堆栈内存
**栈内存** >基本数据类型值是直接存放在栈内存中的 栈内存中的变量一般都是已知大小或者有范围上限的,算作一种简单存储.而堆内存存储的对象类型数据对于大小这方面,一般都是未知的.个人认为,这也 ...