![](https://img-blog.csdnimg.cn/20201014180756724.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
------------数论-------------
努力写题的tyf
这个作者很懒,什么都没留下…
展开
-
求最大公约数的三种方法
最小公倍数=x*y/最大公约数1.辗转相除法:(第一种写法)int gcd(int x,int y){//x>y int z; while(x%y!=0){ z=x%y; x=y; y=z; } return y;}(第二种写法)int gcd(int x,int y){//x>y...原创 2018-03-08 20:02:37 · 526 阅读 · 0 评论 -
牛客OI赛制测试2 B路径数量 - 矩阵快速幂
题目描述给出一个 n * n 的邻接矩阵A.A是一个01矩阵 .A[i][j]=1表示i号点和j号点之间有长度为1的边直接相连.求出从 1 号点 到 n 号点长度为k的路径的数目.输入描述:第1行两个数n,k (20 ≤n ≤ 30,1 ≤ k ≤ 10)第2行至第n+1行,为一个邻接矩阵输出描述:题目中所求的数目 示例1输入复制4 2...原创 2018-09-10 19:30:31 · 185 阅读 · 0 评论 -
ACM-ICPC 2018 焦作赛区网络预赛 - G Give Candies - 规律+欧拉降幂
思路:找规律可以知道是2^(n-1)n的范围很大,所以用欧拉降幂公式降幂:欧拉φ函数:在数论中,对正整数n,欧拉函数是小于或等于n的正整数中与n互质的数的数目。此函数以其首名研究者欧拉命名,它又称为φ函数、欧拉商数等。φ(1)=1m=1e9+7是个质数,所以φ(1e9+7)=1e9+6给质数降幂后就愉快地用快速幂求解好了~代码如下:#include<i...原创 2018-09-16 09:58:47 · 200 阅读 · 0 评论 -
ACM-ICPC 2018 南京赛区网络预赛 J Sum - 线性筛+找规律
思路: 先利用线性筛筛出素数和每个数的最小质因数x然后地推出每个数的f函数值1、若i%(x*x*x)==0 =>02、若i%(x*x)==0 =>cnt[i]=cnt[i/(x*x)]3、若i%(x)==0 =>cnt[i]=2*cnt[i/x]代码如下:#include<iostream>#include<...原创 2018-09-05 20:50:55 · 148 阅读 · 0 评论 -
F - 子序列 - 容斥、组合数、费马小定理(详解)
链接:https://www.nowcoder.com/acm/contest/181/F来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K64bit IO Format: %lld题目描述给出一个长度为n的序列,你需要计算出所有长度为k的子序列中,除最大最小数之外所有数的乘积相乘的结果 输入描述:第...原创 2018-09-02 20:37:14 · 207 阅读 · 0 评论 -
组合数有关的公式及常用求和
转自:https://blog.csdn.net/bigtiao097/article/details/77242624转载 2018-09-02 16:24:20 · 4206 阅读 · 0 评论 -
模运算法则
模运算与基本四则运算有些相似,但是除法例外。其规则如下: (a + b) % p = (a % p + b % p) % p (a - b) % p = (a % p - b % p) % p (a * b) % p = (a % p * b % p) % p (a^b) % p = ((a % p)^b) % p推论: 若a≡b (% p),则对于任意的c,都有(a + c) ≡ ...转载 2018-08-25 20:25:56 · 1440 阅读 · 0 评论 -
Wannafly挑战赛22 - A - 裴蜀定理
链接:https://www.nowcoder.com/acm/contest/160/A来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K64bit IO Format: %lld题目描述有一个计数器,计数器的初始值为0,每次操作你可以把计数器的值加上a1,a2,...,an中的任意一个整数,操作次数不限(可以为...原创 2018-08-24 22:00:47 · 203 阅读 · 0 评论 -
Biorhythms - poj1006 - 中国剩余定理(互质)
BiorhythmsTime Limit:1000MS Memory Limit:10000K Total Submissions:146409 Accepted:47252 DescriptionSome people believe that there are three cycles in a person's life that s...原创 2018-08-27 20:12:46 · 146 阅读 · 0 评论 -
牛客OI赛制测试2 F假的数学游戏 - 斯特林公式+二分
链接:https://www.nowcoder.com/acm/contest/185/F来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K64bit IO Format: %lld题目描述输入一个整数X,求一个整数N,使得N!恰好大于XX。 输入描述:第一行:一个整数X输出描述:第一行:一...原创 2018-09-10 21:09:06 · 123 阅读 · 0 评论 -
牛客国庆集训派对Day2 A - 矩阵分块
题目描述深度学习算法很大程度上基于矩阵运算。例如神经网络中的全连接,本质上是一个矩阵乘法;而卷积运算也通常是用矩阵乘法来完成的。有一些科研工作者为了让神经网络的计算更快捷,提出了二值化网络的方法,就是将网络权重压缩成只用两种值表示的形式,这样就可以用一些 trick 加速计算了。例如两个二进制向量点乘,可以用计算机中的与运算代替,然后统计结果中 1 的个数即可。然而有时候为了降低压缩带来的误...原创 2018-10-03 11:20:25 · 255 阅读 · 0 评论 -
牛客练习赛44 - D小y的盒子 - 欧拉降幂
思路:列举几个就可以发现规律:n=1 4^1-1n=2 4^2-1n=3 4^3-1或者可以这么想,比如n=3时,盒子里有:1 1 3 3四个物品,对于每个物品有“取”或者“不取”两种状态,所以,一共有2^4-1种取法(减1是因为不能全不取),然后对于每种取法,都能保证取出来的数之和是在[1,8]之间的,所以,对应的是2^(2n)-1种取法,也就是4^n-...原创 2019-04-20 10:12:41 · 293 阅读 · 0 评论 -
洛谷P1414 又是毕业季II - 思维
题目背景“叮铃铃铃”,随着高考最后一科结考铃声的敲响,三年青春时光顿时凝固于此刻。毕业的欣喜怎敌那离别的不舍,憧憬着未来仍毋忘逝去的歌。1000多个日夜的欢笑和泪水,全凝聚在毕业晚会上,相信,这一定是一生最难忘的时刻!题目描述彩排了一次,老师不太满意。当然啦,取每位同学的号数来找最大公约数显然不太合理。于是老师给每位同学评了一个能力值。于是现在问题变为,从n个学生中挑出k个人使得他们的...原创 2019-05-04 10:10:33 · 227 阅读 · 0 评论 -
西南民族大学第十届校赛(同步赛)- H - 差分标记
题目描述HJ养了很多花(99999999999999999999999999999999999盆),并且喜欢把它们排成一排,编号0~99999999999999999999999999999999998,每天HJ都会给他的花浇水,但是他很奇怪,他会浇n(1 <= n <= 2 * 105)次水,每次都会选择一个区间[l, r],(0 <= l <= r <= 10...原创 2019-01-01 20:55:06 · 290 阅读 · 0 评论 -
求逆元的3种方法
//扩展欧几里得求逆元//a*x=1(mod m)=>a*x+m*y=1int extgcd(int a,int b,int &x,int &y){ int d=a; if(b!=0){ d=extgcd(b,a%b,y,x); y-=(a/b)*x; } else { x=1;y=0; ...原创 2018-09-22 20:47:32 · 621 阅读 · 0 评论 -
牛客国庆集训派对Day4 - G、I、H
G - 区间权值题目描述小 Bo 有 n 个正整数 a1..an,以及一个权值序列 w1…wn,现在他定义 现在他想知道 的值,需要你来帮帮他你只需要输出答案对 109+7 取模后的值输入描述:第一行一个正整数 n第二行 n 个正整数 a1..an第三行 n 个正整数 w1..wn输出描述:输出答案对 109+7 取模后的值 示例1输入复制...原创 2018-10-04 19:45:47 · 284 阅读 · 3 评论 -
牛客国庆集训派对Day3 - H
题目描述魔方国有n座城市,编号为。城市之间通过n-1条无向道路连接,形成一个树形结构。澜澜打算在魔方国进行m次旅游,每次游览至少一座城市。为了方便,每次旅游游览的城市必须是连通的。此外,澜澜希望游览所有城市恰好一次。澜澜想知道有多少种旅游方案满足条件,两个方案不同当且仅当存在某一次旅游游览了不同的城市。澜澜不会数数,所以只好让你来帮他数方案。输入描述:第一行一个整数t表示数据组...原创 2018-10-04 10:59:46 · 160 阅读 · 0 评论 -
ckoj1119 - 斐波那契数列 - 循环节
1119: 斐波那契数列题目描述斐波那契数列F可以按如下方式表示:F0=0, F1=1, Fi=Fi-1+Fi-2 (i>1)给定T个整数xi和P, 请求出Fxi mod P的值输入数据第一行两个整数T, P如上所述 接下来T行每行一个整数xi 同样如上所述输出数据T行, 每行一个整数表示对应数据的答案样例输入3 7 3 5 233样例输出2...原创 2018-10-04 09:08:01 · 242 阅读 · 0 评论 -
牛客国庆集训派对Day1 - C、J、L
C - Utawarerumono题目描述算术是为数不多的会让Kuon感到棘手的事情。通常她会找Haku帮忙,但是Haku已经被她派去买东西了。于是她向你寻求帮助。给出一个关于变量x,y的不定方程ax+by=c,显然这个方程可能有多个整数解。Kuon想知道如果有解,使得p2*x2+p1*x+q2*y2+q1*y最小的一组整数解是什么。为了方便,你只需要输出p2*x2+p1*x+q2*y2...原创 2018-10-02 19:45:31 · 211 阅读 · 0 评论 -
X问题 - hdu1573 - 中国剩余定理
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1573思路:就是用中国剩余定理啦,再加一点变化。题中让求的是小于等于N的x的个数,并且要是正整数。我们用中国剩余定理把方程组:x%m1=a1x%m2=a2……x%mn=an合并成一个方程x%m=a,即x=a(mod m)我们解得的a是最小值,m*i+a都是满足...原创 2018-08-27 15:49:32 · 164 阅读 · 0 评论 -
Strange Way to Express Integers - poj2891 - 中国剩余定理(不互质)
Strange Way to Express IntegersTime Limit: 1000MS Memory Limit: 131072K Total Submissions: 20313 Accepted: 6858 Description Elina is reading a book written by Rujia Liu, whic...原创 2018-08-27 15:24:12 · 103 阅读 · 0 评论 -
有关素数的基础算法代码总结
1.素性测试给定整数n,判断n是不是素数//素性测试 bool is_prime(int n){ for(int i=2;i*i<=n;i++){ if(n%i==0)return false; } return n!=1;}//约数枚举vector<int> divisor(int n){ vector<int> res; for(i...原创 2018-08-05 11:07:55 · 240 阅读 · 0 评论 -
欧几里得算法、扩展欧几里得算法
欧几里得算法1.引入问题:求线段上格点的个数给定平面两个格点P1=(x1,x2)和P2=(x2,y2),线段P1P2上,除P1,P2外还有几个格点。格点:横纵左边均为整数的点。答案:|x1-x2|与|y1-y2|的最大公约数-1(若点P1,P2重合,答案为0)原因:设d1 =|x1-x2|,d2 = |y1-y2|,d1=a*b,d2=b*c,b是二者的最大公约数...原创 2018-08-05 10:31:52 · 1012 阅读 · 0 评论 -
Wannafly挑战赛13 - B Jxc的军训
链接:https://www.nowcoder.com/acm/contest/80/B来源:牛客网时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K 64bit IO Format: %lld题目描述 在文某路学车中学高一新生军训中,Jxc正站在太阳下站着军姿,对于这样的酷热的阳光,Jxc 表示非常不爽。 Jxc将天空看做一个n...原创 2018-04-07 10:26:16 · 442 阅读 · 0 评论 -
n的阶乘末尾有几个0
链接:https://www.nowcoder.com/acm/contest/93/E来源:牛客网时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K 64bit IO Format: %lld题目描述 这个问题很简单,就是问你n的阶乘末尾有几个0? 输入描述:输入第一行一个整数T(1<=T<=100),代表测试组数接下来T行,每行...原创 2018-04-05 18:57:07 · 991 阅读 · 0 评论 -
牛客练习赛14 A-n的约数(详解)
链接:https://www.nowcoder.com/acm/contest/82/A来源:牛客网时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K 64bit IO Format: %lld题目描述 t次询问,每次给你一个数n,求在[1,n]内约数个数最多的数的约数个数 输入描述:第一行一个正整数t之后t行,每行一个正整数n输出描述:输出t...原创 2018-04-04 11:35:15 · 236 阅读 · 0 评论 -
cmath常用库函数预览
cmath是c++语言中的库函数,其中的c表示函数是来自c标准库的函数,math为数学常用库函数。cmath中常用库函数:int abs(int i);//返回整型参数i的绝对值double fabs(double x);//返回双精度参数x的绝对值long labs(long n);//返回长整型参数n的绝对值double exp(double x);//返回指数函数e^x的值double lo...转载 2018-03-22 10:30:52 · 319 阅读 · 0 评论 -
L1-028. 判断素数 - PAT
本题的目标很简单,就是判断一个给定的正整数是否素数。输入格式:输入在第一行给出一个正整数N(<=10),随后N行,每行给出一个小于231的需要判断的正整数。输出格式:对每个需要判断的正整数,如果它是素数,则在一行中输出“Yes”,否则输出“No”。输入样例:211111输出样例:YesNo通过这道题,我学会了判断素数的几种方法:方法一:bool judge(int n){ ...原创 2018-03-25 13:52:56 · 1006 阅读 · 0 评论 -
矩阵A乘以B - vector+矩阵乘法
题目描述给定两个矩阵A和B,要求你计算它们的乘积矩阵AB。需要注意的是,只有规模匹配的矩阵才可以相乘。即若A有Ra行、Ca列,B有Rb行、Cb列,则只有Ca与Rb相等时,两个矩阵才能相乘。输入输入先后给出两个矩阵A和B。对于每个矩阵,首先在一行中给出其行数R和列数C,随后R行,每行给出C个整数,以1个空格分隔,且行首尾没有多余的空格。输入保证两个矩阵的R和C都是正数,并且所有整数的绝对值不超过10...原创 2018-03-10 19:39:40 · 3108 阅读 · 0 评论 -
牛客训练赛23 - ABCD题
A - 托米的赌球 - 贪心链接:https://www.nowcoder.com/acm/contest/156/A来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K64bit IO Format: %lld题目描述紧张刺激的世界杯正在进行中(在托米的世界线里),欧洲人托米沉迷于赌球无法自拔。托米的口袋里有...原创 2018-07-28 15:25:48 · 389 阅读 · 0 评论 -
Prime Path - poj3126 - 区间筛素数+bfs求最短路
Prime PathTime Limit: 1000MS Memory Limit: 65536K Total Submissions: 27722 Accepted: 15137 DescriptionThe ministers of the cabinet were quite upset by the message from the Chie...原创 2018-08-09 22:13:18 · 381 阅读 · 0 评论 -
Find Integer - hdu6441 - 费马大定理+奇偶数列法则
题目来源:http://acm.hdu.edu.cn/showproblem.php?pid=6441思路:输入n,若n>2由费马大定理知无解。输出-1 -1,n=0时也无解,输出-1 -1 若n=1,则a+b=c,输出1 a+1 若n=2,则a^2+b^2=c^2,由奇偶数列法则得到:a=2n+1 (n=1,2,3…){ b= n^2+(n+1)^2-1 c...原创 2018-08-26 19:09:28 · 161 阅读 · 0 评论 -
直角三角形a^2+b^2=c^2整数解的定a公式直求法
直角三角形a^2+b^2=c^2整数解的定a公式直求法 庄 严 (中国辽阳 111000) 摘要:在直角三角形边长abc关系中,利用ab边条件求得第三边,这是人们的普遍做法。定a公式直求法的发现打破了人们的传统认识。利用定a公式直求法,在abc三边关系中,只要给定一个a值...转载 2018-08-26 19:00:52 · 6774 阅读 · 0 评论 -
牛客小白月赛6 - J洋灰三角形 - 等比数列、逆元
链接:https://www.nowcoder.com/acm/contest/136/J来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K64bit IO Format: %lld题目描述 洋灰是一种建筑材料,常用来筑桥搭建高层建筑,又称,水泥、混凝土。 WHZ有很多铸造成三角形的洋灰块,他想把这...原创 2018-08-26 15:45:55 · 221 阅读 · 0 评论 -
Dead Fraction - poj1930 - 循环小数转分数
Dead FractionTime Limit: 1000MS Memory Limit: 30000K Total Submissions: 3840 Accepted: 1286 DescriptionMike is frantically scrambling to finish his thesis at the last minute. H...原创 2018-08-07 17:25:07 · 164 阅读 · 0 评论 -
X-factor Chains - poj3421 - 分解质因数、多重集排列
X-factor ChainsTime Limit: 1000MS Memory Limit: 65536K Total Submissions: 8730 Accepted: 2775 DescriptionGiven a positive integer X, an X-factor chain of length m is a sequence...原创 2018-08-06 22:55:33 · 214 阅读 · 0 评论 -
牛客练习赛24 - AB
A - 石子阵列 - 排列组合链接:https://www.nowcoder.com/acm/contest/157/A来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 262144K,其他语言524288K64bit IO Format: %lld题目描述xb有m种石子,每种无限个,Ta想从这些石子中取出n个,并按顺序排列起来,为了好看,相邻的石...原创 2018-08-10 23:05:05 · 209 阅读 · 0 评论 -
阶乘之和
思路:可以先输出一下,20!有19位,满足题目要求,判断一个数n是否由阶乘加和得到,可以从第一个小于等于n的阶乘判断起,比如,n=25,25=4!+1!4!是第一个小于等于25的数,25-=4!,n=1,然后1!是第一个满足小于等于n的,n-=1!,此时n等于0,可以用几个数的阶乘表示。(我在想,判断时为什么总是从大到小判断,有没有可能25=3!+2!(打个比方),后来想想还是没可...原创 2018-08-06 11:08:27 · 376 阅读 · 0 评论 -
Semi-prime H-numbers - poj3292 - 筛数
Semi-prime H-numbersTime Limit: 1000MS Memory Limit: 65536K Total Submissions: 10785 Accepted: 4838 DescriptionThis problem is based on an exercise of David Hilbert, who pedago...原创 2018-08-09 22:31:01 · 151 阅读 · 0 评论 -
poj2991 - Crane - 线段树+计算几何
CraneTime Limit: 2000MS Memory Limit: 65536K Total Submissions: 10058 Accepted: 2641 Special Judge DescriptionACM has bought a new crane (crane -- jeřáb) . The crane co...原创 2019-08-20 13:16:09 · 180 阅读 · 0 评论