Steaunk-CSDN博客

目标给定一个字符串 SSS，求出以每一个点为中心的最长的回文子串的长度。预处理首先对于字符串如 abaccaa，将其两两字符之间插入一个#，并在首尾各插入一个不常用的字符（避免边界问题），如变为?#a#b#a#c#c#a#a#$，这样的一个新串在处理时就无需考虑长度为偶数的回文子串，如回文子串 acca 在新...

2019-04-23 13:35:20 230

原创 CodePlus 第五次网络赛掐指会算

失踪人口暂时回归。临近 NOIP 了，退役选手准备打一打 Div.2 来练练手（应该不是天气冷了，没衣服穿了）游戏体验差，OJ 又和第一次一样卡了半天。T1 我有矩阵，你有吗？根据异或的性质，不难发现 AAA 矩阵的每行每列最多只能异或一次。所以我们可以假设 AAA 矩阵的第一行是否被异或了，然后把所有状态递推出来，最后判断一下是否符合题设。#include <bits/st...

2018-11-03 19:10:23 349

原创「Codeforces 632F」Magic Matrix

称一个矩阵为魔法矩阵，当且仅当满足一下三点： 1. ∀1≤i,j≤n,ai,j=aj,i∀1≤i,j≤n,ai,j=aj,i\forall 1 \le i,j \le n, a_{i,j}=a_{j,i} 2. ∀1≤i≤n,ai,i=0∀1≤i≤n,ai,i=0\forall 1 \le i \le n, a_{i,i} = 0 3. ∀1≤i,j,k≤n,ai,j≥ma...

2018-09-14 21:50:58 650

原创「BZOJ 1798」[Ahoi2009]维护序列seq

给定长度为 NNN 的序列，执行以下有三种模 PPP 意义下的操作之一 MMM 次： 1. 区间乘 ccc； 2.区间加 ccc；3.查询区间和。 1≤P≤109,1≤M,N≤105,0≤c≤1091≤P≤109,1≤M,N≤105,0≤c≤1091\le P\le 10^9 ,1\le M,N\le 10^5,0\le c\le 10^9直接线段树，维护两个标记，一个维...

2018-07-14 11:58:01 228

原创「BZOJ 3343」教主的魔法

给一个长度为 NNN 的正整数序列，初始值都小于 100010001000。你有两种操作： 1. 使一个区间内的数都增加 WWW； 2. 查询一段区间内大于等于 CCC 的数有多少个；执行 QQQ 此操作。 N≤1000000N≤1000000N \le 1000000 ， Q≤3000Q≤3000Q \le 3000 ， 1≤W≤10001≤W≤10001...

2018-07-11 21:06:00 257

原创「Codeforces 871E」Restore the Tree

Petya 有一个由 nnn 个点构成的树。但 Petya 只记得这颗树的 kkk 个关键点与所有点 nnn的距离. 你的任务是构造一种树的方案使其满足 Petya 记忆中的树，若无法构成则输出 −1−1-1。 2 ≤n≤ 30000,1 ≤ k ≤ min(200, n)2 ≤n≤ 30000,1 ≤ k ≤ min(200, n)2 \le n \le 30000,1...

2018-06-30 22:47:49 535

原创 K-D Tree 模板

#include &lt;bits/stdc++.h&gt;#define Min(_A, _B) (_A &lt; _B ? _A : _B)#define Max(_A, _B) (_A &gt; _B ? _A : _B)#define Abs(_A) (_A &gt; 0 ? _A : -(_A))struct pair{ int x, y; } a[500010];stru...

2018-06-21 11:31:41 313

原创「BZOJ 5379」Tree

HNOI 2018 省队集训 Tree假设没有操作 111，那显然就是把 dfs 序求出来，子树区间加，子树区间查询，用线段树维护，O(nlogn)O(nlog⁡n)O(n\log n)。考虑有操作 111：以 111 为根，建线段树。假设此时的根变为 rrr。然后发现求以 rrr 为根的 LCA(u,v)LCA(u,v)\text{LCA}(u,v)，相当于求 LCA(u,r...

2018-06-20 22:06:31 378

原创「Codeforces 989D」A Shade of Moonlight

给 nnn 个长度为 lll 且互不相交的开区间 (xi,xi+l)(xi,xi+l)(x_i, x_i+l)，每个区间有一个移动速度 vvv，v∈{1,−1}v∈{1,−1}v \in\{1,-1\}。你可以在某一时刻给所有区间同时增加一个速度 www，要求满足 |w|≤wmax|w|≤wmax|w| \le w_{max}。求有对多少对 (i,j)(i&amp;amp;amp;lt;j)(i,j)(i&amp;amp;amp;...

2018-06-12 10:08:32 508

原创「BZOJ 1972」[SDOI2010]猪国杀

这什么鬼鬼题呀大概就所有东西都封装成函数，不然我根本调不动。恰好两百行！#include &lt;cstdio&gt;#include &lt;cstdlib&gt;#include &lt;list&gt;int n, m, id[12], yes[12] = {1}, HP[12], FP;bool arm[12], master[12];typedef std::lis...

2018-05-23 19:56:10 374 1

原创「Codeforces 280D」k-Maximum Subsequence Sum

给定一个长度为 nnn 的序列 AAA，会进行 mmm 次操作，每次操作为以下两种操作之一： 1. 将 aiaia_i 赋值为 valvalval； 2. 询问区间 [l,r][l,r][l,r] 的最多 kkk 段不相交的序列之和的最大值。（1≤n,m≤1051≤n,m≤1051 \le n,m \le 10^5, |ai|,|val|≤500|ai|,|val|≤50...

2018-05-13 00:21:10 333

原创「51Nod 1824」染色游戏

Description很奇妙，FMT 有趣题。首先可以一眼（或两眼）看出（或推出）式子：f(t)≡∑x=0t(tx)rxbt−xmod2f(t)≡∑x=0t(tx)rxbt−xmod2f(t)\equiv \sum_{x=0}^{t}{t \choose x} r_x b_{t-x} \mod 2这好像是个卷积，FFT？？？然后就想歪了。根据 Lucas 定理，我们有 (t...

2018-04-27 09:12:16 351 1

原创「BZOJ 1975」[SDOI2010]魔法猪学院

Description求 kkk 短路，最短路+可持久化可并堆，O(mlogm+klogk+nlogm)O(mlog⁡m+klog⁡k+nlog⁡m)O(m\log m + k \log k + n\log m)。又是一道论文题。#include &lt;cstdio&gt;#include &lt;cstring&gt;#include &lt;algorithm&gt;...

2018-04-26 15:45:29 339

原创「BZOJ 3181」 [Coci2012]BROJ

求最小质因子等于 ppp 的第 nnn 小的正整数。ppp 一定是质数。若答案超过 10910910^9 则输出 000。(1≤n,q≤1091≤n,q≤1091\le n,q \le 10^9) 论文题。《2014信息学奥林匹克中国国家队候选队员论文——寻找k优解的几种方法》，俞鼎力。 #include <cstdio>#include <cstring...

2018-04-26 07:56:44 306

原创「BZOJ 3511」土地划分

这题好简单好简单呀！不难发现，这是一个最小割模型。首先记一个 sumsumsum = ∑N−1i=2Ai+Bi∑i=2N−1Ai+Bi\sum_{i=2}^{N-1}A_i+B_i，所以 ans=sum−minCutans=sum−minCutans = sum -minCut，我们需要做的就是最小化 minCutminCutminCut（最小割），来最大化 ansansans（答案）。...

2018-04-18 20:03:21 421 3

原创 ZKW费用流模板

以 LOJ 102 为例#include <cstdio>#include <cstring>#define Min(_A, _B) (_A < _B ? _A : _B)#define R registerint F(){ R int x; R char ch; while(ch = getchar(), ch < '0' |...

2018-04-02 09:18:20 244

原创「LOJ 6308」「雅礼国庆 2017 Day1」Mod

LOJ传送门这道题居然让我想到了 jury_2 的集训队论文，其实确实还是有异曲同工之妙的。用线段树维护区间和以及区间最大值。

2018-03-28 21:25:18 347

原创「Codeforces 827D」Best Edge Weight

给定一个点数为 nn，边数为 mm，权值不超过 10910^9 的带权连通图，没有自环与重边。现在要求对于每一条边求出，这条边的边权最大为多少时，它还能出现在所有可能的最小生成树上，如果对于任意边权都出现，则输出 −1-1。 (2 ≤ n ≤ 2×105,n − 1 ≤ m ≤ 2×1052 \le n \le 2\times 10^5, n - 1 ≤ m ≤ 2\times 10^

2018-03-09 23:46:35 504

原创「Codeforces 464D」World of Darkraft - 2

RomaRoma\mathrm{Roma} 在游戏“World of Darkraft”（理论上应该是 World of darkcraft，MineCraft 的一个版本）找到一个新角色。 RomaRoma\mathrm{Roma} 有 kkk 种装备，一开始每种装备各有 111 个，且每种装备的初始等级均为 111。游戏中可以靠打怪来获取新装备，总共有 nnn 只怪兽，每打赢 11...

2018-03-08 19:04:40 328

原创 SAM模板

以 luogu P3804 为例感觉我对 SAM 理解的不是很透彻。#include <cstdio>#include <cstring>#define R register#define Null b#define Max(_A, _B) (_A > _B ? _A : _B)char s[1000010];struct Data{ int len, size, In;

2018-03-02 22:20:30 505

原创「Codeforces 940F」Machine Learning

初始时有一个长度为 nn 数字串。有以下 qq 个操作： 1.询问。定义 cic_i 为区间 [l,r][l,r] 中 ii 出现的次数，求出 Mex{c0,c1,...}\mathrm{Mex} \{c_0,c_1,...\}，其中 Mex\mathrm{Mex} 表示集合中未出现的最小的正整数。 2.修改。将第 pp 个数改成 xx。 (1≤n,q≤100 000)(1\le

2018-02-24 22:26:27 427

原创 Codeforces Round #462 (Div. 2)

我好菜呀，居然还在打 div.2… 比赛传送门A. A Compatible Pair 给两个集合 A,BA,BA,B，Tommy 先从 AAA 中删除一个数，Banban 再从 A,BA,BA,B 中各选一个数 a,ba,ba,b， Tommy 想使 a×ba×ba\times b 尽可能小，Banban 想使 a×ba×ba\times b 尽可能大。求在两者都选择最...

2018-02-15 00:17:30 585

原创「Codeforces 915G」Coprime Arrays

CF915G#include <cstdio>#include <cstring>#define R registerconst int Mod = 1e9 + 7;int n, k, tot, P[149010], Ans, mu[2000010], pow[2000010], d[2000010];bool vis[2000010];int Pow(...

2018-01-16 10:46:52 632

原创「51Nod 1868」彩色树

[算法马拉松31C](http://www.51nod.com/contest/Problem.html#!problemId=1868&contestId=43)[51Nod 1868](http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1868)

2018-01-15 21:06:48 363

原创「Codeforces 809D」Hitchhiking in the Baltic States

CF 809D 题意简述：给你长度为 nnn 的序列，序列中的每个元素 iii 有一个区间限制 [li,ri][li,ri][l_i,r_i]，你从中选出一个子序列，并给它们标号 xixix_i，要求满足 ∀i<j∀i<j\forall i < j，xi<xjxi<xjx_i < x_j，且 ∀i∀i\forall i，xi∈[li,ri]xi∈[li,ri...

2018-01-12 08:08:29 675

原创「THUSC 2016」成绩单

BZOJ-4897 LOJ-2292#include #include #define Min(_A, _B) (_A #define Max(_A, _B) (_A > _B ? _A : _B)#define cmin(_A, _B) (_A > _B ? _A = _B : 0);#define R registerint n, a, b, w[55], c[55], f[

2018-01-11 16:13:20 1115

原创数论三合一

Abstract本文将探讨与下列三个同余方程有关的问题，并进行一些数论相关知识的拓展。ax≡bmodpax≡bmodpax \equiv b \mod p ax≡bmodpax≡bmodpa^x \equiv b \mod p xa≡bmodpxa≡bmodpx^a \equiv b \mod p 分别求出 xxx 的最小值。注：如非特别说明，则本文中所有的 ppp 均表示为素数...

2018-01-06 13:55:56 828

原创欧拉筛法与积性函数

线性筛欧拉筛本质上是在对每一个数找到它最小的质因数，然后把它筛除，复杂度 O(n)O(n) 。本文中所有的 pp 都表示一个质数。PROOF： From the code below we have that pj∣ap_j \mid a. Let n=apk (k>j，pk>pj)n =ap_k\ (k > j，p_k > p_j). ∵pj∣a\because p_j \mi

2017-12-21 08:39:13 664

原创 LCT模板

以【清华集训2016】温暖会指引我们前行为例 UOJ 274

2017-12-20 21:43:38 369

原创 FFT模板

《算法导论》上写的是真不错。FFT递归版，常数有点大呀。#include <cstdio>#include <cmath>#include <cstring>#define R register#define db long double#define ll long longconst int Length = 1 << 1...

2017-12-16 15:41:19 253

原创 FWT模板(loading)

ci=∑j⊗k=iajbk=∑0≤k<2n∑0≤j<2n[j⊗k=i]ajbkci=∑j⊗k=iajbk=∑0≤k<2n∑0≤j<2n[j⊗k=i]ajbkc_i = \sum_{j \otimes k = i} a_j b_k= \sum_{0 \le k < 2^n} \sum_{0\le j < 2^n} [j \otimes k = i] a_j b_k &nbsp...

2017-12-14 16:59:19 348

原创「POJ 3693」Maximum repetition substring

题意：求一个串的最大重复次数的重复子串，要求输出字典序最小。后缀数组，可以看罗穗骞的论文《后缀数组——处理字符串的有力工具》。#include #include #define Max(_A, _B) (_A > _B ? _A : _B)#define Min(_A, _B) (_A #define R registerchar s[100010], t[100010];

2017-11-30 22:04:15 350

原创后缀排序

以 luogu P3809 为例#include #include #define R registerconst int MaxN = 1000010;int SA[MaxN], Rank[MaxN], Height[MaxN], x[MaxN], sum[MaxN], y[MaxN];char s[MaxN];bool cmp(R int i, R int j, R int k

2017-11-28 09:36:22 898

原创 Trie图

Trie图是AC自动机的优化。以 luogu P3808 为例

2017-11-27 13:15:06 297

APIO2016课件

mysql 8.0 修改 root 密码