【数论】--概率与期望第24页

ldu-自测二-A，B（位运算符+位运算线段树（32树））

自测二的题都是数论，推出来题就很简单，推不出就发呆吧。淦！众所周知，int范围内，搞位运算符一般开到32就够了。

蒲公英之殇·2020-10-11 20:03

C++等差数列（数论、欧几里得辗转相除gcd）

数学老师给小明出了一道等差数列求和的题目。但是粗心的小明忘记了一部分的数列，只记得其中N个整数。现在给出这N个整数，小明想知道包含这N个整数的最短的等差数列有几项？输入格式输入的第一行包含一个整数N。第二行包含N个整数A1,A2,⋅⋅⋅,AN。(注意A1∼AN并不一定是按等差数列中的顺序给出)输出格式输出一个整数表示答案。数据范围2≤N≤100000,0≤Ai≤109输入样例：52641020输出

赴星辰大海·2020-10-10 00:34

一些基础数论的知识和证明

算术基本定理N=pα1∗pα2∗...∗pαkN=p^{\alpha_{1}}*p^{\alpha_{2}}*...*p^{\alpha_{k}}N=pα1∗pα2∗...∗pαk约数个数(α1+1)∗(α2+1)...∗(αk+1)(\alpha_{1}+1)*(\alpha_{2}+1)...*(\alpha_{k}+1)(α1+1)∗(α2+1)...∗(αk+1)证明：已知N=pα1∗pα

Dejavu1z·2020-10-07 17:26

概率与期望习题总结

总结概率题一般正着推期望题一般倒着推图上的问题如果是$DAG$可以直接转移否则可能要用到高斯消元$20$的数据范围大概率是装压有些看似无限循环的式子其实可以倒着递推1、骰子基础版题目描述众所周知，骰子是一个六面分别刻有一到六点的立方体，每次投掷骰子，从理论上讲得到一点到六点的概率都是$1/6$。今有骰子一颗，连续投掷$N$次，问点数总和大于等于$X$的概率是多少？输入仅有一行包

·2020-10-02 17:00

【数论】GCD

F.GCD-GCD-GCD[Problem4045]Description给你两个数a,m;问你在[0,m−1]范围内有几个数x满足gcd(a,m)=gcd(a+x,m);Input第一行一个整数T；接下来T行，每行两个整数a,m;1≤ay,则有gcd(x,y)=gcd(x−y,y)x>y,则有gcd(x,y)=gcd(x-y,y)x>y,则有gcd(x,y)=gcd(x−y,y)。gcd(a+x

biu~跃哥冲冲冲·2020-09-28 16:00

数论 —— 整式方程

【概念】整式方程就是方程中所有的未知数均在分子上，分母只是常数且无未知数。通常情况下，常年用字母x、y、z来表示未知数，方程中含有几个不同的未知数就叫做几元，未知数的最高次数是几就叫做几次。例如：ax+b=c就是一个一元一次整式方程【一元一次方程求解】对于方程ax+b=c，有：x=(c-b)/adoublecalculate(doublea,doubleb,doublec){return(c-b)

Alex_McAvoy·2020-09-17 15:37

【基础题】数论系列1-HDU1061、1097、1098

HDU1061RightmostDigit题目HDU1061RightmostDigit题目分析题目大意：输入一个正整数n，计算nn的个位数字。解题思路：刚开始是想n次方的个位数字只与n的个位数有关，所以就计算了n的个位数的n次方，提交了之后果然TLE了。后来对每个数字的幂进行了计算，发现1~9每个数字的幂都是有周期的，1、5、6的任何次幂都是它本身，4、9幂的周期是2，2、3、8幂的周期是4，其

vxiao_shen_longv·2020-09-17 13:05

【从0开始ACM】【LV3】【数学】【 FFT 快速傅里叶变换 && NTT 更高精度的快速数论变换】

DiscreteFouriertransformDFT离散傅里叶变换——一种过程InverseDiscreteFourierTransformIDFT离散傅里叶逆变换FastFourierTransformationFFT快速傅里叶变换——一种算法FFT在ACM中解决啥呢多项式乘法A*B=CA=a0+a1x^1+a2x^2+a（n-1）x^（n-1）B=b0+...C=c0+...+c(n-1)x

under_sky_dxj·2020-09-17 13:58

HDU_4248_A Famous Stone Collector(组合数学+DP)

题型：数论题意：有N堆石子，每堆颜色相同，任意两堆颜色不同。问所有的石子可以排成多少种不同的序列。分析：采用插空法的思想：设dp(i,j)表示用前i堆石头排成了长度为j的序列。

light_starlight·2020-09-17 13:56

数论中四方定理

/*标题：四方定理数论中有著名的四方定理：所有自然数至多只要用四个数的平方和就可以表示。我们可以通过计算机验证其在有限范围的正确性。对于大数，简单的循环嵌套是不适宜的。下面的代码给出了一种分解方案。

mtler·2020-09-17 13:18

[2017HNCPC] Strange Optimization 简单数论

给出正整数n,m≤1e9n,m\leq1e9n,m≤1e9，并且定义函数f(t)=min⁡i,j∈Z∣in−jm+t∣f(t)=\min_{i,j\inZ}|\frac{i}{n}-\frac{j}{m}+t|f(t)=mini,j∈Z∣ni−mj+t∣，求一实数α\alphaα使得f(12+α)f(\frac{1}{2}+\alpha)f(21+α)最大，求出值。由于α\alphaα是任意实数，

SigmaQuadrant·2020-09-17 12:27

数论之路慢慢之GCD性质

题目链接StrangeOptimizationBoboisfacingastrangeoptimizationproblem.Givenn,m,heisgoingtofindarealnumberαsuchthatf(12+α)ismaximized,wheref(t)=mini,j∈Z|in−jm+t|.Helphim!Note:Itcanbeprovedthattheresultisalway

Phoenix丶HN·2020-09-17 12:04

UVa 12683 Odd and Even Zeroes（数论+数位DP）

题意：问小于等于n的数中（0#include#include#include#include#include#include#include#includeusingnamespacestd;typedeflonglongLL;vectordigit;#defineREP(_,a,b)for(int_=(a);_<=(b);_++)LLdp[30][2][2];LLn;LLdfs(intpos,i

mowayao·2020-09-17 12:14

找有最多个约数的数（数论）

题目：一个数的素数因子分解是数论里面的基本问题，我们在课堂上专门讨论过这样的问题，大家也做过实验。

GOTOTHEBAD·2020-09-17 11:43

数论模板

扩展欧几里得：LLextgcd(LLa,LLb,LL&x,LL&y){if(b==0){x=1,y=0;returna;}LLg=extgcd(b,a%b,y,x);y-=(a/b)*x;returng;}快速幂：LLpowmod(LLa,LLb,LLc){LLans=1;while(b){if(b&1)ans*=a,ans%=c;b>>=1;a*=a;a%=c;}returnans;}中国剩余定

GOTOTHEBAD·2020-09-17 11:12

HDU 4254 A Famous Game（概率与期望）

题目描述传送门题目大意：一个口袋里有n个红色或蓝色的球。n+1种颜色分布情况（i个红球n−i个蓝球）的概率是相等的。B从口袋中不放回地摸出了p个球，其中有q个是红色的。求B再摸一个球时，摸出的球是红色的概率。题解设Nk表示n个球中有k个红球的概率。A表示p个球中有q个红球B表示下次摸出的是红球那么P(Nk)=1n+1P(A)=C(k,q)C(n−k,p−q)C(n,p)P(B|ANk)=k−qn−

clover_hxy·2020-09-17 11:25

HDU 5753 Permutation Bo （概率与期望）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5753#includeusingnamespacestd;#definedebugputs("YES");#definerep(x,y,z)for(int(x)=(y);(x)#definemk(x,y)make_pair(x,y)#definefifirst#definesesecondconstin

等我学会后缀自动机·2020-09-17 10:19

《算法导论》笔记汇总

｜发表于2011/08/2118条评论14,136views列表里没有的，或者是纯属理论，不适合写，比如第1~5章；或者是我也不怎么明白的，比如斐波那契堆中抽取最小结点的平摊代价分析；还有没看的，比如数论和

KinneV·2020-09-17 06:19

常见算法：C语言求最小公倍数和最大公约数三种算法

最小公倍数：数论中的一种概念，两个整数公有的倍数成为他们的公倍数，其中一个最小的公倍数是他们的最小公倍数，同样地，若干个整数公有的倍数中最小的正整数称为它们的最小公倍数，维基百科：定义点击打开链接求最小公倍数算法

乐呵1呵·2020-09-17 04:00

[数论]Coprime

题目描述对于两个整数k和m，如果k和m的最大公约数为1，则k和m互质。给出两个正整数n和m(m≤n),定义f(n,m)为1~n！中与m！互质的数的个数。其中n！=1*2*3*..*(n-1)*n。Task：给定n和m，要求计算f(n,m)。Input本题设多组数据。输入文件的第一行有一个整数T(1≤T≤100000)，表示有T组数据。接下来有T行，每行两个整数n和m(2≤n≤100000,2≤m≤

在征途·2020-09-17 03:20

NYOJ 1066 CO-PRIME（数论）

CO-PRIME时间限制：1000ms|内存限制：65535KB难度：3描述Thisproblemissoeasy!Canyousolveit?Youaregivenasequencewhichcontainsnintegersa1,a2……an,yourtaskistofindhowmanypair(ai,aj)(i#include#includeusingnamespacestd;consti

lyhvoyage·2020-09-17 03:15

C语言求最小公倍数和最大公约数三种算法(经典)

最小公倍数：数论中的一种概念，两个整数公有的倍数成为他们的公倍数，其中一个最小的公倍数是他们的最小公倍数，同样地，若干个整数公有的倍数中最小的正整数称为它们的最小公倍数，维基百科：定义点击打开链接求最小公倍数算法

Diiiiiiiiiiiiiiii·2020-09-17 03:14

【莫比乌斯反演最简单的入门题】B.Coprime Integers

提交或者vjudge思路真是个傻逼题，但是当时没有学完莫比乌斯init是在用线性筛获得莫比乌斯函数的前缀和然后用容斥原理得到解注意溢出前置知识莫比乌斯反演线性筛积性函数说明莫比乌斯真是毒瘤，看了几天的《初等数论及其应用

STL_CC·2020-09-17 03:32

数论 - GCD LCM - UVA 11388

数论-GCDLCM-UVA11388题意：T组测试数据，T组测试数据，T组测试数据，每组包括两个正整数，分别表示一个最大公约数G和最小公倍数L。

njuptACMcxk·2020-09-17 02:52

公钥密码（密码学数学基础、RSA、ElGamal、Rabin、椭圆曲线密码体制）

私钥本人私用公钥密码体制也称非对称密码体制公钥密码体制主要用于密钥分发公钥密码体制优势密钥分发：公钥采用公开信道传输密钥管理：在N个用户的系统中，每个用户只需要保管自己的私钥以及其他N-1个用户的公钥，整个系统只需要维护N个公钥密码学数学基础之数论同余类

许永亮·2020-09-16 19:03

数论。gcd求和

https://blog.csdn.net/luojiayu14/article/details/7305030找到了这篇题解，太难了对我来说大概做了四个小时，真是太弱了。下面是最近遇到的问题swap(a,b)a和b的值互换整形，double，float，char都可以互换，但是一个变量a&1是合法的代表的意思是a和1做二进制的且运算即看a的最后边那一位是不是1是1的话返回1否则返回0追问那是不是

天天寒暑假·2020-09-16 19:52

gcd(数论)

题目描述给定整数N，求1#include#include#includeusingnamespacestd;intn;#definelllonglongconstintMAXN=1e7+3;intprime[MAXN],pn;llphi[MAXN];boolvis[MAXN];voidpr(){ for(inti=2;i<=n;i++){ if(!vis[i]){ prime[+

DD(XYX)·2020-09-16 19:14

bzoj 5394: [Ynoi2016]炸脖龙数论+树状数组

给一个长为n的序列，m次操作，每次操作：Input第一行两个整数n,m表示序列长度和操作数接下来一行，n个整数，表示这个序列接下来m行，可能是以下两种操作之一：1lrx表示区间[l,r]加上x2lrp表示对区间[l,r]进行一次查询，模数为pn,musingnamespacestd;#definemaxn500020#defineN20000000#definerep(i,l,r)for(regi

Thomas_ZQQ@Runespoor·2020-09-16 19:01

UVA 10006 - Carmichael Numbers 数论（快速幂取模 + 筛法求素数）

CarmichaelNumbersAnimportanttopicnowadaysincomputerscienceiscryptography.Somepeopleeventhinkthatcryptographyistheonlyimportantfieldincomputerscience,andthatlifewouldnotmatteratallwithoutcryptography.A

weixin_34290000·2020-09-16 19:25

[学习笔记]NTT——快速数论变换

先要学会FFT[学习笔记]FFT——快速傅里叶变换一、简介FFT会爆精度。而且浮点数相乘常数比取模还大。然后NTT横空出世了虽然单位根是个好东西。但是，我们还有更好的东西我们先选择一个模数，$const\spaceint\spacep=998244353$设g为p的单位根。这里就是3那么有：$(\omega_n^1)^n=g^{p-1}=1\spacemod\spacep$那么，假设$x=(\om

weixin_34126557·2020-09-16 19:46

[数论]莫比乌斯反演2

索引莫比乌斯反演1定理莫比乌斯反演2证明莫比乌斯反演3技巧前言本篇内容为定理的证明定理请参考：>传送门<三个性质的证明性质1证明:这个式子是莫比乌斯函数真正的定义式但是我们还是有证明当$n=1$时，显然\[\sum_{d|n}\mu(d)=\mu(1)=1\]根据定义直接得到的结论当$n\neq1$时，\[\sum_{d|n}=\mu(a_1)+\mu(a_2)+\dots+\mu(a_m

weixin_30815427·2020-09-16 19:51

【BZOJ】【P4407】【于神之怒加强版】【题解】【数论】

id=4407这两天刚好在给学弟讲数论，bzoj上就来了一道裸题……答案就是\sum_DF(D)*n/d*m/dF(D=\sum{d|D}d^Kmu(D/d)线性筛，分块求F(p)=p^K-1F(x*p

iamzky·2020-09-16 19:16

数论 —— 莫比乌斯反演

【反演】假设我们手头有个数列F，通过某种变换H，可以得到函数G。，即：但现在只有函数G，需要求F，那么我们就需要寻找一种变换，使得G在经过这种变换后能够获得F，这个过程即为反演，即：【整除分块】对于式子：，其时间复杂度为O(n)，当有多组数据时，O(n)并非正确的时间复杂度，此时有一种时间复杂度为O(√n)的算法：整除分块对于每一个通过打表发现，很多的值是相同的，它们呈一个块状分布，对每一个值相同

Alex_McAvoy·2020-09-16 19:41

bzoj4176 Lucas 的数论

======∑i=1n∑j=1nf(ij)∑i=1n∑j=1n∑d=1n2[dgcd(i,d)|j]∑i=1n∑d=1n2⌊n∗gcd(i,d)d⌋∑d=1n∑i=1⌊nd⌋∑j=1⌊n2d⌋⌊nj⌋e(gcd(i,j))∑d=1n∑i=1⌊nd⌋∑j=1n⌊nj⌋∑x|gcd(i,j)μ(x)∑x=1nμ(x)∑d=1n⌊ndx⌋∑j=1⌊nx⌋⌊njx⌋∑x=1nμ(x)(∑j=1⌊nx⌋⌊n

sdfzyhx·2020-09-16 19:42

数论

（数论那玩意儿是谁发明的给我站出来）感谢黄老师为我打开数论（绝望）的大门进位计数制b进制向十进制转换乘以基数并展开：（按位权）十进制向b进制转换整数部分除以基数并倒取余数，小数部分乘以基数，并顺取整数部分

yingyudou·2020-09-16 19:22

莫比乌斯反演与积性函数部分(转载)

article/details/77888386https://blog.csdn.net/skywalkert/article/details/50500009（偏理论）莫比乌斯与积性函数之前做过不少的数论题

等我学会后缀自动机·2020-09-16 18:03

【BZOJ4407】于神之怒加强版（莫比乌斯反演）

(i,j)k题解根据惯用套路把公约数提出来∑d=1ndk∑i=1n∑j=1m[gcd(i,j)==d]再提一次∑d=1ndk∑i=1n/d∑j=1m/d[gcd(i,j)==1]后面这个东西很显然可以数论分块

小蒟蒻yyb·2020-09-16 18:07

BZOJ3601：一个人的数论（莫比乌斯反演+伯努利数）

题面题意：给出d和n（n以分解质因数给出），问所有与n互质的数的d次幂之和，即∑xxd[gcd(x,n)==1]∑xxd[gcd(x,n)==1]套路推♂倒=∑i|nμ(i)∗id∗∑x=1nixd=∑i|nμ(i)∗id∗∑x=1nixdfa♂现最右边是个幂和，设为h(ni)h(ni)，整个就是一个狄利克雷卷积。根据题目n以分解质因数输入的套路，应该是找到积性函数，然后一个个质因数乘起来μ(i)

KKiseki·2020-09-16 18:07

莫比乌斯反演证明

首先定义几个概念：1，卷积：设是两个数论函数（也就是说，以自然数集为定义域的复数值函数），则卷积运算定义为可以证明，卷积运算满足：1）交换律：由定义显然。

nano_jz·2020-09-16 18:32

浅谈一类积性函数的前缀和

前置技能积性函数的定义若f(n)的定义域为正整数域，值域为复数，即f:Z+→C，则称f(n)为数论函数。

nano_jz·2020-09-16 18:27

bzoj 4174 tty的求助数论莫比乌斯反演

题意：求∑Nn=1∑Mm=1∑m−1k=0⌊nk+xm⌋先膜PoPoQQQ大爷。∑Nn=1∑Mm=1∑m−1k=0⌊nk+xm⌋=∑Nn=1∑Mm=1∑m−1k=0⌊nk%m+xm⌋+nkm−nk%mm然后先不考虑n,m，分三部分算。设gcd(n,m)=d。∑m−1k=0⌊nk%m+xm⌋=∑m−1k=0⌊d(ndk%md)+xm⌋然后因为nd和md互质，所以：=d∑md−1k=0⌊dk+xm⌋=

make_it_for_good·2020-09-16 18:42

bzoj 4176 Lucas的数论莫比乌斯反演

∑ni=1∑nj=1d(i,j)=∑ni=1∑nj=1∑n2k=1[k|ij]=∑ni=1∑nj=1∑n2k=1[kgcd(i,k)|j]=∑ni=1∑n2k=1⌊n∗gcd(i,k)k⌋=∑nd=1∑⌊nd⌋i=1∑⌊n2d⌋k=1⌊nk⌋[gcd(i,k)=1]=∑nd=1∑⌊nd⌋i=1∑nk=1⌊nk⌋[gcd(i,k)=1]=∑nd=1∑⌊nd⌋i=1∑nk=1⌊nk⌋∑t|i,t|kμ

make_it_for_good·2020-09-16 18:42

算法学习FFT系列(2):快速数论变换NTT &&bzoj3992: [SDOI2015]序列统计例题详解

bzoj3992:[SDOI2015]序列统计Description小C有一个集合S，里面的元素都是小于M的非负整数。他用程序编写了一个数列生成器，可以生成一个长度为N的数列，数列中的每个数都属于集合S。小C用这个生成器生成了许多这样的数列。但是小C有一个问题需要你的帮助：给定整数x，求所有可以生成出的，且满足数列中所有数的乘积modM的值等于x的不同的数列的有多少个。小C认为，两个数列{Ai}和

lvzelong2014·2020-09-16 18:02

积性函数与Dirichlet卷积学习小记

先给一波定义：（这里也感谢skywalkert大佬的精心讲解）数论函数的定义若f(n)f(n)的定义域为正整数域，值域为复数，即f:Z+→Cf:Z+→C，则称f(n)f(n)为数论函数。

Felix-Lee·2020-09-16 18:50

[数论数学]莫比乌斯反演定理

近期学习了莫比乌斯反演，算是一个学习笔记吧…在这里首先要说明：1：本文讨论的所有函数为数论函数，即定义域为D=N∗D=N^*D=N∗的函数；2：∑d∣nf(d)\sum\limits_{d|n}f(d)

ix35·2020-09-16 18:43

莫比乌斯反演+数论分块

对于给定的整数a,b和d，有多少正整数对x,y，满足x#definelllonglongusingnamespacestd;constintmaxn=1e5+7;boolvis[maxn];llprime[maxn],mu[maxn];llcnt;voidInit(){llN=maxn;mu[1]=1;cnt=0;for(lli=2;i>T;while(T--){cin>>b>>d>>k;b/=k

henulmh·2020-09-16 18:51

数论二项式反演

反演公式c和d是两个跟n和r有关的函数根据用法不同，c和d是不同的一般数学家会先随便弄c函数然后经过复杂的计算和证明，得到d函数然后公式就可以套用了二项式反演公式(划重点)UVALive7040传送门题意：给n盆花涂色，从m种颜色中选取k种颜色涂，保证正好用上k种颜色，你必须用上这k种颜色去涂满n个相邻的花，并且要求相邻花的颜色不同，求方案数。(1≤n,m≤1e9,1≤k≤1e6,k≤n,m)思路

henucm·2020-09-16 18:19

洛谷 P2257 YY的GCD 莫比乌斯反演

P2257YY的GCD学习数论之莫比乌斯反演、杜教筛推荐：peng-ym推理：令：我们要求的是：令显然F(x)很容易求：我们反演一下：假设n#definelllonglongusingnamespacestd

一只叫橘子的猫·2020-09-16 18:41

[bzoj3601] 一个人的数论 [莫比乌斯反演+高斯消元]

题面传送门思路这题妙啊先把式子摆出来$f_n(d)=\sum_{i=1}^n[gcd(i,n)==1]i^d$这个$gcd$看着碍眼，我们把它反演掉$f_n(d)=\sum_{i=1}^n\sum_{j|i,j|n}\mu(j)i^d=\sum_{j|n}\mu(j)\sum_{i=1}^{\frac{n}{j}}(ij)^d=\sum_{j|n}\mu(j)j^d\sum_{i=1}^{\fra

aiou7071·2020-09-16 18:17

【数论】bzoj4174tty的求助

题目链接（ps：这是一道权限题…）计算∑n=1N∑m=1M∑k=0m−1⌊nk+xm⌋其中N，M#include#defineLLlonglongint#definemod998244353#defineMAXN500005usingnamespacestd;intn,m,prime[MAXN],tot;LLu[MAXN];boolflag[MAXN];doublex;inlinevoidinit

RZXZ·2020-09-16 18:44

推荐频道

【数论】--概率与期望