
数学数论
文章平均质量分 78
oorzhang
这个作者很懒,什么都没留下…
展开
-
逆元详解
今天我们来探讨逆元在ACM-ICPC竞赛中的应用,逆元是一个很重要的概念,必须学会使用它。 对于正整数和,如果有,那么把这个同余方程中的最小正整数解叫做模的逆元。 逆元一般用扩展欧几里得算法来求得,如果为素数,那么还可以根据费马小定理得到逆元为。 推导过程如下 求现在来看转载 2014-07-23 15:18:25 · 707 阅读 · 0 评论 -
FFT 快速傅立叶变换 dhu1402 A*B
A * B Problem PlusTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 14007 Accepted Submission(s): 2538Problem DescriptionCalcul原创 2015-04-23 19:44:02 · 714 阅读 · 0 评论 -
HDU4651 Partition
PartitionTime Limit: 6000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 916 Accepted Submission(s): 527Problem DescriptionHow many ways can the原创 2015-05-14 21:45:35 · 361 阅读 · 0 评论 -
hdu 4658 Integer Partition
Integer PartitionTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 475 Accepted Submission(s): 218Problem DescriptionGiven n, k, ca原创 2015-05-14 22:40:02 · 407 阅读 · 0 评论 -
HDU5490 Simple Matrix
Simple MatrixTime Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 128 Accepted Submission(s): 35Problem DescriptionAs we know, sequenc原创 2015-09-28 22:20:36 · 924 阅读 · 0 评论 -
数列组合数求和( HDU5490 Simple Matrix)
Simple MatrixTime Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 128 Accepted Submission(s): 35Problem Description As we know, sequence in the原创 2016-03-20 00:16:43 · 1870 阅读 · 0 评论 -
FFT快速傅立叶变换 HDU4609 3-idots
3-idiotsTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 2656 Accepted Submission(s): 904Problem DescriptionKing OMeGa catche原创 2015-04-23 21:12:55 · 727 阅读 · 0 评论 -
A^B所有因子和(poj1845 Sumdiv)
题意: 求ABA^B所有因子的和。mod 9901思路: 把A表示成:∑piki\sum{p_i}^{k_i}. 那么所有因子的和就为:(∑pi1)(∑pi2)...(∑pin)(∑p_1^i)(∑p_2^i)...(∑p_n^i); ABA^B就相应让指数变为BB倍。 然后再用等比数列求和公式。#include <iostream>#include <cmath>#include <原创 2016-03-20 14:04:09 · 365 阅读 · 0 评论 -
素数打表,思维(cf300E)
cf300E题目大意:给出k个a1,a2...ak(1<=k<=10^10),求最小的n,使得n!被a1!*a2!*...*ak!整除;原创 2014-05-22 21:33:36 · 910 阅读 · 0 评论 -
∑(Fib(N)%Fib(k))(51nod 1399 Fib(N) mod Fib(K) V2)
题目连接Fib(N)表示斐波那契数列的第N项(F(0) = 0, F(1) = 1),给出N和K,求 ∑i=1N[Fib(i) mod Fib(K)] 由于结果太大,输出Mod 1000000007的结果。 Input 第1行:一个数T,表示后面用作输入测试的数的数量。(1 <= T <= 50000) 第2 - T + 1行:每行2个数N, K(1 <= N, K <= 10^18)原创 2016-03-21 21:57:52 · 1532 阅读 · 0 评论 -
浅谈一类积性函数的前缀和
前置技能积性函数的定义若f(n)的定义域为正整数域,值域为复数,即f:Z+→C,则称f(n)为数论函数。若f(n)为数论函数,且f(1)=1,对于互质的正整数p,q有f(p⋅q)=f(p)⋅f(q),则称其为积性函数。若f(n)为积性函数,且对于任意正整数p,q都有f(p⋅q)=f(p)⋅f(q),则称其为完全积性函数。积性函数的性质与例子常见的积性函数有转载 2016-04-09 16:48:15 · 1484 阅读 · 0 评论 -
约数个数和(数论,莫比乌斯反演)BZOJ3994
题目3994: [SDOI2015]约数个数和Time Limit: 20 Sec Memory Limit: 128 MB Submit: 492 Solved: 319 [Submit][Status][Discuss] Description设d(x)为x的约数个数,给定N、M,求 Input输入文件包含多组测试数据。 第一行,一个整数T,表示测试数据的组数。 接下来的T行,每原创 2016-03-20 12:06:55 · 586 阅读 · 0 评论 -
51nod 1192 Gcd表中的质数
有一个M * N的表格,行与列分别是1 - M和1 - N,格子中间写着行与列的最大公约数Gcd(i, j)(1 <= i <= M, 1 <= j <= N)。 例如:M = 5, n = 4。1 2 3 4 5 1 1 1 1 1 1 2 1 2 1 2 1 3 1 1 3 1 1 4 1 2 1 4 1给出M和N,求这张表中有多少个质数。 Input 第1行:一个数T,表示后面用原创 2016-05-18 13:34:51 · 985 阅读 · 0 评论 -
2016广东省省赛H题
题意:计算:∑i=1n∑j=1mi2j2gcd(i,j)\sum_{i=1}^n\sum_{j=1}^mi^2j^2gcd(i,j)(时限1000ms)input:T(t<=1000),接下来t行,m,n(n,m<1e6).题目难点在于多了i,j的平方,由对称性,设n=====d′=kd===∑i=1n∑j=1mi2j2gcd(i,j)∑i=1n∑j=1m∑d|i,d|ji2j2d[(id,jd)=原创 2016-05-08 23:35:57 · 2642 阅读 · 0 评论 -
莫比乌斯反演证明
首先定义几个概念:1,卷积:设是两个数论函数(也就是说,以自然数集为定义域的复数值函数),则卷积运算定义为可以证明,卷积运算满足:1)交换律:由定义显然。2)结合律:考察两边作用在上,左边是右边是故两边相等。3)存在单位元使得我们需要故不难猜到应该定义为事实上,直接验证可得以上说明数论函数在卷积意义下构成一转载 2015-04-20 00:39:56 · 1130 阅读 · 0 评论 -
2015华农校赛B
B TowerTime Limit:1000MS Memory Limit:65535K题型: 编程题 语言: 无限制描述There is a tower, which have n+1 floors, numbered from 0 to n. Each floor have m distinct doors, one of them i原创 2015-04-07 21:07:16 · 683 阅读 · 0 评论 -
华农校赛E题
E FibonacciLandTime Limit:1000MS Memory Limit:65535K题型: 编程题 语言: 无限制描述Zayhero like travelling very much, and his current stop is FibonacciLand.FibonacciLand consist of tota原创 2015-04-06 23:08:32 · 501 阅读 · 0 评论 -
POJ2151 概率题
Check the difficulty of problemsTime Limit: 2000MS Memory Limit: 65536KTotal Submissions: 4684 Accepted: 2050DescriptionOrganizing a programming contest is not an原创 2014-07-21 23:54:59 · 394 阅读 · 0 评论 -
Catalan数计算及应用
原著:http://blog.youkuaiyun.com/wuzhekai1985 应用1描述:n对括号有多少种匹配方式? 思路:n对括号相当于有2n个符号,n个左括号、n个右括号,可以设问题的解为f(2n)。第0个符号肯定为左括号,与之匹配的右括号必须为第2i+1字符。因为如果是第2i个字符,那么第0个字符与第2i个字符间包含奇数个字符,而奇数个字符是无法构成匹配的。转载 2015-01-20 21:35:33 · 474 阅读 · 0 评论 -
hdu4870 Rating (高斯消元或者dp)
RatingTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 213 Accepted Submission(s): 126Special JudgeProblem DescriptionA little girl lov原创 2014-07-23 11:01:58 · 2026 阅读 · 6 评论 -
Rating ,hdu4870(高校联赛第一场)
有m个扑克牌,初始全部背面朝上有n个操作,每个操作可以同时对xi张牌反转,问最后有多少种情况。原创 2014-07-23 15:58:45 · 551 阅读 · 0 评论 -
hdu4828卡特兰数
hdu4828,卡特兰数原创 2014-05-27 12:14:24 · 570 阅读 · 0 评论 -
poj3292 Semi-prime H-numbers
题意:H-number 是一种4*n+1的数。H-primes 是只可以被1和本身整除的H-numbers(在4n+1的数域内)。H-semi-primes 是有且仅有两个H-numbers因子(除了1和本身)。问给你一个n,求从1到n一共有多少个H-semi-primes数。原创 2015-01-21 14:51:02 · 324 阅读 · 0 评论 -
poj1142 Smith Numbers
题目大意:一个数他的各个数字的和等于他的各个质因子的各个数字的和,这个数就称为smith Number求n开始的下一个Smith Number。思路:暴力搜。#include #include #include #include #include using namespace std;#define LL long long #defin原创 2015-01-21 16:08:42 · 493 阅读 · 0 评论 -
poj3358 Period of an Infinite Binary Expansion
题目大意:给你一个分数,就他二进制的循环节。求开始循环的点和循环长度。思路:已知t*2^phi(r)=t (mod r),先求出phi(r),然后求phi(r)的因子,求一个最小符合条件的。这样就可以大大缩小了查找的范围。#include #include #include #include #include using n原创 2015-01-21 16:35:37 · 382 阅读 · 0 评论 -
poj3682 King Arthur's Birthday Celebration
题目:国王庆祝生日,第i天花费2i-1(千)个硬币,至于哪一天结束,给定一个正整数k,每天抛一次硬币,如果出现第k次正面朝上,就结束庆祝,正面朝上的概率为p。问庆祝的天数和花费的金币的期望。思路:纯概率题,设f[i]为第i天结束的概率,则f[i]=c(i-1,k-1)*p^k*(1-p)^(i-k),然而总概率为 ∑f[i]=1,其中i∈[k,+∞),以下均为此原创 2015-01-21 15:25:24 · 353 阅读 · 0 评论 -
poj1604 Just the Facts
题意:计算n!最后一位不为0的数。思路:模拟一下:#include#includeusing namespace std;int main(){ int n; while(cin>>n) { int tmp = 1; for(int i = 1; i <= n; i++) {原创 2015-01-21 16:52:09 · 453 阅读 · 0 评论 -
poj1150
题意:计算A(n,m)的最后一位非0的数。分析引用下面文章。出处:http://duanple.blog.163.com/blog/static/7097176720081016113033592/http://acm.pku.edu.cn/JudgeOnline/problem?id=1150最最原始的问题,是那个计算n!的结果中的0的个数原创 2015-01-21 17:08:53 · 793 阅读 · 0 评论 -
poj1150 求排列Anm的末尾非0数字
题意:求nPm,即n个元素的m中排列方式结果中末尾非0 的数字。思路: 首先我们知道,2和5相乘末尾非0数字为1,相当于可以相消。我们先把这两个因子抽离出来,考虑末尾为3,7,9的数字。设f(n,x)为:n!因子中,抽离了2,5,后末尾数字为x的因子个数。分析可以知道f(n,x)=f(n/2,x)+g(n,x);g(n,x)表示n!中所有奇因子末尾数字为x的因子个数。可以这原创 2015-03-25 20:19:41 · 430 阅读 · 0 评论 -
ZOJ 3868 GCD Expectation 莫比乌斯反演
ZOJ Problem Set - 3868GCD ExpectationTime Limit: 4 Seconds Memory Limit: 262144 KBEdward has a set of n integers {a1, a2,...,an}. He randomly picks a nonempty subset {x1, x2,…,xm} (原创 2015-04-22 20:26:45 · 435 阅读 · 0 评论 -
lcm/gcd
A(n)=∑ni=1lcm(i,n)gcd(i,n)A(n)=\sum_{i=1}^n\frac{lcm(i,n)}{gcd(i,n)} ,F(n)=∑ni=1A(i) F(n)=\sum_{i=1}^nA(i),求F(n)F(n)%1000000007A(n)=∑i=1nlcm(i,n)gcd(i,n)=∑i=1nn∗i(i,n)2=∑i=1n∑d|n[(i,n)=d]n∗id2=∑d|n∑i=1原创 2017-04-23 22:27:03 · 949 阅读 · 0 评论