数论

网络安全是一门涉及计算机、网络、通讯、密码、信息安全、应用数学、数论、信息论等多种学科的综合性学科，涉及面极广，而且不断更新和发展... 13-03-20

　　从根本上来说，性交实际上是一个性能量的“释放”过程。性交时间短、质量差的时候，身体积聚的性冲动能量并没有得到完全“释放”，身体就有将“剩余”的能量排出的需求。所以，此时的不应期是相对较短的;而在性

　　大多数男性认为在一次性生活中重复性交次数多就是性能力强的表现。所以在一些纯男性聚会的场合，经常可以听到某些男人吹嘘一晚能进行多少多少次，以满足自己的虚荣心。一些性知识网站也为此介绍一些能够“二次攻

题意 墨墨突然对等式很感兴趣，他正在研究a1x1+a2y2+…+anxn=B存在非负整数解的条件，他要求你编写一个程序，给定N、{an}、以及B的取值范围，求出有多少B可以使等式存在非负整数解。 Sol maya神仙题啊，感觉自己做题难度跨度太大了qwq。 这里有一篇讲的非常好的博客https:// ...

Greatest Common Divisor(GCD) 欧几里得算法据说是最早的算法，用于计算最大公约数，也是数论的基础算法之一。 1.欧几里德算法的思想： 欧几里德算法的思想基于辗转相除法的原理，辗转相除法是欧几里德算法的核心思想，欧几里德算法说白了其实就是辗转相除法的计算机算法的实现而已。下面 ...

题意 "题目链接" 求满足$i^2 + j^2 \% M = 0$的数对$(i, j)$的个数，$1 \leqslant i, j \leqslant 10^9, M \leqslant 1000$ Sol 发这篇博客的目的就是为了证明一下我到底有多菜。 mdzz小学组水题我想了40min都没想出来 ...

题意 "题目链接" Sol 啊啊这题好恶心啊，推的时候一堆细节qwq $a \% i = a \frac{a}{i} i$ 把所有的都展开，直接分块。关键是那个$i \not= j$的地方需要减。。。。 然后就慢慢写就好了 cpp include define Pair pair define MP ...

　　性爱出了问题后果可大可小，小的话顶多也就是不尽如人意，但如果大了，那么势必会严重降低夫妻感情和生活指数。因此，如果我们在性生活中发现问题，就应该及时解决。

题意 "题目链接" Sol 这题是来搞笑的吧。。 考虑一个数的贡献是$O(\frac{N}{i})$ 直接数论分块。 cpp include define Pair pair define MP(x, y) make_pair(x, y) define fi first define se seco ...

在我们需要判断一个数是否是素数的时候，最容易想到的就是那个熟悉的O(√n)的算法。那个算法非常的简单易懂，但如果我们仔细想想，当n这个数字很大的时候，这个算法其实是不够用的，时间复杂度会相对比较高。 怎么解决呢？我们先来了解一下“费马小定理”。假设我们有一个素数p，且另一个数a和p互素，就可以得到a ...

主要讨论三个比较常见的博弈游戏 Bash Game,Nim Game和Wythoff Game，较为领人惊叹的是，他们最后都是通过数论或者自然数性质完美解决： Bash Game:同余理论 Nim Game:异或理论 Wythoff Game:黄金分割 （1）Bash Game：一堆n个物品，两人轮 ...

又是几天过去，最近在刷数论基础和dp基础。kuangbin数论专题刷了差不多一大半吧。深深感jio到自己的菜。唉，温故而知新，所以决定在踩一遍坑点。 Bi-shoe and Phi-shoe https://vjudge.net/problem/LightOJ-1370 题目大意：给出一些数字，对于 ...

[ bzoj2820] YY的GCD Time Limit : 3000 ms Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种傻×必然不会了，于是向你来请教…… 多组输入 Inp ...

来自Perm排列计数的悲伤 lucas说过 C(n,m)%p=C(n%p,m%p)*C(n/p,m/p)%p 于是乎 在好记的情况下没有搞原理证明 导致PermWA了一上午 1 ll pow(ll a,ll b){ 2 ll ans=1; 3 while(b){ 4 if(b&1)ans=(ans* ...

费马小定理新手入门\+总结 纵有疾风起，人生不言弃。 前言 最近新手的我做了几个和快速幂有关的题目，发现他们还经常和费马小定理联系在一起，所以有必要写一篇文章来总结一下费马小定理，以便后面更好的学习。 内容介绍 费马小定理是数论中的一个重要定理，再1636年提出。 ​核心：如果p是一个质数，并且整数 ...

wdnmd什么垃圾网站，把我两道题分搞没了，实际得分50/300，应该得分150/300； 无话可说； T1暴力枚举； T2数论+DP； T3几何等； T1先枚举行的情况，再枚举列的情况； #include<iostream> #include<cstdio> #include<algorithm> ...

[TOC] 题目 "P1965 转圈游戏" 思路 数论。题目就是让你判断一下第$x$个人走了$m \times 10^k$步之后到了哪里。总共走的步数为$m \times 10^k$直接走肯定会 ，对$n$取模之后再走就好了。 $Code$ cpp include include include i ...

T1 一道数论神题 题目 【题目描述】 LYK有一张无向图G={V,E}，这张无向图有n个点m条边组成。并且这是一张带权图，只有点权。 LYK想把这个图删干净，它的方法是这样的。每次选择一个点，将它删掉，但删这个点是需要代价的。 假设与这个点相连的还没被删掉的点是u1,u2,…,uk。LYK将会增加 ...

写在前面 这是蒟蒻第一次写这么长的博文 $gyh\ nb$，$\text{OI Wiki}\ nb$ 如果觉得写得凑合就点个支持吧$qwq$ 前置知识 "积性函数" 、 "狄利克雷卷积" 、 "数论分块（这一篇去找gyh吧我讲也讲不好）" ~~(有空慢慢补)~~ Mobius函数 定义 莫比乌斯函数 ...

思路：先按照正常构造自然数：S=2×3×…×(n−Δx)×(n−Δx+1)×…×nS=2\times 3\times \ldots \times \left( n-\Delta x\right) \times \left( n-\Delta x+1\right) \times \ldots \times nS=2×3×…×(n−Δx)×(n−Δx+1)×…×n一般会多出一个数Δx\Delta xΔx (Δx≤n\Delta x\leq nΔx≤n)S=2×3×…×(n−Δx)×(n−Δx+2)×….

线性求逆元快速版#include<bits/stdc++.h>#define MOD mod#define MAX maxnusing namespace std;typedef pair<int, int> pii;typedef long long ll;const double eps = 1e-6;const int maxn = 2e3 + 10, inf = 0x3f3f3f3f, mod = 1e9 + 7;ll quickpow(ll x, ll k

题目题目链接：https://www.luogu.com.cn/problem/P1516思路好麻烦啊不想写呜呜呜数论怎么这么麻烦代码#include<cstdio>#include<cstring>#include<cmath>#include<cstdlib>#include<cctype>#include<ctime>#include<iostream>#include<string&g

依旧是这个生成函数。。怎么这么多水题啊。。#include<iostream>#include<cstring>#include<algorithm>#include<vector>using namespace std;typedef long long ll;#define ls (o<<1)#define rs (o<<1|1)#define pb push_backconst int M = 40+7;

依旧是最基础的生成函数。(1+x^i)*……直接用lst每次上界最大值优化就够了。（因为最坏情况是尽量跑满，所以没必要记录已有的数值）#include<iostream>#include<cstring>#include<algorithm>#include<vector>using namespace std;typedef long long ll;#define ls (o<<1)#define rs (o<

还是水题。。这题每个种类数目起点不为0；#include<iostream>#include<cstring>#include<algorithm>#include<vector>using namespace std;typedef long long ll;#define ls (o<<1)#define rs (o<<1|1)#define pb push_backconst int M = 110+

ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a;}ll quick_mult(ll a, ll b, ll mod) { ll ans = 0; while(b) { if(b & 1) ans = (ans + a) % mod; a = (a + a) % mod; b >>= 1; } return ans;}ll quick_po

还是经典生成函数：(1+x^i)#include<iostream>#include<cstring>#include<algorithm>#include<vector>using namespace std;typedef long long ll;#define ls (o<<1)#define rs (o<<1|1)#define pb push_backconst int M = 300+7;/

(1+X^i)*……最经典的生成函数。#include<iostream>#include<cstring>#include<algorithm>#include<vector>using namespace std;typedef long long ll;#define ls (o<<1)#define rs (o<<1|1)#define pb push_backconst int M = 1e4+7

这些题都可以暴力直接做，但用于练习生成函数。（虽然并没有。。。马上做几道难题练练）#include<iostream>#include<cstring>#include<algorithm>#include<vector>using namespace std;typedef long long ll;#define ls (o<<1)#define rs (o<<1|1)#define pb push_back

Java实现求解两个正整数的最大公约数和最小公倍数公约数 ：公因数，显示着若干个整数之间的数论关系。如果一个数同时是几个数的约数，称这个数为它们的“公因数”；公约数中最大一个的称为最大公因数。公倍数 ：公倍数，显示着若干个整数之间的数论关系。如果一个数同时是几个数的倍数，称这个数为它们的“公倍数”；公倍数中最小一个的称为最小公倍数。普通思路 ：①求解最大公约数：显然，最大公约数小于等于两正整数中较小的一个，同时满足公约(m % i == 0 && n % i == 0)，故使用fo