【十二省联考2019】字符串问题【后缀自动机】【拓扑排序】

题意:给一个字符串 SSS,以子串的形式给出一些 A 类串和 B 类串以及 mmm 对 A 类串支配 B 类串的关系。求一个总长度最长的 A 类串序列,使得每个串都存在一个 B 类串前缀被后一个串支配。无穷输出 −1-11

∣S∣,m≤2×105|S|,m\leq 2\times 10^5S,m2×105

显然需要先把反串 SAM 建出来,然后树上倍增定位每个子串。

对于每一步,相当于是可以不断丢掉最后一个字符变成前缀,如果是个 B 类串,可以走到一个被支配的 A 类串。

所以 SAM 上每个点向 fail 树上的父亲连一条 000 的边, uuu 支配 vvv 就从 vvvuuu∣u∣|u|u 的边。然后拓扑排序后做最长路即可。因为没有 000 环,也没有不联通等奇奇怪怪的东西,所以有环就是 −1-11

然后这道题就做完了……

才怪。你的 A,B 类串定位的是一类子串,而限制是死的。具体来说,一个 B 类串定位到了某个结点的一坨子串上,然后一个 A 类串定位也定位到了这个结点,而且还比 B 类串短,那么它是不能走到这个 B 的。但是 SAM 上它们是同一个结点,所以上面的算法就会出错。

怎么办?拆了呗

当然不用真的拆了,把这些定位的点放结点的 vector 里,然后在内部排序建一下就可以了。

不算难写,单纯的码量大。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <vector>
#include <utility>
#include <queue>
#include <algorithm>
#define MAXN 800005
#define MAXM 1200005
using namespace std;
typedef long long ll;
inline int read()
{int ans=0;char c=getchar();while (!isdigit(c)) c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return ans;
}
char s[MAXN];
int ch[MAXN][26],len[MAXN],fa[MAXN],las=1,tot=1;
void insert(int c)
{int cur=++tot;len[cur]=len[las]+1;int p=las;for (;p&&!ch[p][c];p=fa[p]) ch[p][c]=cur;if (!p) fa[cur]=1;else{int q=ch[p][c];if (len[q]==len[p]+1) fa[cur]=q;else{int _q=++tot;len[_q]=len[p]+1;fa[_q]=fa[q],fa[cur]=fa[q]=_q;memcpy(ch[_q],ch[q],sizeof(ch[_q]));for (;p&&ch[p][c]==q;p=fa[p]) ch[p][c]=_q;}}las=cur;
}
int up[MAXN][20];
inline void build()
{for (int i=1;i<=tot;i++) up[i][0]=fa[i];for (int j=1;j<20;j++)for (int i=1;i<=tot;i++)up[i][j]=up[up[i][j-1]][j-1];
}
inline int find(int x,int l)
{for (int i=19;i>=0;i--)if (len[up[x][i]]>=l)x=up[x][i];return x;
}
int na,nb,la[MAXN],lb[MAXN],ra[MAXN],rb[MAXN];
struct node
{	int idx,type;inline int len()const{return type? ra[idx]-la[idx]+1:rb[idx]-lb[idx]+1;}inline int pos()const{if (type==1) return tot+idx;if (type==0) return tot+na+idx;return idx;}
};
inline bool operator <(const node& a,const node& b){return a.len()<b.len()||(a.len()==b.len()&&a.type<b.type);}
vector<node> lis[MAXN];
struct edge{int u,v,w;}e[MAXM];
int head[MAXN],nxt[MAXM],deg[MAXN],cnt;
inline void addnode(int u,int v,int w)
{e[++cnt]=(edge){u,v,w};nxt[cnt]=head[u];head[u]=cnt;++deg[v];
}
int pos[MAXN],dfn[MAXN],tim;
ll dp[MAXN];
int q[MAXN],qs,qt;
void clear()
{for (int i=1;i<=tot;i++) {memset(ch[i],0,sizeof(ch[i])),fa[i]=len[i]=0;memset(up[i],0,sizeof(up[i]));vector<node>().swap(lis[i]);}for (int i=1;i<=tot+na+nb+1;i++) head[i]=deg[i]=dp[i]=0;for (int i=1;i<=cnt;i++) nxt[i]=0;tim=cnt=0,las=tot=1;
}
int main()
{for (int T=read();T;T--){scanf("%s",s+1);int n=strlen(s+1);for (int i=n;i>=1;i--) insert(s[i]-'a'),pos[i]=las;build();	na=read();for (int i=1;i<=na;i++){la[i]=read(),ra[i]=read();lis[find(pos[la[i]],ra[i]-la[i]+1)].push_back((node){i,1});}nb=read();for (int i=1;i<=nb;i++){lb[i]=read(),rb[i]=read();lis[find(pos[lb[i]],rb[i]-lb[i]+1)].push_back((node){i,0});}for (int i=1;i<=tot;i++) {sort(lis[i].begin(),lis[i].end());lis[i].push_back((node){i,-1});if (fa[i]) addnode(lis[i].front().pos(),fa[i],0);for (int j=0;j<(int)lis[i].size()-1;j++) addnode(lis[i][j+1].pos(),lis[i][j].pos(),0);}for (int m=read();m;m--){int u,v;u=read(),v=read();addnode(tot+na+v,tot+u,ra[u]-la[u]+1);}int S=tot+na+nb+1;for (int i=1;i<=na;i++) addnode(S,tot+i,ra[i]-la[i]+1);qs=1,qt=0;for (int i=1;i<=S;i++) if (!deg[i]) q[++qt]=i;while (qs<=qt){int u=dfn[++tim]=q[qs++];for (int i=head[u];i;i=nxt[i])if (!(--deg[e[i].v]))q[++qt]=e[i].v;}if (tim<S) puts("-1");else{for (int k=S;k>=1;k--){int u=dfn[k];for (int i=head[u];i;i=nxt[i])dp[u]=max(dp[u],dp[e[i].v]+e[i].w);}printf("%lld\n",dp[S]);}clear();}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314799.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Codeforces Round #586 (Div. 1 + Div. 2) B. Multiplication Table 思维 + 公式

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个n∗nn*nn∗n的矩阵&#xff0c;每个位置由ai∗aja_i*a_jai​∗aj​得来&#xff0c;主对角线为000&#xff0c;让你求出来aia_iai​。 n≤1e3n\le1e3n≤1e3 思路&#xff1a; 由公式ai,j∗ai,kaj,…

不好意思,这么久没有更新《从零开始掌握ASP.NET Core 》

点击上方蓝字&#xff0c;关注「我们」等了快个月了&#xff0c;终于开始更新了。因为感冒&#xff0c;弄的嗓子有点沙哑。所以停了半个月才是更新&#xff0c;目前一口气更新了12个章节&#xff0c;大家可以耐心观看内容了。《从零开始学ASP.NET Core 》-- 更新通知视频课程更…

【NOI2018】你的名字【后缀自动机】【可持久化线段树合并】【乱搞】

题意&#xff1a;给一个串 SSS&#xff0c;qqq 次询问&#xff0c;每次给定串 TTT 和 l,rl,rl,r &#xff0c;求有多少个本质不同的串是 TTT 的子串而不是 Sl…rS_{l\dots r}Sl…r​ 的子串。 ∣S∣≤5105,q≤105,∑∣T∣≤106|S|\leq 5\times 10^5,q\leq 10^5,\sum|T|\leq 10^…

2021牛客暑期多校训练营1 A.Alice and Bob 博弈 SG函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 有两堆石子&#xff0c;两个人每次可以进行如下操作&#xff1a;从某一堆狮子中拿出x(x>0)x(x>0)x(x>0)个&#xff0c;从另一堆石子中拿出s∗x(s>0)s*x(s>0)s∗x(s>0)个。谁不能操作谁输&…

【CTSC2010】珠宝商【后缀自动机】【点分治】【根号分治】

题意&#xff1a;给一棵 nnn 个点的树&#xff0c;每个点有个字符&#xff0c;另给一个长度为 mmm 的特征串&#xff0c;求树上 n2n^2n2 条有向路径在特征串中出现的次数之和。 n,m≤5104n,m\leq 5\times 10^4n,m≤5104 看到母串先建 SAM &#xff08;bushi 树上路径统计问题…

使用Azure云原生构建博客是怎样一种体验?(下篇)

点击上方蓝字关注“汪宇杰博客”接上篇《使用Azure云原生构建博客是怎样一种体验&#xff1f;&#xff08;上篇&#xff09;》DNSAzure DNS 是一套分布全球的域名解析服务。具有超高可用性和接近实时的记录更新及生效速度。我的博客也使用了这项服务。Azure 现在可以提供域名注…

2021牛客暑期多校训练营1 G Game of Swapping Numbers 思维 + 巧妙的转换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个数组A,BA,BA,B&#xff0c;你可以选择AAA的两个位置i,j,i<ji,j,i<ji,j,i<j交换Ai,AjA_i,A_jAi​,Aj​&#xff0c;需要交换正好kkk次&#xff0c;问你最大的∑i1n∣Ai−Bi∣\sum_{i1}^n|A_…

.NET Core 3.0之深入源码理解HttpClientFactory(一)

写在前面创建HttpClient实例的时候&#xff0c;在内部会创建HttpMessageHandler链&#xff0c;我们知道HttpMessageHandler是负责建立连接的抽象处理程序&#xff0c;所以HttpClient的维护实际上就是维护HttpMessageHandler的使用&#xff0c;释放HttpClient并不会及时释放连接…

WTM 构建DotNetCore开源生态,坐而论道不如起而行之

作为一个8岁开始学习编程&#xff0c;至今40岁的老程序员&#xff0c;这辈子使用过无数种语言&#xff0c;从basic开始&#xff0c;到pascal, C, C&#xff0c;到后来的 java, c#,perl,php,再到现在流行的python。小时候的我总觉得多掌握一门语言&#xff0c;我的技术能力就又前…

【WC2014】时空穿梭【组合数】【莫比乌斯反演】【整除分块】【暴力多项式】

题意&#xff1a;TTT 组数据&#xff0c;给一个 nnn 维空间&#xff0c;第 iii 维大小为 [1,mi]∩Z[1,m_i]\cap \Z[1,mi​]∩Z&#xff0c;求大小为 ccc 的严格偏序上升的共线点集个数。答案模 100071000710007。 T≤100,n≤11,m≤105,c≤20T\leq 100,n\leq 11,m\leq 10^5,c\le…

2021牛客暑期多校训练营1 H Hash Function FFT\NTT

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个数组aaa&#xff0c;你需要找一个最小的模数xxx&#xff0c;使得aaa中每个数都模上xxx之后互不相同。 n≤5e5,ai≤5e5,ai!ajn\le5e5,a_i\le5e5,a_i!a_jn≤5e5,ai​≤5e5,ai​!aj​ 思路&#xff1a…

架构杂谈《六》

超时处理模式在服务化或者微服务架构里&#xff0c;传统的整体应用拆分成多个职责单一的微服务&#xff0c;微服务之间通过某种网络通信协议互相通信和交互&#xff0c;完成特定的功能&#xff0c;然而由于网络通信的不稳定&#xff0c;在设计系统时必须考虑到对网络通信的容错…

【BZOJ4543】Hotel加强版【神仙树形dp】【长链剖分】

题意&#xff1a;给一棵 nnn 个点的树&#xff0c;求两两距离相等的三元组个数。 n≤105n\leq 10^5n≤105 显然相当于是找一个点到这三个点距离相等。子树内和子树外到当前点的距离为某个值的点的个数可以长链剖分快速得到&#xff0c;但统计答案非常棘手。 接下来是个鬼才想…

hdu 6962 I love tree 线段树维护二次函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个点的一颗树&#xff0c;有mmm次询问&#xff0c;每次询问有两个操作&#xff1a; (1)(1)(1)将[a,b][a,b][a,b]路径上的点依次加上12,22,32,...,len2,lenpath(a,b)1^2,2^2,3^2,...,len^2,lenpath(a,…

基于surging 的stage组件设计,谈谈我眼中的微服务

一、前言surging 开源地址&#xff1a;https://github.com/dotnetcore/surging随着业务的发展&#xff0c;并发量的增多&#xff0c;业务的复杂度越来越大&#xff0c;对于系统架构能力要求越来越高&#xff0c;这时候微服务的设计思想应运而生&#xff0c;但是对于微服务需要引…

【PKUSC2018】星际穿越【结论】【倍增dp】

题意&#xff1a;有一张边权为 111 的无向图&#xff0c;对 i∈[2,n]i\in [2,n]i∈[2,n]&#xff0c;iii 与 [li,i−1][l_i,i-1][li​,i−1] 间有边。 qqq 次询问 l,r,xl,r,xl,r,x&#xff0c;表示 xxx 与 [l,r][l,r][l,r] 中的所有点的最短路长度的平均值&#xff0c;其中 l&l…

HDU - 6971 K - I love max and multiply sosdp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 直接求i&j>ki\And j>ki&j>k不是很好求&#xff0c;所以转换成i&jki\And jki&jk的情况。 考虑对a,ba,ba,b求一遍超集&#xff0c;让后从[0,n−1][0,n-1][0,n−1]扫…

推荐10个技术圈优质的公众号大号

公众号有很多但需要什么只有自己知道本次筛选了一批技术圈优质的公众号&#xff0c;主要与python、人工智能、机器学习、技术人生相关希望对你有所帮助!▼★长按二维码&#xff0c;选择“识别二维码”进行关注。▲长按二维码&#xff0c;识别关注简介&#xff1a;Python爱好者社…

【SDOI2013】项链【莫比乌斯反演】【Polya定理】【递推式求通项】【数论】

题意&#xff1a;TTT 组数据&#xff0c;每组给定 n,an,an,a&#xff0c;求满足下列条件的项链数量&#xff1a; 有 nnn 个珠子。每个珠子上有三个 [1,a]∩Z[1,a]\cap \Z[1,a]∩Z 的数&#xff0c;且三个数 gcd⁡\gcdgcd 为 111。相邻两个珠子不同。 珠子旋转、翻转同构&…

HDU - 6955 Xor sum tire树 + 贪心

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个数列aaa&#xff0c;你需要找出来一个长度最小且左端点最靠前的区间&#xff0c;使其异或和≥k\ge k≥k。 n≤1e5,0≤ai,k<230n\le1e5,0\le a_i,k<2^{30}n≤1e5,0≤ai​,k<230 思路&#…