数论：容斥第33页

中国剩余定理(CRT)及其扩展(EXCRT)详解

ailanxier·2020-09-14 05:19

【BZOJ3616】War，KD树+bitset压位

思路：一开始想到顺序无关、轮数独立啥的，想用DP来做，但发现怎么定义状态都有后效性，想套容斥发现也没有什么用，后来聪哥提供思路(baochuzhengjie)，对于一些概率期望的题目，除DP外还可以考虑每一个元素的贡献沿着这种思路想一下

iamxym·2020-09-14 05:23

买不到的数目（数论）

买不到的数目小明开了一家糖果店。他别出心裁：把水果糖包成4颗一包和7颗一包的两种。糖果不能拆包卖。小朋友来买糖的时候，他就用这两种包装来组合。当然有些糖果数目是无法组合出来的，比如要买10颗糖。你可以用计算机测试一下，在这种包装情况下，最大不能买到的数量是17。大于17的任何数字都可以用4和7组合出来。本题的要求就是在已知两个包装的数量时，求最大不能组合出的数字。输入：两个正整数，表示每种包装中糖

smiling~·2020-09-14 05:03

蓝桥杯 - 买不到的数目（动态规划 + 数论）

历届试题买不到的数目时间限制：1.0s内存限制：256.0MB问题描述小明开了一家糖果店。他别出心裁：把水果糖包成4颗一包和7颗一包的两种。糖果不能拆包卖。小朋友来买糖的时候，他就用这两种包装来组合。当然有些糖果数目是无法组合出来的，比如要买10颗糖。你可以用计算机测试一下，在这种包装情况下，最大不能买到的数量是17。大于17的任何数字都可以用4和7组合出来。本题的要求就是在已知两个包装的数量时，

zzuspy·2020-09-14 04:43

埃拉托斯特尼筛法算法复杂度（n*lnlnn）的证明（[欧拉数学]素数倒数之和）

上一篇文章我通过欧拉数学的方式简单地讲了数论中的“黎曼ζ函数”和“金钥匙”。

Self-Discipline·2020-09-14 04:06

快速幂 + 快速幂取模

HDU5363KeySet（二项式定理）AcWing875.快速幂（模板）POJ1995RaisingModuloNumbers洛谷P1226快速幂||取余运算（模板）AcWing1289.序列的第k个数洛谷P3197越狱（容斥

WA-Accepted·2020-09-14 04:35

ACM数论----欧几里得算法与拓展欧几里得

一.首先介绍一下什么叫欧几里得算法欧几里得算法又称为辗转相除法，用于求两个自然数数的最大公约数，若有负数，全变为正数再运算，这里直接给出代码：（1）非递归版：intgcd(inta,intb){//returnb==0?a:gcd(b,a%b);if(a当然有关系了，ax+by=c是不一定有解的，但是由于拓展欧几里得定律：对于不完全为0的非负整数a，b，gcd(a,b)表示a,b的最大公约数，必定

阿阿阿安·2020-09-14 04:53

【数论】[HNOI2008]越狱

题目一开始没思路后来看题解学会了先考虑所有情况每个人都有m种可能的宗教所以总方案数为mn然后是考虑不能越狱的因为这个情况好像比较好考虑因为只要相邻两个不一样就可以所以假如第一个人可以是m个宗教那第2个人到第n个人一定只能有m-1种所以这种情况的方案数就是m∗(m−1)(n−1)然后两边做差就可以了需要用到快速幂边搞边取膜注意最后一定要有一个操作以防减出来负数因为都是取膜后的数还不一定谁大谁小代码如

ars4me·2020-09-14 04:30

SPOJ FACVSPOW(数论+二分)

传送门：http://www.spoj.com/problems/FACVSPOW/en/题意：给出一个数a，要求找到最小的n使得n!>a^n题解：先将不等式两边取对数得到ln(1)+ln(2)+...+ln(n)>nln(a)可以发现ausingnamespacestd;constintMX=3e6+5;doublesum[MX];voidinit(){sum[0]=0;for(inti=1;i

丿残念灬·2020-09-14 04:15

hdu 4704 Sum（费马小定理）解题报告

ProblemDescriptionSampleInput2SampleOutput2Hint1.ForN=2,S(1)=S(2)=1.2.Theinputfileconsistsofmultipletestcases.费马小定理是数论四大定理之一

小夥·2020-09-14 04:08

快速幂取余

小余最近迷恋上了数论，他认为大素数的判定是非常迷人的，但是他的程序总是数一大就死循环==，结果找了一晚上，发现原来是求a*b%c这个环节上出现了错误，聪明的你能帮助可怜的小余求出这个问题吗？

ITAK·2020-09-14 04:56

数论初步之快速幂取模

快速幂的写法完全是我自己完成的哦,你们不要跟我强功,呵呵,其实是自己找不到,呵呵;没事自己写的感觉还不错呢.快速幂取模就是用到了线性取模,呵呵.很简单的,.现在贴出我的代码:/**输入正整数a,n和m,输出a^n%m的值,a,n,m#include#include#include/**运用二分,也就是分治法,快速求幂;*/usingnamespacestd;longlongx=1;longlong

AledaLee·2020-09-14 04:56

luogu P3197 [HNOI2008]越狱题解（数论）

题目链接：P3197[HNOI2008]越狱思路：这题第一眼看见以为是一个dp，后来发现好像也是道数学题。我们发现这道题可以来求补集和全集，补集就是没有越狱的情况，第一间房子的情况有种，那么第二个房间就有种情况，第三个房间也有种情况，所以根据乘法原理没有越狱的情况就有种。还可以计算出全集，第一个房间有种情况，第二个也有种情况，以此类推还是根据乘法原理我们就会发现全集是种。所以答案就是全集减去补集就

阿嚏_·2020-09-14 04:49

【luogu2759】奇怪的函数（二分，数论）

一个换底公式解决的事https://www.luogu.org/problemnew/show/P2759题摘使得x^x达到或超过n位数字的最小正整数x是多少？求x^x次方的位数,可用x*log10(x)+1#include#include#definelllonglongusingnamespacestd;lln;intpd(llx){return(longlong)(x*log10(x)+1)

fl_334·2020-09-14 03:55

第十五届中北大学算法与程序设计竞赛（公开赛）

文章目录A-俄罗斯方块（模拟）B-签到题C-港口（差分）*E-简单的线性代数（数论）-找规律蒙的G数位操作1（模拟）K-签个到（思维）-忘记特判比赛入口A-俄罗斯方块（模拟）用数组c[]表示每一列能放方块的最大值

coder370·2020-09-14 03:27

完全平方数(数论)

完全平方数题目概述:从1−N中找一些数乘起来使得答案是一个完全平方数,求这个完全平方数最大可能是多少(答案取模100000007)(数据说明:对于20%的数据,1≤N≤100.对于50%的数据.1≤N≤5000.对于70%的数据,1≤N≤10^5.对于100%的数据,1≤N≤5×10^6.)题目分析:首先考虑什么是合法解——完全平方数,那么对其进行质因数分解之后得到的每种质因数的个数应当是偶数个.

The_useless·2020-09-14 03:27

数论变换 - 莫比乌斯反演篇

2020.8.17今天第一次屁股坐在椅子上学莫比乌斯反演的一天。本来是想让队友学的，现在队友不知道换了多少人了，也没几个真的能靠得住的，还都得看自己。其实学到现在这个程度，除了一些极其吃天赋的问题，基本只有没认真学和学会两种内容，自己原先觉得不会，但是只要肯坐在椅子上学那是能学个十有八九的，所以还是得好好学的。这个莫比乌斯函数就是用来加速求gcd，lcm，约数之类问题的问题。对于区间[1,i]和[

tiany7·2020-09-14 03:50

数论->快速幂

数论基础快速幂算法介绍算法利用了二分的思想，可以达到O(logn)。

_无尽夏_·2020-09-14 03:16

数论——快速幂算法快速计算a^b mod c的值

//快速计算(a^p)%m的值__int64FastM(__int64a,__int64p,__int64m){if(p==0)return1;__int64r=a%m;__int64k=1;while(p>1){if((p&1)!=0){k=(k*r)%m;}r=(r*r)%m;p>>=1;}return(r*k)%m;}

BrotherDongDong·2020-09-14 03:45

luogup1226取余运算||快速幂（数论）

输入b，p，k的值，求b^pmodk的值。其中b，p，k*k为长整型数。varb,p,k,t:int64;beginreadln(b,p,k);write(b,'^',p,'mod',k,'=');b:=bmodk;t:=1;while(p<>0)dobeginif(pmod2=1)thent:=t*bmodk;p:=pdiv2;b:=(b*b)modk;end;write(t);end.

fl_334·2020-09-14 03:36

Luogu P1226 取余运算||快速幂(数论,分治)

P1226取余运算||快速幂题目描述输入b，p，k的值，求b^pmodk的值。其中b，p，k*k为长整型数。输入输出格式输入格式：三个整数b,p,k.输出格式：输出“b^pmodk=s”s为运算结果输入输出样例输入样例#1：2109输出样例#1：2^10mod9=7这是一道很有趣的水题，如果知道公式。一般求解会溢出，导致答案错误。这里介绍取模的一个公式:a*b%k=(a%k)*(b%k)%k.在我

banban8290·2020-09-14 03:25

初等数论笔记第二讲数论函数（2,3,4节）

Mathcha传送门https://www.mathcha.io/editor/L0XOTnjTGxs6mCDm

随便写写~·2020-09-14 01:55

总结（10.27）

第一题的数论要尽量去想但是一个小时是最多的了，想不出来就拿了暴力分赶紧走，不要觉得自己想得出来。。。2.DP题不是每道题都很难。。其实

v1Ncen70·2020-09-14 00:10

Bzoj2839:集合计数:组合数学+容斥

题目链接:集合计数答案是含有至少k个的-至少k+1个的+至少k+2个的……从n个数中选出k个作为交集中的数,是C(n,k)，这样的集合共有2^(2^(n-k))-1个2^(n-k)是包含选定的k个数的可选集合的数量，选取方案有2^(2^(n-k))-1个(不能有空集否则无法保证k个元素)所以ans=C(n,k)*C(k,k)*(2^(2^(n-k))-1)-C(n,k+1)*C(k+1,k)*2^

TheWolfWhistlingSong·2020-09-14 00:18

bzoj 3812 状压dp 容斥原理

题意：一个n个点m条边的有向强连通图，去掉一些边使其仍然强连通，求方案数。以前做的题，现在看已经不知道自己在写什么了。写一点题解。如果一个图缩点后变成一个有多个点的DAG，那么这玩意一定不连通。设f[i]表示拆边使集合i强连通的方案数，g[i]表示i集合的点缩点后成为奇数个彼此没有边的点的方案数，p[i]表示缩成偶数个彼此没有边的点的方案数。对于g[i]和p[i]，枚举所在集合序号最大的点的所属强

make_it_for_good·2020-09-14 00:44

【BZOJ】4804 欧拉心算莫比乌斯函数+欧拉函数+数论分块

题目传送门来来来，推式子啦：∑i=1n∑j=1nϕ(gcd(i,j))=∑i=1n∑j=1n∑d=1n[gcd(i,j)=d]×ϕ(d)=∑d=1n(ϕ(d)×∑i=1⌊nd⌋∑j=1⌊nd⌋[gcd(i,j)=1])然后我们来证明一个结论：∑ni=1∑nj=1[gcd(i,j)=1]=2×∑ni=1ϕ(i)−1证明：首先我们考虑∑ni=1∑ij=1[gcd(i,j)=1]=∑ni=1ϕ(i)，这

Chester_King·2020-09-14 00:37

墨托卡投影(Mercator projection)

由于数论中直线有无数的点构成，墨托卡投影也不能表示出极点。它是一个用于地球的平

Flueg·2020-09-14 00:10

Learning：数论（三）欧拉函数

定义：对正整数，欧拉函数是小于等于的数中与互质的数的数目，又称函数。例如。引理：如果为某个素数，则。如果为某个素数的幂次。函数为积性函数。设为正整数的素数幂乘积表达式，则：。证明都不难，自己推一推吧。由引理1,2,3，我们不难可以想到怎么对欧拉函数进行线性筛。代码如下：intn,m,cnt,ans,x,y,phi[N],prm[N],vis[N];voidinit(){phi[1]=1;for(i

ezoiHQM·2020-09-14 00:09

[题解] NOIP2012 国王游戏（数论+高精）

题目分析：我们对于国王身后的两个点来分析队列可能是这样的：*LeftRightking:a0a0b0b0p1a1a1b1b1p2a2a2b2b2那么我们计算可得ans1ans1=max(a0b1,a0∗a1b2)max(a0b1,a0∗a1b2)队列也有可能是这样的*LeftRightking:a0a0b0b0p2a2a2b2b2p1a1a1b1b1那么我们计算可得ans2ans2=max(a0b

Chlience·2020-09-14 00:02

Miller-Rabin

对素数的研究属于数论范畴，你可以看到许多数学家没事就想出一些符合某种性质的素数并称它为某某某素数。整个数论几乎就围绕着整除和素数之类的词转过去转过来。

weixin_30566149·2020-09-13 23:48

[SinGuLaRiTy-1002] Miller Rabin Prime Judge 米勒·罗宾素数判定法

ByWenJian【背景】数论学家利用费马小定理研究出了多种素数测试办法，Miller-Rabbin素数测试算法是其中较快的一种。

SinGuLaRiTy-CQWJ·2020-09-13 22:45

miller_rabin学习笔记数论

首先介绍一下miller_rabin算法。miller_rabin是一种素性测试算法，用来判断一个大数是否是一个质数。miller_rabin是一种随机算法，它有一定概率出错，设测试次数为ss，那么出错的概率是4−s4−s，至于为什么我也不会证明。我觉得它的复杂度是O(slog2n)O(slog2n)，因为你要进行ss次，每次要进行一次快速幂，每次快速幂要lognlogn次快速乘，每次快速乘又是l

forever_shi·2020-09-13 21:58

数论原根及其求法

定义2：若G为n的原根，则当gcd(i,φ(n))==1，Gi也为n的原根可以通过先求出n的最小原根，来枚举得到其他的合法原根。#include#include#include#definelllonglongusingnamespacestd;constintmaxn=1e6+10;boolvis[maxn]={0};intpri[maxn/10]={0};intphi[maxn]={0};in

Glory_g·2020-09-13 20:51

ACM-ICPC 2018 南京赛区网络预赛--J Sum（数论）

题意：6=1*6=2*3=3*2=6*1，12=2*6=6*2（12可以被2^2整除），这样f(6)=4，f(12)=2，求f(1)...f(n)的和。题解：其实我们把数字拆成质数相乘就会发现只要这个质数的个数不超过2个就有值，比如6的有2,3，即6=1*2^1*3^1,24=1*2^3*3^1。而且我们可以发现数字如果是质数则f(n)=2。我们假设n=p^i*x，如果i==2则f(n)=f(x)

曾英俊·2020-09-13 20:15

中国剩余定理

数论太弱啦，还是来点基础的吧..终于不用担心小学奥数题不会做啦……中国剩余定理是解决线性模方程组的东西x≡a1(modp1)x≡a2(modp2).....考虑将两个方程合并x=k1∗p1+a1x=k2

ws_yzy·2020-09-13 20:47

2020年SDUTACM暑假集训阶段总结

2020年SDUTACM暑假集训阶段总结9#501计科1902杜广优一、学习内容–暑假新学了哪些专题，每个专题刷了多少题(1)背包进阶6(2)字典树3(3)线段树3(4)线段树扫描线1(5)数论基础5(

SDUT_you·2020-09-13 20:54

ACM-ICPC 2018 南京赛区网络预赛 J-Sum

ACM-ICPC2018南京赛区网络预赛J-Sum题意思路一（欧拉筛）代码一思路二（数论）代码二题目链接：ACM-ICPC2018南京赛区网络预赛J-Sum题意 \;\;\;\;f[i]f

nero_my·2020-09-13 20:34

一些算法学习的推荐博文阅读（数论居多，图论没有）

weixin_30652491·2020-09-13 20:04

多项式之快速傅里叶变换(FFT)/数论变换(NTT)/常用套路【入门】

原文链接https://www.cnblogs.com/zhouzhendong/p/Fast-Fourier-Transform.html多项式之快速傅里叶变换(FFT)/数论变换(NTT)/例题与常用套路

weixin_30522183·2020-09-13 20:28

扩展中国剩余定理

扩展中国剩余定理算法作用对于给定的一个同余方程组：x≡c[i](modm[i]){x\equivc[i]\pmod{m[i]}}x≡c[i](modm[i])初始技巧1、如何使用扩展欧几里得算法求逆元2、简单数论解决方法第一眼

EasternCountry·2020-09-13 19:25

浅析中国剩余定理（CRT）

中国剩余定理用来求解同余方程组的最小非负整数解，其中都互质首先让M等于所有的最小公倍数，对于求解每一个的方程先设一个，再求解其逆元则会有一组最小解其通解就是如果没有看懂，可以看详细求解同余方程这一篇博客代码数论大法好

C20201018·2020-09-13 19:44

ACM数论----中国剩余定理与拓展中国剩余定理

一.问题引入：在《孙子算经》中有这样一个问题：“今有物不知其数，三三数之剩二（除以3余2），五五数之剩三（除以5余3），七七数之剩二（除以7余2），问物几何？”这个问题称为“孙子问题”，该问题的一般解法国际上称为“中国剩余定理”。（壮哉我大中华）这个问题如何求解呢？来看一我一步步推导：1.假设n1%3=2,n2%5=3,n3%7=2;那么如果要使,(n1+n2)%3=2,那么就要求n2也是3的倍数

阿阿阿安·2020-09-13 19:17

模逆矩阵

整数a对同余n之乘法模逆元是指满足以下公式的整数b乘法模逆元又称为数论倒数，其实可以看作是普通倒数在模算术中的推广。同理，乘法模逆矩阵可以看作是普通逆矩阵在模算术中的推广。

天雨龙马·2020-09-13 19:17

HDU5838 G.Mountain(容斥+dfs+状压)

链接nnn和mmm范围很小,考虑状压但是因为谷底的存在不方便状压所以我们定义dp[i][j]dp[i][j]dp[i][j]为从小到大填了iii个数,谷底填过的状态是jjj这样当我们需要把数字i+1i+1i+1填上去的话Ⅰ.可以填充在一个没有被填过的谷底Ⅰ.可以填充在一个没有被填过的谷底Ⅰ.可以填充在一个没有被填过的谷底Ⅱ.可以填充在一个没有被填过的非谷底Ⅱ.可以填充在一个没有被填过的非谷底Ⅱ.可

issue敲腻害·2020-09-13 19:30

数论之原根

简介：原根是数论中一个非常重要的概念，它在密码学中有着很广泛的应用。原根从直观上非常好理解，数g对与p是原根，则(g^i)%p的结果互不相同，其中，i∈[1,p-1],g∈[2,p-1]。

vufw_795·2020-09-13 18:20

【bzoj3129】【SDOI2013】方程【exLucas】【容斥】

洛谷链接bzoj链接题目大意给一个不定方程∑i=1nxi=m\sum_{i=1}^{n}x_i=m∑i=1nxi=m(xi>1x_i>1xi>1，xi∈N+x_i\inN^+xi∈N+)，求方程的解的数量，并满足：当0<i≤n10<i\leqn10#include#include#include#include#include#include#include#include

ezoixx118·2020-09-13 18:30

中国剩余定理（孙子定理）

具体解释过程可看：数论四大定理代码：voidex_gcd(lla,llb,ll&x,ll&y,ll&d){if(!

baodream·2020-09-13 18:49

P4149 距离为K的点对(最少边数) n=200000 点分治

这题数据范围变成了200000n^2就过不了同时要求求的是最少的边数不能容斥#includeusingnamespacestd;constintMAXN=2e5+5;constintMAXM=2e5+5

baichuan9723·2020-09-13 18:43

数论概论读书笔记 39.斐波那契与线性递归序列

斐波那契与线性递归序列比内公式斐波那契序列FnFn用递归公式描述如下：F1=F2=1，Fn=Fn−1+Fn−2，n=3,4,5,...F1=F2=1，Fn=Fn−1+Fn−2，n=3,4,5,...则斐波那契序列的第nn项可用公式Fn=15–√{(1+5–√2)n−(1−5–√2)n}Fn=15{(1+52)n−(1−52)n}给出。斐波那契序列模m对m=2,3,4,5,6...m=2,3,4,5

Feynman1999·2020-09-13 18:14

求原根（模板）

数论令人头秃，原理就看看别人博客吧：数论之原根#include#include#include#include#include#include#include#include#include#include

别动我的白羊毛·2020-09-13 18:45

推荐频道

数论：容斥