
数 论
文章平均质量分 81
数 论
「已注销」
这个作者很懒,什么都没留下…
展开
-
Wannafly挑战赛6 A(二分)B(dfs)C(公式推导)D(思维)E(技巧)
A题:题目描述多次查询[l,r]范围内的完全平方数个数 定义整数x为完全平方数当且仅当可以找到整数y使得y*y=x 输入描述:第一行一个数n表示查询次数之后n行每行两个数l,r输出描述:对于每个查询,输出一个数表示答案示例1输入51 31 42 44 41 1000000000输出121原创 2017-12-23 14:32:00 · 718 阅读 · 0 评论 -
hdu 6069 统计区间约数的个数 2017 Multi-University Training Contest - Team 4
Counting DivisorsTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others)Total Submission(s): 1720 Accepted Submission(s): 629Problem DescriptionIn mathema原创 2017-08-04 13:21:28 · 421 阅读 · 0 评论 -
hdu 6044 组合数+分治+模拟元 2017 Multi-University Training Contest - Team 1
Limited PermutationTime Limit: 4000/2000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others)Total Submission(s): 1240 Accepted Submission(s): 313Problem DescriptionAs to a p原创 2017-07-27 09:55:00 · 508 阅读 · 0 评论 -
hdu 6053 莫比乌斯反演函数的利用 2017 Multi-University Training Contest - Team 2
TrickGCDTime Limit: 5000/2500 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Total Submission(s): 622 Accepted Submission(s): 236Problem DescriptionYou are given an arra原创 2017-07-27 22:20:02 · 482 阅读 · 0 评论 -
hdu 2256 公式推导,小数取模问题
Problem of PrecisionTime Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 1547 Accepted Submission(s): 940Problem Description I转载 2017-08-01 11:34:04 · 2660 阅读 · 0 评论 -
hdu 2604 矩阵快速幂
QueuingTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5951 Accepted Submission(s): 2594Problem DescriptionQueues and Priority Q原创 2017-08-01 10:31:29 · 334 阅读 · 0 评论 -
最容易理解的莫比乌斯反演
对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数1≤n≤50000,1≤a≤b≤50000,1≤c≤d≤50000,1≤k≤50000Sample Input 22 5 1 5 11 5 1 5 2Sample Output 143原创 2017-07-28 16:39:56 · 734 阅读 · 0 评论 -
csuoj X - Rectangles (公式推导或者二分)
如果限定矩形的边长必须为整数,且周长为定值L,那么面积在[A, B]范围内不同的矩形一共有多少个呢?在这个问题中,当且仅当两个矩形面积不同时,视作是两个不同的矩形。Input 输入数据的第一行包含一个整数T (1 对于每组测试数据,包含三个整数L (1 Output 对于每组测试数据,用一行输出一个整数,代表上述问题的答案。Sample Inpu原创 2017-06-29 16:36:28 · 412 阅读 · 0 评论 -
FFT 模板 大数相乘
大数乘法基准时间限制:2 秒 空间限制:131072 KB 分值: 80 难度:5级算法题 收藏 关注给出2个大整数A,B,计算A*B的结果。Input第1行:大数A第2行:大数B(A,B的长度 = 0)Output输出A * BInput示例123456234567Output示例2895870原创 2017-06-06 05:26:45 · 976 阅读 · 0 评论 -
Sheldon Numbers 暴力枚举
题意:求在区间内Sheldon数字的个数题解:枚举n m,然后判断是否可以组成相应的位数然后再判断是否这个数字是在这个区间内#include#include#includeusing namespace std;#define LL long longint a[100];unsigned long long x,y;bool deal(int le原创 2017-05-19 16:57:34 · 442 阅读 · 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互质的数,且小于等于原创 2017-05-17 08:50:35 · 329 阅读 · 0 评论 -
51nod 1040 欧拉定理
题解:我们很容易可以想到去统计n的因子在1~n中出现的次数假设此因子为 i ,那么就是求满足gcd(num,n)== i 的num有多少个,最后乘以 i化为 gcd(num / i ,n / i)==1,那么就是求与 i 互质的个数#include#include#includeusing namespace std;int euler(int原创 2017-05-17 20:01:11 · 353 阅读 · 0 评论 -
费马小定理
费马小定理(Fermat Theory)是数论中的一个重要定理,其内容为:假如p是质数,且gcd(a,p)=1,那么 a(p-1)≡1(mod p)。即:假如a是整数,p是质数,且a,p互质(即两者只有一个公约数1),那么a的(p-1)次方除以p的余数恒等于1。引理1. 若a,b,c为任意3个整数,m为正整数,且(m,c)=1,则当ac≡bc(modm)时,有a≡b(原创 2017-05-17 10:25:13 · 937 阅读 · 0 评论 -
51nod 1179 求数组中一对数的最大gcd
1179 最大的最大公约数题目来源: SGU 基准时间限制:1 秒 空间限制:131072 KB 分值: 40 给出N个正整数,找出N个数两两之间最大公约数的最大值。例如:N = 4,4个数为:9 15 25 16,两两之间最大公约数的最大值是15同25的最大公约数5。Input第1行:一个数N,表示输入正整数的数量。(2 <= N <= 50000)原创 2017-05-17 08:31:29 · 638 阅读 · 0 评论 -
扩展欧几里得
扩展欧几里德算法 谁是欧几里德?自己百度去 先介绍什么叫做欧几里德算法 有两个数 a b,现在,我们要求 a b 的最大公约数,怎么求?枚举他们的因子?不现实,当 a b 很大的时候,枚举显得那么的naïve ,那怎么做? 欧几里德有个十分又用的定理: gcd(a, b) = gcd(b , a%b) ,这样,我们就可以在几乎是 log 的时间复杂度里求解出转载 2017-05-16 21:30:33 · 317 阅读 · 0 评论 -
poj 3904 莫比乌斯函数灵活运用
Stancu likes space travels but he is a poor software developer and will never be able to buy his own spacecraft. That is why he is preparing to steal the spacecraft of Petru. There is only one problem原创 2017-08-19 09:45:28 · 529 阅读 · 0 评论 -
hdu 4609 FFT 2013年多校
3-idiotsTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5646 Accepted Submission(s): 1960Problem DescriptionKing OMeGa catched t原创 2017-08-05 11:08:08 · 348 阅读 · 0 评论 -
牛客练习赛 8 (1~n约数和)(前缀和好题)(加边得连通图)(容斥原理)
A题目描述给个n,求1到n的所有数的约数个数的和~ 输入描述:第一行一个正整数n输出描述:输出一个整数,表示答案示例1输入3输出5说明样例解释:1有1个约数12有2个约数1,23有2个约数1,3备注:n <= 100000000#incl原创 2017-12-16 13:53:03 · 986 阅读 · 0 评论 -
hdu 5446 Lucas定理,中国剩余定理,处理爆long long 的乘法取模运算
Unknown TreasureTime Limit: 1500/1000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others)Total Submission(s): 3056 Accepted Submission(s): 1128Problem DescriptionOn the way原创 2017-12-01 18:59:47 · 785 阅读 · 0 评论 -
反素数入门 zoj 2562 and codefore
学习资料来源传送们反素数的定义:对于任何正整数,其约数个数记为,例如,如果某个正整数满足:对任意的正整 数,都有,那么称为反素数定义理解:素数的约数只有两个 1 和 本身,而反素数的约数是尽可能多,比这个数小到正数的约数都要多在ACM竞赛中,最常见的问题如下:(1)给定一个数,求一个最小的正整数,使得的约数个数为(2)求转载 2017-11-20 16:52:18 · 318 阅读 · 0 评论 -
bzoj 3288 欧拉函数线性筛法
3288: Mato矩阵Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 348 Solved: 269[Submit][Status][Discuss]DescriptionMato同学最近正在研究一种矩阵,这种矩阵有n行n列第i行第j列的数为gcd(i,j)。例如n=5时,矩阵如下:1 1 1 1 11原创 2017-11-21 21:16:54 · 406 阅读 · 0 评论 -
bzoj 1013 高斯消元
1013: [JSOI2008]球形空间产生器sphereTime Limit: 1 Sec Memory Limit: 162 MBSubmit: 6149 Solved: 3198[Submit][Status][Discuss]Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在,你被困在了这个n维球体中,你只知道球面上n原创 2017-11-21 13:05:06 · 343 阅读 · 0 评论 -
bzoj 1407 扩展欧几里得
1407: [Noi2002]SavageTime Limit: 5 Sec Memory Limit: 64 MBSubmit: 2336 Solved: 1045[Submit][Status][Discuss]DescriptionInput第1行为一个整数N(1第2行到第N+1每行为三个整数Ci, Pi, Li表示每个野人所住的初始洞穴编原创 2017-11-21 10:54:05 · 342 阅读 · 0 评论 -
uva 12716 约数关系的打表优化方法 lrj-P318
题意:给出一个区间,问这个区间内有多少个数对满足 a^b==gcd(a,b) 题解:打表发现 gcd ==a-b ( a > b)然后就打表就行了,打表的时候有优化,枚举 c ,然后枚举 c 的倍数作为 a然后计算 b ,再判断即可,最后通过递推得到最终答案#include#define MAXN 30000005int ans原创 2017-09-19 17:39:15 · 549 阅读 · 0 评论 -
uva 10375 质数筛选 质数分解 lrj-P316
题意:求解C(p,q)/ C(r,s),保留五位小数题解:刘汝佳老师的代码真的很好,思路很清晰,好好品味,好好学习#include#include#include#includeusing namespace std;#define MAXN 10010#define LL long longconst LL mod=1000000007原创 2017-09-18 17:19:47 · 518 阅读 · 0 评论 -
uva 12169 ex_gcd lrj—P316
题意:给出一个递推式,通过这个递推式可以得到一系列数,但是题目只给出位置为奇数的数字(从1开始),求偶数位置的数题解:枚举+扩展欧几里得通过学习,发现刘汝佳的这本书简直是神书,以前在网上找了很久的资料才看懂的算法他居然三言两语就解决了我疑惑的地方,以后就会用这个算法了扩展欧几里得算法是解决 求解ax+by=c 这个函数的整数解 x y当 c 是gc原创 2017-09-18 16:51:52 · 627 阅读 · 0 评论 -
uva 11582 lrj-P316 斐波那契大数
给出无符号64位的a b两个数,要求斐波那契数列中第a的b次方个数是多少?(此数对n取模)取模运算a^b%n=(a%n)^b,运用快速幂搞一手就好了然后这里需要注意的是无符号输入的格式问题,见代码,以及n为1的情况,深入思考一下就好了#include#define LL unsigned long longLL mod;int f[100000+5];LL原创 2017-09-18 12:12:08 · 584 阅读 · 0 评论 -
spoj 7001 zoj 3435 莫比乌斯反演
Consider a N*N*N lattice. One corner is at (0,0,0) and the opposite one is at (N,N,N). How many lattice points are visible from corner at (0,0,0) ? A point X is visible from point Y iff no other latti原创 2017-08-07 11:43:47 · 335 阅读 · 0 评论 -
spoj Primes in GCD Table 莫比乌斯反演
Johnny has created a table which encodes the results of some operation -- a function of two arguments. But instead of a boring multiplication table of the sort you learn by heart at prep-school, he ha原创 2017-08-07 10:24:12 · 348 阅读 · 0 评论 -
xtu oj 1236 小数转分数
题目传送门这里可以找到规律就是分子等于分母乘以n的四舍五入四舍五入的计算方法见下面注意题目的数据里面包括分子为0的情况,也就是说n可能为0.00000#includedouble fbs(double x){return x<0?-x:x;}int gcd(int a,int b){ while(b) { int t原创 2017-04-17 15:26:48 · 712 阅读 · 0 评论 -
poj 3685 二分套二分 叉姐的魔法训练
MatrixTime Limit: 6000MS Memory Limit: 65536KTotal Submissions: 6251 Accepted: 1791DescriptionGiven a N × N matrix A, whose element in the i-th row and j-th col原创 2016-08-17 21:16:31 · 545 阅读 · 0 评论 -
扩展欧几里得算法
用类似辗转相除法,求二元一次不定方程 { 47x+30y=1} 的整数解47 = 30 * 1 + 1730 = 17 * 1 + 1317 = 13 * 1 + 413 = 4 * 3 + 1 //这里的 1 表示最小公约数然后把它们改写成“余数等于”的形式17 = 47 * 1 + 30 * (-1) //式113原创 2016-08-11 15:13:34 · 431 阅读 · 0 评论 -
poj 3070 矩阵快速幂
FibonacciTime Limit: 1000MS Memory Limit: 65536KTotal Submissions: 13073 Accepted: 9295DescriptionIn the Fibonacci integer sequence, F0 = 0, F1 = 1, and Fn =原创 2016-08-09 11:27:31 · 361 阅读 · 0 评论 -
矩阵乘法证明
理解矩阵乘法作者: 阮一峰日期: 2015年9月 1日大多数人在高中,或者大学低年级,都上过一门课《线性代数》。这门课其实是教矩阵。刚学的时候,还蛮简单的,矩阵加法就是相同位置的数字加一下。矩阵减法也类似。矩阵乘以一个常数,就是所有位置都乘以这个数。但是,等到矩阵乘以矩阵的时候,转载 2016-08-09 09:47:37 · 4552 阅读 · 0 评论 -
poj 3761 反逆序数+快速取模+阶乘预处理
Bubble SortTime Limit: 1000MS Memory Limit: 65536KTotal Submissions: 2724 Accepted: 927DescriptionBubble sort is a simple sorting algorithm. It works by repeate转载 2016-08-09 09:20:53 · 899 阅读 · 0 评论 -
poj1220(多种进制转换) poj3191(负数进制)
The Moronic CowmpouterTime Limit: 1000MS Memory Limit: 65536KTotal Submissions: 3767 Accepted: 1962DescriptionInexperienced in the digital arts, the cows tried原创 2016-02-01 00:11:14 · 953 阅读 · 0 评论 -
hdu 2035 (快速乘方)快速乘方取模
人见人爱A^BTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 35507 Accepted Submission(s): 24132Problem Description求A^B的最后三位数表示的整数。原创 2016-08-08 18:08:06 · 1486 阅读 · 0 评论 -
poj 1996多项式计算
The Highest ProfitsTime Limit: 1000MS Memory Limit: 30000KTotal Submissions: 656 Accepted: 286DescriptionThis program will not be part of the KOKODáKH collectio原创 2016-07-25 19:33:07 · 722 阅读 · 2 评论 -
poj 2635 素数打表+同余求模定理+千进制
The Embarrassed CryptographerTime Limit: 2000MS Memory Limit: 65536KTotal Submissions: 13831 Accepted: 3762DescriptionThe young and very promising cryptographer原创 2016-07-24 21:58:40 · 626 阅读 · 0 评论 -
中国剩余定理 poj 1006 poj2891
BiorhythmsTime Limit: 1000MS Memory Limit: 10000KTotal Submissions: 129686 Accepted: 41277DescriptionSome people believe that there are three cycles in a person原创 2016-08-19 10:20:20 · 799 阅读 · 0 评论