迷之数论
文章平均质量分 72
luer9
假如我年少有为。
展开
-
PAT 我要通过
“答案正确”是自动判题系统给出的最令人欢喜的回复。本题属于PAT的“答案正确”大派送 —— 只要读入的字符串满足下列条件,系统就输出“答案正确”,否则输出“答案错误”。得到“答案正确”的条件是:1. 字符串中必须仅有P, A, T这三种字符,不可以包含其它字符;2. 任意形如 xPATx 的字符串都可以获得“答案正确”,其中 x 或者是空字符串,或者是仅由字母 A 组成的字符串;3. 如果 aPb...原创 2018-03-18 17:31:32 · 1620 阅读 · 1 评论 -
HDU-1717 小数转换成分数 运用GCD
Ray 在数学课上听老师说,任何小数都能表示成分数的形式,他开始了化了起来,很快他就完成了,但他又想到一个问题,如何把一个循环小数化成分数呢? 请你写一个程序不但可以将普通小数化成最简分数,也可以把循环小数化成最简分数。 Input第一行是一个整数N,表示有多少组数据。 每组数据只有一个纯小数,也就是整数部分为0。小数的位数不超过9位,循环部分用()括起来。Output...原创 2018-01-22 19:19:55 · 382 阅读 · 0 评论 -
HDU 1568 Fibonacci
FibonacciTime Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 6214 Accepted Submission(s): 2962 Problem Description2007年到来了。经过2006年一年的修炼,数...原创 2018-08-07 16:59:21 · 161 阅读 · 0 评论 -
HDU 5187 zhx's contest(快速幂+快速乘法)
zhx's contestTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 3730 Accepted Submission(s): 1207 Problem DescriptionAs one of the most ...原创 2018-08-07 18:19:19 · 195 阅读 · 0 评论 -
LightOJ - 1336 Sigma Function
Sigma function is an interesting function in Number Theory. It is denoted by the Greek letter Sigma (σ). This function actually denotes the sum of all divisors of a number. For example σ(24) = 1+2+3+4...原创 2018-08-08 15:41:28 · 156 阅读 · 0 评论 -
HDU - 4828 Grids
度度熊最近很喜欢玩游戏。这一天他在纸上画了一个2行N列的长方形格子。他想把1到2N这些数依次放进去,但是为了使格子看起来优美,他想找到使每行每列都递增的方案。不过画了很久,他发现方案数实在是太多了。度度熊想知道,有多少种放数字的方法能满足上面的条件?Input 第一行为数据组数T(1<=T<=100000)。 然后T行,每行为一个数N(1<=N<=1000...原创 2018-08-08 15:50:11 · 183 阅读 · 0 评论 -
数论之费马小定理
费马小定理:假如p是素数,且(a,p)=1,那么a^(p-1)≡1(mod p)ps:a≡b(modm)表示a,b对模m的余数相同,如3三5(mod2)等证明略注意:1、费马小定理只能在 gcd(a,p)=1 条件成立时使用2、费马定理是,已知素数p,得到 。但是已知 并不能确定p是素数。3、 若 ,则p一定为合数(费马定理的逆反命题)。费马小定理在acm中的应...原创 2018-08-08 17:40:53 · 1544 阅读 · 0 评论 -
HDU 6441 Find Integer(费马大定理)
Find IntegerTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 1199 Accepted Submission(s): 340Special Judge Problem Descriptionpeople i...原创 2018-08-27 09:34:23 · 465 阅读 · 0 评论 -
LightOJ - 1370 Bi-shoe and Phi-shoe(欧拉函数)
Bamboo Pole-vault is a massively popular sport in Xzhiland. And Master Phi-shoe is a very popular coach for his success. He needs some bamboos for his students, so he asked his assistant Bi-Shoe to go...原创 2018-08-20 11:37:10 · 222 阅读 · 0 评论 -
卡特兰数
Catalenhttps://blog.csdn.net/wu_tongtong/article/details/78161211求法:Cn=C(2n,n)/(n+1);h(n)=h(n-1)*(4*n-2)/(n+1) 高精度用数组来解决,,,java大数也可以 但是我不会 hhhhhvoid catalan() //求卡特兰数{ int i, j, len, ...原创 2018-12-22 18:16:31 · 181 阅读 · 0 评论 -
HDU Big Number(斯特林公式,求位数)
Big NumberTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 45155 Accepted Submission(s): 22070 Problem DescriptionIn many applications...原创 2018-12-17 23:20:44 · 180 阅读 · 1 评论 -
康托展开和逆康托展开
康托展开:已知一个排列,求这个排列在全排列中是第几个。。。康托展开逆运算:已知在全排列中排第几,求这个排列。。。 康托展开:对于{1,2,3,...,n}生成的已经从小到大排序好的全排列x=a[n]*(n-1)!+a[n-1]*(n-2)!+...a[1]*0! a[i] 指的是位于位置i后面的数小于a[i]值的个数,后面乘的就是后面还有多少个数的阶乘 说明 :这个算...原创 2019-01-20 09:53:04 · 194 阅读 · 2 评论 -
HDU 1043 Eight (bfs 康托展开 )
EightTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 32621 Accepted Submission(s): 8488Special JudgeProblem DescriptionThe 15-puzzle ha...原创 2019-01-20 13:28:04 · 159 阅读 · 0 评论 -
Codeforces Global Round 1 C. Meaningless Operations
题目链接题目大意:给定一个a(2 ≤ a ≤2^25−1),找到一个b(0<b<a),求得最大GCD(a^b,a&b);看了别人的题解好久才明白过来,果然……这个世界本来就应该是二进制的世界。①、当 a≠ (1<<k)-1 ,我们总能找到一个 b 使得 a^b==(1<<k)-1 a&b=0 ,此时gcd为最大②、当a==...原创 2019-02-09 23:17:06 · 163 阅读 · 0 评论 -
Codeforces Global Round 1 A. Parity
题目链接题意:根据题目中的那个公式,判断最后的和是奇数还是偶数思路:(ps:ly大佬说可以…暴力…)……奇数+奇数=偶数 偶数+偶数=偶数 奇数+偶数=奇数 奇数×奇数=偶数只需判断 公式(题目里有哦) 中的 a 和 b代码:#include <bits/stdc++.h>#define ll long long#define Max 100003u...原创 2019-02-09 23:27:40 · 302 阅读 · 2 评论 -
POJ 2407 Relatives (欧拉函数 模板题)
Given n, a positive integer, how many positive integers less than n are relatively prime to n? Two integers a and b are relatively prime if there are no integers x > 1, y > 0, z > 0 such that...原创 2018-08-06 14:59:09 · 277 阅读 · 0 评论 -
继续(3n+1)猜想
卡拉兹(Callatz)猜想已经在1001中给出了描述。在这个题目里,情况稍微有些复杂。当我们验证卡拉兹猜想的时候,为了避免重复计算,可以记录下递推过程中遇到的每一个数。例如对n=3进行验证的时候,我们需要计算3、5、8、4、2、1,则当我们对n=5、8、4、2进行验证的时候,就可以直接判定卡拉兹猜想的真伪,而不需要重复计算,因为这4个数已经在验证3的时候遇到过了,我们称5、8、4、2是被...转载 2018-03-18 19:11:04 · 464 阅读 · 0 评论 -
PAT 1051. 复数乘法
复数可以写成(A + Bi)的常规形式,其中A是实部,B是虚部,i是虚数单位,满足i2 = -1;也可以写成极坐标下的指数形式(R*e(Pi)),其中R是复数模,P是辐角,i是虚数单位,其等价于三角形式 R(cos(P) + isin(P))。现给定两个复数的R和P,要求输出两数乘积的常规形式。输入格式:输入在一行中依次给出两个复数的R1, P1, R2, P2,数字间以空格分隔。输出格式:在一行...转载 2018-04-08 18:38:39 · 187 阅读 · 1 评论 -
PAT 1026. 程序运行时间(数据处理 四舍五入)
要获得一个C语言程序的运行时间,常用的方法是调用头文件time.h,其中提供了clock()函数,可以捕捉从程序开始运行到clock()被调用时所耗费的时间。这个时间单位是clock tick,即“时钟打点”。同时还有一个常数CLK_TCK,给出了机器时钟每秒所走的时钟打点数。于是为了获得一个函数f的运行时间,我们只要在调用f之前先调用clock(),获得一个时钟打点数C1;在f执行完成后再调用c...原创 2018-04-09 23:10:34 · 154 阅读 · 0 评论 -
(╥╯^╰╥) 一元三次方程求解
Description 有形如:ax3+bx2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d 均为实数),并约定该方程存在三个不同实根(根的范围在-100至100之间),且根与根之差的绝对值>=1。要求三个实根。。Input输入描述: 四个实数:a,b,c,d输入样例:1 -5 -4 20Output输出描述: 由小到大依次在同一行输出这三个实根(根与根...原创 2018-04-30 21:55:09 · 1004 阅读 · 0 评论 -
(╥╯^╰╥)求和
Description请你找一个适当的C,使得上式结果最小(2<n<10000,ai<1e9,n,ai均为整数)Input第一行一个数n 表示有几个数第二行n个数用空格隔开分别代表a1……anOutput请输出一个结果Sample Input 1 21 3Sample Output 12Sample Input 2 41 2 3 4Sample Output 24数学问题。C...原创 2018-05-01 17:08:22 · 234 阅读 · 0 评论 -
NYOJ 多边形重心问题.
多边形重心问题时间限制:3000 ms | 内存限制:65535 KB难度:5描述在某个多边形上,取n个点,这n个点顺序给出,按照给出顺序将相邻的点用直线连接, (第一个和最后一个连接),所有线段不和其他线段相交,但是可以重合,可得到一个多边形或一条线段或一个多边形和一个线段的连接后的图形; 如果是一条线段,我们定义面积为0,重心坐标为(0,0).现在求给出的点集组成的图形的面积和重心横纵坐标...转载 2018-06-03 20:13:38 · 194 阅读 · 0 评论 -
NYOJ 7 街区最短路径问题(曼哈顿距离)
街区最短路径问题时间限制:3000 ms | 内存限制:65535 KB难度:4描述一个街区有很多住户,街区的街道只能为东西、南北两种方向。住户只可以沿着街道行走。各个街道之间的间隔相等。用(x,y)来表示住户坐在的街区。例如(4,20),表示用户在东西方向第4个街道,南北方向第20个街道。现在要建一个邮局,使得各个住户到邮局的距离之和最少。求现在这个邮局应该建在那个地方使得所有住户距离之和最...原创 2018-06-05 11:29:10 · 643 阅读 · 0 评论 -
圆面分割问题
链接:https://www.nowcoder.com/acm/contest/135/F来源:牛客网 题目描述 签到题来了,送你们一个Python秒的题。 Apojacsleam来到了OI大陆,经过了连年征战,成为了一方国王。 Apojacsleam把他的王国命名为“Apo国”,Apo国的领土是一个标准的圆形。 Apojacsleam现在想封赏他的...原创 2018-07-24 09:50:06 · 1260 阅读 · 0 评论 -
HDU 1014 Uniform Generator
Uniform GeneratorTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 33113 Accepted Submission(s): 13135Problem DescriptionComputer simulations o...原创 2018-07-14 16:17:52 · 136 阅读 · 0 评论 -
数论之斐波那契
斐波那契数列(Fibonacci sequence),又称 黄金分割数列、因数学家列昂纳多·斐波那契(Leonardoda Fibonacci)以兔子繁殖为例子而引入,故又称为“ 兔子数列”,指的是这样一个数列:1、1、2、3、5、8、13、21、34、……在数学上,斐波纳契数列以如下被以 递归的方法定义:F(0)=0,F(1)=1, F(n)=F(n-1)+F(n-2)(n>=2,n∈N*...原创 2018-08-09 09:15:37 · 588 阅读 · 0 评论 -
POJ 1061 青蛙的约会
青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 127188 Accepted: 27545 Description 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上...原创 2018-07-27 18:19:38 · 178 阅读 · 0 评论 -
数论之扩展欧几里德
欧几里德求整数a,b的最小公约数gcd(a,b)的算法。递归代码:#include <iostream>using namespace std;#define ll long longll gcd(ll a,ll b){ if(b==0)return a; else return gcd(b,a%b);}int main(){ ll ...原创 2018-08-09 16:34:04 · 230 阅读 · 0 评论 -
数论之Lucas
Lucas 用来简化 上图 的 计算。例:求C(n, m) % p , n<=1e18,m<=1e18,p<=1e5 这样就会严重超时诶。。。。。。。BUT:卢卡斯说:C(n, m) % p = C(n / p, m / p) * C(n%p, m%p) % p对于C(n / p, m / p),如果n / p 还是很大,可以递归下去,一直到...原创 2018-08-09 17:02:08 · 182 阅读 · 0 评论 -
数论之欧拉函数
请思考以下问题: 任意给定正整数n,请问在小于等于n的正整数之中,有多少个与n构成互质关系?(比如,在1到8之中,有多少个数与8构成互质关系?)计算这个值的方法就叫做欧拉函数,以φ(n)表示。在1到8之中,与8形成互质关系的是1、3、5、7,所以 φ(n) = 4。φ(n) 的计算方法并不复杂,但是为了得到最后那个公式,需要一步步讨论。第一种情况如果n=1,则 φ(1) ...原创 2018-08-09 17:30:24 · 620 阅读 · 0 评论 -
快速幂详解
#include <bits/stdc++.h>using namespace std;/*快速幂*/int power(int a,int b){ int ans=1,base=a; while(b!=0) { if(b&1) ans*=base; base*=base; ...转载 2018-02-27 21:05:55 · 546 阅读 · 2 评论 -
Codeforces Global Round 1E. Magic Stones
题目链接题意,分别有两个数组 c 和 t将c数组中的元素c[i],进行 c[i]=c[i-1]+c[i+1]-c[i] 操作,进行几次操作后,判断c是否可以跟t一样。思路,显然 开头和结尾的元素是显然不会被改变的其次对于:7 2 4 127 15 10 12差分数组都是 -5,2,8。。。。改变前后,差分数组是不变的。所以最后只需要判断 开头和结尾的元素是否相同和两个数组...原创 2019-02-10 15:38:47 · 352 阅读 · 5 评论