计数
文章平均质量分 67
YiPeng_Deng
这个作者很懒,什么都没留下…
展开
-
CF1528F AmShZ Farm
CF1528F我们要求∑j≤icj≥i\sum_{j\le i}c_j\ge i∑j≤icj≥i的n!∏ci!\frac{n!}{\prod c_i!}∏ci!n!的和。打表发现答案为(n+1)n−1(n+1)^{n-1}(n+1)n−1,考虑这是一个n+1n+1n+1个点的无根树的方案,考虑怎么将答案与一棵树构造一个双射。有一个很神奇的构造,我们有一个以1为根的大小为n+1n+1n+1的有标号无根树,我们在上面dfs,先走编号小的儿子,将走过每一个点的儿子个数写成一个序列,可以得原创 2021-07-10 15:40:10 · 181 阅读 · 0 评论 -
arc120简要题解
arc120A:每一次加的都会变成最大值,因此简单计算即可。B:最后的对角线上都相同。C:对于ai+i=aj+ja_i+i=a_j+jai+i=aj+j的一一匹配即可。D:考虑111到nnn的和n+1n+1n+1到2n2n2n的两两匹配最大,设1−n1-n1−n为0,n+1−2nn+1-2nn+1−2n为1,贪心匹配即可。E:首先二分答案,然后考虑从前往后DP,设f[i],g[i]f[i],g[i]f[i],g[i]分别表示iii先往左走并与i−1i-1i−1匹配,回头的时间,以及iii先往原创 2021-05-27 22:50:28 · 353 阅读 · 0 评论 -
JZOJ6912. 【2020.12.01提高组模拟】数论(math)
Descriptionn,m≤1e7n,m\le1e7n,m≤1e7Solution后面那个互质数之和由于对称性可以得到是ϕ(n)n/2\phi(n)n/2ϕ(n)n/2,需要特判n=1n=1n=1然后把ϕ(ij)\phi(ij)ϕ(ij)拆成ϕ(i)ϕ(j)(i,j)ϕ((i,j))\phi(i)\phi(j)\frac{(i,j)}{\phi((i,j))}ϕ(i)ϕ(j)ϕ((i,j))(i,j),相当于是一个容斥掉1p−1\frac{1}{p-1}p−11然后直接枚举(i,j)(i原创 2020-12-01 21:37:48 · 260 阅读 · 0 评论 -
【牛客挑战赛45】G.致两千年后的你
DescriptionSolution首先一个集合SSS可以把gcd(a1,a2,a3,...,P)gcd(a_1,a_2,a_3,...,P)gcd(a1,a2,a3,...,P)的倍数拼出。因此我们把P,x,aiP,x,a_iP,x,ai分解质因数,把aia_iai根据每一个质因数的次数是否大于xxx的次数作为二进制数,然后高位后缀和一下,就可以知道某个TTT的倍数的个数。考虑转移kkk次对于f0(S)f_0(S)f0(S)的贡献实际上是knk∣S∣\frac{k^n}{k^{原创 2020-11-15 21:42:53 · 236 阅读 · 0 评论 -
agc039F Min Product Sum
agc039F Min Product SumSolution考虑行和列分别按照最小值排序,那么Ai,jA_{i,j}Ai,j应该是这样的:111111111111111111111122222111122333为若干个LLL型。每个位置乘在一起就是贡献。考虑方案数为这样的:111122333111122333222222333333333333表示的是那个位置填的数要>=上面的矩阵中那个位置的数。考虑按照1−K1-K1−K从左上角往下加入一行一列,但是并原创 2020-10-15 21:04:12 · 292 阅读 · 0 评论 -
JZOJ6702. 【2020.06.07省选模拟】仙人掌(cactus)
Descriptionn<=1e5,m<=2e5n<=1e5,m<=2e5n<=1e5,m<=2e5Solution一开始对于矩阵的行列式啥的以为跟矩阵树有关系,但是又不好套,就以为是一些我所不知道的神仙知识点,于是就没有直接往行列式的定义上去想了。考虑行列式的定义来计算——选择一个排列,计算它的逆序对个数,加在一起(因为所有权值都是1)。那么再考虑排列在图上的贡献,相当于是选了n条边,i连向p[i],这样会形成若干个环。但是由于这些p在原矩阵中是交叉的,原创 2020-06-06 19:34:13 · 401 阅读 · 0 评论 -
JZOJ6693. 【2020.06.05省选模拟】紫色彼岸樱推迟绽放
Description幽幽子饿了,妖梦需要给幽幽子准备食物。有 T 天,每天幽幽子划分成了 k 个时段,妖梦需要安排每一天的日程。第 i 天妖梦准备了 D+i-1 道菜,每道菜有无数个。第 1 个时段是早餐,幽幽子会选择 L 道不同的菜吃。接下来 k-1 个时段,每个时段可以选择 D+i-1 道菜中的一道吃或者选择 A 个活动中的一个参加,但是出于健康考虑,幽幽子不能连续两个时段同时吃菜。k 和 A 是幽幽子事先决定好的,她给出了 Q 个询问,每次询问给出 L,D,T,问这 T 天每一天的安排的方原创 2020-06-06 08:12:12 · 236 阅读 · 0 评论 -
Min_25筛学习小计
简单介绍由Min_25在他的博客中介绍的做法,又称Min_25筛。用于求积性函数f(n)f(n)f(n)的前缀和,其中要求f(p)f(p)f(p)可以表示成多项式,并且f(pk)f(p^k)f(pk)可以快速算出。能够在O(n0.75log)O(\frac{n^{0.75}}{log})O(logn0.75)的时间内算出。基本思路使用DPDPDP先求出所有素数的f(p)f(p)f(p)的前缀和,然后再通过枚举最小的质因子计算出合数的前缀和。求素数的前缀和下面的所有除法都为整除。我原创 2020-06-03 19:16:59 · 298 阅读 · 0 评论 -
JZOJ6675. 【2020.05.30省选模拟】交通网络
Descriptionn<=5e5n<=5e5n<=5e5SolutionO(n^4)直接放弃思考矩阵树即可过n<=80n<=80n<=80,获得23分。O(n^2)考虑先枚举一些特殊边,将n个点分成m个联通块,再用prufer序列计算这个完全图的方案数:nm−2∏i=1msin^{m-2}\prod_{i=1}^ms_inm−2∏i=1msi考虑prufer序列中的一个值代表了一条边,且长度为m−2m-2m−2,那么连边的时候就可以有n个顶点原创 2020-06-03 11:38:10 · 222 阅读 · 0 评论 -
JZOJ6679. 【2020.06.02省选模拟】数
Descriptionn<=1e6,ai<=1e9n<=1e6,a_i<=1e9n<=1e6,ai<=1e9,输出∑i=0nfi(n)∗(1e9+7)i(mod 232)\sum_{i=0}^{n}f_{i}(n)*(1e9+7)^i (mod\ 2^{32})∑i=0nfi(n)∗(1e9+7)i(mod 232)Solution完全没有想到可以直接将组合数用多项式表示,来推式子,经典姿势++原式aia_iai的系数显然是(1−x原创 2020-06-03 08:59:23 · 253 阅读 · 0 评论 -
burnside定理学习小计
基本概念置换:对一个集合的映射,简单来说就是重排列。一个集合SSS经过映射a[1..n]a[1..n]a[1..n]后得到S′S'S′的即S′[1]=S[a[1]],S′[2]=S[a[2]]....S'[1]=S[a[1]],S'[2]=S[a[2]]....S′[1]=S[a[1]],S′[2]=S[a[2]]....不动点:如果一个集合SSS在通过置换aaa后生成的S′S'S′与SS...原创 2019-11-30 16:58:18 · 620 阅读 · 1 评论 -
JZOJ6423. 【NOIP2019模拟11.11】画
DescriptionSolutionNOIP前做NOI模拟???(神奇的题目)DP+DP+DP???首先肯定是考虑如果没有边的限制怎么做啦。有一种显然的做法是数位DP,记一个状态S表示哪些顶住了上限。但是这样最优是只能做到64*3n的。我们没有运用上异或的性质。还有一种更优的做法,枚举二进制下的一个位置表示在这个位置以上的所有位置所有点都顶住了上限(显然要先满足前面这些位置异...原创 2019-11-12 11:52:18 · 245 阅读 · 0 评论 -
【简单计数知识2】JZOJ6405. 【NOIP2019模拟11.04】c
Descriptionn<=1e6,m<1e9+7n<=1e6,m<1e9+7n<=1e6,m<1e9+7Solution刚开始看到矩阵求逆后发现连裸的矩阵求逆我都不会(其实我去年应该是学过的。。。),只会一发n6的暴力对n2个点进行高斯消元。实际上,矩阵求逆有一种很好做的n3高斯消元。对于矩阵AAA,将它变为单元矩阵III,同时根据AAA的操作...原创 2019-11-06 22:25:30 · 157 阅读 · 0 评论 -
反演学习小计
参考资料炫酷反演魔术,打开了我新世界的大门,原来这些反演都是一个东西!!!二项式反演f(n)=∑i=0nCnig(i)f(n)=\sum_{i=0}^{n}C_n^ig(i)f(n)=i=0∑nCnig(i)g(n)=∑i=0n(−1)n−iCnif(i)g(n)=\sum_{i=0}^{n}(-1)^{n-i}C_n^if(i)g(n)=i=0∑n(−1)n−iCnif(i)推...原创 2019-11-06 21:59:20 · 399 阅读 · 0 评论 -
JZOJ6407. 【NOIP2019模拟11.05】小 D 与随机
Descriptionn<=5000n<=5000n<=5000Solution神仙计数题。考虑树形DP。为了DP的方便,有一种很骚的操作,因为是对于每一个恰好为c个好点的状态求kck^ckc,我们为了复杂度不好记录c。但是我们可以将它转化为至少c个好点。kc=∑i=0cCci∗(k−1)ik^c=\sum_{i=0}^{c}C_c^i*(k-1)^ikc=i=...原创 2019-11-05 22:45:07 · 618 阅读 · 0 评论 -
【CSP-S2019模拟】11.01比赛模拟
日均挂50分比赛思路传送门T1(森林):如果是加的话就倒着变成删。排序后贪心就好了。T2(玩具):理解题意之后发现并不是很好做,先做了T3。推了一波DP,发现中间的标号会算重,又加了一维,原本n3的就变成n4了。T3(飘雪圣城):求边数。刚开始以为是可持久化线段树。打的时候发现直接用树状数组就可以维护了,然后就又短又水了。赛后消化T1因为会加入一些房子,数组要开大两倍。。。还好...原创 2019-11-01 16:26:34 · 361 阅读 · 0 评论 -
【简单计数知识】JZOJ6395. 【NOIP2019模拟2019.10.28】消失的序列
Description传送门求有多少个n排列可以被一下的栈排序排成升序,其中有一个元素已经确定每一次可以对序列进行入栈和出栈两种操作。入栈操作:如果序列非空,那么三叶会取出序列的第一个元素,把它扔进栈的顶部。出栈操作:如果栈非空,那么三叶可以取出栈的顶部,把这个元素接在结果序列的末尾。注意结果序列和原序列是不一样的,也就是出栈操作不是把元素接在原序列后。如果所有元素都经过了一次入栈...原创 2019-10-28 22:46:28 · 288 阅读 · 0 评论 -
JZOJ6024. 【GDOI2019模拟2019.2.16】网格
Description多组数据,1<=n,T<=5e5Solution通过找规律 我们发现我们可以将答案分为左上到右下,右上到左下两种吗,并且既不重复,也不遗漏。例如:1100000110000000100000001110000111000011100000001这个状态,最后必然合法。对于这个左上到右下的状态里面的每一个1构成的正方形,又分别是一个右上到...原创 2019-02-16 20:23:48 · 316 阅读 · 0 评论 -
JZOJ6021. 【GDOI2019模拟2019.2.15】车
DescriptionData ConstraintSolution容斥!首先如果没有对角线的限制,方案显然是n!容斥1:将不能放改为必须放,枚举集合,容斥系数为(-1)size容斥2:两条对角线必须有,改为总方案减去一条没有的方案加上两条没有的方案。容斥3:一条对角线没有的方案,改为总方案减去有若干个的方案,容斥系数(-1)cnt,组合数计算一下。容斥4:两条对角线没有...原创 2019-02-15 22:36:03 · 193 阅读 · 0 评论 -
JZOJ3206. 【HNOI模拟题】Axonometric Projection
Description让我们来考虑一个单位立方体建成的模型。这个建筑的底是一个n m的单位正方形网格。在每个正方形上面,堆着若干个(可能是0)个单位立方体。每个立方体属于其中一个立方体堆。给出了一个建筑的左视图和正视图。请计算有多少种建筑,符合给出的左视图和正视图。答案可能很大,只要返回它除以10^9 + 9的余数即可。Input第一行是整数n。第二行描述了建筑的左视图。第i个数表示了由上...原创 2019-01-29 08:04:03 · 302 阅读 · 0 评论 -
【THUWC2019模拟2019.1.18】Counting
Description羽月最近发现,她发动能力的过程是这样的:构建一个 V (V<=100)个点的有向图 G,初始为没有任何边,接下来羽月在脑中构建出一个长度为 E 的边的序列,序列中元素两两不同,然后羽月将这些边依次加入图中,每次加入之后计算当前图的强连通分量个数并记下来,最后得到一个长度为E 的序列,这个序列就是能力的效果。注意到,可能存在边的序列不同而能力效果相同的情况,所以羽月...原创 2019-01-18 22:55:05 · 734 阅读 · 0 评论 -
【PKUWC2019模拟2019.1.15】Permutation
DescriptionData Constraint一个简单的计数。稍微化简一下发现我们需要求(i-1)/2*(i-2)!/(i-y)!考场上的时候想到了这里,也想到了NTT,就是不会转换!!!??我们可以这样卷:设f[i]=(i-1)/2*(i-2)!,g[i]=1/(n-i)!那么那个式子就是sigma(f[i]*g[j])[i+j==n+y]妙啊妙…但...原创 2019-01-16 22:48:15 · 180 阅读 · 0 评论 -
JZOJ5968.【NOIP提高组11.24】电竞选手
PROBLEMSOLUTION不难发现如下的性质: 1.必须从小到大按顺序消除所有数 2.每一块(这里表示所有数值相等的集合)最后消剩一个数之前的二元组操作可以在之前的任意一个时刻进行。先考虑所有的a相等时的情况:设答案为f[n],则有f[n]=f[n-1]*C(n,2) 对于每一块的答案我们可以单独考虑:我们枚举放在前面的这一块的前面的数并在...原创 2018-11-24 16:58:40 · 201 阅读 · 0 评论 -
JZOJ5803.【2018.8.12省选模拟】girls(三元环)
PROBLEM有0至n-1个元素,给出m对元素的冲突,给出A,B,C,定义一个满足i&amp;lt;j&amp;lt;k的三元组(i,j,k)的贡献为A∗i+B∗j+C∗k,求所有没有冲突的三元组的贡献和SOLUTION考虑容斥,设F(i)表示每个三元组考虑i组冲突的总和Ans=F(0)−F(1)+F(2)−F(3)没有冲突的每个元素分别为A,B,C时独立计算次数一条冲突的枚举冲突,讨论这两个端点的前...原创 2018-11-22 17:30:35 · 197 阅读 · 0 评论 -
【分治NTT/多项式求逆】JZOJ3303. 城市规划
Description求出n 个点的简单(无重边无自环)无向连通图数目.n <= 130000Solution设 f[i] 表示大小为i的答案,g[i]为2C(i,2)表示大小为i的任意无向图个数。考虑运用容斥。f[i] = g[i]-sigma( f[j] * g[i-j] * C(i-1,j-1) )用所有方案减去不连通的方案。枚举1所在的连通块(保证不重复不遗漏),再考...原创 2019-07-06 17:01:54 · 338 阅读 · 0 评论 -
JZOJ6272.【NOIP2019提高组A】整除
Description传送门Solutionn|xm-x等价于:任意p,满足 p|xm-x。p为质因子。将n的质因子p分开考虑。不难得出,在模p意义下满足以上条件的解y(y<p)。那么我们要求最终解x同时满足所有质因子条件,每一个质因子都可以选择一些解y中的一个并且要求x与y在模p意义下同余,那么这个条件肯定满足了。如果cnti表示对于第i个质因子y的个数,那么共有cnt...原创 2019-08-04 16:41:12 · 301 阅读 · 0 评论 -
JZOJ1898. 【2010集训队出题】密码系统
DescriptionLambda受任于某情报站,他的工作是获取敌人情报。一次他在破解密码系统时,得到了一个N位B进制数φ,满足φ≡V (mod M)。他发现组成φ的数字很奇特。为了验证φ的特殊性,他将所有模M为V的N位B进制数,按照各数位构成的集合分类,并想知道每一类数各有多少个。Solution不难想到矩乘,枚举所选集合再进行矩乘,f[i]表示余数为i的方案数,但是这样并不能保证...原创 2019-08-04 17:05:32 · 335 阅读 · 1 评论 -
JZOJ3402. 【GDOI2014模拟】Pty的字符串
咕了半个暑假的题目Description有一棵树有n个节点,深度浅的点往深度深的点连有向边,边上带有字符ch(a~c)给定字符串S,设一条从任意点出发往下的树上路径对应边上字符连接起来的字符串,求S的所有子串对应树上字符串的总匹配数和。|S|<=8e6,n<=8e5Solution关于子串问题,不难想到后缀自动机,刚开始我想以S建一个SAM,结果发现对于从任意点出发并...原创 2019-08-01 19:19:28 · 187 阅读 · 0 评论 -
JZOJ6343. 【NOIP2019模拟2019.9.7】Medium Counting
Description1≤n≤50,1≤L≤201\leq n \leq 50,1\leq L \leq 201≤n≤50,1≤L≤20Solution简单的计数,但是我的考虑方向并不正确,设的状态也不好做,导致正解其实并不难,但是我在考场上没有做出来。既然已经给定了顺序,那么一段满足条件一定是前几个字符完全一样,后面再区分开来。如果我们按照从后往前的顺序做,就很好转移,有点类似区间...原创 2019-09-07 20:27:02 · 347 阅读 · 0 评论 -
JZOJ6365. 【NOIP2019模拟2019.9.20】膜拜大会(fake)
Description定义一次对于i位置的操作为a[l]+=a[i],a[r]+=a[i],a[i]=0,其中l和r为i的左边和右边位置(1的左为n,n的右为1),每一个位置只可以操作一次给定一个长度n的序列a,求长度为m的操作序列的个数,使得最终a[1]>=K。模998244353T<=5,n<=1e5,m<=n-2,ai<=1e9Solution由...原创 2019-09-22 07:42:20 · 526 阅读 · 0 评论 -
JZOJ6383. 【NOIP2019模拟2019.10.07】果实摘取
Description小 D 的家门口有一片果树林,果树上果实成熟了,小 D 想要摘下它们。为了便于描述问题,我们假设小 D 的家在二维平面上的 (0, 0) 点,所有坐标范围的绝对值不超过 N 的整点坐标上都种着一棵果树。((0, 0) 这个点没有果树)小 D 先站在 (0, 0) 处,正对着 (1, 0) 的方向。每次摘果实时,小 D 会逆时针选择他能看到的第 K 棵还未摘取果实的果...原创 2019-10-11 12:30:26 · 333 阅读 · 0 评论 -
JZOJ4683. 【GDOI2017模拟8.11】矩阵
Description给定一个仅含有大写字母的n*m(n,m<=110)的矩阵,求不同的矩阵个数。Solution考虑枚举按照宽度将矩阵分类,分别统计。将每一行的每个位置开头,且长度该宽度的字串压一个编号,同一列的该宽度字串上下拼接即可形成任意一个该宽度下的矩阵。将每一列的编号排列成一个字符串,设为Si将所有列的S排列成:S1#S2#S3#…Sk我们只需要用后缀自动机求出...原创 2019-05-31 20:44:02 · 259 阅读 · 0 评论