[AHOI2013] 差异 题解

合集 - 题解(82)
1.[JSOI2008]火星人 题解2.[BZOJ3786] 星系探索 题解3.[ZJOI2015]幻想乡战略游戏 题解4.[SDOI2015] 寻宝游戏5.[HNOI2011] 数学作业 题解6.[HDU5396] Expression 题解7.[NOIP2018] 旅行 题解8.[POI2012] Rendezvous 题解9.[BZOJ3037] 创世纪 题解10.[HEOI2014]大工程 题解11.[ZJOI2019] 语言 题解12.[国家集训队] 矩阵乘法 题解13.[THUPC2017] 天天爱射击 题解14.[SCOI2007] 蜥蜴 题解15.[SDOI2015] 星际战争 题解16.[国家集训队] happiness 题解17.[SDOI2016] 数字配对 题解18.[SDOI2009] 晨跑 题解19.[SCOI2007] 修车20.acwing329 围栏障碍训练场 题解21.无限之环 题解22.[BZOJ4350] 括号序列再战猪猪侠 题解23.[TJOI2015] 弦论 题解24.[SDOI2008] Sandy的卡片 题解
25.[AHOI2013] 差异 题解
26.[Ynoi2016] 镜中的昆虫 题解27.[JOISC 2023 Day3] Tourism 题解28.[PA2021] Od deski do deski 题解29.[TJOI2019] 甲苯先生的字符串 题解30.[POI2014] HOT-Hotels 加强版题解31.[Ynoi2015] 盼君勿忘 题解32.[luogu2123] 皇后游戏33.[JXOI2017] 加法 题解34.[luogu1248] 加工生产调度 题解35.[NOIP2018] 赛道修建36.[Ynoi2015] 我回来了 题解37.[HNOI2009] 图的同构计数38.[BZOJ3811] 玛里苟斯 题解39.[BZOJ3569] DZY Loves Chinese II 题解40.[SCOI2016] 幸运数字 题解41.[HDU5603] the soldier of love 题解42.[QOJ8672][PKUSC2024] 排队43.[BZOJ3489] A simple rmq problem44.[HAOI2010] 软件安装 题解45.[IOI2020] 连接擎天树 题解46.[LOJ6669] Nauuo and Binary Tree 题解47.[JOISC2019] 聚会 题解48.[CERC2014] Parades 题解49.[SHOI2017] 摧毁“树状图”50.[NOI2014] 购票 题解51.[WC2018] 通道 题解52.[BZOJ2741][FOTILE模拟赛] L 题解53.[THUSC2015] 异或运算 题解54.[BZOJ4771] 七彩树 题解55.[BZOJ4605] 崂山白花蛇草水 题解56.[BZOJ3600] 没有人的算术 题解57.[WC2014] 紫荆花之恋 题解58.[SDOI2008] 洞穴勘测 题解59.[POJ3237] 树的维护 题解60.[国家集训队] Tree2 题解61.[WC2006] 水管局长 题解62.[BZOJ3514] [Codechef MARCH14] GERALD07加强版 题解63.[luogu4114] Qtree1 题解64.[THUWC2017] 在美妙的数学王国中畅游 题解（内附求导小技巧）65.[BZOJ3159] 决战 题解66.[BZOJ2194] 快速傅立叶之二 题解67.[ZJOI2014] 力 题解68.[BZOJ3451] Normal 题解69.[BZOJ3771] Triple 题解70.[BZOJ3160] 万径人踪灭 题解71.[SDOI2015] 序列统计 题解72.[联合省选 2020A] 组合数问题 题解73.[HDU4625] JZPTREE+[国家集训队] Crash 的文明世界 题解74.[TJOI/HEOI2016] 求和 题解75.[BZOJ5093] 图的价值 题解76.[FJOI2016] 建筑师 题解77.[BZOJ4665] 小w的喜糖 题解78.[BZOJ3622] 已经没有什么好害怕的了 题解79.[BZOJ4671] 异或图 题解80.[BZOJ4833] 最小公倍佩尔数 题解81.[BalticOI 2022] Uplifting Excursion (Day1) 题解82.[JLOI2016] 成绩比较 题解
收起

后缀自动机维护子串公共后缀方便一点，所以直接倒序插入字符串即可。

我们给所有前缀打上标记，然后跑树形 \(dp\)，设 \(sum_i\) 表示第 \(i\) 个点的子树内有多少个前缀，\(ans\) 统计 \(\sum \text{LCP}(T_i,T_j)\)，则有：

\[ans=\sum\limits_{i=1}^{id}\sum\limits_{j\in ison} {sum}_j({sum}_i-{sum}_j) \]

简化求解式发现就是 \(\dfrac{n(n-1)(n+1)}{2}-ans\)，时间复杂度 \(O(n)\)。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e6+5;
struct SAM{
	int id,sm[N],ln[N];
	int tl,tr[N][26],pr[N];
	vector<int>g[N];ll ans;
	SAM(){pr[0]=-1;}
	void cpy(int x,int y){
		for(int i=0;i<26;i++)
			tr[x][i]=tr[y][i];
		pr[x]=pr[y];ln[x]=ln[y];
	}void add(int x){
		ln[++id]=ln[tl]+1;
		int p=tl;tl=id;
		while(~p&&!tr[p][x])
			tr[p][x]=id,p=pr[p];
		if(p<0){
			sm[tl]++;
			return;
		}int u=p,v=tr[p][x];
		if(ln[u]+1==ln[v]){
			pr[id]=v;
			sm[tl]++;
			return;
		}cpy(++id,v);
		ln[id]=ln[u]+1;
		pr[v]=pr[id-1]=id;
		while(~p&&tr[p][x]==v)
			tr[p][x]=id,p=pr[p];
		sm[tl]++;
	}void jb(){
		for(int i=1;i<=id;i++)
			g[pr[i]].push_back(i);
	}void dfs(int x){
		ll sum=sm[x];
		for(int i=0;i<g[x].size();i++){
			int y=g[x][i];dfs(y);
			ans+=sum*sm[y]*ln[x];
			sum+=sm[y];
		}sm[x]=sum;
	}
}sam;char s[N];int n;
int main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	cin>>s;n=strlen(s);
	for(int i=n-1;~i;i--)
		sam.add(s[i]-'a');
	sam.jb();sam.dfs(0);
	cout<<(ll)(n-1)*n*(n+1)/2-sam.ans*2;
	return 0;
}//man!what can I say!
posted @   长安一片月_22  阅读(8)  评论(0)  编辑  收藏  举报
编辑推荐：
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列：向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
· Linux系列：如何用heaptrack跟踪.NET程序的非托管内存泄露
· 开发者必知的日志记录最佳实践
阅读排行：
· TypeScript + Deepseek 打造卜卦网站：技术与玄学的结合
· Manus的开源复刻OpenManus初探
· 写一个简单的SQL生成工具
· AI 智能体引爆开源社区「GitHub 热点速览」
· C#/.NET/.NET Core技术前沿周刊 | 第 29 期（2025年3.1-3.9）
点击右上角即可分享
微信分享提示