Nickwzk_莫队,FFT,数论-CSDN博客

原创 2019ICPC南京网络赛 E.K Sum 反演+杜教筛

E.K SumProblemDef. fn(k)=∑l1=1n∑l2=1n...∑lk=1n(gcd(l1,l2,...,lk))2 f_n(k) = \sum_{l_1=1}^n \sum_{l_2 = 1}^n ...\sum_{l_k=1}^n (gcd(l_1, l_2, ...,l_k))^2fn(k)=l1=1∑nl2=1∑n...lk=1∑n(gcd(l1,l2...

2019-09-22 16:29:59 210

原创 2019牛客暑期多校训练营（第一场） D.Parity of Tuples (FWT)

D.Parity of Tuples题目大意:给一个n*m矩阵，行向量为viv_ivi, 求⊕x=02k(count(x)∗3x mod(109+7))\oplus_{x=0}^{2^k}(count(x) * 3^x\space mod(10^9+7))⊕x=02k(count(x)∗3x mod(109+7))其中count(x)count(x)count(x)表...

2019-07-19 02:21:47 376

原创 NOIP2017

//Cease to struggle and you cease to live.Day0从学校赶到酒店，入住已经是晚上9点了。感冒加重了，有些头疼。Day1配好vim发现gdb有点问题，开考10分钟后工作人员告诉我说DEV的MinGW有问题，换成codeblocks的MinGW就好了。折腾半天才开始看题，感觉药丸的节奏。 T1 竟然是数学题，画风有点不对。 T2 竟然是模拟题，画风果

2017-11-25 21:09:37 452 3

原创 SCOI2017

Day0有点紧张呢。繁花终是落，迷途人将返。 SCOI加油，希望这不是这条路的终点吧。 ……未完待续

2017-04-07 23:11:36 1241

原创数论函数变换

对于各种积性函数，都可以通过两种方法进行计算：分解质因数线性筛那么，我们就可以得到一些常用的积性函数值。typedef pair<int, int> P;int prime[MAXN], cnt, phi[MAXN], mu[MAXN], mnsum[MAXN], a[MAXN];bool isp[MAXN];P f[MAXN];//约数和void Euler(int n) {

2017-03-27 11:09:31 681

原创 bzoj4503两个串快速傅里叶变换（FFT）

题目传送门：bzoj4503题目大意：给出两个串S和T，T中有通配符（可以与任意字符匹配），求T在S中匹配的次数及位置。题目一看，怎么做啊，kmp？AC自动机？好像很困难啊。但是，我们转换一下思路。我们定义ck=Σki=0(Si−Tk−i)2∗Tk−jc_k = \Sigma_{i = 0}^{k} (S_i - T_{k - i} ) ^ 2*T_{k - j}。当且仅当ck==0c_k ==

2017-02-22 19:40:58 1120

原创快速傅里叶变换（FFT）

快速傅里叶变换FFT是用来计算离散傅里叶变换（DFT）及其逆变换（IDFT）的快速算法。两个n次多项式直接相乘所需的时间为O（n2{n}^{2}），而FFT可以将其复杂度降低为O（nlogn{n}logn）。令A(x) = ∑n−1j=0ajxj\sum_{j = 0}^{n-1}a_jx^j B(x) = ∑n−1j=0bjxj\sum_{j = 0}^{n-1}b_jx^j

2017-02-21 22:04:48 1495

原创 bzoj3052糖果公园树上莫队

题目传送门: http://www.lydsy.com/JudgeOnline/problem.php?id=3052最近不太想写博客了，退役前还是写点学的东西吧。这个坑从我会莫队开始就挖在这里了。。莫队： http://blog.csdn.net/nickwzk/article/details/52954097 可是在树上怎么莫队呢？这有两种非常妙的方法，dfs序和王室联邦分块（bz

2017-02-16 21:22:47 984

原创 bzoj2982: combination 组合数卢卡斯定理

Time Limit: 1 Sec Memory Limit: 128 MBDescription LMZ有n个不同的基友，他每天晚上要选m个进行[河蟹]，而且要求每天晚上的选择都不一样。那么LMZ能够持续多少个这样的夜晚呢？当然，LMZ的一年有10007天，所以他想知道答案mod 10007的值。(1<=m<=n<=200,000,000)Input 第一行一个整数t，表示有t组数据。(

2016-11-14 10:24:08 712

原创 2127: happiness 最小割

Time Limit: 51 Sec Memory Limit: 259 MBDescription 高一一班的座位表是个n*m的矩阵，经过一个学期的相处，每个同学和前后左右相邻的同学互相成为了好朋友。这学期要分文理科了，每个同学对于选择文科与理科有着自己的喜悦值，而一对好朋友如果能同时选文科或者理科，那么他们又将收获一些喜悦值。作为计算机竞赛教练的scp大老板，想知道如何分配可以使得全班的喜悦

2016-11-10 10:50:38 729

原创 bzoj1411: [ZJOI2009]硬币游戏找规律神题

Time Limit: 10 Sec Memory Limit: 162 MBDescription Orez很喜欢玩游戏，他最近发明了一款硬币游戏。他在桌子的边缘上划分出2*n个位置并按顺时针把它们标号为1，2，……，2n，然后把n个硬币放在标号为奇数的位置上。接下来每次按如下操作：在任意两个硬币之间放上一个硬币，然后将原来的硬币拿走；所放硬币的正反面由它两边的两个硬币决定，若两个硬币均为正面

2016-10-29 17:36:28 1529 1

原创 bzoj1251: 序列终结者 splay

Time Limit: 20 Sec Memory Limit: 162 MB Submit: 3646 Solved: 1527 [Submit][Status][Discuss] Description网上有许多题，就是给定一个序列，要你支持几种操作：A、B、C、D。一看另一道题，又是一个序列要支持几种操作：D、C、B、A。尤其是我们这里的某人，出模拟试题，居然还出了一道这样的，真是

2016-10-28 13:10:09 1110

原创 bzoj2565: 最长双回文串 manacher

Time Limit: 10 Sec Memory Limit: 128 MBDescription 顺序和逆序读起来完全一样的串叫做回文串。比如acbca是回文串，而abc不是（abc的顺序为“abc”，逆序为“cba”，不相同）。输入长度为n的串S，求S的最长双回文子串T,即可将T分为两部分X，Y，（|X|,|Y|≥1）且X和Y都是回文串。Input 一行由小写英文字母组成的字符串S。

2016-10-28 11:00:40 1081

原创 bzoj2038: [2009国家集训队]小Z的袜子(hose) 莫队

Time Limit: 20 Sec Memory Limit: 259 MB Description 作为一个生活散漫的人，小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天，小Z再也无法忍受这恼人的找袜子过程，于是他决定听天由命…… 具体来说，小Z把这N只袜子从1到N编号，然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双，甚至不在意两只袜子是否一左一右，

2016-10-28 10:36:12 789

Nickwzk的博客