数论
文章平均质量分 92
Sqwlly
有自己热爱的东西,真好。
展开
-
HDU 4549(矩阵快速幂+费马小定理)
M斐波那契数列Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others)Problem DescriptionM斐波那契数列F[n]是一种整数数列,它的定义如下:F[0] = a F[1] = b F[n] = F[n-1] * F[n-2] ( n > 1 )...原创 2018-03-10 19:17:26 · 277 阅读 · 0 评论 -
LightOJ 1259 Goldbach`s Conjecture(素数筛)
Goldbach`s Conjecture题解:素数打表+暴力枚举。#include<bits/stdc++.h>using namespace std;const int N = 1e7+10, tot = N/log(N)+ 100000;int prime[tot];bool vis[N];int init(){ int k = 0; for(int i ...原创 2018-12-03 22:08:03 · 172 阅读 · 0 评论 -
LightOJ - 1197 Help Hanzo(区间素数筛)
Help Hanzo题解:因为区间太大,所以我们需要考虑再仔细看一下题,求的是[a,b][a,b][a,b]区间内的素数,并且b−a≤100000b - a ≤ 100000b−a≤100000,因此我们可以借用素数筛的思想,将[a,b][a,b][a,b]区间内的合数筛去。具体做法就是先找到区间[a,b]内的第一个合数int t = prime[i] * ceil(1.0 * a / prim...原创 2018-12-03 22:19:19 · 215 阅读 · 0 评论 -
LightOJ 1236 Pairs Forming LCM(唯一分解定理+素数筛)
Pairs Forming LCM题解:考虑唯一分解定理:a=p1n1⋅p2n2...pknka = {p_1}^{n_1}\cdot {p_2}^{n_2}...{p_k}^{n_k}a=p1n1⋅p2n2...pknk,b=p1m1⋅p2m2...pkmkb = {p_1}^{m_1}\cdot {p_2}^{m_2}...{p_k}^{m_k}b=p1m1⋅p2m2.....原创 2018-12-03 23:18:40 · 178 阅读 · 0 评论 -
LightOJ 1336 - Sigma Function(算术基本定理)
Sigma Function题解:首先,分析题目公式(1)σ(n)=p1e1+1−1p2−1⋅p2e2+1−1p2−1⋅......⋅pkek+1−1pk−1\sigma(n) = \frac{p_1^{e_1+1}-1}{p_2-1}\cdot \frac{p_2^{e_2+1}-1}{p_2-1}\cdot......\cdot\frac{p_k^{e_k+1}-1}{p_k-1} \tag...原创 2018-12-10 23:28:13 · 216 阅读 · 0 评论 -
POJ 1061 青蛙的约会(扩展欧几里德)
青蛙的约会题解:显然可以列出方程(x+k⋅m)%L=(y+k⋅n)%L(x+k\cdot m)\%L=(y+k\cdot n)\%L(x+k⋅m)%L=(y+k⋅n)%L然后可以转换为(m−n)⋅a+L⋅b=y−x(m-n)\cdot a+L\cdot b=y-x(m−n)⋅a+L⋅b=y−x,最后求最小正整数解。exgcdexgcdexgcd了呀~代码#include<iostre...原创 2018-12-10 23:37:39 · 187 阅读 · 0 评论 -
POJ 2142 The Balance(扩展欧几里德)
The Balance由题意可得:{ax1−by1=dbx2−ay2=d\begin{cases}ax_1 - by_1 = d \\\\ bx_2 - ay_2 = d\end{cases}⎩⎪⎨⎪⎧ax1−by1=dbx2−ay2=d用exgcdexgcdexgcd解出两个最小解,然后输出最小的那一对即可。注意如果x2+y2x_2 + y_2x2+y2更小要先输出y2y_...原创 2018-12-11 01:11:05 · 196 阅读 · 0 评论 -
HDOJ 2669 Romantic(扩展欧里几得)
Romantic题解:太裸了。。求exgcdexgcdexgcd最小正整数解代码#include<bits/stdc++.h>using namespace std;typedef long long LL;void exgcd(LL a,LL b,LL &x,LL &y, LL &d) { if(b == 0) { x = 1, y =...原创 2018-12-11 01:20:41 · 198 阅读 · 0 评论 -
LightOJ 1341 - Aladdin and the Flying Carpet(算术基本定理唯一分解定理)
Aladdin and the Flying Carpet题解:首先必须知道一个数n=p1e1⋅p2e2⋅p3e3⋅...pkekn = p_1^{e_1}\cdot p_2^{e_2}\cdot p_3^{e_3}\cdot ...p_k^{e_k}n=p1e1⋅p2e2⋅p3e3⋅...pkek的因子个数ϕ(n)=∏i=1k(ei+1)\phi(n) = \prod_{i =...原创 2018-12-11 01:31:43 · 220 阅读 · 0 评论 -
HDU 5685 (前缀+逆元)
Problem A题意:给出你哈希值的计算方式,然后多次询问子串的哈希值。题解:我们通过观察哈希值的计算式子就可以发现是连乘,又是多次询问,因此我们可以想到打表的方式。前缀积即可。ans[a,b]=dp[b]dp[a−1]ans[a,b] = \frac{dp[b]}{dp[a - 1]}ans[a,b]=dp[a−1]dp[b]然后要注意到取模,所以需要乘法逆元。用扩欧或者费马小定理都可...原创 2018-12-22 16:38:08 · 330 阅读 · 0 评论 -
CCPC-Wannafly Winter Camp Day3 div2 F. 小清新数论* 莫比乌斯反演
小清新数论心情:蒻蒻的第一道莫比乌斯反演!!看了好几个小时QAQ,终于看懂些了!开心!^_^题解:(1)∑i=1n∑j=1nμ(gcd(i,j)) \sum_{i = 1}^n\sum_{j = 1}^n \mu(gcd(i,j)) \tag 1i=1∑nj=1∑nμ(gcd(i,j))(1)(2)∑d=1n∑i=1n∑j=1nμ(d)[gcd(i,j)==d]\sum_{d = 1}^...原创 2019-02-02 22:57:33 · 398 阅读 · 0 评论 -
Codeforces Round #552 (Div. 3) G. Minimum Possible LCM(埃氏筛法枚举GCD)
G. Minimum Possible LCM题意:求nnn个数中最小公倍数数值最小的两个数的下标。题解:参考于https://blog.csdn.net/qq_41157137/article/details/89353527,因为LCM(x,y)=x×ygcd(x,y)LCM(x,y) = \frac{x\times y}{gcd(x,y)}LCM(x,y)=gcd(x,y)x×y对于...原创 2019-04-17 23:50:33 · 264 阅读 · 0 评论 -
CodeForces - 476C-Dreamoon and Sums(数学思维)
C. Dreamoon and Sums题解:设(1)x=q⋅b+rx = q\cdot b + r \tag{1}x=q⋅b+r(1)则有(2)div(x,b)=q,&amp;MediumSpace;mod(x,b)=r∈[1,b−1]div(x,b) = q ,\: mod(x,b) = r\in[1,b-1]\tag{2}div(x,b)=q,mod(x,b)=r∈[1,b−1](2)...原创 2018-11-07 00:48:15 · 202 阅读 · 0 评论 -
牛客国庆集训派对Day5 L数论之神(CCPC_CAMP)除法分块
数论之神题解:除法分块。如果是k≤nk\leq \sqrt nk≤n的话第二个数值就是⌊ak⌋\lfloor \frac{a}{k} \rfloor⌊ka⌋,否则就是n+1\sqrt n+1n+1对应的⌊ni⌋\lfloor \frac{n}{i} \rfloor⌊in⌋减掉kkk超过n\sqrt nn的部分然后在加上111。代码#include&lt;bits/stdc++.h...原创 2018-10-08 20:47:03 · 193 阅读 · 0 评论 -
Wannafly挑战赛11 B 白兔的式子(卢卡斯定理+费马小定理求逆元)
B 白兔的式子时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524288K 64bit IO Format: %lld题目描述已知f[1][1]=1,f[i][j]=a*f[i-1][j]+b*f[i-1][j-1](i >= 2,1<=j<=i)。 对于其他情况f[i][j]=0有T组询问,每次给出a,b,n,...原创 2018-03-11 17:52:51 · 376 阅读 · 0 评论 -
BAPC 2014 Preliminary A
Choosing Ice Cream问题分析求一个最小的ii\,i\,满足kimodn=0kimodn=0k^{i}\:mod\:n = 0 因为数据最大只有1e91e9\,1e9\,,所以我们用快速幂,每次最多枚举29次就够了。#include &lt;cstdio&gt;#include &lt;iostream&gt;#include &lt;algorithm&gt;...原创 2018-07-07 22:33:24 · 309 阅读 · 0 评论 -
Benelux Algorithm Programming Contest 2014 Final I Interesting Integers(扩展欧里几得)
Interesting Integers问题分析我们稍微列一下式子就可以发现 aF(k−1)+bF(k−2)=naF(k−1)+bF(k−2)=na F(k-1) + bF(k-2) = n 因此,最笨的办法就是直接枚举a和b,这里nnn最大只有1e91e91e9,所以可以刚好卡过,但是不要傻傻的枚举,需要优化一下式子,去掉一层forforfor#include&lt;iost...原创 2018-07-15 16:24:09 · 165 阅读 · 0 评论 -
关于欧拉定理的一些知识
欧拉定理aφ(n)≡1(modn),gcd(a,n)=1aφ(n)≡1(modn),gcd(a,n)=1a^{\varphi(n)}\equiv 1(mod\:n),\:gcd(a,n)=1 对于正整数nnn,代表小于等于nnn的与nnn互质的数的个数,记作φ(n)φ(n)\varphi(n)。 比如φ(6)=2φ(6)=2\varphi(6)=2,因为与666互质并且小于等于666的正整...原创 2018-08-07 12:59:05 · 793 阅读 · 0 评论 -
牛客多校第八场 G-Counting regions(组合数学+几何欧拉公式)
Counting regions这道题呢,其实是一道组合数学的题,再加上欧拉公式。题意:问边数为奇数的正多边形所有的顶点两两相连,直线所分割的区域数量为多少,答案对1e9+71e9+71e9+7取模。我们会发现题目中的图形,除了顶点是不会有三线交于一点的。所以分割区域的区域数量就只取决于点的数量。现在我们未知: 1.1.1.\:一共有多少条直线? 2.2.2.\:这些直线在正多...原创 2018-08-15 16:14:06 · 486 阅读 · 0 评论 -
牛客多校第四场 A-Ternary String(欧拉降幂)
Ternary String问题分析题意:问多长时间后可以将整个序列完全消除,每秒序列会在111后面插入一个000,在222后面插入一个111,作为序列的第一个数字将会消失。最后答案对1e9+71e9+71e9+7取模。 emmm,通过打表就可以发现 假设现有一个长度为nnn的序列, 每次序列插入一个000后,我们需要n+1n+1n+1秒才能将这个序列消除; 每次序列插入一个1...原创 2018-08-08 22:27:48 · 320 阅读 · 8 评论 -
2018ccpc网络预选赛1010(hdu6441) 费马大定理
Find Integer问题分析题意:已知an+bn=cnan+bn=cna^n+b^n=c^n,给出nnn和aaa,求bbb,ccc,如果无解输出−1−1-1。 费马大定理 1. an+bn=cnan+bn=cna^n+b^n=c^n,n>2n>2n>2时无解。 2. 当aaa为奇数时, a=2⋅k+1a=2⋅k+1a = 2 \cdot k + 1c=k2+...原创 2018-08-26 15:15:25 · 289 阅读 · 0 评论 -
bzoj 2299 向量(裴蜀定理)
向量描述给你一对数a,b,你可以任意使用(a,b),(a,−b),(−a,b),(−a,−b),(b,a),(b,−a),(−b,a),(−b,−a)(a,b),(a,−b),(−a,b),(−a,−b),(b,a),(b,−a),(−b,a),(−b,−a)(a,b), (a,-b), (-a,b), (-a,-b), (b,a), (b,-a), (-b,a), (-b,-a)这些...原创 2018-08-28 13:56:45 · 178 阅读 · 0 评论 -
bzoj 2257 瓶子和燃料(裴蜀定理)
瓶子和燃料描述jyy就一直想着尽快回地球,可惜他飞船的燃料不够了。 有一天他又去向火星人要燃料,这次火星人答应了,要jyy用飞船上的瓶子来换。jyy 的飞船上共有 N个瓶子(1<=N<=1000) ,经过协商,火星人只要其中的K 个 。 jyy 将 K个瓶子交给火星人之后,火星人用它们装一些燃料给 jyy。所有的瓶子都没有刻度,只 在瓶口标注了容量,第i个瓶子的容...原创 2018-08-28 15:04:46 · 205 阅读 · 0 评论 -
洛谷P3951 小凯的疑惑(找规律|扩展欧里几得)
P3951 小凯的疑惑题解:扩展殴里几得的做法请点我蒻的做法。。。首先,根据数据范围可以肯定是有规律的题。(遇事先暴力打个表,既增加对题目的理解,有时还能直接ac,比如近日的南京网赛A)发现3⋅7−(3+7)=113\cdot7 - (3+7) = 113⋅7−(3+7)=11,然后又测了几组数据发现同样满足a∗b−(a+b)a*b-(a+b)a∗b−(a+b),然后就莽夫般的交了,意外的...原创 2018-09-17 22:17:04 · 452 阅读 · 0 评论 -
牛客国庆集训派对Day4 I 连通块计数(排列组合)
连通块计数题意:有nnn条链和一个点作为中心组成的树,求非空连通子树数目。题解:包含中心的连通子树数目:∏i=1n(ai+1)\prod_{i=1}^n(a_i+1)∏i=1n(ai+1)(这里画个图就知道了)不包含中心的连通子树数目:∑i=1nC(ai+1,2)\sum_{i=1}^nC(a_i+1,2)∑i=1nC(ai+1,2)代码#include<bits/std...原创 2018-10-04 22:06:17 · 206 阅读 · 0 评论 -
The Preliminary Contest for ICPC Asia Nanjing 2019 B super_log(欧拉降幂)
super_log题意求aaa...a^{a^a} ...aaa...,一共有bbb个aaa。题解由扩展欧拉定理可得ab(modc)={ab%ϕ(c)+ϕ(c)b ≥ ϕ(c) abb < ϕ(c)a^b\pmod c=\begin{cases}a^{b\%\phi(c)+\phi(c)}& \text{b...原创 2019-09-02 11:32:40 · 283 阅读 · 0 评论