数论
文章平均质量分 52
algzjh
这个作者很懒,什么都没留下…
展开
-
莫比乌斯反演
莫比乌斯(Mobius)函数定义:μ(n)={1n=1(−1)kn=p1p2⋯pk0其他\mu(n)=\left\{\begin{array}{}1 & & n=1 \\(-1)^k & & n=p_1p_2\cdots p_k \\0 & & 其他\end{array}\righ...原创 2018-09-30 14:24:49 · 247 阅读 · 0 评论 -
中南大学2016年校队选拔赛第一场Problem A
Problem A: 小Z的黑白棋Time Limit: 1 Sec Memory Limit:128 MBSubmit: 18 Solved: 6[Submit][Status][Ask]Description小Z有一些黑白棋,他觉得黑白混杂在一起极具美感,所以他总喜欢将这些棋子排成一排序列S1,但是小Y就喜欢跟小Z作对,她会趁小Z不注意偷偷将小Z最右边的棋子拿走原创 2016-08-23 19:22:23 · 773 阅读 · 0 评论 -
CodeForces 27E
E. Number With The Given Amount Of Divisorstime limit per test2 secondsmemory limit per test256 megabytesinputstandard inputoutputstandard outputGiven the number n, find原创 2016-08-21 16:22:00 · 728 阅读 · 2 评论 -
求两个数的最大公约数的四种方法
关于两个数的最大公约数的求法探究 2016年9月12日求最大公约数是一个很基本的问题,在一些数论题中也有广泛的应用~一些公司的面试题也会涉及这方面的知识,所以有必要好好掌握~解法一:辗转相除法(1)递归int gcd(int x,int y)//辗转相原创 2016-09-16 21:26:51 · 1861 阅读 · 0 评论 -
CF837E-Vasya's Function
E. Vasya’s Functiontime limit per test1 second memory limit per test256 megabytes inputstandard input outputstandard output Vasya is studying number theory. He has denoted a function f(a, b) such t原创 2017-08-06 14:32:10 · 912 阅读 · 0 评论 -
2017"百度之星"程序设计大赛 - 初赛- 小C的倍数问题
小C的倍数问题Accepts: 1990 Submissions: 4931 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Problem Description根据小学数学的知识,我们知道一个正整数x是3的倍数的条件是x每一位加起来的和是3的倍数。反之,如果一个数每一位加起来原创 2017-08-13 00:54:10 · 714 阅读 · 0 评论 -
线性筛法求解积性函数(莫比乌斯函数)
莫比乌斯函数μ\muμ(n)=⎧⎩⎨1n=1(−1)kn=p1p2…pk0其余情况\mu(n)= \begin{cases}1\qquad \qquad n=1 \\(-1)^k \quad n=p_1p_2\dots p_k \\0 \qquad \qquad 其余情况 \end{cases}#include<iostream>#include<cstdio>#include<cs原创 2017-08-04 17:02:29 · 1688 阅读 · 0 评论 -
线性筛(同时得到欧拉函数和素数表)
#include<iostream>#include<cstring>#include<cstdio>#include<string>#include<cmath>using namespace std;const int INF=0x3f3f3f3f;const int MAXN=1e7+5;const int MOD=1e9+7;bool check[MAXN];int ph转载 2017-08-18 16:22:00 · 225 阅读 · 0 评论 -
HDU6033-Add More Zero
Add More ZeroTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others) Total Submission(s): 849 Accepted Submission(s): 585Problem Description There is a youngster know原创 2017-07-27 08:49:33 · 436 阅读 · 0 评论 -
递推总结
卡特兰数: f(n)=Cn2nn+1=(2n)!n!(n+1)!f(n)=\dfrac{C_{2n}^n}{n+1}=\dfrac{(2n)!}{n!(n+1)!} f(n+1)=4n−6nf(n)f(n+1)=\dfrac{4n-6}{n}f(n)原创 2017-07-30 11:22:39 · 512 阅读 · 0 评论 -
HDU6124-Euler theorem
Euler theoremTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others) Total Submission(s): 157 Accepted Submission(s): 132Problem Description HazelFan is given two pos原创 2017-08-15 20:59:31 · 406 阅读 · 0 评论 -
Bell数
将基数为nn的集合划分为任意个非空集的方法的数目称为BellBell数#include<iostream>#include<cstdio>#include<algorithm>#include<queue>#include<cstring>using namespace std;typedef long long LL;const int INF=0x3f3f3f3f;const i原创 2017-08-26 13:10:59 · 374 阅读 · 0 评论 -
CF983A
/*判断p/q在b进制下是否是无限小数如果1/q是在b进制下是有限小数,那么1/q能够写成a1/b+a2/b^2+...+an/b^n也就是q的所有素因子也是b的素因子*/#include<iostream>#include<cstdio>#include<vector>#include<algorithm>#include<...原创 2018-05-16 14:34:04 · 803 阅读 · 0 评论 -
HourRank 25-Maximum Palindromes
题目链接:Maximum Palindromes知识点:逆元、阶乘、排列、前缀和、乘法原理、预处理#include <cmath>#include <cstdio>#include <vector>#include <cstring>#include <iostream>#include <algorithm&...原创 2019-03-28 15:19:06 · 394 阅读 · 0 评论 -
有关逆元的求法
详见Modular multiplicative inverse转载 2018-01-03 15:31:00 · 582 阅读 · 0 评论 -
Lucas定理
#include <iostream>#include <vector>#include <algorithm>#include <cstdio>using namespace std;typedef long long LL;const int MAXN=1e4+5;const int MAX_P=1e4+5;const double eps=1e-8;int gcd(int a原创 2017-08-31 16:39:56 · 296 阅读 · 0 评论 -
阶乘取模预处理
int fact[MAX_P];//预处理n! mod p的表 O(p)//分解n!=ap^e,返回a mod p O(log_p n)int mod_fact(int n,int p,int& e){ e=0; if(n==0) return 1; //计算p的倍数的部分 int res=mod_fact(n/p,p,e); e+=n/p; //由于原创 2017-08-31 16:34:06 · 1261 阅读 · 0 评论 -
快速幂
LL Q_pow(LL x,LL n,LL MOD){ LL res=1; while(n>0) { if(n&1) res=res*x%MOD;//二进制最低位为1,则乘上x^(2^i) x=x*x%MOD; n>>=1; } return res;}原创 2017-08-31 10:06:05 · 234 阅读 · 0 评论 -
素数筛选
#include<bits/stdc++.h>using namespace std;const int MAXN=1000010;bool prime[MAXN];//true表示是素数//判断<MAXN的数是不是素数void init(){ memset(prime,true,sizeof(prime)); prime[0]=prime[1]=false; fo原创 2017-03-15 12:01:52 · 407 阅读 · 0 评论 -
素性测试
//假设输入的都是正数//素性测试O(sqrt(n))bool is_prime(int n){ for(int i=2;i*i<=n;i++) { if(n%i==0) return false; } return n!=1;//1是例外}//约数枚举vector<int> divisor(int n){ vector<int>原创 2017-08-31 09:35:12 · 392 阅读 · 0 评论 -
CF850B-Arpa and a list of numbers
B. Arpa and a list of numberstime limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard output Arpa has found a list containing n numbers. He calls a list ba原创 2017-09-05 21:54:32 · 477 阅读 · 0 评论 -
HDU6069-Counting Divisors
Counting DivisorsTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others) Total Submission(s): 442 Accepted Submission(s): 132Problem Description In mathematics, the原创 2017-08-03 19:12:30 · 359 阅读 · 0 评论 -
快速指数算法
一般地,求ama^m可如下进行(其中,a、m是正整数): 将m表示为二进制形式bk2k+bk−12k−1+⋯+b12+b0b_k2^k+b_{k-1}2^{k-1}+\cdots+b_12+b_0 因此 am=(⋯(((abk)2abk−1)2abk−2)2⋯ab1)2ab0a^m=(\cdots(((a^{b_k})^2a^{b_{k-1}})^2a^{b_{k-2}})^2\cdots a原创 2017-03-15 09:45:27 · 4685 阅读 · 0 评论 -
单变元模线性方程
复杂度O(logn) 使用范例:POJ2115#include<iostream>#include<vector>using namespace std;typedef long long LL;//ax+by=gcd(a,b)LL extend_gcd(LL a,LL b,LL &x,LL &y){ if(b==0) { x=1;y=0;原创 2017-04-03 16:18:51 · 801 阅读 · 0 评论 -
POJ2115-C Looooops
C Looooops Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 25950 Accepted: 7352 DescriptionA Compiler Mystery: We are given a C-language style for loop of type for (variable原创 2017-04-03 16:15:20 · 543 阅读 · 0 评论 -
欧几里得算法与扩展欧几里得算法与逆元
欧几里得算法与扩展欧几里得算法预备知识1. [a,b]=|ab|(a,b) [a,b] = \dfrac{|ab|}{(a,b)} 实际上,计算最好使用lcm(a,b)=a/gcd(a,b)∗blcm(a,b)=a/gcd(a,b)*b 2. 设 a,b,c∈Z\quad\ a,b,c\in\mathbb Z (1)∀x∈Z,(ax+b,a)=(a,b)\quad\forall x\in原创 2017-01-31 21:56:45 · 627 阅读 · 0 评论 -
2017湖南多校第五场-A(1043): Biorhythms
A(1043): Biorhythms Submit Page Summary Time Limit: 1 Sec Memory Limit: 256 Mb Submitted: 50 Solved: 27 Description Some people believe that there are three cycles in a person’s原创 2017-04-03 13:48:27 · 512 阅读 · 0 评论 -
质因数分解
#include<bits/stdc++.h>using namespace std;typedef long long LL;LL factor[100];int fatcnt;void getfactors(LL x){ memset(factor,0,sizeof(factor)); fatcnt=0; for(int i=2;i<=sqrt(x);i++)原创 2017-03-15 17:36:25 · 1294 阅读 · 0 评论 -
中国剩余定理
设m1,m2,⋯,mk∈Z+m_1,m_2,\cdots,m_k\in Z^+,两两互素,k≥2k\ge 2,则对任意k个整数a1,a2,⋯,aka_1,a_2,\cdots,a_k,同余方程组: ⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪x≡a1(mod m1)x≡a2(mod m2)⋅⋅⋅x≡ak(mod mk)\begin{equation}\begin{cases}x\equiv a_1 (m原创 2017-03-11 17:31:37 · 642 阅读 · 0 评论 -
剩余类与剩余系
1、模m的两两不同的剩余类有m个2、同余方程的解就是一个或几个剩余类3、一个完全剩余系(完系)中,与m互质的整数的个数叫作欧拉函数,记作φ(m)\varphi(m)。 这φ(m)\varphi(m)个与m互质的整数组成模m的一个简化剩余系(缩系)。4、φ(m)\varphi(m)是偶数。5、当p为素数时,φ(p)\varphi(p)=p-1;若q也为素数,m=pq时,φ(m)\varphi(m)=原创 2017-03-11 15:27:04 · 7097 阅读 · 0 评论 -
一次同余方程
1、ax≡b(modm)ax\equiv b(mod m) 一次同余方程有解的充要条件是(a,m)|b。若有解,解的个数为d=(a,m),它们是 x≡x0+mt/d(modm),t=0,1,⋯,d−1x\equiv x_0+mt/d(mod m), t=0,1,\cdots, d-1,其中x0x_0为方程的一个特解。2、逆元若m,a∈Z,m≠0,(a,m)=1m,a\in \mathbf{Z}原创 2017-03-11 15:01:21 · 2206 阅读 · 0 评论 -
HDU2669 Romantic
RomanticTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5625 Accepted Submission(s): 2383Problem DescriptionThe Sky is Sprite原创 2017-01-31 23:42:40 · 293 阅读 · 0 评论 -
中国剩余定理
任务:求出方程组x≡ai (mod mi) (0<=i<n) 的解 x其中m_0, m_1,m_2, ... , m_{n-1}两两互质使用范例:[POJ1006](http://poj.org/problem?id=1006)#include#includeusing namespace std;typedef long long LL;LL a[5],m原创 2017-04-03 16:45:27 · 367 阅读 · 0 评论 -
HDU1286找新朋友(求欧拉函数)
找新朋友Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 12787 Accepted Submission(s): 6785Problem Description 新年快到了,“猪头帮协会”准备搞一个聚会,已经知道现有会员N人,把原创 2017-05-17 15:13:33 · 662 阅读 · 0 评论 -
HDU6063-RXD and math
RXD and mathTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others) Total Submission(s): 448 Accepted Submission(s): 230Problem Description RXD is a good mathematicia原创 2017-08-02 10:03:48 · 604 阅读 · 0 评论 -
HDU6053-TrickGCD
#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>using namespace std;typedef long long LL;const int INF=0x3f3f3f3f;const LL MOD=1e9+7;const int MAXN=1e5+7;原创 2017-07-31 11:41:11 · 336 阅读 · 0 评论 -
欧拉定理与费马小定理
1、欧拉(Euler)定理 设m∈Z+,m≥2m\in Z^+,m\ge2,则∀a∈Z,(m,a)=1\forall a \in Z,(m,a)=1,有 aφ(m)≡1(modm)a^{\varphi(m)}\equiv 1(mod m)2、费马(Fermat)小定理 设p是素数,则∀a∈Z\forall a\in Z有 ap≡a(modp)a^p\equiv a(mod p)原创 2017-03-11 17:02:40 · 431 阅读 · 0 评论 -
HackerRank-Mehta and the Typical Supermarket
Mehta and the Typical SupermarketMehta is a very rich guy. He has types of coins, and each type of coin is available in an unlimited supply.So Mehta goes to a supermarket to buy monthly groceries. The转载 2017-07-18 11:44:40 · 423 阅读 · 0 评论 -
AtCoder Grand Contest 017-A - Biscuits
A - BiscuitsTime limit : 2sec / Memory limit : 256MBScore : 200 pointsProblem Statement There are N bags of biscuits. The i-th bag contains Ai biscuits.Takaki will select some of these bags and eat al原创 2017-07-09 22:21:21 · 423 阅读 · 0 评论 -
AtCoder Regular Contest 077-D - 11
D - 11Time limit : 2sec / Memory limit : 256MBScore : 600 pointsProblem Statement You are given an integer sequence of length n+1, a1,a2,…,an+1, which consists of the n integers 1,…,n. It is known tha原创 2017-07-02 01:51:19 · 575 阅读 · 0 评论