数论
★漂浮ing☆
奋斗ing
展开
-
矩阵乘法
问题描述 输入两个矩阵,分别是m*s,s*n大小。输出两个矩阵相乘的结果。输入格式 第一行,空格隔开的三个正整数m,s,n(均不超过200)。 接下来m行,每行s个空格隔开的整数,表示矩阵A(i,j)。 接下来s行,每行n个空格隔开的整数,表示矩阵B(i,j)。输出格式 m行,每行n个空格隔开的整数,输出相乘後的矩阵C(i,j)的值。样例输入2 3 21 0 -11 1 -30 31 23 1样例输出-3 2-8 2提示矩阵C应该是m行n列,其原创 2020-07-24 18:00:43 · 290 阅读 · 0 评论 -
最大最小公倍数
点这里【思路】三个数的最大最小公倍数,就是找三个数要两两互质,,如果这个数小于2,直接就是它本身了;如果为奇数,就是,,如果为偶数,还要看是否能被3整除,如果不能,,n和(n-2)都是偶数,往前推一个n-3,就是,,如果能被3整除,,6就是一个很好例子,3和6有公约数3,不是1,,这时就选择不开longlong错误。。#include <bits/stdc++.h>using namespace std;typedef long long ll;int main(){原创 2020-07-24 17:28:24 · 324 阅读 · 0 评论 -
矩阵快速幂(好用的板子)
试题 基础练习 矩阵乘法提交此题资源限制时间限制:1.0s 内存限制:512.0MB问题描述 给定一个N阶矩阵A,输出A的M次幂(M是非负整数) 例如: A = 1 2 3 4 A的2次幂 7 10 15 22输入格式 第一行是一个正整数N、M(1<=N<=30, 0<=M<=5),表示矩阵A的阶数和要求的...原创 2020-04-16 12:06:34 · 155 阅读 · 0 评论 -
斐波那契数列找循环节
首先贴一个大佬的广义方法,,数学真的diao。。广义斐波那契数列的循环节暴力求法,一般情况下这样求long long 范围内的也就够了。#include<bits/stdc++.h>#define ll long longusing namespace std;const int maxn=1e5+10;// 自调ll f[maxn];int ...原创 2020-04-14 11:24:11 · 423 阅读 · 0 评论 -
筛选区间内素数的个数
问题描述 给定区间[L, R] , 请计算区间中素数的个数。输入格式 两个数L和R。输出格式 一行,区间中素数的个数。样例输入2 11样例输出5数据规模和约定 2 <= L <= R <= 2147483647 R-L <= 1000000核心代码:for(ll j = max(2LL, (a+i-1)/i) *...原创 2020-04-10 15:28:45 · 858 阅读 · 1 评论 -
数字游戏
问题描述 栋栋正在和同学们玩一个数字游戏。 游戏的规则是这样的:栋栋和同学们一共n个人围坐在一圈。栋栋首先说出数字1。接下来,坐在栋栋左手边的同学要说下一个数字2。再下面的一个同学要从上一个同学说的数字往下数两个数说出来,也就是说4。下一个同学要往下数三个数,说7。依次类推。 为了使数字不至于太大,栋栋和同学们约定,当在心中数到 k-1 时,下一个数字从0开始数。例如,当k=13时,...原创 2020-04-04 17:02:48 · 122 阅读 · 0 评论 -
大数阶乘
问题描述 输入一个正整数n,输出n!的值。 其中n!=1*2*3*…*n。算法描述 n!可能很大,而计算机能表示的整数范围有限,需要使用高精度计算的方法。使用一个数组A来表示一个大整数a,A[0]表示a的个位,A[1]表示a的十位,依次类推。 将a乘以一个整数k变为将数组A的每一个元素都乘以k,请注意处理相应的进位。 首先将a设为1,然后乘2,乘3,当乘到n时,即得到了...原创 2020-03-19 12:37:09 · 107 阅读 · 0 评论 -
codeforces1312 D(组合数)
D. Count the Arraystime limit per test2 secondsmemory limit per test512 megabytesinputstandard inputoutputstandard outputYour task is to calculate the number of arrays such that:ea...原创 2020-03-10 23:36:14 · 272 阅读 · 0 评论 -
2038: [2009国家集训队]小Z的袜子(hose)(莫队算法+分块)
Description作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L尽管小Z并不在意两只袜子是不是完整的一双,甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬。你的任务便是告诉小Z,他有多...原创 2019-08-21 16:17:25 · 128 阅读 · 0 评论 -
AtCoder Beginner Contest 156D - Bouquet(Lacus)
Time Limit: 2 sec / Memory Limit: 1024 MBScore :400400pointsProblem StatementAkari hasnnkinds of flowers, one of each kind.She is going to choose one or more of these flowers to make a bou...原创 2020-02-24 20:24:19 · 473 阅读 · 0 评论 -
组合数(卢卡斯定理)
组合数从n个不同元素中,任取m(m≤n)个元素并成一组,叫做从n个不同元素中取出m个元素的一个组合;从n个不同元素中取出m(m≤n)个元素的所有组合的个数,叫做从n个不同元素中取出m个元素的组合数。公式:C(n,m)=n!/((n-m)!*m!)(m≤n)性质1:C(n,m)= C(n,n-m)性质2:C(n,m)=C(n-1,m-1)+C(n-1,m)定义组合是数学...原创 2020-02-24 20:08:59 · 4738 阅读 · 0 评论 -
矩阵快速幂
矩阵快速幂,看着听高深,其实就是字面意思,矩阵的快速幂,再说清楚点就是矩阵的幂运算,很好理解,但这涉及到矩阵的乘法,先来矩阵的乘法:上模板:const int N=2;int tmp[N][N];void multi(int a[][N],int b[][N],int n){ memset(tmp,0,sizeof tmp); for(int i=0;i<n;i+...原创 2020-02-22 16:40:53 · 148 阅读 · 0 评论 -
完美平方数
We say that x is a perfect square if, for some integer b, x = b 2 . Similarly, x is a perfect cube if, for some integer b, x = b 3 . More generally, x is a perfect pth power if, for some integer b, x ...原创 2020-02-18 17:19:17 · 509 阅读 · 0 评论 -
区间内约数最多的书的约数个数
Mathematicians love all sorts of odd properties of numbers. For instance, they consider 945 to be an interesting number, since it is the first odd number for which the sum of its divisors is larger th...原创 2020-02-15 13:54:44 · 271 阅读 · 0 评论 -
*佩尔方程*
A computer programmer lives in a street with houses numbered consecutively (from 1) down one side of the street. Every evening she walks her dog by leaving her house and randomly turning left or right...原创 2020-02-11 19:09:29 · 168 阅读 · 0 评论 -
最大公约数(打表)
Well, here is another math class task. In mathematics, GCD is the greatest common divisor, and it's an easy task to calculate the GCD between two positive integers.A common divisor for two positive ...原创 2020-02-01 23:43:58 · 264 阅读 · 0 评论 -
扩展欧拉降幂
求 b个a的次幂解法:对于,我们可以直接使用欧拉降幂去求,但这个式子是a次幂,又该怎么办?那么他的下一步一定还是递归的形式已知往下,直到q == 1 时,向上返回结果,所以只需在欧拉降幂的基础上加一个dfs模板:#include<bits/stdc++.h>using namespace std;typedef long long LL;...原创 2020-02-01 16:09:42 · 86 阅读 · 0 评论 -
欧拉降幂模板
欧拉降幂公式:第一个是a和p互质的时候,后两个是a和p不互质的时候。来一道裸题,稍微有点变化,两层的。一般情况下数据比较大的时候,a和b都定义为字符数组型。Given A,B,C, You should quickly calculate the result of A^B mod C. (1<=A,C<=1000000000,1<=B<...原创 2020-02-01 16:02:18 · 228 阅读 · 0 评论 -
扩展欧几里得算法
发现很多以前学过的都没有整理。。扩展欧几里得算法用于求解二元一次不定式的通解,因而被广泛使用。设定一常数为gcd=(A,B)令AX0+BY0=gcd(A,B)(1),如果完成欧几里得的算法后,发现:当B=0时,gcd(A,B)=A。要满足(1)式,显然x=1,y=0。此时,把递推过程完成,得到了一组结果。但是要求得出AX+BY=gcd(A,B)这种式子的满足的根,这样还不够,因为此时...原创 2020-01-31 15:07:12 · 156 阅读 · 0 评论 -
小规律(阶乘)
Mr. Santa asks all the great programmers of the world to solve a trivial problem. He gives them an integermand asks for the number of positive integersn, such that the factorial ofnends with exac...原创 2020-01-29 16:23:27 · 355 阅读 · 0 评论 -
数学思维小结论
It is so boring in the summer holiday, isn't it? So Alice and Bob have invented a new game to play. The rules are as follows. First, they get a set ofndistinct integers. And then they take turns to ...原创 2020-01-27 23:25:06 · 110 阅读 · 0 评论 -
费马大定理+勾股数奇偶规律
people in USSS love math very much, and there is a famous math problem .give you two integersnn,aa,you are required to find22integersbb,ccsuch thatanan+bn=cnbn=cn.Inputone line contains on...原创 2020-01-21 23:28:58 · 646 阅读 · 0 评论 -
积性函数
积性函数定义:对于正整数n的一个算术函数 f(n),若f(1)=1,且当a,b互质时f(ab)=f(a)f(b),在数论上就称它为积性函数。若对于某积性函数 f(n) ,就算a, b不互质,也有f(ab)=f(a)f(b),则称它为完全积性的常见积性函数:φ(n) -欧拉函数,计算与n互质的正整数之数目 μ(n) -莫比乌斯函数,关于非平方数的质因子数目gcd(n,k) -最大公因...原创 2020-01-14 15:58:44 · 231 阅读 · 0 评论 -
容斥原理
在组合数学中,容斥是常常被用到的,用容斥求解一些带有条件的组合数。容斥原理:具有性质A和性质B的元素个数等同于具有性质A的个数和具有性质B的个数的和再减去同时具有性质A和性质B的元素的个数。数学公式表示为 |A∪B|=|A|+|B|-|A∩B|。同样以此类推对于三个性质来说其数学公式为|A∪B∪C|=|A|+|B|+|C|-|A∩B|-|A∩C|-|B∩C|+|A∩B∩C|为什么要...原创 2020-01-14 15:19:02 · 722 阅读 · 0 评论 -
欧拉函数
The Euler function phi is an important kind of function in number theory, (n) represents the amount of the numbers which are smaller than n and coprime to n, and this function has a lot of beautiful c...原创 2020-01-13 22:21:40 · 188 阅读 · 0 评论 -
欧拉函数模板
欧拉函数的定义:在数论中,对于正整数N,少于或等于N ([1,N]),且与N互质的正整数(包括1)的个数,记作φ(n)。φ函数的值:φ(x)=x(1-1/p(1))(1-1/p(2))(1-1/p(3))(1-1/p(4))……(1-1/p(n)) 其中p(1),p(2)…p(n)为x的所有质因数;x是正整数; φ(1)=1(唯一和1互质的数,且小于等于1)。注意:每种质因数只有一个。例...原创 2020-01-13 22:19:35 · 75 阅读 · 0 评论 -
原根
原根的定义:原根,是一个数学符号。设m是正整数,a是整数,若a模m的阶等于φ(m),则称a为模m的一个原根。1. 原根的定义设m是正整数,a是整数,若a模m的阶等于φ(m),则称a为模m的一个原根。(其中φ(m)表示m的欧拉函数)假设一个数g对于P来说是原根,那么g^i mod P的结果两两不同,且有 1<g<P, 0<i<P,那么g可以称为是P的一个原根,归...原创 2019-12-22 09:08:53 · 1115 阅读 · 0 评论 -
打表找规律(筛素数)
This time I need you to calculate the f(n) . (3<=n<=1000000)f(n)= Gcd(3)+Gcd(4)+…+Gcd(i)+…+Gcd(n).Gcd(n)=gcd(C[n][1],C[n][2],……,C[n][n-1])C[n][k] means the number of way to choose k things...原创 2019-12-21 17:45:45 · 93 阅读 · 0 评论 -
第K大的因子
In mathematics, the greatest common divisor (gcd), also known as the greatest common factor (gcf), highest common factor (hcf), or greatest common measure (gcm), of two or more integers (when at least...原创 2019-12-20 22:19:28 · 124 阅读 · 0 评论 -
快速幂+gcd变形定理
To think of a beautiful problem description is so hard for me that let's just drop them off. :)Given four integers a,m,n,k,and S = gcd(a^m-1,a^n-1)%k,calculate the S.InputThe first line conta...原创 2019-12-19 16:32:22 · 131 阅读 · 0 评论 -
数论中的互质数的最大不能组合数
题意:给定A和B,A和B互质,求最大不能组合数,和不能组合数的个数。基础知识:Gcd(A, B) = 1 → Lcm(A, B) = AB剩余类,把所有整数划分成m个等价类,每个等价类由相互同余的整数组成任何数分成m个剩余类,分别为 mk,mk+1,mk+2,……,mk+(m-1)分别记为{0(mod m)},{1(mod m)}……而n的倍数肯定分布在这m个剩余类中因为Gcd(m,n)=1,所以...原创 2019-12-18 22:32:42 · 590 阅读 · 0 评论 -
逆元
对于正整数a和m,如果有,那么把这个同余方程中的最小正整数解叫做a模m的逆元。逆元一般用扩展欧几里得算法来求得,如果为素数,那么还可以根据费马小定理得到逆元为。推导过程如下求现在来看一个逆元最常见问题,求如下表达式的值(已知)当然这个经典的问题有很多方法,最常见的就是扩...原创 2019-12-17 19:31:27 · 449 阅读 · 0 评论 -
变同余求最小公倍数
我知道部分同学最近在看中国剩余定理,就这个定理本身,还是比较简单的:假设m1,m2,…,mk两两互素,则下面同余方程组:x≡a1(mod m1)x≡a2(mod m2)…x≡ak(mod mk)在0<=<m1m2…mk内有唯一解。记Mi=M/mi(1<=i<=k),因为(Mi,mi)=1,故有二个整数pi,qi满足Mipi+miqi=1,如果记e...原创 2019-12-17 11:56:03 · 437 阅读 · 0 评论 -
求素数的位置
Everybody knows any number can be combined by the prime number.Now, your task is telling me what position of the largest prime factor.The position of prime 2 is 1, prime 3 is 2, and prime 5 is 3, ...原创 2019-12-16 21:38:29 · 116 阅读 · 0 评论 -
找一个数的因子个数,因子和
1.所有因子个数如果一个数是因数,就不断除这个数,保存这个因子次方的数 temp++;运用所有因子个数计算公式(见上图),保存因子个数的 ans不断乘( temp+1 )。注意 :当最后,在 x 不断除因数得到的值有两种情况:x == 1,这说明 x 没有其他因子了。x != 1, 这时 x 为其一个素数因子(且这个因子大于 根号x ),所以最后再乘(1+1)。代码...原创 2019-12-16 21:36:55 · 384 阅读 · 0 评论 -
佩尔方程与连分数
【概述】连分数是一种特殊的繁分数,其形式为:,通常记为:,其中和 称为连分数多项式,对于任意的 a 均为一次式,它们的比值称为第 n 个渐进值渐进分数。佩尔(Pell)方程是一种不定二次方程,其与连分数,二次型,代数论等有着重要的联系。其形式为:,其中 d不为非平方数【佩尔方程迭代公式】定义:设 p、q 为整数,且满足,则称给出该方程的解推论:设给出方程...原创 2019-12-10 17:51:07 · 397 阅读 · 0 评论 -
扩展欧几里得求逆元
什么叫乘法逆元?这里,我们称 x 是 a 关于 m 的乘法逆元这怎么求?可以等价于这样的表达式: a*x + m*y = 1看出什么来了吗?没错,当gcd(a , m) != 1 的时候是没有解的这也是 a*x + b*y = c 有解的充要条件: c % gcd(a , b) == 0接着乘法逆元讲,一般,我们能够找到无数组解满足条件,...原创 2019-12-07 10:28:27 · 420 阅读 · 0 评论 -
快速幂/快速乘/矩阵快速幂模板
一、快速乘ll mul(ll m,ll k,ll p){ ll ans=0; while(k) { if(k&1) ans=(ans+m)%p; m=m*2%p; k>>=1; } return ans;}二、快速幂ll qmi(ll m,ll k,ll ...原创 2019-12-05 16:53:56 · 111 阅读 · 0 评论 -
求区间内素数个数
题目背景区间查询素数个数题目描述输入输出格式输入格式:一行两个整数 询问次数n,范围m接下来n行,每行两个整数 l,r 表示区间输出格式:对于每次询问输出个数 t,如l或r∉[1,m]输出 Crossing the line输入输出样例输入样例#1:复制2 51 32 6输出样例#1:复制2Crossing the ...原创 2019-12-04 20:58:38 · 1879 阅读 · 0 评论 -
筛素数
筛法一#include<iostream>#include<cstring>#include<cstdio>#define mem(a,x) memset(a,x,sizeof(a))#define inf (1<<29)using namespace std;typedef long long ll;const int N = ...原创 2019-12-04 20:43:33 · 66 阅读 · 0 评论