数论
dance_in_the_dark
这个作者很懒,什么都没留下…
展开
-
[Gym102135][G - Digital characteristic]2017-2018 8th BSUIR Open Programming Contest
description给定a,b,m,k,产生长度为k的数组,a1=a%m,ai=(ai-1+b)%m(i>1),将这k个数首尾相接形成一个新数,问这个数的f函数的值‘’‘solution’’’首先f(x)可以等价于询问x%9的值,由于10的幂模9都为1,x%9=(∑(a+bi)%m)(\sum{(a+bi)\%m})%9(∑(a+bi)%m)=(∑a+bi−⌊a+bim⌋∗m)=(\sum{a+bi-\lfloor{\dfrac{a+bi}{m}}\rfloor*m})%9=(∑a+bi−⌊原创 2020-07-07 22:46:22 · 277 阅读 · 0 评论 -
第四次icpc集训(Saturday)(待修改)
总结第三次训练感觉我的状态比上一次又好了一些,但还不够熟练,主要是小错误不断犯,不能一次性打出bug较少的代码。A - oval-and-rectangle (签到)description给定一个椭圆的半长轴a,b,求在椭圆中[0.b]内任意选一点后画一条线后勾勒出的矩形的周长的期望solution积一下分code#include<bits/stdc++.h>#d...原创 2020-03-26 11:23:20 · 188 阅读 · 0 评论 -
【51nod1239】 欧拉函数之和
Description对正整数n,欧拉函数是小于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名,它又称为Euler’s totient function、φ函数、欧拉商数等。例如:φ(8) = 4(Phi(8) = 4),因为1,3,5,7均和8互质。 S(n) = Phi(1) + Phi(2) + …… Phi(n),给出n,求S(n),例如:n = 5,S(n) = 1 +原创 2017-03-08 21:54:11 · 1714 阅读 · 2 评论 -
【51nod 1244】 莫比乌斯函数之和
Description莫比乌斯函数,由德国数学家和天文学家莫比乌斯提出。梅滕斯(Mertens)首先使用μ(n)(miu(n))作为莫比乌斯函数的记号。具体定义如下: 如果一个数包含平方因子,那么miu(n) = 0。例如:miu(4), miu(12), miu(18) = 0。 如果一个数不包含平方因子,并且有k个不同的质因子,那么miu(n) = (-1)^k。例如:miu(2), miu原创 2017-03-09 15:33:56 · 667 阅读 · 0 评论 -
【51nod1363】最小公倍数之和
Description给出一个n,求1-n这n个数,同n的最小公倍数的和。 例如:n = 6,1,2,3,4,5,6 同6的最小公倍数分别为6,6,6,12,30,6,加在一起 = 66。 由于结果很大,输出Mod 1000000007的结果。Solution这道题化简到一半,差点废了,后来经高人指点才明白…… ans=∑i=1ni∗ngcd(i,n)=n∗∑d|n∑i=1n/di[原创 2017-03-07 16:46:16 · 1167 阅读 · 0 评论 -
【51nod 1253】Kundu and Tree
Description树包含N个点和N-1条边。树的边有2中颜色红色(‘r’)和黑色(’b’)。给出这N-1条边的颜色,求有多少节点的三元组(a,b,c)满足:节点a到节点b、节点b到节点c、节点c到节点a的路径上,每条路径都至少有一条边是红色的。 注意(a,b,c), (b,a,c)以及所有其他排列被认为是相同的三元组。输出结果对1000000007取余的结果。 Solution这道题看一下就原创 2017-03-29 12:04:45 · 563 阅读 · 0 评论 -
【51nod 1190】最小公倍数之和 V2
Description给出2个数a, b,求LCM(a,b) + LCM(a+1,b) + .. + LCM(b,b)。 例如:a = 1, b = 6,1,2,3,4,5,6 同6的最小公倍数分别为6,6,6,12,30,6,加在一起 = 66。 由于结果可能很大,输出Mod 10^9 + 7的结果。(测试数据为随机数据,没有构造特别坑人的Test)Solutionans=∑i=abi∗bgc原创 2017-03-14 17:17:59 · 638 阅读 · 0 评论 -
【51nod 1237】 最大公约数之和 V3
description给出一个数N(<= 10^10),输出小于等于N的所有数,两两之间的最大公约数之和。相当于计算这段程序(程序中的gcd(i,j)表示i与j的最大公约数): 由于结果很大,输出Mod 1000000007的结果。G=0; for(i=1;iSolution这道题想出来后被人告知是杜教筛…… 前面的都很顺利,非常快的推出式子。 ans=∑i=1n∑j=1ngcd(i,j)=原创 2017-03-14 16:14:27 · 1957 阅读 · 2 评论 -
【JZOJ5069】【GDSOI2017第二轮模拟】蛋糕
DescriptionCJY很喜欢吃蛋糕,于是YJC弄到了一块蛋糕,现在YJC决定和CJY分享蛋糕。 这块蛋糕上有n^2颗葡萄干,排成了一个n*n的点阵,每颗葡萄干互不相同且被编号为1~n^2。YJC决定沿着一条直线把蛋糕切成两份。YJC和CJY都很喜欢吃葡萄干,所以切出的两份蛋糕必须都包含至少一颗葡萄干。同时他们都不希望吃到不完整的葡萄干,所以切的时候不能经过任意一颗葡萄干。CJY喜欢1号葡萄干原创 2017-04-18 17:11:02 · 693 阅读 · 0 评论 -
【JZOJ5105】【GDSOI2017】魔兽争霸 x
DescriptionData ConstraintSolution我们发现答案最多只有两条直线构成,自己简单画一下就成。 现在问题成了如何求两条直线构成的最优方案。 我们设两条向量分别为(a,b),(c,d),价值分别为e,f,那么设(a,b)的时间为x,我们可以发现在n−a∗xc<m−b∗xd{n-a*x\over c}<{m-b*x\over d}时,f(x)=ex+f∗n−a∗xc=(e原创 2017-05-07 13:32:47 · 724 阅读 · 0 评论 -
【JZOJ5250】【GDOI2018模拟8.11】质数
DescriptionData ConstraintSolution我们发现2f(n)2^{f(n)}就等于∑i|n[gcd(i,n/i)==1]\sum_{i|n}[gcd(i,n/i)==1] 所以题目就是求:∑i=1n∑j|i[gcd(j,i/j)==1]\sum_{i=1}^{n}\sum_{j|i}[gcd(j,i/j)==1]=∑j=1n∑i=1⌊n/j⌋[gcd(i,j)==1]=\原创 2017-08-20 22:14:30 · 357 阅读 · 0 评论 -
【JZOJ5317】【清华集训2017模拟8.19】func
DescriptionData ConstraintSolution根据性质我们可以发现,对于一个f(x),f(x-1),若x为奇数,那么f(x)=f(x/2+1)+f(x/2),f(x-1)=f(x/2),若x为偶数,那么f(x)=f(x/2),f(x-1)=f(x/2)+f(x/2-1)。而且转移出来的数仍相邻,这也意味着里面有一个奇数有一个偶数。因为就对奇数满足f(x)=f(x-1)+f(x+原创 2017-08-24 08:00:05 · 426 阅读 · 0 评论 -
【JZOJ5224】【GDOI2018模拟7.12】C
DescriptionData ConstraintSolution首先必须讲讲自然数幂求和。 我们设Sk(n)=∑i=1nikSk(n)=\sum_{i=1}^ni^k 我们用第一类斯特林数来计算。 由第一类斯特林数的定义可知Cpn=Ppnp!=S(p,p)np−S(p,p−1)np−1+……S(P,0)n0p!C^p_n={P^p_n\over p!}={S(p,p)n^p-S(p,p-1原创 2017-07-12 22:43:44 · 437 阅读 · 0 评论 -
【JZOJ5413】【NOIP2017提高A组集训10.22】清兰
DescriptionData Constraint对于30%的数据,满足n<=500,m<=200 对于70%的数据,满足n<=20000,m<=100000 对于100%的数据,满足 1<=n<=50000,1<=m<=100000000,-100<=L<=100,|si|<=100 对于每个部分,均有50%的数据满足L=0Solution我们发现对于一个间距若分成k个间隔,一定是均分更原创 2017-10-25 21:18:05 · 421 阅读 · 0 评论 -
【51nod1238】 最小公倍数之和 V3
Description出个数N,输出小于等于N的所有数,两两之间的最小公倍数之和。相当于计算这段程序(程序中的lcm(i,j)表示i与j的最小公倍数): 由于结果很大,输出Mod 1000000007的结果。G=0; for(i=1;i< N;i++) for(j=1;j<=N;j++) { G = (G + lcm(i,j)) % 1000000007; }Solutionan原创 2017-03-19 10:48:40 · 925 阅读 · 0 评论 -
【JZOJ4965】【GDKOI2017模拟1.21】Equation
Description听着自己美妙的曲子,小Z进入了梦乡。在梦中,小Z仿佛又回到了自己纵横考场的年代。在梦中,小Z参加了一场考试,这场考试一共有n道题,每道题的最终得分都是一个大于等于0的整数。然而醒来后,小Z忘记了自己每道题的得分。他只记得自己计算过m次一些题目的分数和,每道题都被计算过,并且只被计算过一次。除此之外他还记得其中t道题的满分分别是多少(一道题的得分不会超过满分)。现在小Z想知道他这原创 2017-01-21 20:58:15 · 567 阅读 · 0 评论 -
【JZOJ3892】【NOIP2014模拟10.25A组】放棋子
DescriptionData ConstraintSolution这又是一道组合数的题。我们枚举现在至少有x行y列没有放至少一个棋子,那么要选择这x行y列的方案显然为CxN∗CyMC^{x}_{N}*C^{y}_{M},由于剩下的N-x行和M-y列需要放棋子,这些棋子有c+1种颜色(我们把不放棋子视为第c+1种颜色),且除了第c+1种颜色,其他颜色不能为空。这就让人想到了第二类Stirling,我原创 2017-01-21 07:47:06 · 421 阅读 · 0 评论 -
JZOJ4866【NOIP2016提高A组集训第8场11.5】禅与园林艺术
Description上了大学之后,小W和小Z一起报了一门水课,在做作业时遇到了问题。有一个长度为 n 的数列{ai},为一列树木的美观值。现在有m 次询问,每次给出三个数l,r和P,询问对于所有的l <= l’ <= r’ <= r (a[l’] + a[l’ + 1] + … + a[r’]) mod P的最小值。Data Constraint对于20%的数据 1 <= n, m <= 100原创 2016-11-05 14:33:54 · 510 阅读 · 0 评论 -
JZOJ4870. 【NOIP2016提高A组集训第9场11.7】涂色游戏
DescriptionData ConstraintSolution我们设f[i][j]表示当前有i个格子恰好放了j种颜色的方案数,那么f[i][j]=f[i−1][j−1]∗(p−(j−1))+f[i−1][j]∗jf[i][j]=f[i-1][j-1]*(p-(j-1))+f[i-1][j]*j。我们设g[j]表示n个格子恰好放了j种颜色的方案数,那么g[j]=f[n][j]。对于假如上一列放了原创 2016-11-07 21:09:53 · 661 阅读 · 3 评论 -
NOIP2016提高组复赛day2 组合数问题
DescriptionData ConstraintSolution由于k很小,最多也就分解成两个质因数,所以我们先将k分解成两个质因数x,y。然后统计出1~n的每个数是x或y的几次幂,设为a[i]。一次组合数运算相当于a[i]-a[j]-a[i-j]。判断一下即可。Code#include<iostream>#include<cmath>#include<cstring>#include<c原创 2016-11-28 22:04:36 · 1234 阅读 · 0 评论 -
JZOJ4888【NOIP2016提高A组集训第14场11.12】最近公共祖先
DescriptionYJC最近在学习树的有关知识。今天,他遇到了这么一个概念:最近公共祖先。对于有根树T的两个结点u、v,最近公共祖先LCA(T,u,v)表示一个结点x,满足x是u、v的祖先且x的深度尽可能大。YJC很聪明,他很快就学会了如何求最近公共祖先。他现在想寻找最近公共祖先有什么性质,于是他提出了这样的一个问题:n层的满k叉树T,求对于每一对(i,j)(1≤i,j≤T的点数),LCA(T,原创 2016-11-12 19:21:03 · 443 阅读 · 0 评论 -
JZOJ4877. 【NOIP2016提高A组集训第10场11.8】力场护盾
DescriptionZMiG成功粉碎了707的基因突变计划,为了人类的安全,他决定向707的科学实验室发起进攻!707并没有想到有人敢攻击她的实验室,一时间不知所措,决定牺牲电力来换取自己实验室的平安。 在实验室周围瞬间产生了一个无限大的力场护盾,它看上去无懈可击!不过ZMiG拥有惊人的双向观察能力,经过他的反复观察,找到了这个护盾的N个弱点,他本想逐一击破,却发现一股神秘力量阻止了他的原创 2016-11-08 22:36:43 · 661 阅读 · 0 评论 -
JZOJ 4921. 【NOIP2017提高组模拟12.10】幻魔皇
Description幻魔皇拉比艾尔很喜欢斐波那契树,他想找到神奇的节点对。所谓斐波那契树,根是一个白色节点,每个白色节点都有一个黑色节点儿子,而每个黑色节点则有一个白色和一个黑色节点儿子。神奇的节点对则是指白色节点对。请问对于深度为n的斐波那契树,其中距离为i的神奇节点对有多少个?拉比艾尔需要你对于1<=i<=2n的所有i都求出答案。Solution我们设f[i]表示第i层有多少白色节点,g1[i原创 2016-12-12 20:44:24 · 602 阅读 · 0 评论 -
JZOJ 4919. 【NOIP2017提高组模拟12.10】神炎皇
Description神炎皇乌利亚很喜欢数对,他想找到神奇的数对。对于一个整数对(a,b),若满足a+b<=n且a+b是ab的因子,则成为神奇的数对。请问这样的数对共有多少呢?Solution我们设一个数对(a,b)的最大公因数为d,a’*d=a,b’ *d=b,gcd(a’,b’)=1,那式子一定满足d(a′+b′)|d2a′b′d(a'+b')|d^2a'b'。因为gcd(a’,b’)=1,所以原创 2016-12-12 14:45:04 · 655 阅读 · 0 评论 -
【51 Nod】 1120 机器人走方格 V3
descriptionN * N的方格,从左上到右下画一条线。一个机器人从左上走到右下,只能向右或向下走。并要求只能在这条线的上面或下面走,不能穿越这条线,有多少种不同的走法?由于方法数量可能很大,只需要输出Mod 10007的结果。2 <= N <= 10^9Solution这道题找一下规律就发现了……(虽然我并没有发现这是卡特兰数)所以对于一个n,ans=Cn−12(n−1)∗2nmodpans原创 2016-12-03 20:10:43 · 510 阅读 · 0 评论 -
【51nod1403】 有趣的堆栈
Description大家都熟悉堆栈操作。一个堆栈一般有两种操作,push和pop。假设所有操作都是合法的并且最终堆栈为空。我们可以有很多方法记录堆栈的操作, (1) 对每个pop操作,我们记录它之前一共有多少个push操作。 (2) 对每个pop操作,我们记录这个被Pop的元素曾经被压上了几个。 例如:操作push, push, pop, push, push, pop, push, pop原创 2016-12-16 22:44:52 · 392 阅读 · 0 评论 -
【JZOJ 4924】 【NOIP2017提高组模拟12.17】向再见说再见
DescriptionData ConstraintSolution我们用dp来处理这道题。我们想将第一对的人从小到大排序,d[i]表示第二队有多少人的能力比第一队中第i人的能力低。设f[i][j]表示前i个人赢j场的方案数,但先不考虑第i个人被击败时被谁击败,只考虑赢的状况。那么f[i+1][j+1]+=f[i][j]*(d[i+1]-j),f[i+1][j]=f[i][j]。对于一个f[n][i原创 2016-12-17 21:07:10 · 551 阅读 · 0 评论 -
【JZOJ4932】【NOIP2017提高组模拟12.24】B
Description现在你有N个数,分别为A1,A2,…,AN,现在有M组询问需要你回答。每个询问将会给你一个L和R(L<=R),保证Max{Ai}-Min{Ai}<=R-L,你需要找出并输出最小的K(1<=K<=N,不存在输出-1)满足以下两个条件: ①能够在原来的N个数中选出不重复(下标不重复)的K个数,使得这K个数的和在区间[L,R]内。 ②能够在原来的N个数中选出不重复(下标不重复)的原创 2016-12-24 17:20:56 · 513 阅读 · 0 评论 -
【JZOJ3871】【NOIP2014八校联考第4场第1试10.19】无聊的游戏(game)
Description学校的运动会开始了,体能很菜的小可可没报任何比赛项目,于是和同学们玩一个十分无聊的游戏。 游戏在一个由n*n个方格组成的正方形棋盘上进行,首先在每个方格上均匀随机地填入1到m之间的正整数(每个方格填的数均不同),然后小可可均匀随机地选出k个1到m的数字(可能选的数不在棋盘上),把它们出现在棋盘上的方格涂黑,设有R行被整行涂黑,有C列被整列涂黑,小可可便可以得到2^(R+C)分原创 2017-01-18 07:42:19 · 687 阅读 · 0 评论 -
【JZOJ3874】【NOIP2014八校联考第4场第2试10.20】准备复赛(exam)
Description今年的NOIP初赛真是简单,小可可不用吹灰之力就考进了复赛,但是复赛可没有那么简单了,小可可想要好好准备复赛,争取复赛拿个省一。今天小可可在复习树和图的最大匹配时就碰到这样的一个难题:n个节点满足以下性质的不同的树有多少种。 1、树是有标号的,每个节点被标上1到n之间的整数; 2、每个节点最多和其他3个节点相连,但是1号节点最多和其他2个节点相连; 3、这棵树的最大匹配(原创 2017-01-18 20:39:08 · 783 阅读 · 0 评论 -
CF 702D Road to Post Office
Description这道题目是讲想在A想去一个城市,这个城市距离A所处的地方距离为d。A有一辆车,但这辆车每开k的距离就会熄火,熄火后A要用t的时间去维修。值得注意的是,这辆车一开始时维修好的。A的车每开1的距离需要的时间为a。另外,A也可以选择在任何一个地方下车,然后徒步走向终点,这徒步走1的距离为b(b>a)。问A到这个城市的最小时间为多少。Data Constraintd,k<=101210原创 2016-09-06 16:34:17 · 342 阅读 · 0 评论