数论
文章平均质量分 59
henuzxy
已经成为刷题过千的男人,即将成为cf 1900分的男人。
展开
-
leetcode 1739. 放置盒子 (数学,二分)
放置盒子有一个立方体房间,其长度、宽度和高度都等于 n 个单位。请你在房间里放置 n 个盒子,每个盒子都是一个单位边长的立方体。放置规则如下:你可以把盒子放在地板上的任何地方。如果盒子 x 需要放置在盒子 y 的顶部,那么盒子 y 竖直的四个侧面都 必须 与另一个盒子或墙相邻。给你一个整数 n ,返回接触地面的盒子的 最少 可能数量。示例 1:输入:n = 3输出:3解释:上图是 3 个盒子的摆放位置。这些盒子放在房间的一角,对应左侧位置。示例 2:输入:n = 4输出:3解释.原创 2021-01-26 01:24:23 · 705 阅读 · 3 评论 -
一些数学几何知识和小技巧
发现自己几何sb,还是总结一下已知三点求圆心,但三点不能共线Point Getcir(Point A,Point B,Point C){//给予三个点,求圆心。 double a = 2*(B.x - A.x); double b = 2*(B.y - A.y); double c = (B.x*B.x+B.y*B.y) - (A.x*A.x+A.y*A.y);...原创 2017-11-12 14:40:29 · 535 阅读 · 0 评论 -
D - Noldbach problem CodeForces - 17A
Nick is interested in prime numbers. Once he read about Goldbach problem. It states that every even integer greater than 2 can be expressed as the sum of two primes. That got Nick’s attention and he d原创 2017-02-08 23:16:24 · 675 阅读 · 0 评论 -
ZOJ-3609(Modular Inverse )(拓展欧几里得)
The modular modular multiplicative inverse of an integer a modulo m is an integer x such that a-1≡x (mod m). This is equivalent to ax≡1 (mod m).Input There are multiple test cases. The first line o...原创 2017-04-16 22:49:41 · 431 阅读 · 0 评论 -
快速幂的两种算法
快速幂是一种在O(logN)的时间复杂度下求出幂的一种方法。 第一种算方法就是对于一个X^N我们可以把N转化成2进制的形式,如X^22 = X^16+X^4+X^2;因为22(二进制:10110) = 16 + 4 + 2; 所以我们可以推到出来这种算法,完整代码如下:typedef long long ll;ll mod_pow(ll x,ll n,ll mod){ ll ...原创 2017-04-20 12:25:23 · 548 阅读 · 0 评论 -
2017广东工业大学程序设计竞赛决赛-E(倒水)
Problem E: 倒水(Water) Description一天,CC买了N个容量可以认为是无限大的瓶子,开始时每个瓶子里有1升水。接着~~CC发现瓶子实在太多了,于是他决定保留不超过K个瓶子。每次他选择两个当前含水量相同的瓶子,把一个瓶子的水全部倒进另一个里,然后把空瓶丢弃。(不能丢弃有水的瓶子) 显然在某些情况下CC无法达到目标,比如N=3,K=1。此时CC会重新买一些新的瓶子(新...原创 2017-03-28 21:29:18 · 559 阅读 · 0 评论 -
“景驰科技杯”2018年华南理工大学程序设计竞赛 H 对称与反对称(逆元)
链接:https://www.nowcoder.com/acm/contest/94/H 来源:牛客网题目描述 给出一个N*N的方阵A。构造方阵B,C: 使得A = B + C.其中 B为对称矩阵,C为反对称矩阵。 对于方阵S中的任意元素,若(S)ij = (S)ji,则称S为对称矩阵 对于方阵T中的任意元素,若(T)ij = -(T)ji,则称T为反对称矩阵 注意,所有运算在模...原创 2018-04-07 21:31:30 · 267 阅读 · 0 评论 -
codeforces 964C. Alternating Sum(等比数列)
题意简单,给你一个序列s0,s1,s2…sn,共n+1个数,这个序列是满足k循环,即s[i] = s[i+k],且这些值只可能是-1,或1,且(n+1) % k == 0。问你上图式子的和是多少。思路:因为这个序列是满足k循环的,所以我们应该主要考虑a,b的变化,如果把这个求和公式展开后我们可以发现如果我们按k个隔开,发现他们是符合等比序列的,公比q = (b/a)^k 。所以我们可以利用...原创 2018-04-18 22:08:08 · 654 阅读 · 2 评论 -
康托展开与逆康托展开
康托展开可以理解为把一个全排列映射到一个数上面,因为全排列如果按照从小到大或者从大到小,肯定是有一个确定的序列的。一般是从小到大的序列个数。我们就是要求出这个序列的位置。,想法很简答,就是求出前面比他小的个数就可以了。 理解为一个每位都是阶乘进位的数转化为10进制的数。思路如下: 先准备求每一位的阶乘,然后从高位开始统计后面有多少个数比他小记录这个个位数,然后乘以后面个数的阶乘,再把它...原创 2018-07-18 10:30:50 · 736 阅读 · 0 评论 -
牛客网暑期ACM多校训练营(第二场)I car (规律)
题目链接题意非常绕,最开始直接就理解错了,发了公告才明白。题目的意思是给你一个N*N的地图,然后又M个障碍,每个车都只能放到边上,然后行驶到另一边,边角可以有两个方向的选择。然后如果两个车在一个地方相遇会发生碰撞,每个车速度相同,问你在不发生碰撞的情况下,地图最多能放多少辆车。自己放一下就会发现个规律,每次放的车都是按某个方向旋转,这样是最多的,不过当地图是奇数时,可以在中间多放一个。 ...原创 2018-07-21 18:32:09 · 190 阅读 · 0 评论 -
BZOJ 3884: 上帝与集合的正确用法 欧拉降幂
根据一些书上的记载,上帝的一次失败的创世经历是这样的: 第一天, 上帝创造了一个世界的基本元素,称做“元”。 第二天, 上帝创造了一个新的元素,称作“α”。“α”被定义为“元”构成的集合。容易发现,一共有两种不同的“α”。 第三天, 上帝又创造了一个新的元素,称作“β”。“β”被定义为“α”构成的集合。容易发现,一共有四种不同的“β”。 第四天, 上帝创造了新的元素...原创 2018-08-01 23:37:00 · 190 阅读 · 0 评论 -
FZU 1759 Super A^B mod C (欧拉降幂)
Given A,B,C, You should quickly calculate the result of A^B mod C. (1<=A,C<=1000000000,1<=B<=10^1000000).Input There are multiply testcases. Each testcase, there is one line contains th...原创 2018-08-02 18:02:04 · 159 阅读 · 0 评论 -
Codeforces Round #554 (Div. 2) C. Neko does Maths(数论)
C. Neko does Mathstime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputNeko loves divisors. During the latest number theory lesson, he got an inter...原创 2019-04-25 23:20:46 · 519 阅读 · 2 评论 -
HDU 5976 Detachment (数学规律+逆元)
In a highly developed alien society, the habitats are almost infinite dimensional space. In the history of this planet,there is an old puzzle. You have a line segment with x units’ length represent原创 2017-12-07 12:51:04 · 294 阅读 · 0 评论 -
haut oj 1299 fireworks(杨辉三角规律)
题目链接1299: fireworks 时间限制: 1 秒 内存限制: 64 MB 提交: 42 解决: 16 提交 状态 题目描述Hmz likes to play fireworks, especially when they are put regularly. Now he puts some fireworks in a line. This time he put a t原创 2017-10-02 15:07:33 · 375 阅读 · 0 评论 -
HDu 6216 && 2017 ACM/ICPC Asia Regional Qingdao Online 1011
A Cubic number and A Cubic NumberTime Limit: 1000/1000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others) Total Submission(s): 89 Accepted Submission(s): 52Problem Description A cubic n原创 2017-09-17 19:30:30 · 299 阅读 · 0 评论 -
SCU 4437 Carries(数论)
Carriesfrog has (n) integers (a_1, a_2, \dots, a_n), and she wants to add them pairwise.Unfortunately, frog is somehow afraid of carries (进位). She defines hardness (h(x, y)) for adding (x) and (y) the原创 2017-03-18 22:07:03 · 646 阅读 · 0 评论 -
zstu oj 4274 约素(求一个数约数的个数)(暴力法)
Problem I: 约素 Time Limit: 1 Sec Memory Limit: 128 MB Submit: 1808 Solved: 468 Description 判断一个正整数n的约数个数是否为p,其中p是素数。 Input 第一行给测试总数T(T <= 10000)。 接下来有T行,每行有两个数字n(1 <= n <= 1000000000)和p(2 < p <原创 2017-03-22 14:45:46 · 876 阅读 · 0 评论 -
求最大公约数和最小公倍数(利用递归实现)
今天看了学长一段代码真是惊为天人啊。原来他们最大公约数gcd都是这样求的:#include<stdio.h>#include<string.h>#include<math.h>int gcd(int x,int y){ return y? gcd(y,x%y):x; } /*这里可以看出如果是x比y小时,x对y取余还是x,x则变为y,只比较一次就相当于交换了位置可以看出,因为这原创 2017-03-22 15:08:59 · 1586 阅读 · 0 评论 -
HDU-1576(A/B)拓展欧几里得
Problem Description 要求(A/B)%9973,但由于A很大,我们只给出n(n=A%9973)(我们给定的A必能被B整除,且gcd(B,9973) = 1)。Input 数据的第一行是一个T,表示有T组数据。 每组数据有两个数n(0 <= n < 9973)和B(1 <= B <= 10^9)。Output 对应每组数据输出(A/B)%99...原创 2017-04-17 23:09:49 · 349 阅读 · 0 评论 -
HDU-1211(RSA)-拓展欧几里得
Problem Description RSA is one of the most powerful methods to encrypt data. The RSA algorithm is described as follow: choose two large prime integer p, q calculate n = p × q, calculate F(n) = (p原创 2017-04-19 13:28:54 · 395 阅读 · 0 评论 -
河南工业“玲珑杯”最大字段和
这是一道类似DP的题,不过非dp的方法也能做,比赛时总是无法确定区域的扫描,真是菜啊。 题目描述一个大小为n的数组a1到an(?10^4≤ai≤10^4)。请你找出一个连续子段,使子段长度为奇数,且子段和最大。输入第一行为T(1≤T≤5),代表数据组数。 之后每组数据,第一行为n(1≤n≤10^5),代表数组长度。 之后一行n个数,代表a1到an。输出每组数据输出一行,表示满足要求的子段和最大原创 2017-05-06 21:33:27 · 245 阅读 · 0 评论 -
河南工业大学“玲珑杯”地狱飞龙(高数上微积分~~)
一道微积分题目啊。链接如下: http://218.28.220.249:50015/JudgeOnline/problem.php?id=1261 被我舍友一下就把微积分求出来了,年级第一果然强啊,看来是时候拉他入坑了题目描述最近clover迷上了皇室战争,他抽到了一种地狱飞龙,很开心。假设地域飞龙会对距离为d的敌人每秒造成k/d2伤害。假设地域飞龙位于坐标轴原点,以每秒v1的速度向y轴正方向移原创 2017-05-07 13:45:01 · 948 阅读 · 0 评论 -
阶乘中因子k的个数
参考例题 zzuli 2177传送门题意就是求n!里面因子2的个数,改怎么去做呢?我们可以思考2的倍数是每隔一个就会出现的,例如1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 是16!里面的全部的数,我们可以发现2的倍数的数是哪些,我们可以先对这些数都除以2,然后他们2的因子的数就是变为1,2,3,4,5,6,7,8,这里由于8个数都除以了2,所以我们已经提取出来8个2原创 2017-08-18 21:52:14 · 1417 阅读 · 0 评论 -
hpu 1410: QAQ & 火星情报局(哥德巴赫猜想定理)
题目传送门哥德巴赫猜想传送利用了哥德巴赫的两个定理 任何一个大于2的偶数都可以写成两个素数之和, 任何一个大于5的整数都可以写成三个素数之和。以及一个推论 存在大于5的整数可以写为2和一个素数之和(这是推论)代码如下#include<cstdio>#include<cmath>#include<iostream>#include<cstring>#include<strin原创 2017-08-12 12:04:50 · 593 阅读 · 0 评论 -
费马小定理,拓展欧几里得求逆元。zzuli 2180
逆元:在mod m的运算下,像满足 ax%m=1ax\%m = 1这样的求解出的x,我们称作a的逆元即a−1a^{-1} 首先简介费马小定理 传送门 对于p是素数时,任意正整数都有(xp)%p=x(x^p)\%p= x 然后可以推出(xp−1)%p=1(x^{p-1} ) \%p = 1最后可以变形得到x−1%p=xp−2x^{-1} \% p = x^{p-2} 所以我们对于一个数求原创 2017-08-20 22:54:42 · 466 阅读 · 0 评论 -
codeforce 850A Five Dimensional Points(特殊判别)
传送You are given set of n points in 5-dimensional space. The points are labeled from 1 to n. No two points coincide.We will call point a bad if there are different points b and c, not equal to a, from t原创 2017-09-05 19:14:05 · 816 阅读 · 0 评论 -
POJ 3233(矩阵快速幂的转化法)
传送门 Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak.Input The input contains exactly one test case. The first line of input contains three positive integers n (原创 2017-08-28 11:23:36 · 227 阅读 · 0 评论 -
将一个数分解质因数(初级篇)
如何将一个数分解质因数呢? 就是直接一个for循环下标i从2开始如果这个数能整除i就不断整除i然后记录,直到这个数为1时停止这样就可以达到分解质因数的目的了。 代码如下:#include<stdio.h>#include<string.h>#include<math.h>int main(void){ int n,i,j,k,a; int ans[1000],top;原创 2017-02-27 22:46:41 · 2876 阅读 · 0 评论