【数论】--扩展欧几里得第45页

【bzoj4872】[Shoi2017]分手是祝愿数论+期望dp

题目描述ZeitundRaumtrennendichundmich.时空将你我分开。B君在玩一个游戏，这个游戏由n个灯和n个开关组成，给定这n个灯的初始状态，下标为从1到n的正整数。每个灯有两个状态亮和灭，我们用1来表示这个灯是亮的，用0表示这个灯是灭的，游戏的目标是使所有灯都灭掉。但是当操作第i个开关时，所有编号为i的约数（包括1和i）的灯的状态都会被改变，即从亮变成灭，或者是从灭变成亮。B君发

weixin_30795127·2020-08-17 10:24

【数论求逆元】东华大学2020年程序设计竞赛（同步赛）B题： A Number Theoretical Problem

题目传送门简单说一下题目意思，就是给你两个数，一个整数yyy和一个质数ppp，然后让你再另外找一个数xxx满足(x∗y(x*y(x∗y)%p==1p==1p==1，如果这个xxx不存在的话，就输出−1-1−1;存在的话就输出xxx%ppp，好题意就是这么个意思，接下来展开思路的分析：首先对于x不存在的情况就只有一种，即：若yyy%p==0p==0p==0，则xxx一定不存在，因为余数一直就是000

biu~跃哥冲冲冲·2020-08-17 10:17

ACM数论----高斯消元法解方程组（半成品）

#include#include#include#definemaxn100intnum[maxn][maxn+1];//增广矩阵doubleanswer[maxn];//解集voiddisplay(intN)//打印结果{inti;for(i=0;iabs(num[max_r][l]))//找到该列中最大的元素max_r=k;}if(max_r!=r){for(k=r;k=0;k--)//否则有

阿阿阿安·2020-08-17 10:04

Double Deep Q Network 算法

（Y值通过Q评估值、Q目标值的双网计算）符号：：表示Q目标值神经网络（target_net）的网络参数：表示Q评估值神经网络(eval_net)的网络参数论文：

段智华·2020-08-17 10:43

除法取模(费马小定理+扩展欧几里得)

除法取模以下只说这两个方法1）费马小定理a^(p-1)==1(mod)pa*a^(p-1)=1%Ma/b%mod=a*b^(mod-2)%mod;只有p,mod为素数时，而一般都为素数。b^(mod-2)一般用快速幂typedeflonglongll;llquick_pow(lla,lln){llans=1;while(n){if(n&1){ans*=a;ans%=mod;//mod自己设置}a=

DeathYmz·2020-08-17 10:11

牛客练习赛48 B.小w的a=b问题(差分+数论)

题目一个a[]数组，长度为n，第i个元素ai，另一个是b[]数组，长度为m，第j个元素bj。先对每个元素求阶乘，再对数组求乘积，问两个数组最终所得结果是否相等1usingnamespacestd;typedeflonglongll;constintmaxn=1e5+10;constintmod=1e9+7;intt,n,m,v,dp[maxn];bitsetnp;boolok(){for(inti

Code92007·2020-08-17 10:59

除法取模与逆元/费马小定理

逆元一般用扩展欧几里得算法来求得，如果为素数，那

tony_bfx·2020-08-17 10:23

RSA

rsa算法原理根据数论，寻求两个大素数比较

ruigeling·2020-08-17 10:17

混合个人训练第四十六场（A,B题解） A 数论线代全家桶 (all） B: 看星星 (stars)

A数论线代全家桶(all)题意如图。

ldu 王跃·2020-08-17 10:44

Jumping Frog（时间限制: 1 Sec 内存限制: 128 MB）数论&思维

题目描述Afrogislocatedatthecoordinate(x1,y1).Hewantstogotothecoordinate(x2,y2).Hewillperformoneormorejumpstoreachhisdestination.Theruleofthejumpingisasfollows:Supposethefrogislocatedatthecoordinate(x,y);t

ldu 王跃·2020-08-17 10:39

Java实现求解两个正整数的最大公约数和最小公倍数

Java实现求解两个正整数的最大公约数和最小公倍数公约数：公因数，显示着若干个整数之间的数论关系。如果一个数同时是几个数的约数，称这个数为它们的“公因数”；公约数中最大一个的称为最大公因数。

阿浩喜欢敲代码·2020-08-17 07:38

幸运的袋子(穷举法、数学、数论)

1.题目来源链接：计算日期到天数转换来源：牛客网2.题目说明一个袋子里面有n个球，每个球上面都有一个号码(拥有相同号码的球是无区别的)。如果一个袋子是幸运的当且仅当所有球的号码的和大于所有球的号码的积。例如：如果袋子里面的球的号码是{1,1,2,3}，这个袋子就是幸运的，因为1+1+2+3>1*1*2*3你可以适当从袋子里移除一些球(可以移除0个,但是别移除完)，要使移除后的袋子是幸运的。现在让你

Y_puyu·2020-08-17 07:57

计算扩展欧几里得算法求1024位大整数a和b的系数x，y 执行时间

要点：大整数a、b均匀分布在0-2^1024次方内使用扩展欧几里得算法实现求a和b的最大公因子的同时求解ax+by=gcd（a,b）的系数x和y同时观察一下该递归算法的运行时间直接上代码：importjava.math.BigInteger

yky__xukai·2020-08-17 07:56

欧几里得算法和扩展欧几里得算法-java递归实现

定义已知整数a,b,扩展欧几里得算法可以在求得a、b的最大公约数的同时，能找到整数x、y，满足ax+by=gcd（a，b）分析设a>b（1）显然当b=0，gcd（a，b）=a。

yky__xukai·2020-08-17 07:56

数论入门篇

欧几里得定理gcd(a,b)=gcd(a%b,b)gcd(a,b)=gcd(a\%b,b)gcd(a,b)=gcd(a%b,b)(a>=b)给出简单证明:可以知道gcd(a,b)=b,gcd(a,b)=gcd(a%b,b)gcd(a,b)=gcd(a\%b,b)gcd(a,b)=gcd(a%b,b)。递归求解,当b=0时返回a即可。拓展性质:对于a>=b,gcd(a,b)=gcd(a−b,b)\\

悭吝韶华·2020-08-17 07:40

学习记录----2020.7.24

学习记录这几天看了一些关于驱动的东西，然后学了一些WinDbg的双机调试以及IDA使用WinDbg插件调试驱动，踩了很多坑然后之前还做了一个比赛的crackme，然后数学题太顶了，不懂数论就没做出来WinDbg

0x菜鸡·2020-08-17 06:49

【数论】Codeforces1027G X-mouse in the Campus

分析：显然，i->ix(modm)连一条边，则最终一定会形成若干个环，并且，环上每个点与m的gcd值必定相同。并且，gcd值相同的环大小也一定相同。所以，如果能算出对于所有数中，与m的gcd为d的个数f(d)f(d)f(d)，并算出相应的当gcd为d时的每个环的大小l(d)l(d)l(d)，那么答案就是∑f(d)l(d)\sum\frac{f(d)}{l(d)}∑l(d)f(d)很容易发现，f(d

616156·2020-08-17 06:49

Codeforces 1027G X-mouse in the Campus 数论+Pollard_rho

题意给定mm和xx，满足gcd(m,x)=1gcd(m,x)=1。现在把每个小于mm的整数都看作一个点，然后ii向ixix连边，问最后最少需要选出多少个点使得每个点的后继中至少有一个点被选。m≤1014m≤1014分析感谢sam队长教我做这题。首先因为gcd(m,x)=1gcd(m,x)=1，所以最后形成的图一定是若干个环，显然某个环上每个点与m的gcd都相等，且每个和m的gcd相等的环的大小都一

SFN1036·2020-08-17 06:42

SP5971 LCMSUM 数论

题面题目要我们求这个：\[\sum_{i=1}^nlcm(i,n)\]开始化式子：\[\sum_{i=1}^{n}\frac{i*n}{gcd(i,n)}\]\[\sum_{d|n}\sum_{i=1}^{\frac{n}{d}}i*n[gcd(i,\frac{n}{d})=1]\]\[n*\sum_{d|n}\sum_{i=1}^{d}i[gcd(i,d)=1]\]注意那个\(\sum_{i=1

baobi9069·2020-08-17 05:49

ZOJ 3593 One Person Game 【带简单处理的扩展欧几里得】

题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3593题意：给定一维坐标的出发点的终点，每次可以向左或向右走a,b或a+b的距离，问是否可达，如可达求最小步数分析：这种每步有多种决策，初始感觉可能是DP，然而这个数据大小已经覆盖了整个int的大小，dp显然不现实。考虑到最终总是要求最小步数，某个对于a,b这两种步

babao9492·2020-08-17 05:16

Sum of Digits(数论+思维)

如果要让答案最小，显然我们是希望不进位的，因为一旦进位，会让对答案的贡献减少9*进位次数(这样明显会让答案变大).所以我们优先考虑不进位的。其实只要(n-(1+k)*k/2)%(k+1)即可。((1+k)*k/2)是后面k个数个位上比答案大的差值。那如果不满足就要考虑是否存在进位情况产生了,为了让答案最小,结论是:每个数最多只进一次位。那么如果答案是两位数以上,让十位最大为8即可。#include

悭吝韶华·2020-08-17 05:02

牛客练习赛23（题解）

这题刚看到的时候还以为是DP，或者是数论（前几天做这类题做疯了），但一看数据范围……果断排除。重新看了一下托米有的硬币的面值，这不就是贪心吗，每次能用面值大的就用面值大的，然后就很开心地去写代码了。

lahlah_·2020-08-17 04:12

数论基础知识补充

一、阶数论术语其在数论中的定义为：设a，p是整数，有：anΞ1(modp)可以使上式成立的最小正整数n叫做a模p的阶。

黎明前的沉睡forever·2020-08-17 04:15

ACM数论基础之一_______质因数分解

也就是每一个数都可以分解成素数相乘的结果只是简单的代码，相信大家都看得懂吧，于是我就不啰嗦了#include#include#include#include#include#include#includeusingnamespacestd;vectorv;intmain(){intn;while(scanf("%d",&n)!=EOF){v.clear();intok=n;for(inti=2;o

zhao9642·2020-08-17 04:50

codeforces 1295D Same GCDs 数论欧拉函数

https://vjudge.net/problem/CodeForces-1295D题目大意：给定n、mn、mn、m，对于任意x(0#defineINF0x3f3f3f3fusingnamespacestd;typedeflonglongll;lleuler(lln)//得到φ(n){llret=n,i;for(i=2;i*i1)//素数情况ret-=ret/n;returnret;}llgcd

csu_xiji·2020-08-17 04:53

ACM数论总结

断断续续的学习数论已经有一段时间了，学得也很杂，现在进行一些简单的回顾和总结。学过的东西不能忘啊。。。

xieshimao·2020-08-17 04:53

HDU - 5685 Problem A 求逆元（扩展欧几里得方法）

这里采用扩展欧几里得的方法求逆元预处理，求得每个点的Hash值，然后计算x到y段的答案的时候只要求（H[y]/H[x-1]%mod）就好了这里用到了除法但是在模运算世界世界里不能用除法，所以要求H[x-

冰冰的小宝贝·2020-08-17 04:15

20181220第二周周训思路整理

ProblemListBZOJ5027:数学题Clickme扩展欧几里得算法BZOJ1977:次小生成树Clickme次小生成树HDU1536:S-NimClickme博弈论BZOJ4919:大根堆Clickme

wlxsq·2020-08-17 04:17

最大公约数和最小公倍数的经典求法（C++）及例题实战练习

最大公约数和最小公倍数的经典求法（C++）最小公倍数：数论中的一种概念，两个整数公有的倍数成为他们的公倍数，其中一个最小的公倍数是他们的最小公倍数，同样地，若干个整数公有的倍数中最小的正整数称为它们的最小公倍数求最小公倍数算法

笑着的程序员·2020-08-17 04:56

逆元模板

(lla,llb){llres=1;while(b>0){if(b&1)res=res*a%mod;a=a*a%mod;b>>=1;}returnres;}//a的逆元=qpow(a,mod-2)2.扩展欧几里得

weixin_ajls·2020-08-17 03:18

快速求组合数

1同余同余是数论中的重要概念。给定一个正整数m，如果两个整数a和b满足a-b能够被m整除，即(a-b)/m得到一个整数，那么就称整数a与b对模m

-Coca·2020-08-17 03:25

辣鸡(ljh) NOIP模拟赛模拟平面几何数论化学相关（雾）

【题目描述】辣鸡ljhNOI之后就退役了，然后就滚去学文化课了。然而在上化学课的时候，数学和化学都不好的ljh却被一道简单题难住了，受到了大佬的嘲笑。题目描述是这样的：在一个二维平面上有一层水分子，请问形成了多少个氢键？这个二维平面可以看做一个类似棋盘的东西，每个格子可以容纳一个水分子，左下角的格子为(0,0)，这个格子右边的格子为(1,0)，上方格子为(0,1)，以此类推。辣鸡ljh当然不会做了

weixin_34376986·2020-08-17 03:30

weixin_30402343·2020-08-17 03:19

【数论】勒让德定理(笔记)

勒让德定理是用来求n!的素因子分解的一个定理,表述为:记素数p的指数为Lp,Lp(n!)=∑[n/pk](k>=1).维基百科有详细的证明一个要用到该定理的问题:uvalive5828求一个最大的数d,满足:1.能分解成不同的因子,这些因子不能超过n;2.是一个完全平方数.输入n输出d%1000000007分析:为满足完全平方数只需素分解后所有素数p的指数是偶数,因此这个最容易满足的条件我们放在最

weixin_30367169·2020-08-17 03:15

数论常用内容——数根

对于数根可能接触的不多，但是我觉得这里还是应该做一下简单的总结和介绍数根数根(又称数字根Digitalroot)是自然数的一种性质，每个自然数都有一个数根。数根是将一正整数的各个位数相加(即横向相加)，若加完后的值大于等于10的话，则继续将各位数进行横向相加直到其值小于十为止，最后得到的数字就是该数的数根数根的性质root(a+b)=root(root(a)+root(b))root(a*b)=r

tick_tokc97·2020-08-17 02:41

【ACM训练六】数论入门

素数问题素数判定枚举：2-n（时间复杂度：O（n））优化：2-[√n]（时间复杂度：O（√n））素数筛如何筛出2到n中所有素数？枚举：时间复杂度O（n√n）素数筛：通过预处理，做到每次O（1）判断一个数（某个范围内）是否为素数。埃氏筛：初始将所有数(从2开始)默认为素数，从小到大遍历每个数x，如果x是素数，那么将2x,3x,4x……划掉，即标记成非素数优化：每次改为划掉x2，x2+x，x2+2x…

三分天注定·2020-08-17 02:00

[kuangbin带你飞]专题十四数论基础

A-Bi-shoeandPhi-shoe——筛素数题意一个竹竿长度为p,它的score值就是比p长度小且与且与p互质的数字总数，比如9有1，2，4，5，7，8这六个数那它的score就是6。给你T组数据，每组n个学生，每个学生都有一个幸运数字，求出要求买nnn个竹子每个竹子的score都要大于或等于该学生的幸运数字，每个竹竿长度就是花费，求最小花费。思路ppp长度小且与且与ppp互质的数字总数就是

Miserable_ccf·2020-08-17 02:27

数论板子

returnans;}//求a与b的最大公约数，其中b>aintgcd(inta,intb){if(a>b)swap(a,b);while(a){intt=b%a;b=a;a=t;}returnb;}//扩展欧几里得

rain_falls·2020-08-17 02:14

hdu4828

个数放到这些格子中,每个格子放一个数,要求每一行后一个总比前一个大,每一列后面的比前面的大,即aij这本书课后习题的第二个习题.因为这里模数M=109+7M=10^9+7M=109+7是一个素数,可以用扩展欧几里得计算出逆元

2018slgys·2020-08-17 02:48

数论的基础入门（初读数论概论有感）（acm知识储备）

在寒假自己对自己的硬核知识进行了充电，这本书是学长极力推荐的，是专门给数学系的学生进行对数论进行全面认识和理解的入门书籍，虽然自己是计算机的，但是并不影响去阅读它，读的过程中就发现了，书中不少数学系的知识映入眼帘

洋-葱·2020-08-17 02:39

Two Divisors【GCD数论】

Youaregivennnintegersa1,a2,…,ana1,a2,…,an.Foreachaiaifinditstwodivisorsd1>1d1>1andd2>1d2>1suchthatgcd(d1+d2,ai)=1gcd(d1+d2,ai)=1(wheregcd(a,b)gcd(a,b)isthegreatestcommondivisorofaaandbb)orsaythatthere

Andres_Lionel·2020-08-17 02:00

数论基础 1

1.欧几里得算法（辗转相除法）递归写法：intgcd(inta,intb){if(b==0)returna;elsereturngcd(b,a%b);}更简洁的写法：intgcd(inta,intb){returnb?gcd(b,a%b):a;}扩展：求最大公倍数lcmlcm(a,b)=(a*b)/gcd(a,b);//为了防止a*b超longlong溢出//所以建议这样写lcm(a,b)=a/g

不想悲伤到天明·2020-08-17 02:57

ACM 数论入门题（附代码解释）

目录51Nod-1119机器人走方格V2（费马小定理）HDU2710MaxFactor（素数筛选）POJ2142TheBalance（扩展欧几里得）POJ1061青蛙的约会（扩展欧几里得）洛谷P1069

SpongeBob_Y·2020-08-17 02:24

【天梯】数论入门-1212 最大公约数

题目描述Description求两个数A和B的最大公约数。1usingnamespacestd;intgcd(inta,intb)//最大公约数{inti;for(i=a;i>=1;i--){if(a%i!=0)continue;elseif(b%i==0)returni;}}intmain(){inta,b,ans;cin>>a>>b;ans=gcd(a,b);cout<

dyyz_wlx·2020-08-17 01:52

数论板子(持续更新

boolissqr(__int128_tx)//开方{__int128_ty=(__int128_t)ceil(sqrt((longdouble)x));for(;y*yx;--y);returny*y==x;}booliscub(__int128_tx)//三次方{__int128_ty=(__int128_t)ceil(pow((longdouble)x,1.0/3));for(;y*y*yx

twh233·2020-08-17 01:20

拓展欧几里得算法

根据数论中的定理：gcd(a,b)=gcd(b,a%b)，我们可以得到ax+by=gcd(a,b)=gcd(b,a%b)=bx1+(a%b)y1而取模公式为：a%b=a-[a/b]*b([a

Racyevick·2020-08-17 01:37

数论入门1——最大公约数与最小公倍数

gcd(a,b)，就是求a和b的最大公约数lcm(a,b)，就是求a和b的最小公倍数欧几里德算法又称辗转相除法，是指用于计算两个正整数a，b的最大公约数。应用领域有数学和计算机两个方面。。下面记几个公式，证明以后再补。1.gcd(a,b)=gcd(b,a)=gcd(-a,b)=gcd(|a|,|b|)2.gcd(a,b)=gcd(b,amodb);3.a，b互质等价于gcd(a,b)=1;4.a*

小玖玖要加油吖·2020-08-17 01:28

数论题目入门基础积攒

基本性质若p|(a-b)，则a≡b(%p)。例如11≡4(%7)，18≡4(%7)(a%p)=(b%p)意味a≡b(%p)对称性：a≡b(%p)等价于b≡a(%p)传递性：若a≡b(%p)且b≡c(%p)，则a≡c(%p)运算规则模运算与基本四则运算有些相似，但是除法例外。其规则如下：(a+b)%p=(a%p+b%p)%p（1）(a-b)%p=(a%p-b%p)%p（2）(a*b)%p=(a%p*

折腾怪·2020-08-17 01:54

hdu5914 Triangle _数论

TriangleTimeLimit:2000/1000MS(Java/Others)MemoryLimit:65536/65536K(Java/Others)TotalSubmission(s):937AcceptedSubmission(s):596ProblemDescriptionMr.Froghasnsticks,whoselengthsare1,2,3⋯nrespectively.Wal

15wylu·2020-08-17 01:45

hdu 4828 - Grids

分析：组合，计数，卡塔兰数，数论。结果是卡塔兰数。这里要计算第1000000个卡特兰数，这个数字十分巨大，第10

小白菜又菜·2020-08-17 01:19

推荐频道

【数论】--扩展欧几里得

【bzoj4872】[Shoi2017]分手是祝愿 数论+期望dp