- 2024-11-2124.11.21
A怎么只有我一个写这种唐诗做法啊/kk当括号匹配时会对若干区间造成贡献。如果我们考虑每个右括号作为右端点统计贡献区间的话,左侧所有和它同一括号范围内(或最外层)的同层的左括号作为左端点和其构成一个贡献区间。举例子来说\(({\color{blue}(}))({\color{yellow}(}){\color{
- 2024-11-19BSGS
给定\(a,b,p\)。求最小非负整数\(x\)使得\(a^x\equivb\pmodp\),或报告无解。保证\((a,p)=1\)。首先根据欧拉定理,\(a^x\equiva^{x\bmod\varphi(p)}\bmodp\)。所以最优的\(x\)一定不大于\(\varphi(p)\)。换一个比较松上限\(p\)。不妨先随便找一个数\(k\)
- 2024-11-15(ex)CRT / (ex)Lucas
以下所有分数默认取整。CRT有\(\begin{cases}x\equiva_1\pmod{p_1}\\x\equiva_2\pmod{p_2}\\\cdots\\x\equiva_n\pmod{p_n}\end{cases}\),且所有\(p\)互质。令\(p'_i=\dfrac{\prodp}{p_i}\),且\(p'_i\)在\(\bmodp_i\)意义下逆元为\(q_i\),则
- 2024-11-12多项式板子
一、数组版本数组版本和poly版本都只涵盖目录中第\(4\sim10\)部分。namespacePoly{intp[maxn],q[maxn],r[maxn],w[maxn];intinum[maxn];intqpow(inta,intk){intres=1;for(;k;a=1ll*a*a%mod,k>>=1)if(k&1)res=1ll*res*a%
- 2024-11-11Pollard-rho & Miller Rabin
Pollard-rho找到\(n\)的一个非平凡因子。暴力发现\(n\)的因子数\(\omega(n)\)实际很少,我们考虑随机一个数,判断是否和\(n\)有公因子,显然很劣。生日悖论:随机\(k\)个值域大小为\(n\)的数,当\(k\ge\sqrtn\)时,\(k\)个数两两不同的几率几乎为\(0\)。以下忽
- 2024-11-0211.02
A.故障机器人天生具备大常熟,劳资就爱写递归用vector写唐怎么你了,复杂度对了凭什么不让过,时间卡这么紧有意思吗?贡献可以拆为识别为↑的字符与识别为→的字符间的贡献,而字符间的贡献又互相独立,所以可以先预处理\(val[x][y]\)代表字符\(x\)识别为↑,字符\(y\)识别为→
- 2024-11-02基础数论
基础数论Update:2024/11/02。素数素数和合数定义若\(p\in\Zeta\),且\(p\not=0,\pm1\),其约数集合中的元素只有\(1\)和\(p\)本身,那么称\(p\)为素数。若\(a\in\Zeta\),且\(a\not=0,\pm1\),\(a\)不为素数,则为合数。素数一般指正的素数。素数计数\(\pi(x)
- 2024-11-02多项式
多项式的表示系数表示法即\(F(x)=a_0x^0+a_1x^1+...+a_nx^n\)点值表示法一个\(n\)次多项式可以被\(n+1\)个点唯一确定可以用这\(n+1\)个点表示该多项式多项式卷积\[(f*g)(x)=\sum_{i=0}^{n}\sum^n_{j=0}a_ib_{j}x^{i+j}\]说白了就是暴力展开快速傅里叶变换(FFT)单位根对
- 2024-10-24多校A层冲刺NOIP2024模拟赛12
多校A层冲刺NOIP2024模拟赛12\(T1\)A.Alice和璀璨花\(65pts/65pts/65pts\)部分分测试点\(1\sim10\):设\(f_{i,j}\)表示前\(i\)位中生长趋势子序列长度为\(j\)时的末尾最小元素,然后进行暴力转移。测试点\(11\sim13\):观察到至多选择\(\left\lceil\log_
- 2024-10-24RSA算法详解及相关数学原理解析
RSA算法详解及相关数学原理解析前言为了记录自己学习密码学的过程,也是为了便于个人应付相关课程的考核,故写此博客。本博客总结了怎么用C++手搓一个RSA算法,以及补补欠缺的一些数学知识和可能欠缺的一些其他算法的实现。参考了其他人的相关博客,用便于我自己理解的话和方式和
- 2024-10-222024 BUPT Programming Contest F
简要题意给定一个\(n\timesn\)矩阵,矩阵中的每一个元素的计算方式如下:随机从\([0,n-1]\)中选出两个整数(可以重复)\(a,b\),矩阵的元素为\(a\timesb\bmodn\)求矩阵中元素\(m\)出现次数的期望。\(0\lem<n\le10^{12}\)题解对于矩阵中的任意一个元素是独
- 2024-10-2124.10.19
A数学题,不会。随便取一数\(v\),询问得到\(t\equiv\log_gv\pmodp\)。我们希望找到\(x\)使得\(v^x\equivg\pmodp\),即\(g^{tx}\equivg\pmodp\Leftrightarrowtx\equiv1\pmod{p-1}\)。那么只要\(t\)与\(p-1\)互质即可求得逆元。有原根相关知识可以知
- 2024-10-21[ARC185A] mod M Game 2 Solution
ARC185A-modMGame2简要题意Alice和Bob玩卡牌游戏。每个人都有一副\(N\)张卡牌,分别标上数字\(1\simN\)。现从Alice开始,两人轮流出牌放入牌堆,每人每局恰好出一张牌,出过的牌不能再出;如果在某一时刻,牌堆里所有牌的数字总和是\(M(N<M)\)的倍数,则刚刚出牌的玩家输,
- 2024-10-192024.10.19总结
本文于github博客同步更新。A:考虑随便取一个数\(v\),用一次询问问出\(t=\log_gv\)。我们希望找到一个\(x\)使得\(v^x\equivg\pmodp\),也即\(g^{tx}\equivg\pmodp\ifftx\equiv1\pmod{p-1}\)。于是,我们希望找到的\(v\)使得\(t\)与\(p-1\)互质即可。由原根的
- 2024-10-1910.19
别样的\(\text{NOI}\)模拟赛。\(A\)十几分钟能写完的随机化都放过去了,\(B\)题面的代码\(CE\)了,\(C\)边分治的思路仅闪过一瞬就忘了。A.离散猜数你说得对,但是若答案正确,且你的代码使用的询问次数为\(x\),std使用的询问次数为\(y\),计算\(c=\dfrac{x}{y}\)。若\(c\l
- 2024-10-17多项式全家桶
每次复习完下一次都会忘,这次下定决心一定要记下来!!!FFT和NTT板子,直接拿过来(包括了其他的定义):intn,m,rev[maxn];intqpow(intx,intk,intp){ intres=1; while(k){ if(k&1) res=res*x%p; k>>=1,x=x*x%p; } returnres;}voidprepare(
- 2024-10-16Montgomery Curves and Weierstrass Curves
WeierstrassCurvesWeierstrassCurves形如\(y^2=x^3+Ax+B\)其中\(4A^3+27B^2≠0\),这种形式称为WeierstrassForm。WeierstrassCurves上的运算在椭圆曲线(此处即为WeierstrassCurves)上,可以定义点之间的加法运算,其满足:单位元\(O\)为无穷远点对于曲线上的两点\(P\)和\(Q\)
- 2024-10-1010.10 爆零赛(2023 炼石计划 NOIP 模拟赛 #2)
炼石计划9月10日NOIP模拟赛#2【补题】-比赛-梦熊联盟(mna.wang)模拟赛恒等式:\(0+0+0+0=0\)。复盘T1好像可做。有个显然的\(n^2\)DP。推式子的时候猜到了\(\gcd=1\)的做法。进一步尝试正解未果。T2一眼只会爆搜。想到了\(b\timesb!\)的做法,应该能过\(
- 2024-10-04国庆 CF 做题记录
CF2002F2考虑F1,当\(n=m\)时,我们默认\(l\gef\)。此时我们可以发现一个比较正确的策略:先从\((0,0)\)跳到满足\(p\)是质数的\((p,0)\)处,然后再跳到满足\(q\)是小于\(p\)的质数的\((p,q)\)处,然后再暴力BFS。不会证明,可以达标找出这样的结论。当\(n>m\)时,注
- 2024-09-25从CF1920C看同余的一个性质
https://codeforces.com/problemset/problem/1920/C同余的一个性质:证明很显然,但是想不到这个性质题意给你\(n\)个数,划分\(k\)段,每段在对\(m(m\ge2)\)取模之后相等即为一个合法方案,问有多少个合法方案。断点//check是能O(n)的//问题在于怎么check//经验证,m=2
- 2024-09-24【信息安全数学基础】二次剩余(Quadratic residue)
什么是二次剩余呢?小小定义设m是大于1的整数,a是与m互素的整数,若x2≡a
- 2024-09-21乘法逆元
逆元定义如果\(ax\equiv1\pmodp\),则称\(x\)为\(a\)在模\(p\)意义下的乘法逆元。逆元存在当且仅当\(a\perpp\),即\(\gcd(a,p)=1\)。将\(ax\equiv1\pmodp\)转化可得\(x\equiv\dfrac{1}{a}\pmodp\),那么模意义下\(t\diva\)就相当于\(t\timesx\)。快速求单
- 2024-09-20数论学习笔记(一)(2024.7.25)
一、最大公约数定义不全为\(0\)的整数\(a,b\)的最大公约数是指能够同时整除\(a\)和\(b\)的最大整数。欧几里得算法(gcd)gcd是用来求解两个整数的最大公约数定理1.2.1对于整数\(a,b,m,n\),若\(c\mida,c\midb\),则\(c\mid(ma+nb)\)证:\(\becausec\mida
- 2024-09-18yzh 的石子
题意简述yzh在一个特殊的日子里,送给你了\(n\)颗具有魔法的石子!你当然知道这是她为你准备的惊喜。她偷偷告诉你,为了使它们的魔法被释放出来,只有把这些石子分成若干堆,并且每一堆石子个数都能表示成\(3x^2-3x+1\)的形式,其中yzh希望让\(x>1\)。以你的聪明才智,不久就想到了