数论
文章平均质量分 69
_beginend
这个作者很懒,什么都没留下…
展开
-
【51nod 1986 Jason曾不想做的数论题】【数论】
题意给定n,mn, mn,m,求∏X∈[1,m]nlcm(X1,⋯ ,Xn)gcd(X1,⋯ ,Xn)\prod_{X\in [1,m]^n}\mathrm{lcm}(X_1,\cdots,X_n)^{\gcd(X_1,\cdots,X_n)}X∈[1,m]n∏lcm(X1,⋯,Xn)gcd(X1,⋯,Xn)X∈[1,m]nX\in [1,m]^nX∈[1,m]n表示XXX取遍所有长度为nnn的序列,且XiX_iXi为111到mmm之间的整数。n≤109,m≤108n\le 10^9,m原创 2020-06-30 22:47:20 · 383 阅读 · 2 评论 -
【luogu P5655 基础数论函数练习题】【分治+数论】
题意给出nnn个数,有qqq次询问,每次询问一个区间内的lcm对1e9+71e9+71e9+7取模后的值。n,q,T≤300,1≤ai≤260n,q,T\le 300,1\le a_i\le 2^{60}n,q,T≤300,1≤ai≤260分析在求一个区间的lcm的时候,可以把lcm表示成∏bi\prod b_i∏bi的形式,其中bi∣aib_i|a_ibi∣ai。在新加入一个ai...原创 2019-11-16 00:17:29 · 383 阅读 · 2 评论 -
二次剩余Cipolla算法学习小记
前言之前遇到过一些关于二次剩余的题目,因为姿势不够都跳了。由于最近觉得自己的数论姿势严重不足,便觉得有必要来学习一下二次剩余相关的算法。 学完后感觉这个算法作者的脑洞真的是比较大,居然能想出这么巧妙的构造。 关于这个算法的介绍,我比较推荐czy大爷的博文。前置技能(以下内容均在模数pp为奇素数的前提下讨论)二次剩余首先要搞懂什么是二次剩余。若方程x2≡a(modp)原创 2018-01-22 08:01:09 · 3705 阅读 · 1 评论 -
Codeforces 235E Number Challenge 莫比乌斯反演+数论
题意设d(s)表示s的约数个数,给出a,b,c,求∑i=1a∑j=1b∑k=1cd(ijk)∑i=1a∑j=1b∑k=1cd(ijk)\sum_{i=1}^a\sum_{j=1}^b\sum_{k=1}^cd(ijk) a,b,c<=2000分析题解貌似是一个很鬼畜的做法。。。 约数个数函数σ0(d)σ0(d)\sigma_0(d)有一个小结论,就是σ0(ij)=∑p|i...原创 2018-02-28 20:08:40 · 419 阅读 · 0 评论 -
Codeforces 906D Power Tower 拓展欧拉定理
题意定义一个序列w[1..n]的价值为 给出一个长度为n的序列w,每次询问一个区间[l,r]的价值模m。 n≤200000,1≤m≤109n≤200000,1≤m≤109n\le200000,1\le m\le10^9分析首先注意到,若序列中某个数为1,那么从这一位开始后面的元素全部是没用的。那么我们可以先把这部分去掉,这样序列中的全部元素都大于1。 定义函数solve(l,r,...原创 2018-04-07 20:59:46 · 304 阅读 · 0 评论 -
bzoj 4635: 数论小测验 数论+容斥原理
题意有一个长度为N的数组Ai,每个元素可以取1~M中的一个正整数。那么一共有M^N种可能的数组。因为 SHUXK 对数 论有特殊的爱好,所以他立刻想到了下面两个问题: 1. 对于给定的正整数K,有多少个数组Ai满足GCD(A1,A2…An) = K 2. 对于给定的正整数K,有多少个数组Ai满足K|Lcm(A1,A2…An) (我相信机智的你在看到这道题的英文名称时就已经猜得八九不离十...原创 2018-03-25 19:44:13 · 344 阅读 · 0 评论 -
bzoj 4833: [Lydsy1704月赛]最小公倍佩尔数 数论
题意令(1+2–√)n=e(n)+f(n)∗2–√(1+2)n=e(n)+f(n)∗2(1+\sqrt 2)^n=e(n)+f(n)*\sqrt 2,其中e(n),f(n)e(n),f(n)e(n),f(n)都是整数,显然有(1−2–√)n=e(n)−f(n)∗2–√(1−2)n=e(n)−f(n)∗2(1-\sqrt2)^n=e(n)-f(n)*\sqrt2。令g(n)表示f(1),f(2)...原创 2018-04-03 17:14:56 · 247 阅读 · 0 评论 -
hdu 5780 gcd 线性筛+数论
题意给出n个x,求∑1≤a,b≤ngcd(xa−1,xb−1)∑1≤a,b≤ngcd(xa−1,xb−1)\sum_{1\le a,b\le n}gcd(x^a-1,x^b-1) 1≤n,x≤1000001≤n,x≤1000001\le n,x\le100000,300组数据。分析有个结论就是gcd(xa−1,xb−1)=xgcd(a,b)−1gcd(xa−1,xb−1)=xgc...原创 2018-04-03 17:38:27 · 301 阅读 · 0 评论 -
bzoj 5028: 小Z的加油店 线段树+差分+数论
题意小Z经营一家加油店。小Z加油的方式非常奇怪。他有一排瓶子,每个瓶子有一个容量vi。每次别人来加油,他会让 别人选连续一段的瓶子。他可以用这些瓶子装汽油,但他只有三种操作: 1.把一个瓶子完全加满; 2.把一个瓶子完全倒空; 3.把一个瓶子里的汽油倒进另一个瓶子,直到倒出瓶子空了或者倒进的瓶子满了。 当然,为了回馈用户,小Z会时不时选择连续一段瓶子,给每个瓶子容积都增加x。 为了...原创 2018-05-11 23:32:41 · 614 阅读 · 0 评论 -
LibreOJ #2523.「HAOI2018」奇怪的背包 动态规划+数论
题意小CCC非常擅长背包问题,他有一个奇怪的背包,这个背包有一个参数PPP,当他向这个背包内放入若干个物品后,背包的重量是物品总体积对PPP取模后的结果.现在小CCC有nnn种体积不同的物品,第iii种占用体积为ViViV_i,每种物品都有无限个.他会进行qqq次询问,每次询问给出重量wiwiw_i,你需要回答有多少种放入物品的方案,能将一个初始为空的背包的重量变为wiwiw_i.注意,两种方...原创 2018-05-16 16:13:52 · 367 阅读 · 0 评论 -
LibreOJ 6247&bzoj 3328: PYXFIB 数论
题意给出n,k,pn,k,pn,k,p,设FiFiF_i表示斐波那契数列第iii项,求∑i=0nCinFi[k|i]∑i=0nCniFi[k|i]\sum_{i=0}^nC_n^iF_i[k|i] 满足n≤1018,k≤104,p≤109,p是素数,pmodk=1n≤1018,k≤104,p≤109,p是素数,pmodk=1n\le 10^{18},k\le10^4,p\le10^9,p是素...原创 2018-06-11 18:31:31 · 184 阅读 · 0 评论 -
Codeforces 819D Mister B and Astronomers 数论+exgcd
题意有nnn个人在观察天空,第imodn+1imodn+1i\bmod n+1个人会在第iii个人观察完后的第xi+1xi+1x_{i+1}个时刻观察一个时刻的天空,其中第一个人在第000个时刻观察天空。现在天空中有一颗星星,只知道其每隔TTT秒会闪烁一次。对于每一个人,询问在[0,T)[0,T)[0,T)中有多少个整数时刻满足若星星在这个时刻闪烁,则这个人是第一个观察到星星的人。 n≤20...原创 2018-06-13 11:43:14 · 262 阅读 · 0 评论 -
Codeforces 1027G X-mouse in the Campus 数论+Pollard_rho
题意给定mmm和xxx,满足gcd(m,x)=1gcd(m,x)=1gcd(m,x)=1。现在把每个小于mmm的整数都看作一个点,然后iii向ixixix连边,问最后最少需要选出多少个点使得每个点的后继中至少有一个点被选。 m≤1014m≤1014m\le10^{14}分析感谢sam队长教我做这题。 首先因为gcd(m,x)=1gcd(m,x)=1gcd(m,x)=1,所以最后...原创 2018-08-19 14:10:04 · 737 阅读 · 2 评论 -
牛客网Wannafly挑战赛23 F-计数 矩阵树定理+拉格朗日插值法
题意给出一个n个点m条边的带权无向图,问有多少棵生成树满足边权和模k等于0,答案模p输出。 n,k≤100,p≤109,p为质数且满足k≡1(modp)n,k≤100,p≤109,p为质数且满足k≡1(modp)n,k\le100,p\le10^9,p为质数且满足k\equiv1\pmod p。分析一种直观的套路就是把每条边的边权看成是xwxwx^w,然后做矩阵树定理,最后得出来...原创 2018-09-09 16:13:42 · 408 阅读 · 0 评论 -
LibreOJ #525.「LibreOJ β Round #4」多项式 拓展欧拉定理
题意构造一个整多项式f(x)f(x)f(x)使得对任意的整数xxx都有f(x)≡0(modk)f(x)\equiv 0\pmod kf(x)≡0(modk)。k≤30000k\le30000k≤30000,多项式的次数不超过600006000060000分析根据拓展欧拉定理,对任意的x,px,px,p,若n≥φ(p)n\ge\varphi(p)n≥φ(p)则有xn≡xn+φ(p)(modp...原创 2019-07-08 08:42:12 · 196 阅读 · 0 评论 -
51nod 2026 Gcd and Lcm 杜教筛+数论
题意已知f(x)=∑d|xμ(d)∗df(x)=\sum_{d|x}\mu(d)*d 现在请求出下面式子的值 ∑i=1n∑j=1nf(gcd(i,j))∗f(lcm(i,j))\sum\limits_{i=1}^n\sum\limits_{j=1}^nf(gcd(i,j))*f(lcm(i,j)) 答案模1e9+7。 n<=109n<=10^9分析设x=pa11∗...∗pakkx=p_1^原创 2017-12-28 15:06:38 · 375 阅读 · 0 评论 -
bzoj 4522: [Cqoi2016]密钥破解 pollard_rho+欧拉定理
题意一种非对称加密算法的密钥生成过程如下: 1. 任选两个不同的质数 p ,q 2. 计算 N=pq , r=(p-1)(q-1) 3. 选取小于r ,且与 r 互质的整数 e 4. 计算整数 d ,使得 ed≡1 mod r 5. 二元组 (N,e) 称为公钥,二元组 (N,d) 称为私钥当需要加密消息 n 时(假设 n 是一个小于 N 整数,因为任何格式的消息都可转为整数表示),使用原创 2017-12-27 11:25:48 · 304 阅读 · 0 评论 -
CodePlus 2017 12 月赛 可做题2 矩阵乘法+exgcd
题意 分析注意到a1a_1对aka_k贡献的系数为斐波那契数列第k-2项,a2a_2对aka_k贡献的系数为第k-1项,那么可以把这两项求出来,那么a2a_2就一定是某个不定方程ax+by=c的解x。 设d为gcd(a,b),先解出该方程一个解x0,y0,有个定理就是该方程的全部解是x0+(b/d)n,y0-(a/d)n,那么就很容易求出区间解的数量了。代码#include<iostream>原创 2017-12-25 21:41:25 · 503 阅读 · 0 评论 -
51nod 1355 斐波那契的最小公倍数 莫比乌斯反演+数学
题意F(0) = 0 F(1) = 1 F(n) = F(n-1) + F(n-2) 给出n个正整数a1, a2,…… an,求对应的斐波那契数的最小公倍数,由于数字很大,输出Mod 1000000007的结果即可。 2 <= N <= 50000,1 <= ai <= 1000000。分析跪zyz大爷写的题解。首先有个结论就是 gcd(fn,fm)=fgcd(n,m)gcd(f_n,f_m原创 2017-12-24 15:44:17 · 502 阅读 · 0 评论 -
51nod 1318 最大公约数与最小公倍数方程组 2-SAT+数学
题意求解一个有趣的方程组,方程组有N个未知正整数,x[0],x[1],x[2],…,x[N-1]。 方程组由M个方程组成,方程只有两种类型: 1)GCD(x[i],x[j]) = G ;其中,i != j,且GCD(a,b)为正整数a与b的最大公约数的值 2)LCM(x[i],x[j]) = L ;其中,i != j,且LCM(c,d)为正整数c与d的最小公倍数的值 你需要判断这样一个方程原创 2017-12-23 13:13:41 · 605 阅读 · 0 评论 -
bzoj 1406: [AHOI2007]密码箱 数学
题意在一次偶然的情况下,小可可得到了一个密码箱,听说里面藏着一份古代流传下来的藏宝图,只要能破解密码就能打开箱子,而箱子背面刻着的古代图标,就是对密码的提示。经过艰苦的破译,小可可发现,这些图标表示一个数以及这个数与密码的关系。假设这个数是n,密码为x,那么可以得到如下表述: 密码x大于等于0,且小于n,而x的平方除以n,得到的余数为1。 小可可知道满足上述条件的x可能不止一个,所以一定要把所有满足原创 2017-10-22 21:41:41 · 311 阅读 · 0 评论 -
如何线性求[1,p-1]的逆元
今天做了一下bzoj 4011,顺便学习了一波如何线性递推求[1,p-1]在模p下的逆元。注意p为素数。 设aa的逆元为a−1a^{-1},根据逆元的定义显然有a∗a−1≡1(modp)a*a^{-1}\equiv1\pmod p首先有1的逆元就是1. 假设现在已经球出了[1,a−1][1,a-1]的逆元,现在考虑如何求aa的逆元。设k=⌊pa⌋,r=pmodak=\lfloor\frac{p}原创 2017-06-13 20:11:54 · 742 阅读 · 0 评论 -
bzoj 3884: 上帝与集合的正确用法 欧拉定理+数学
题意给出p,求2222....2^{2^{2^2....}}无限个2 mod p的值。 p<=107p<=10^7分析PoPoQQQ好强啊!!! 复制一波题解: 代码#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#include<cmath>#define原创 2017-02-16 19:56:13 · 281 阅读 · 0 评论 -
bzoj 2142: 礼物(组合数取模终极版) 组合数学+中国剩余定理+exgcd
题意一年一度的圣诞节快要来到了。每年的圣诞节小E都会收到许多礼物,当然他也会送出许多礼物。不同的人物在小E心目中的重要性不同,在小E心中分量越重的人,收到的礼物会越多。小E从商店中购买了n件礼物,打算送给m个人,其中送给第i个人礼物数量为wi。请你帮忙计算出送礼物的方案数(两个方案被认为是不同的,当且仅当存在某个人在这两种方案中收到的礼物不同)。由于方案数可能会很大,你只需要输出模P后的结果。 m原创 2017-04-10 21:51:04 · 725 阅读 · 0 评论