★数学
墓华
这个作者很懒,什么都没留下…
展开
-
AtCoder Beginner Contest 172 E - NEQ(二项式反演)
AtCoder Beginner Contest 172 E - NEQ题意:求满足下列条件的长度为 NNN且包含 [1,M][1, M][1,M]范围内整数的序列 A1,A2,⋯ ,ANA_1, A_2, \cdots, A_NA1,A2,⋯,AN及 B1,B2,⋯ ,BNB_1, B_2, \cdots, B_NB1,B2,⋯,BN组成的序列对数量。∀1≤i<j≤N, Ai≠Aj∧Bi≠Bj\forall 1 \le i \lt j \le N, \; A_i \neq A_原创 2020-06-30 16:34:34 · 607 阅读 · 0 评论 -
AtCoder Beginner Contest 172 F - Unfair Nim(Nim博弈,位运算)
AtCoder Beginner Contest 172 F - Unfair Nim题意有 N(2≤N≤300)N(2\le N \le300)N(2≤N≤300)堆石子,其中第 iii堆有 Ai(1≤Ai≤1012)A_i(1\le A_i\le 10^{12})Ai(1≤Ai≤1012)个石子,两人进行尼姆(Nim)游戏,问至少要将多少石子从第111堆移到第222堆(不可全部移走),才能使得先手必败。不存在可行解时输出−1-1−1。分析Nim博弈的必败局面要求:A1⊕A2⊕A3⊕⋯⊕A原创 2020-06-30 10:59:13 · 601 阅读 · 0 评论 -
牛客挑战赛34 - D 拉普兰德的愿望(曼哈顿距离转切比雪夫距离)
链接:牛客挑战赛34 - D 拉普兰德的愿望题意:给出平面上N(≤100,000)N(\le100,000)N(≤100,000)个点,坐标绝对值不超过L(≤50,000)L(\le50,000)L(≤50,000),求曼哈顿距离不小于d(≤10,000,000)d(\le10,000,000)d(≤10,000,000)的点对共有多少?对于点对 (x1,y1)(x_1,y_1)(x1,y...原创 2019-11-23 15:24:34 · 364 阅读 · 0 评论 -
2019多校第二场 HDU6601 Keen On Everything But Triangle(静态主席树,斐波那契数列)
链接: HDU6601 Keen On Everything But Triangle题意:给出一段长度为N的序列a[1]、a[2]、… 、a[N],每次询问[L,原创 2019-07-28 19:01:45 · 237 阅读 · 0 评论 -
LightOJ - 1236 Pairs Forming LCMPairs Forming LCM(唯一分解定理+LCM)
链接:LightOJ - 1236 Pairs Forming LCMPairs Forming LCM题意:给出正整数nnn,问有多少对(a,b)(a,b)(a,b)满足 lcm(a,b)=nlcm(a,b)=nlcm(a,b)=n 且 a≤b≤na\le b\le na≤b≤n分析:这就涉及到gcd和lcm的另一种表示方法了,对于a,ba,ba,b,根据唯一分解定理得到标准分解式:...原创 2019-08-05 19:59:41 · 258 阅读 · 0 评论 -
LightOJ - 1245 Harmonic Number (II)(数论分块)
链接:LightOJ - 1245 Harmonic Number (II)题意:给出T(≤1000)T(\le1000)T(≤1000)个n(1≤n≤231)n(1\le n\le2^{31})n(1≤n≤231),每次求解H(n)=[n1]+[n2]+⋯+[nn]H(n)=[\frac{n}{1}]+[\frac{n}{2}]+\dots+[\frac{n}{n}]H(n)=[1n]+[...原创 2019-08-05 09:42:28 · 257 阅读 · 0 评论 -
2019多校第三场 HDU6608 Fansblog(威尔逊定理,逆元,质数间隔)
链接:HDU6608 Fansblog题意:给出一个质数PPP(109≤P≤101410^9\le P\le 10^{14}109≤P≤1014),找出最大的质数QQQ(Q<PQ \lt PQ<P),并求出Q!mod  PQ!\mod PQ!modP分析:对于101410^{14}1014以内,已知相邻质数的最...原创 2019-08-01 10:02:14 · 309 阅读 · 0 评论 -
LightOJ - 1282 Leading and Trailin(快速幂,对数求大数前3位)
链接:LightOJ - 1282 Leading and Trailin题意:给出n  (2≤n<231)n\;(2\le n\lt 2^{31})n(2≤n<231) 和 k  (1≤k≤107)k\;(1\le k\le10^7)k(1≤k≤107),问nkn^knk的前三位和后三位分别为多少?(题目保证nkn...原创 2019-08-04 18:29:44 · 165 阅读 · 0 评论 -
快速幂,快速乘 及 欧拉降幂
【快速幂】abmod  ca^b\mod cabmodclong long qpow(long long a,long long b,long long c) //(a^b)%c{ long long res=1; while(b) { if(b&1) //若b为奇数 ...原创 2019-08-04 16:23:34 · 238 阅读 · 0 评论 -
欧拉函数φ(x)相关性质及计算
定义:对 正整数nnn,欧拉函数是 小于nnn的正整数 中 与nnn互质的数的数目。(φ(1)=1\varphi(1)=1φ(1)=1)互质:aaa与bbb互质,即gcd(a,b)=1\gcd(a,b)=1gcd(a,b)=1性质:对于质数ppp,φ(p)=p−1\varphi(p)=p-1φ(p)=p−1欧拉函数φ(x)\varphi(x)φ(x)为(不完全)积性...原创 2019-08-04 14:31:03 · 5090 阅读 · 0 评论 -
素数表:埃氏素筛 及 欧拉筛
一、埃氏素数筛选法时间复杂度O(N⋅loglogN)O(N\cdot\log{\log N})O(N⋅loglogN)int prime[maxn],tot;bool vis[maxn];void get_prime() //埃氏素筛{ memset(vis,0,sizeof(vis)); tot=0; for(int i=2;i<=maxn;...原创 2019-08-04 11:16:52 · 352 阅读 · 0 评论 -
2019多校第五场 HDU6627 equation(分数,解绝对值方程)
链接:HDU6627 equation题意:解方程:∑i=1N∣ai⋅x+bi∣=C\displaystyle\sum_{i=1}^{N}\left|a_i\cdot x+b_i\right|=Ci=1∑N∣ai⋅x+bi∣=C其中1≤N≤1051\le N\le 10^51≤N≤1051≤ai≤10001\le a_i\le 10001≤ai≤1000−1000≤...原创 2019-08-06 14:21:51 · 247 阅读 · 0 评论 -
辗转求分数中间值
已知正整数a,b,c,da,b,c,da,b,c,d,有:(均是最简分数形式)ab<xy<cd\frac{a}{b}\lt\frac{x}{y}\lt\frac{c}{d}ba<yx<dc要令xxx和yyy都尽可能的小,问x,yx,yx,y最小是多少?首先,因为xy\frac{x}{y}yx是被限定在区间内的,所以要令xxx小,就要yyy...原创 2019-08-06 15:37:07 · 385 阅读 · 0 评论 -
2019多校第五场 HDU6624 fraction(辗转求分数中间值)
链接:HDU6624 fraction题意:求出最小的正整数bbb,使得存在正整数a<ba\lt ba<b,满足a≡bx(modp)a\equiv bx\pmod pa≡bx(modp)p,xp,xp,x已知,且ppp为质数,1<x<p≤10151\lt x\lt p\le10^{15}1<x<p≤1015分析:和解同余方程...原创 2019-08-06 16:03:53 · 333 阅读 · 0 评论 -
Gym101964 - B Broken Watch(排列组合)
链接: Gym101964题意:给出一个分成N (2≤N≤232)N\,(2\le N\le 2^{32})N(2≤N≤232)个刻度的钟,钟上有333个长度分别为A,B,CA,B,CA,B,C的指针,指针都指向刻度,以针尖为三角形顶点,问一共能形成多少个包含圆心的三角形?(可以圆心在边上)答案对2642^{64}264取余。分析:分析当三指针长度不同时的情况:若NNN为偶数,答案...原创 2019-09-25 17:58:55 · 342 阅读 · 0 评论 -
中国剩余定理 及 拓展中国剩余定理模板
求解同余方程组:{x≡r1(modm1)x≡r2(modm2)⋯x≡rn(modmn)\left\{\begin{matrix}x\equiv r_1 \pmod {m_1}\\ x\equiv r_2 \pmod {m_2}\\ \cdots\\ x\equiv r_n \pmod {m_n}\\ \end{matrix}\right.⎩⎪⎪⎨⎪⎪⎧x≡r1(modm1)x≡r...原创 2019-08-26 20:07:22 · 244 阅读 · 0 评论 -
CodeForces - 559C Gerald and Giant Chess(组合数,递推求阶乘逆元)
链接:CodeForces - 559C Gerald and Giant Chess题意:给出一个hhh行www列的棋盘(1≤h,w≤1051\le h,w\le10^51≤h,w≤105),其中有n  (1≤n≤2000)n\;(1\le n\le 2000)n(1≤n≤2000)个点(xi,yi)(x_i,y_i)(xi,yi)不可走,每次只能往右或者往下...原创 2019-08-12 10:42:49 · 166 阅读 · 0 评论 -
LightOJ - 1336 Sigma Function(唯一分解定理,n以内平方数个数)
链接:LightOJ - 1336 Sigma Function题意:σ(x)\sigma(x)σ(x)表示xxx的所有因数之和,给出n  (1≤n≤1012)n\;(1\le n\le10^{12})n(1≤n≤1012),问111 ~ nnn中有多少数的σ\sigmaσ值是偶数?分析:根据唯一分解定理,有x=p1a1⋅p2a2⋯pkakx=p_1^{a_...原创 2019-08-11 15:38:48 · 253 阅读 · 0 评论 -
POJ - 2115 C Looooops(拓展欧几里得,解同余方程)
链接:POJ - 2115 C Looooops题意:一个C语言的for循环:for(int i=A;i!=B;i+=C)st;其中i表示变量,AAA、BBB和CCC分别表示初值、终值和步长,ststst表示循环体。要求计算当循环变量i运算在k位无符号整数体系下(运算值范围:000 ~ 2k−12^k-12k−1,此时,iii+=CCC溢出后截断),循环体ststst会执行多少次?...原创 2019-08-11 13:35:51 · 209 阅读 · 0 评论 -
LightOJ - 1138 Trailing Zeroes (III)(二分,阶乘质因数分解)
链接:LightOJ - 1138 Trailing Zeroes (III)题意:共T  (≤10000)T\;(\le10000)T(≤10000)组测试数据,每组数据给出一个Q  (1≤Q≤108)Q\;(1\le Q\le 10^8)Q(1≤Q≤108),要求出最小的正整数NNN,使得N!N!N!的末尾恰有QQQ个000分析...原创 2019-08-11 13:21:10 · 182 阅读 · 0 评论 -
LightOJ - 1220 Mysterious Bacteria(唯一分解定理)
链接:LightOJ - 1220 Mysterious Bacteria题意:给出xxx,求最大的整数eee满足:存在整数yyy,使得x=yex=y^ex=yexxx是一个32位有符号整数,且保证∣x∣≥2|x|\ge2∣x∣≥2分析:先假设xxx为正,根据唯一分解定理,得到x=p1a1⋅p2a2⋯pnanx=p_1^{a_1}\cdot p_2^{a_2}\cdots p_n^{...原创 2019-08-07 19:58:08 · 234 阅读 · 0 评论 -
LightOJ - 1197 Help Hanzo(埃氏素筛在线求区间素数个数)
链接:LightOJ - 1197 Help Hanzo题意:共T(≤200)T(\le 200)T(≤200)组数据,每组数据给出a,b  (1≤a≤b<231,b−a≤100000)a,b\;(1\le a\le b\lt2^{31},b-a\le 100000)a,b(1≤a≤b<231,b−a≤100000),问区间[a,&Thi...原创 2019-08-07 11:05:34 · 201 阅读 · 0 评论 -
LightOJ - 1234 Harmonic Number(调和级数,分块打表)
链接:LightOJ - 1234 Harmonic Number题意:给出T(≤10000)T(\le 10000)T(≤10000)个n(1≤n≤108)n(1\le n\le10^8)n(1≤n≤108),求调和级数H(n)=∑k=1n1k        &ThickS...原创 2019-08-06 18:54:39 · 177 阅读 · 0 评论 -
质因数分解(N, N!, 大数) 及 唯一分解定理
质因数分解把一个合数分解成若干个质因数的乘积的形式,即求质因数的过程叫做分解质因数。①短除法分解质因数只针对合数。(分解质因数也称分解素因数)求一个数分解质因数,要从最小的质数除起,一直除到结果为质数为止。暴力解法,...原创 2019-08-04 10:50:14 · 3425 阅读 · 1 评论 -
数论四大定理
注:(modp)\pmod p(modp),表示 恒等式两边同时对ppp取模,即恒等式两边同余一、威尔逊定理若 ppp为质数,则有:(p−1)!≡−1(modp)(p-1)!\equiv -1\pmod p(p−1)!≡−1(modp)或写成:(p−1)!≡p−1(modp)(p-1)!\equiv p-1\pmod p(p−1)!≡p−1(modp)同时这是 ppp为质数的充...原创 2019-07-31 11:15:03 · 578 阅读 · 0 评论 -
线性同余方程的求解 及 模板题
定义:ax≡c(modb)ax\equiv c\pmod bax≡c(modb)解的结构:求解ax≡c(modb)ax\equiv c\pmod bax≡c(modb),等价于求解ax+by=cax+by=cax+by=c,设d=gcd(a,b)d=\gcd(a,b)d=gcd(a,b)该方程有整数解的充分必要条件是d∣cd|cd∣c(即ccc被gcd(a,b)\gcd(a,...原创 2019-08-03 20:31:37 · 725 阅读 · 0 评论 -
[牛客竞赛] 处女座的约会(二进制)
链接:https://ac.nowcoder.com/acm/contest/329/I来源:牛客网题目描述处女座放完了&amp;quot;高利贷&amp;quot;,拿到了不少的资金,又可以和小姐姐约会啦!(之前不还是攒钱打比赛的吗)现在处女座拿到了一份宁波市旅游地图决定和小姐姐一起去玩耍。他们来到了动物园,去参观里面的动物。但是很不幸的是,他们在游玩的途中遇到了一只恶龙。恶龙长有n个头,但经过了处女座的调教,恶龙变得...原创 2019-01-27 15:38:05 · 215 阅读 · 0 评论 -
[牛客竞赛] 处女座的测验(一)(质数,构造)
链接:https://ac.nowcoder.com/acm/contest/327/H来源:牛客网题目描述处女座进行了一场c语言的考试,要求很简单,输出2000个正整数,并且满足以下条件:任意两个数互质任意两个数x,y,满足τ(x·y)&amp;amp;amp;amp;gt;10,其中τ(n)为n的因子的个数举例:6的因子有1,2,3,6,所以τ(6)=4输入描述:本题没有输入输出描述:...原创 2019-01-25 13:38:30 · 253 阅读 · 1 评论 -
[牛客竞赛] 求第K大的三角形面积(数学,排序)
链接:https://ac.nowcoder.com/acm/contest/327/A来源:牛客网题目描述平面上有n个点,问:平面上所有三角形面积第k大的三角形的面积是多少?输入描述:第一行T,表示样例的个数。对于每一组样例,第一行两个整数n和k,接下来n行,每行两个整数x,y表示点的坐标T&amp;amp;amp;amp;lt;=803&amp;amp;amp;amp;lt;=n&amp;amp;amp;am原创 2019-01-25 12:33:18 · 435 阅读 · 0 评论 -
PAT Advanced1024 Palindromic Number(高精度运算)
链接:PAT Advanced1024A number that will be the same when it is written forwards or backwards is known as a Palindromic Number. For example, 1234321 is a palindromic number. All single digit numbers ar...原创 2019-01-12 16:07:12 · 186 阅读 · 0 评论 -
PAT Advanced1023 Have Fun with Numbers(高精度运算)
链接:PAT Advanced1023Notice that the number 123456789 is a 9-digit number consisting exactly the numbers from 1 to 9, with no duplication. Double it we will obtain 246913578, which happens to be anoth...原创 2019-01-12 15:00:42 · 370 阅读 · 0 评论 -
PAT Advanced1049 Counting Ones(数学问题)
链接:PAT Advanced1049The task is simple: given any positive integer N, you are supposed to count the total number of 1’s in the decimal form of the integers from 1 to N. For example, given N being 12,...原创 2019-01-03 21:00:39 · 304 阅读 · 0 评论 -
PAT Advanced1059 Prime Factors(建立素数表)
链接:PAT Advanced1059 Given any positive integer N, you are supposed to find all of its prime factors, and write them in the format N = p1k₁ × p2k₂ × ⋯ × pm km .Input Specification:Each in...原创 2018-12-30 17:21:55 · 158 阅读 · 0 评论 -
PAT Basic1034 有理数四则运算
本题要求编写程序,计算 2 个有理数的和、差、积、商。输入格式:输入在一行中按照 a1/b1 a2/b2 的格式给出两个分数形式的有理数,其中分子和分母全是整型范围内的整数,负号只可能出现在分子前,分母不为 0。输出格式:分别在 4 行中按照 有理数1 运算符 有理数2 = 结果 的格式顺序输出 2 个有理数的和、差、积、商。注意输出的每个有理数必须是该有理数的最简形式 k a/...原创 2018-11-28 09:19:52 · 152 阅读 · 0 评论 -
PAT Basic1022 D进制的A+B
链接:PAT Basic1022输入两个非负 10 进制整数 A 和 B (≤230-1),输出 A+B 的 D (1&amp;amp;lt;D≤10)进制数。输入格式:输入在一行中依次给出 3 个整数 A、B 和 D。输出格式:输出 A+B 的 D 进制数。输入样例:123 456 8输出样例:1103就是模拟10进制转其他进制的转换方法不断用sum=sum/D,并且取...原创 2018-11-25 15:49:42 · 165 阅读 · 0 评论 -
PAT Basic1017 A除以B
链接:PAT Basic1017本题要求计算 A/B,其中 A 是不超过 1000 位的正整数,B 是 1 位正整数。你需要输出商数 Q 和余数 R,使得 A=B×Q+R成立。输入格式:输入在一行中依次给出 A 和 B,中间以 1 空格分隔。输出格式:在一行中依次输出 Q 和 R,中间以 1 空格分隔。输入样例:123456789050987654321 7输出样例...原创 2018-11-25 13:53:21 · 144 阅读 · 0 评论 -
[牛客竞赛] 炫酷数字(简单数学,素数)
链接:https://ac.nowcoder.com/acm/contest/331/G来源:牛客网题目描述小希希望你构造一个最小的正整数,使得其有n个因子。输入描述:第一行一个整数T表示数据组数每组数据第一行输入一个正整数n,表示其因子数。n≤1,000,000T≤1,000,000输出描述:输出一行一个整数,表示你构造出的这个数。注意:你需要保证你构造的数≤...原创 2019-02-01 15:44:43 · 555 阅读 · 2 评论 -
[牛客竞赛] 出题(简单数学)
链接:https://ac.nowcoder.com/acm/contest/332/A来源:牛客网题目描述小B准备出模拟赛。她把题目按难度分为四等,分值分别为6,7,8,9。已知小B共出了m道题,共n分。求小B最少出了多少道6分题。输入描述:两个正整数n,m输出描述:一个数,表示答案。若无解,输出"jgzjgzjgz"。示例1输入34 5输出1...原创 2019-02-03 12:47:46 · 400 阅读 · 0 评论 -
欧几里得定理 及 拓展欧几里得 小结
欧几里得定理定义:gcd(a,b)=gcd(b,amod  b)\gcd(a,b)=\gcd(b,a\mod b)gcd(a,b)=gcd(b,amodb)应用:故得到 求解gcd的辗转相除法模板:int gcd(int a,int b){ return b?gcd(b,a%b):a;}int lcm(int...原创 2019-08-03 18:49:18 · 873 阅读 · 0 评论 -
2019多校第四场 HDU6623 Minimal Power of Prime(质因数分解,思维)
链接:HDU6623 Minimal Power of Prime题意:给出T(≤50000)T(\le50000)T(≤50000)个n(≤1018)n(\le10^{18})n(≤1018),对nnn进行质因数分解,问分解后质因数的幂最小的是多少?例如:108=22∗33108=2^2*3^3108=22∗33,所以质因数的幂最小是222分析:直接对nnn进行质因数分解的话时间复...原创 2019-08-03 11:42:26 · 241 阅读 · 0 评论