【乱搞】某 SCOI 模拟赛 T3 串（string）【PAM 乱搞】

最新推荐文章于 2020-07-21 15:50:04 发布
最新推荐文章于 2020-07-21 15:50:04 发布
阅读量312 收藏 1
点赞数
分类专栏： # 来源-模拟赛 # 字符串-回文自动机 # 其它-乱搞
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_39854734/article/details/106650457
版权

目录

题意

对于字符串 S S S 的每个前缀，求其所以回文子串的前缀个数和，本质相同的前缀不重复计数。 N ≤ 3 × 1 0 5 N\leq3\times 10^5 N≤3×105。

乱搞

考虑在末尾新加入一个字符时会产生多少贡献。

假如新出现前缀，它肯定是新产生的回文串（记为 c u r cur cur）的前缀，也就是 PAM 上新加的结点的前缀。

这些前缀有一些是在之前的回文串里面出现过的，因此我们尝试找到 c u r cur cur 与之前任意一个回文串的 LCP 的最大值，则不长于这个 LCP（记为 t t t）的前缀都是出现过的。

t t t 至少有 c u r . f a i l cur.fail cur.fail 那么长。以 t t t 为前缀的回文串要么就是 c u r . f a i l cur.fail cur.fail；要么比 c u r . f a i l cur.fail cur.fail 长，那么这个回文串（记为 u u u）一定是 c u r . f a i l cur.fail cur.fail 在 fail 树上的后代。

实际上第三种情况只用枚举 u . f a i l u.fail u.fail 等于 c u r . f a i l cur.fail cur.fail 的 u u u。原因：假设 u u u 有一个长于 c u r . f a i l cur.fail cur.fail 的回文前缀 v = u . f a i l v=u.fail v=u.fail。若 v v v 是 t t t 的前缀，则 v v v 也应该是 c u r cur cur 的前缀， c u r . f a i l cur.fail cur.fail 应当是 v v v，矛盾；因此 t t t 是 v v v 的前缀，则枚举 v v v 就能统计到 t t t。

因此每次在 PAM 中新加入节点 c u r cur cur 时，暴力枚举 c u r cur cur 在 fail 树上的兄弟 u u u，并判断 c u r cur cur 与 u u u 的最长公共前缀。答案的增量即为 c u r . l e n − max ⁡ LCP cur.len-\max \text{LCP} cur.len−maxLCP。

时间复杂度： O ( 没被卡 ) O(\text{没被卡}) O(没被卡)

代码（代码后还有内容）：

#include<bits/stdc++.h>
using namespace std;
int getint(){
	int ans=0,f=1;
	char c=getchar();
	while(c<'0'||c>'9'){
		if(c=='-')f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		ans=ans*10+c-'0';
		c=getchar();
	}
	return ans*f;
}
const int N=3e5+10;
char s[N];
int n;
namespace PAM{

int ch[N][26],fail[N],len[N],cnt,lst,po[N];
vector<int>fch[N];
string str[N];
long long ans=0;
void init(){
	cnt=2;
	fail[0]=1;len[0]=0;fch[1].push_back(0);
	fail[1]=1;len[1]=-1;
}
void extend(char c,int pos){
	int cur=lst;
	while(s[pos-len[cur]-1]!=c)cur=fail[cur];
	if(!ch[cur][c-'a']){
		lst=cnt++;
		len[lst]=len[cur]+2;
		po[lst]=pos-len[lst]+1;
		fail[lst]=fail[cur];
		while(s[pos-len[fail[lst]]-1]!=c)fail[lst]=fail[fail[lst]];
		fail[lst]=ch[fail[lst]][c-'a'];
		fch[fail[lst]].push_back(lst);
		ch[cur][c-'a']=lst;
		if(0);else{
			int fa=fail[lst];
			int l=len[fa];
			for(int i=0;i<fch[fa].size();i++){
				if(fch[fa][i]==lst)continue;
				int ll=len[fa];
				while(s[po[lst]+ll]==s[po[fch[fa][i]]+ll])++ll;//脑抽没有用 SA 判 LCP
				l=max(l,ll);
			}
			ans+=max(0,len[lst]-l);
		}
	}
	lst=ch[cur][c-'a'];
}
void print(){
	for(int i=0;i<cnt;i++){
		for(int j=0;j<26;j++){
			if(ch[i][j])cerr<<i<<" "<<ch[i][j]<<" "<<char('a'+j)<<endl;
		}
		cerr<<"fail "<<i<<" "<<fail[i]<<endl;
	}
}

}//namespace PAM

int main(){
	freopen("string.in","r",stdin);
	freopen("string.out","w",stdout);
	n=getint();
	scanf("%s",s+1);
	PAM::init();
	for(int i=1;i<=n;i++){
		PAM::extend(s[i],i);
		printf("%d\n",PAM::ans);
	}
	//print();
	return 0;
}

Hack

想到斐波那契字符串很强，于是乱 Hack 一波。

f=open("string.in","w")
f.write("300000\n")
sb="mnbvcfdssdfcvbnm"
sa="qwertyuioplkjhgazxxzaghjklpoiuytrewq"
while len(sa)<=300000:
    sc=sa+sb
    sa=sb
    sb=sc
# f.write(sa[0:300000])
f.write('a'*300000)
f.close()

能将上述乱搞卡到 2s 左右。

Anti-Hack

除了像我这样没有脑子的人以外，大家应该都想得到用 SA 来求 LCP，于是上述 Hack 只能把它卡到 0.5s 左右。

代码（在出题人的数据上跑得更慢了）（代码后还有内容）：

#include<bits/stdc++.h>
using namespace std;
int getint(){
	int ans=0,f=1;
	char c=getchar();
	while(c<'0'||c>'9'){
		if(c=='-')f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		ans=ans*10+c-'0';
		c=getchar();
	}
	return ans*f;
}
const int N=3e5+10;
char s[N];
int n;

namespace SA{

int sa[N],rk[N],tp[N],h[N];
int c[N];
void rsort(int m){
	for(int i=1;i<=m;i++)c[i]=0;
	for(int i=1;i<=n;i++)c[rk[tp[i]]]++;
	for(int i=1;i<=m;i++)c[i]+=c[i-1];
	for(int i=n;i>=1;i--)sa[c[rk[tp[i]]]--]=tp[i];
}
bool cmp(int x,int y,int k){
	return tp[x]==tp[y]&&tp[x+k]==tp[y+k];
}
void get_sa(){
	for(int i=1;i<=n;i++)rk[i]=s[i],tp[i]=i;
	int m='z';
	rsort(m);
	for(int p=0,k=1;p<n;k<<=1,m=p){
		p=0;
		for(int i=n-k+1;i<=n;i++)tp[++p]=i;
		for(int i=1;i<=n;i++)
			if(sa[i]>k)tp[++p]=sa[i]-k;
		rsort(m);
		memcpy(tp,rk,sizeof(tp));
		rk[sa[1]]=p=1;
		for(int i=2;i<=n;i++)
			rk[sa[i]]=cmp(sa[i],sa[i-1],k)?p:++p;
	}
	//for(int i=1;i<=n;i++)cerr<<sa[i]<<" ";cerr<<endl;
	//for(int i=1;i<=n;i++)cerr<<rk[i]<<" ";cerr<<endl;
	int t=0;
	for(int i=1;i<=n;i++){
		if(t)--t;
		while(s[sa[rk[i]-1]+t]==s[i+t])++t;
		h[rk[i]]=t;
	}
	//for(int i=1;i<=n;i++)cerr<<h[i]<<" ";cerr<<endl;
}
const int L=20;
int st[L][N],l2[N];
void init_st(){
	l2[0]=-1;for(int i=1;i<=n;i++)st[0][i]=h[i+1],l2[i]=l2[i>>1]+1;
	for(int i=1;i<L;i++){
		for(int j=1;j<=n-(1<<i-1);j++){
			st[i][j]=min(st[i-1][j],st[i-1][j+(1<<i-1)]);
		}
	}
}
int get_min(int l,int r){
	if(l>r)swap(l,r);
	if(l==r)return n-sa[l]+1;
	int t=l2[r-l];
	return min(st[t][l],st[t][r-(1<<t)]);
}

}

namespace PAM{

int ch[N][26],fail[N],len[N],cnt,lst,po[N];
vector<int>fch[N];
string str[N];
long long ans=0;
void init(){
	cnt=2;
	fail[0]=1;len[0]=0;fch[1].push_back(0);
	fail[1]=1;len[1]=-1;
}
long long qaq=0;
void extend(char c,int pos){
	int cur=lst;
	while(s[pos-len[cur]-1]!=c)cur=fail[cur];
	if(!ch[cur][c-'a']){
		lst=cnt++;
		len[lst]=len[cur]+2;
		po[lst]=pos-len[lst]+1;
		fail[lst]=fail[cur];
		while(s[pos-len[fail[lst]]-1]!=c)fail[lst]=fail[fail[lst]];
		fail[lst]=ch[fail[lst]][c-'a'];
		fch[fail[lst]].push_back(lst);
		ch[cur][c-'a']=lst;
		if(0);else{
			int fa=fail[lst];
			int l=len[fa];
			for(int i=0;i<fch[fa].size();i++){
				if(fch[fa][i]==lst)continue;
				if(s[po[lst]+len[fa]]!=s[po[fch[fa][i]]+len[fa]])continue;
				//int ll=len[fa];
				//while(s[po[lst]+ll]==s[po[fch[fa][i]]+ll])++ll;
				int ll=SA::get_min(SA::rk[po[lst]],SA::rk[po[fch[fa][i]]]);
				l=max(l,ll);
				//++qaq;
			}
			ans+=max(0,len[lst]-l);
		}
	}
	lst=ch[cur][c-'a'];
}
void print(){
	for(int i=0;i<cnt;i++){
		for(int j=0;j<26;j++){
			if(ch[i][j])cerr<<i<<" "<<ch[i][j]<<" "<<char('a'+j)<<endl;
		}
		cerr<<"fail "<<i<<" "<<fail[i]<<endl;
	}
}

}//namespace PAM

int main(){
	freopen("string.in","r",stdin);
	freopen("string.out","w",stdout);
	n=getint();
	scanf("%s",s+1);
	PAM::init();
	SA::get_sa();
	SA::init_st();
	for(int i=1;i<=n;i++){
		PAM::extend(s[i],i);
		printf("%d\n",PAM::ans);
	}
	//print();
	//cerr<<PAM::qaq;
	return 0;
}

Anti-Anti-Hack

尽可能让足够多的回文串的 f a i l fail fail 指到一个节点上。于是把它 Hack 掉了。（3~4s）

f=open("string.in","w")
f.write("300000\n")
a="qwertyuiopasdfghjklzxcvbnm"
s=""
i=0
while len(s)<=300000:
  i=(i+23333)%300000
  c0=a[i%26]
  c1=a[i//26%26]
  c2=a[i//26//26%26]
  c3=a[i//26//26//26%26]
  s=s+'a'+c0+c1+c2+c3+c2+c1+c0+'a'
f.write(s[0:300000])
f.close()

Anti-Anti-Anti-Hack（正确的复杂度）

见 某 SCOI 模拟赛 T3 串（string）【PAM】（在出题人的数据上跑得更慢了）

破壁人五号
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
[SCOI2007]压缩 个人理解分析外加部分程序
07-14
这是一个演示文档，内容是本人对SCOI2007的压缩这一题的个人理解、分析以及部分代码，近似于题解。
某 SCOI 模拟赛 T3 s3mple【生成函数 拉格朗日插值】
破壁人五号的博客
06-07 316
题意 对于序列 aaa，记 viv_ivi​ 为位置距离 aia_iai​ 最近的、比 aia_iai​ 大的数与它的距离（假设 a0a_0a0​ 和 an+1a_{n+1}an+1​ 都为无穷大）。 TTT 次询问，给定 n,xn,xn,x，求有多少长度为 nnn 的排列 ppp 使得 ∑i=1nvi=x\sum\limits_{i=1}^nv_i=xi=1∑n​vi​=x。答案对给定的质数 PPP 取模。n≤200n\leq 200n≤200，T≤10T\leq 10T≤10。 题解 记 FiF_iFi
20200721 T2 s2mple【后缀自动机妙用】
Master.Yi的博客
07-21 316
题目描述 n,Q≤4∗105n,Q\le4*10^5n,Q≤4∗105 题目分析 建 SSS 的 SAM。记询问串为 TTT。 要数 SSS 中本质不同的串中 TTT 的出现次数之和 相当于给 TTT 后面接上一个串 RRR，前面接上一个串 LLL，LTRLTRLTR 出现在 SSS 中，求 L,RL,RL,R 的对数。 假设已经确定了 TR=pTR=pTR=p，考虑怎么数 LLL：以 ppp 为后缀的串，就是在 ppp 对应 fail 树上节点 qqq 子树中的串，加上 qqq 剩下长度对应的串。可以表
某 SCOI 模拟赛 T3 串（string）【PAM】
破壁人五号的博客
06-10 141
题意 对于字符串 SSS 的每个前缀，求其所以回文子串的前缀个数和，本质相同的前缀不重复计数。N≤3×105N\leq3\times 10^5N≤3×105。 题解 建议先阅读 【乱搞】某 SCOI 模拟赛 T3 串（string）【PAM 乱搞】 在 Anti-Hack 的基础上，把 vector 维护每个儿子改成 set 维护每个儿子在 SA 里的 rankrankrank。于是每次出现新节点只需要判断与它 rk 最接近的两个兄弟与它的 LCP，复杂度至多是 O(log⁡n)O(\log n)O(log
某 SCOI 模拟赛 T3 c
破壁人五号的博客
06-12 177
题意 用 nnn 个积木拼一个塔，第 iii 个积木长 aia_iai​。给定 ddd，要求积木 iii 能放在 jjj 上面当且仅当 ai−aj≤da_i-a_j\leq dai​−aj​≤d。求合法的搭积木方案数。答案模 998244353998244353998244353。n≤7×105n\leq 7\times10^5n≤7×105 题解 先排序。 考虑搭好前 i−1i-1i−1 小的积木后，第 iii 个积木能插入进什么地方：它上面是哪个积木不重要，只需要它下面的积木够宽或者它下面没有积木。于是
scoi2006动态最值minmax测试数据
04-25
scoi2006动态最值minmax测试数据
hehelego#WhyNotMarkdown#SCOI2009-windy数-数位DP经典题1
07-25
[SCOI2009] windy数求十进制数$x$的计数, 满足如下约束无前导0相邻数位上的数码之差至少为2首先, 后两条约束和$a,b$没有关系, 于是, 可
soci-backend
04-09
这个小程序运行SCOI Web App的后端API。要求您需要设置一些环境变量才能运行go应用。 如果没有这些设置，它将做出一些（合理的）默认假设： APP_PORT=9000DB_HOST=localhostDB_PORT=3306DB_DATABASE=socidbDB_USER=...
快速 gcd 变换（FGT）
破壁人五号的博客
05-27 511
题目：定义 gcd⁡\gcdgcd 卷积如下：A∗BA*BA∗B 的结果为 Ci=∑[gcd⁡(j,k)=i]AjBkC_i=\sum[\gcd(j,k)=i]A_jB_kCi​=∑[gcd(j,k)=i]Aj​Bk​。求 555555 个长为 N=106N=10^6N=106 的数列卷积的结果。时限 0.6s。 记 (A′)i=∑[i∣j]Aj(A')_i=\sum[i|j]A_j(A′)i​=∑[i∣j]Aj​。则 gcd⁡\gcdgcd 卷积相当于变换后点值相乘。 直接枚举的话复杂度是调和级数 O(
某 SCOI 模拟赛 T3 画图（draw）（无根无标号平面树计数）【卡特兰数 Burnside引理】
破壁人五号的博客
05-31 473
题意 平面上有 nnn 个点，n−1n-1n−1 条只在各自端点处相交的线段将它们连成一棵树。如果两颗树能通过移动顶点位置重合且移动过程中线段只在各自端点处相交，两棵树等价。问：有多少棵两两不等价的有 nnn 个点的树，答案模给定的质数 ppp。多组询问。n≤5×105n\leq 5\times 10^5n≤5×105，108≤p≤10910^8\leq p\leq 10^9108≤p≤109，T≤3T\leq 3T≤3。 题解 首先每一棵顶点数为 n+1n+1n+1 的树可以转化为：圆上有等距的 2n2n
某 SCOI 模拟赛 T4 构造题（construct）（构造 n 元有限域）【群论】
破壁人五号的博客
05-28 421
（20200527） 题意 构造一个 nnn 元有限域。（1<n≤10001<n\leq 10001<n≤1000） 题解 首先只有在 nnn 为质数的整数次幂时才能构造出来（不会证）。记 n=pdn=p^dn=pd。 d=1d=1d=1 时，显然所有运算直接  mod p\bmod pmodp 就行。 d>1d>1d>1 时，考虑把每个数写成 d−1d-1d−1 次多项式 f(x)f(x)f(x) 的形式，系数就对应这个数在 ppp 进制下的每一位，令 x=px=px=
题解 CF1007E Mini Metro（某 SCOI 模拟赛 T3）【DP】
破壁人五号的博客
06-17 343
题意 有 nnn 座火车站（原题是地铁，不过不重要）在一条轨道上依次排列，编号为 1 到 nnn。每个车站在 0:000:000:00 时（初始状态）有 aia_iai​ 个人，每小时增加 bib_ibi​ 个人，容量为 cic_ici​ 人。你还有足够多的火车，每辆车的大小为 kkk。 每个小时的流程大致如下： i:00i:00i:00：第 iii 小时开始； i:30i:30i:30：你可以派出任意多辆火车，它们会从 1 到 nnn 经过每一个站，并在每一个站装尽可能多的人（要么把站台清空、要么把火车
某 SCOI 模拟赛 T1 集合划分（divide）【生成函数 NTT 分治】
破壁人五号的博客
06-26 331
题意 有长为 2n2n2n 的序列 A,BA,BA,B，求有多少个单调不减的序列 CCC，要求 CCC 的每一位均为 AAA 或 BBB 中对应位置上的数，且恰有 nnn 个数来自 AAA。n≤5×104n\leq 5\times 10^4n≤5×104。3s。 题解 设 f[i][j=0/1][k]f[i][j=0/1][k]f[i][j=0/1][k] 表示 DP 确定了前 iii 个数，第 iii 个数是 Ai/BiA_i/B_iAi​/Bi​ ，已经有 kkk 个数来自 AAA 的方案数。直接 DP
某 CSP 模拟赛 T1 金币（coin）
破壁人五号的博客
11-14 277
题意 有一个 n×nn\times nn×n 的网格，每个格子可能是宝藏（T，只会有 kkk 个）、起点（S，只会有一个）、红/绿/蓝色（R/G/B）。玩家从起点出发，当经过红、绿、蓝色格子时，分别花费 r,g,br,g,br,g,b 金币，首次经过一个宝藏格子得到 ttt 金币，收集到 [LX,RX][LX,RX][LX,RX] 个宝藏后，便可以结束游戏。问结束游戏时玩家至多可以有多少金币？（若...
某NOIP模拟赛 T1 字符串（string）
破壁人五号的博客
08-22 257
题目 【描述】 给定两个字符串 s,t，其中s只包含小写字母以及*，t只包含小写字母。 你可以进行任意多次操作，每次选择 s中的一个*，将它修改为任意多个（可以是 0个）它的前一个字符。问是否能将 s修改为t。 有多组数据。 【输入】 第一行一个整数T表示数据组数。 每组数据两行，第一行一个字符串 s，第二行一个字符串 t。 【输出】 每组数据输出一行，如果...
某 SCOI 模拟赛 T1 s1mple【集合幂级数】
破壁人五号的博客
06-07 256
题意 给定 n×nn\times nn×n 的 01 矩阵 MMM，qqq 次询问：给定长度 n−1n-1n−1 的 01 序列 aaa，求有多少 1 到 nnn 的排列 ppp 使得 ∀i∈[1,n−1]Mpi,pi+1=ai\forall i\in[1,n-1]M_{p_i,p_{i+1}}=a_i∀i∈[1,n−1]Mpi​,pi+1​​=ai​。n≤17n\leq 17n≤17，q≤105q\leq10^5q≤105。 题解 先将 MMM 看作一个有向图，1 代表有边。 将 aaa 压缩为 SSS。
某 SCOI 模拟赛 T2 s2mple【后缀自动机】
破壁人五号的博客
06-08 239
题意 定义在确定了模式串 strstrstr 的情况下，字符串 SSS 的权值为 strstrstr 在 SSS 中出现的次数。现给定字符串 SSS，QQQ 次询问：指定 SSS 的某个子串为模式串，求 SSS 的所有本质不同子串的权值和。∣S∣,Q≤4×105|S|,Q\leq 4\times 10^5∣S∣,Q≤4×105。时限 3s 1s。 题解 （官方题解被某神仙 D 了，说是官方题解数据结构学傻了） 建 SSS 的 SAM。记询问串为 TTT。 首先换个思路：如果 TTT 后面接上一个字符串成为
某 SCOI 模拟赛 T3 贪吃的史莱姆（slime）【树形DP】
破壁人五号的博客
06-11 232
题意 给定一棵树，问多少个排列使得每个长为 mmm 的子区间内的点的导出子图是连通的，答案对 109+910^9+9109+9 取模。2≤m≤n≤4002\leq m\leq n\leq 4002≤m≤n≤400。 题解 考虑两种情况： 2m≤n2m\leq n2m≤n； 2m>n2m>n2m>n。 2m≤n2m\leq n2m≤n 此时排列的最前 mmm 个元素一定是一棵子树，中间部分在一条链上，最后 mmm 个元素在另一个子树里。显然以重心为根 DFS 时，如果恰有两棵大小为 mm
[SCOI2011] 糖果
最新发布
08-11
- *1* *3* [【差分约束系统】【SCOI2011】糖果 candy](https://blog.csdn.net/jiangzh7/article/details/8872699)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值