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

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

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

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

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

然后这道题就做完了……

才怪。你的 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;
}

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

  1. [十二省联考2019]字符串问题 后缀自动机 + 拓扑排序 + 最长路 + 倍增

    题目描述: 给定一个长串 $S$,给定若干 $S$ 的子串 $a_{i}$, $b_{i}$,再给出一些 $a$ 串和 $b$ 串的支配关系. 构造一个长度最长的字符串,使得: 字符串只由 $a_{i ...

  2. 【洛谷5284】[十二省联考2019] 字符串问题(后缀树优化建边)

    题目: 洛谷 5284 分析: 首先不要问我标题里的「后缀树」是什么,我也不会,瞎写的 -- (传说就是反串后缀自动机的 fa 树?) 前置技能:[知识总结]后缀自动机的构建 首先有一个很 naive ...

  3. 十二省联考 2019 题解

    [十二省联考2019]异或粽子 首先异或转前缀和,类似超级钢琴,将三元组 ( l , r , p ) (l,r,p) (l,r,p) 插入堆,表示 s u m [ p ] sum[p] sum[p] ...

  4. 【BZOJ5498】[十二省联考2019]皮配(动态规划)

    [BZOJ5498][十二省联考2019]皮配(动态规划) 题面 BZOJ 洛谷 题解 先考虑暴力\(dp\),设\(f[i][j][k]\)表示前\(i\)所学校,有\(j\)人在某个阵营,有\(k ...

  5. [十二省联考2019]春节十二响——长链剖分+堆

    题目链接: [十二省联考2019]春节十二响 可以发现每条链上的所有点都要放在不同的段里,那么最多只需要树的深度这么多段就够了. 因为这样可以保证每条链上的点可以放在不同的段中而且一个点放在这些段中一 ...

  6. 「十二省联考 2019」皮配——dp

    题目 [题目描述] #### 题目背景 一年一度的综艺节目<中国好码农>又开始了.本季度,好码农由 Yazid.Zayid.小 R.大 R 四位梦想导师坐镇,他们都将组建自己的梦想战队,并 ...

  7. 【十二省联考2019】春节十二响

    题面 https://www.luogu.org/problem/P5290 题解 真的是我傻逼,十二省联考$day2$至今还是我的噩梦.$day1$起码一直在调可持久化$trie$树,$day2$真 ...

  8. HAOI(十二省联考)2019 qwq记

    \(\large{Day\ -1}:\) 放假了,白天大概是抱着最后一次在机房的心态复习着板子过去的.看着机房里的各位神仙丝毫不慌的颓倒是有点慌了,敲了一下多项式的板子感觉写的相当自闭,感觉AFO应该 ...

  9. 十二省联考2019酱油记

    在中考前去省选玩一趟. Day -1 对于一个还没有学会所有省选内容的初三Oier来说,这一趟真的是去打酱油的啊.但还是要认真复习. 最近几天在字符串的路上越走越远-晚上才开始复习图论.还有一大堆没有 ...

最新文章

  1. 关于js中的时间处理
  2. WebRTC / Jitsi / 多人视频通讯常用架构 Mesh / MCU / SFU
  3. Visual Entity 手册(十一)代码生成设置
  4. Java package 与 import
  5. AMOS问卷数据建模前传【SPSS 052期】
  6. JavaScript常见的运算符优先级面试题
  7. IAR在写结构体时不提示_U盘被写保护了以后怎么办?
  8. AB触摸屏2711P-T12W22A9P、2711P-T10C4A9,2711P-T10C22D8S-B、2711P-T10C22D9P-B、2711P-T12C10D2
  9. python画椭圆形_如何绘制Python中的半椭圆?
  10. 顶峰Android手机视频转换器,顶峰FLV视频转换器
  11. 罗永浩“吓尿”大会上,我裤子真湿了
  12. nltk文件下载以及word_tokenize运行问题解决
  13. 拒酒词,好难找哟,留到有用
  14. JAVA基于局域网的聊天室系统(源代码+论文)
  15. 单点登录(SSO)解决方案介绍
  16. PCIe扫盲——热插拔简要介绍
  17. 杂谈 之 闲来无事 (一)
  18. c语言漩涡图形,我看大家都在讨论漩涡矩阵,网上搜了一个大家讨论讨论。
  19. 饮用水铁离子超标,各种溶液铁离子超标去除工艺
  20. 公众号榜单 | 2020·6月公众号地区排行榜重磅发布

热门文章

  1. 我怀疑对象做了什么对不起我的事......
  2. 老娘不就是没化妆吗?你几个意思?
  3. 有些图,只要看错一眼就再也回不去了!
  4. 2019优质公众号大盘点,果断收藏了慢慢看吧 ~
  5. 她只用1个方法,就把英语拿下了!
  6. 要成为年薪五十万的数据分析师,除了技术还需要什么?
  7. 【直观详解】支持向量机SVM
  8. 仿真模拟,需要注意这几点!
  9. mysql enum_MySQL数据库中关于ENUM类型的详细解释
  10. docker 查看容器_Docker介绍