------数论------
文章平均质量分 95
Cwolf9
莯
展开
-
2020_CCPC长春L题 Coordinate Paper Exgcd贪心
2020 CCPC Changchun gym102832L题意链接:点我点我你需要构造出一个长度为nnn的序列满足以下三个条件:ai≥0a_i\ge 0ai≥0(∑i=1nai)=s(\sum_{i=1}^na_i)=s(∑i=1nai)=s对于任意i≥2, ai=ai−1+1i\ge 2,\; a_i=a_{i-1}+1i≥2,ai=ai−1+1或者ai=ai−1−ka_i=a_{i-1}-kai=ai−1−k1≤n,k≤1e5,1≤s≤1e181\le n,k\le原创 2020-11-12 21:09:04 · 2656 阅读 · 0 评论 -
记一些数学引理和结论等
一些结论gcd(xa−1,xb−1)=xgcd(a,b)−1gcd(x^a-1,x^b-1)=x^{gcd(a,b)}-1gcd(xa−1,xb−1)=xgcd(a,b)−1gcd(fib[x],fib[y])=fib[gcd(x,y)]gcd(fib[x],fib[y])=fib[gcd(x,y)]gcd(fib[x],fib[y])=fib[gcd(x,y)]费马大定理定理:an+bn...原创 2019-03-01 12:55:48 · 487 阅读 · 0 评论 -
杜教筛学习日记
文章目录前置知识杜教筛上:莫比乌斯函数求和杜教筛上:欧拉函数求和杜教筛上:小结杜教筛下:前置:迪利克雷卷积杜教筛下:求$\sum_{i=1}^{n}\phi(i)\times i$杜教筛下:一般方法P4213 【模板】杜教筛(Sum)求$\sum_{i=1}^n\mu*\mu(i)$的值题目参考博客:jk_chen_acmer前置知识莫比乌斯反演学习日记∑d∣nμ(d)=[n=1]\su...原创 2020-07-16 18:50:26 · 851 阅读 · 0 评论 -
牛客练习赛40 E题 小D的Lemon 莫比乌斯反演
题目描述あの日の悲しみさえ那一天的悲伤也好あの日の苦しみさえ那一天的痛苦也好そのすべてを愛してた あなたとともに我深爱着和你在一起的点点滴滴胸に残り離れない残留在心中久久不离苦いレモンの匂い苦柠檬的香气雨が降り止むまでは帰れない雨停为止都无法回去切り分けた果実の片方の様に如同切开的果实的一面那样今でもあなたはわたしの光至今为止你依然是我的光米津玄师—《Lemon》...原创 2019-02-16 11:37:00 · 428 阅读 · 0 评论 -
莫比乌斯反演学习日记
∑d∣nu(d)=[n=1]\sum_{d|n} u(d) = [n = 1]∑d∣nu(d)=[n=1]原创 2019-02-10 15:53:06 · 682 阅读 · 4 评论 -
Wannafly Winter Camp Day8(Div1, onsite) I题 岸边露伴的人生经验
Wannafly Winter Camp Day8(Div1, onsite)I题 岸边露伴的人生经验有点无聊,写点博客打发时间吧。就按照题解所说的,把它转化成一个20位的二进制数,然后套个fwt的板子,算一下每两位的贡献,注意一下四元组是和顺序有关系的,本题就没有了。AC_code#include<bits/stdc++.h>#define clr(a, b) mem...原创 2019-02-02 10:37:42 · 322 阅读 · 0 评论 -
有关取模,同余,逆元等的一些东西
刚看到一些东西,先记下来a^-1表示a的逆元a^2 ≡ b^2 (mod p), so a ≡ b (mod p)a ≡ b (mod p), so b ≡ a (mod p)a ≡ b (mod p), so 1/b ≡ a^-1 (mod p)a ≡ b (mod p), so a/c ≡ b*c^-1 (mod p) a ≡ b (mod p), so -a ≡ mod-b (m...原创 2019-01-15 20:03:58 · 1538 阅读 · 2 评论 -
SPOJ LCMSUM (BZOJ 2226) 求n以内所有数与n的lcm之和(保证所有人都能看懂)
SPOJ LCMSUM BZOJ 2226#include<bits/stdc++.h>#include <cstdio>#include <cmath>#include &原创 2019-01-20 22:43:31 · 458 阅读 · 2 评论 -
HDU6440 Dream 2018CCPC网络赛-费马小定理
(有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 CatalogCatalogProblem:Portal传送门Solution:AC_Code:Problem Description:Problem:Portal传送门 原题目描述在最下面。 给定一个素数p,要求定义一个加法运算表和乘法运算表,使的(m+n)p=mp...原创 2018-08-25 22:16:16 · 520 阅读 · 0 评论 -
FZU2278 YYS: Java+概率+数论-第八届福建省大学生程序设计竞赛
(有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 CatalogCatalogProblem:Portal传送门Solution:AC_Code:Problem Description:Problem:Portal传送门 原题目描述在最下面。 (n-1)!天可以抽一次卡,共有n张卡,问抽到n张卡的期望天数是多少。...原创 2018-08-20 19:58:43 · 369 阅读 · 0 评论 -
Gym 101466K: Random Numbers-线段树+DFS序
(有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 CatalogCatalogProblem:Portal传送门Solution:AC_Code:Problem Description:Problem:Portal传送门 类似:HDU3974 原题目描述在最下面。 给你一颗树,每个节点有权值。每次询问一个子树...原创 2018-08-19 21:01:56 · 461 阅读 · 0 评论 -
Gym 101652P:Fear Factoring 数论 除法分块 or 素数筛
(有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 CatalogCatalogProblem:PortalSolution:AC_Code:Problem Description:Problem:Portal PDF链接:here 原题目描述在最下面。 函数f(X)f(X)f(X)表示XXX所有约数的和。例如...原创 2018-08-18 12:20:06 · 870 阅读 · 0 评论 -
Codeforces837F-Prefix Sums-二分加组合数(爆longlong了)
(有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 题意:传送门 原题目描述在最下面。 对于一段序列ai,每次操作后ai'等于前缀和ai。 问多少次操作后序列最大值大于aim。 之前在牛客上做了一道类似的简单题,牛客的题是输出k次变化后的矩阵。牛客简单在答案输出mod。本题不能取模,然后爆longlong了???不想用大整数,懒得敲java。没想到居然用lo...原创 2018-07-21 12:58:54 · 233 阅读 · 0 评论 -
牛客练习赛17-C操作数-组合数
(有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 目录目录题目:传送门思路:AC代码:题目:传送门给定长度为n的数组a,定义一次操作为: 1. 算出长度为n的数组s,使得si= (a[1] + a[2] + … + a[i]) mod 1,000,000,007; 2. 执行a = s; 现在问k次操作以后a长什...原创 2018-08-15 23:11:18 · 185 阅读 · 0 评论 -
扩展BSGS-传送门
很好的讲解:ZigZagK 好的讲解:mjtcn 某个模板:here 模板题: BSGS:ZigZagK的poj2417exBSGS:ZigZagK的poj3243 — AC_Gibson一般的板子过不了这个题!GYM101853G:here exLucas:here...原创 2018-08-15 16:47:47 · 346 阅读 · 0 评论 -
HDU3037-卢卡斯(Lucas)定理模板 组合数模板
前言: 写上一个题的时候用到组合数了,就去百度了下,在此更新一下。 typedef long long LL;const LL md = 1e9 + 7;const int maxn = 1e5 + 10;LL ar[maxn],sum[maxn];const int N = 100007;LL f[N];LL quickpow(LL a,LL b,LL m){ LL...原创 2018-07-21 13:06:01 · 357 阅读 · 0 评论 -
求欧拉函数
/*3.欧拉函数对正整数n,欧拉函数是少于或等于n的数中与n互质的数的数目显然对素数n,phi(n)=n-1通式:φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数,x是不为0的整数。φ(1)=1(唯一和1互质的数就是1本身)很简单,就不写代码什么的了*/降幂公式当x≥ϕ(p)时,有ax≡axm...原创 2018-07-23 16:29:42 · 527 阅读 · 0 评论 -
唯一质因数分解定理
唯一质因数分解定理: 任意一个合数a仅能以一种方式,写成如下的乘积形式: a = p1^e1*p2^e2*…*pr^erconst int N = (int)2e5 + 7;int noprime[N], pcnt, p[N / 2];int nump[N / 2], yinzi[N / 2];int n, m, top;void getprime(){ pcnt = 0...原创 2018-07-23 15:22:41 · 1268 阅读 · 0 评论