蓝书数论
zsyz_lb2003
这个作者很懒,什么都没留下…
展开
-
球形空间产生器[BZOJ1013]
题面描述传送门思路因为一个球体的球面上的所有点到圆心的距离相等,因此只需求出一个点(x1,x2,x3,......,xnx_1,x_2,x_3,......,x_nx1,x2,x3,......,xn),使得:∑j=1n(ai,j−xj)2=C(i∈[1,n+1])\large\sum_{j=1}^{n}(a_{i,j}-x_j)^2=C (i\in[1,n+1])j=1∑n(a...原创 2019-04-26 08:50:29 · 83 阅读 · 0 评论 -
Zap[BZOJ1101]\[POI2007]
题面描述传送门思考题目要求等价于求有多少二元组(x,y)(x,y)(x,y)满足x≤a/k,y≤b/kx\le a/k,y\le b/kx≤a/k,y≤b/k并且x,yx,yx,y互质(因为gcd(x,y)=k\operatorname{gcd}(x,y)=kgcd(x,y)=k)。设D(a,b,k)D(a,b,k)D(a,b,k)表示满足x≤a,y≤bx\le a,y\le bx≤a,...原创 2019-04-29 11:55:13 · 179 阅读 · 0 评论 -
Hankson的趣味题[NOIP2009]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路因为lcm(x,b0)=b1\operatorname{lcm}(x,b_0)=b_1lcm(x,b0)=b1,所以xxx是b1b_1b1的因数,所以xxx的质因子一定也是时b1b_1b1的质因子。我们可以对b1b_1b1的每个质因子p,计算xxx可能包含多少个ppp。设a0,b0,a1,b1,xa_0,b_0...原创 2019-04-24 18:36:54 · 349 阅读 · 0 评论 -
Visible Lattice Points[POJ3090]
题面描述传送门思路详细参考#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<cstdlib>using namespace std;const int N=1010;const int inf=1000;int p...原创 2019-04-24 19:36:08 · 106 阅读 · 0 评论 -
莫比乌斯函数
设正整数N按照算数基本定理分解质因数为N=p1c1p2c2⋯pmcmN=p_1^{c_1}p_2^{c_2}\cdots p_m^{c_m}N=p1c1p2c2⋯pmcm,定义函数μ(N)={0∃i∈[1,m],ci>11m≡0(mod2),∀i∈[1,m],ci=1−1m≡1(mod2),∀i∈[1,m],ci=1\mu(N)=\begin{cases}0 &am...原创 2019-04-29 12:55:26 · 463 阅读 · 1 评论 -
The Luckiest Number[POJ3696]
题面描述传送门思路一句话题意:求至少多少个888连在一次组成的正整数是L\operatorname{L}L的倍数。∵x个8\because x个8∵x个8实际上就是(10x−1)/9∗8(10^x-1)/9*8(10x−1)/9∗8,现在我们可以转化一下题意:求最小的xxx,满足L∣(10x−1)/9∗8\operatorname{L} \mid (10^x-1)/9*8L∣(10...原创 2019-04-25 07:38:03 · 118 阅读 · 0 评论 -
Sumdiv[POJ1845]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路一句话题意:求ABA^BAB的所有约数之和mod9901(1≤A,B≤5∗107)\operatorname{mod}9901(1\le A,B\le 5*10^7)mod9901(1≤A,B≤5∗107)提要:等比数列求和公式:Sn=a1∗1−qn1−qS_n=a_1*\frac{1-q^n}{1-q}Sn=a1...原创 2019-04-25 09:23:32 · 102 阅读 · 0 评论 -
同余方程组[caioj1155]
题面描述同余方程是这样的:已知a,b,n,求x的最小正整数解,使得ax=b(mod m)同余方程组是这样:也是求x的最小正整数解,但已知b数组和m数组的情况下,x=b[1](mod m[1]),x=b[2](mod m[2]),x=b[3](mod m[3]),……x=b[n](mod m[n])【输入格式】一行一个整数 n(1<=n<=?)下来n行每行两个整数b[...原创 2019-04-16 20:45:54 · 305 阅读 · 0 评论 -
同余方程[NOIP2012/CH3301]
题面描述传送门思路详细参考#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<cstdlib>#define ll long longusing namespace std;ll exgcd(ll a,ll b,ll ...原创 2019-04-25 10:33:13 · 132 阅读 · 0 评论 -
Strange Way to Express Integers[POJ2891]
欢迎大家访问我的老师的OJ———caioj.cn题目描述传送门思路详细参考这道题输入有坑,且多组数据。#include<cstdio>#include<cstring>#include<algorithm>#include<cstdlib>#include<cmath>#define ll long longusi...原创 2019-04-25 10:43:14 · 121 阅读 · 0 评论 -
矩阵乘法原理
欢迎大家访问我的老师的OJ———caioj.cn一个n∗mn*mn∗m的矩阵可看作一个n∗mn*mn∗m的二维数组。矩阵的加法和减法,仅满足于两个同样大小(同为n∗mn*mn∗m)的矩阵,就是把矩阵对应位置上的数相加减,即C=A+B↔∀i∈[1,n],∀j∈[1,m],Ci,j=Ai,j+Bi,jC=A+B \leftrightarrow\forall i \in[1,n],\forall...原创 2019-04-25 14:21:20 · 2281 阅读 · 0 评论 -
Fibonacci[POJ3070]
题面描述传送门思路原创 2019-04-25 14:22:07 · 129 阅读 · 0 评论 -
GCD4[caioj1283]\[Zoj3435]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路算是仪仗队的进阶版吧。因为这题的左下角为(1,1,1)(1,1,1)(1,1,1),我们不妨将其变为(0,0,0)(0,0,0)(0,0,0),右上角为(a−1,b−1,c−1)(a-1,b-1,c-1)(a−1,b−1,c−1)方便讨论。首先我们可以观察到,三条棱上的(1,0,0),(0,1,0),(0,0,1)(1,...原创 2019-04-30 10:57:06 · 145 阅读 · 0 评论 -
余数之和[BZOJ1257]
题面描述传送门思路∵k mod i=k−⌊ k/i⌋∗i\because k ~\operatorname{mod}~{i}=k-\left\lfloor\ k/i \right\rfloor*i∵k mod i=k−⌊ k/i⌋∗i,原问题转化为:n∗k−∑i=1n⌊ k/i⌋∗i\large n*k-\sum_{i...原创 2019-04-24 15:17:48 · 286 阅读 · 0 评论 -
Devu and Flowes[Codeforces451E]
题面描述传送门(由于CF要翻墙,所以推荐用Luogu)思路经过观察,就是一道多重集的组合数的模板题CN+M−1N−1−∑i=1NCN+M−ni−2+∑1≤i<j≤NCN+M−ni−nj−3N−1−⋯+(−1)NCN+M−∑i=1Nni−(N+1)NC_{N+M-1}^{N-1}-\sum_{i=1}^NC_{N+M-n_i-2}+\sum_{1\le i<j\l...原创 2019-04-29 08:52:44 · 94 阅读 · 0 评论 -
容斥原理
容斥原理设S_1,S_2,…,S_n为有限集合,∣S∣\left\vert\\{S}\right\vert∣S∣表示集合S的大小,则:∣⋃i=1nSi∣=∑i=1n∣Si∣−∑1≤i<j≤n∣Si∩Sj∣+∑1≤i<j<k≤n∣Si∩Sj∩Sk∣+⋯\left\vert\bigcup_{i=1}^nS_i\right\vert=\sum_{i=1}^n...原创 2019-04-29 07:33:39 · 133 阅读 · 0 评论 -
开关问题[POJ1830]
题面描述传送门思路设xix_ixi表示第i个开关的操作情况,xi=1x_i=1xi=1表示按了这个开关,xi=0x_i=0xi=0表示没按。再统计ai,ja_{i,j}ai,j表示第i个开关和第j个开关的联系情况,ai,j=1a_{i,j}=1ai,j=1表示按下j会影响i的状态,ai,j=0a_{i,j}=0ai,j=0表示不会影响,特别地,令ai,i=1a_{i,i}=1ai...原创 2019-04-26 11:07:02 · 240 阅读 · 0 评论 -
装备购买[BZOJ4004]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路提要:做此题时,请务必了解线性空间的基以及矩阵的秩。把nnn件装备看作nnn个长度为mmm的向量,根据题意,购买的装备对应的向量应该是线性无关的(要满足向量不能被其他向量表出)。要买下最多数量的装备,其实就是求线性空间的任意一个基的个数。把ai.j(1≤i≤n,1≤j≤m)a_{i.j}(1\le i\le n,1...原创 2019-04-26 15:14:30 · 123 阅读 · 0 评论 -
XOR[HDOJ3949]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路用高斯消元求出a1,a2,......ana_1,a_2,......a_na1,a2,......an构成的异或空间的基。不妨设这个基由整数b1,b2,...btb_1,b_2,...b_tb1,b2,...bt构成,其中b1>b2>...>btb_1>b...原创 2019-04-26 16:47:54 · 145 阅读 · 0 评论 -
组合计数
加法原理若完成一件事的方法有n类,其中第i类方法包括a_i种不同的方法,且这些方法互补重合,则完成这件事共有a1+a2+...+ana_1+a_2+...+a_na1+a2+...+an种不同的方法。乘法原理若完成一件事需要n个步骤,其中第i个步骤有a_i种不同的完成方法,且这些步骤互不干扰,则完成这件事共有a1∗a2∗...∗ana_1*a_2*...*a_na1∗a2∗...∗...原创 2019-04-28 07:43:13 · 179 阅读 · 0 评论 -
计算系数[NOIP2011]\[CH3601]
题面描述传送门思路根据二项式定理原创 2019-04-28 07:54:52 · 154 阅读 · 0 评论 -
Counting Swaps[CH3602]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路对于一个排列p1,p2,...pnp_1,p_2,...p_np1,p2,...pn,如果从每个iii向pip_ipi连一条边,那么可以得到nnn个点nnn条边的图,并且这张图由若干个环构成。列入排列2,4,6,1,5,32,4,6,1,5,32,4,6,1,5,3对应下图,由1−2−4,3−61-2-4,3-61...原创 2019-04-28 11:00:05 · 178 阅读 · 0 评论 -
Lucas定理
若p是质数,则对于任意整数1≤m≤n1\le m\le n1≤m≤n,有: Cnm≡Cnmodpmmodp∗C⌊np⌋⌊mp⌋(modp)\huge\ C_n^m\equiv C_{n \operatorname{mod}p}^{m \operatorname{mod} p}*C_{\left\lfloor\frac{n}{p}\right\rfloor}^{\left\lflo...原创 2019-04-28 13:55:55 · 121 阅读 · 0 评论 -
中国剩余定理
设m1,m2,...,mnm_1,m_2,...,m_nm1,m2,...,mn是两两互质的整数,m=∏i=1nmi,tim=\prod_{i=1}^nm_i,t_im=∏i=1nmi,ti是线性同余方程Miti≡1(modmi)M_it_i\equiv1(\operatorname{mod} m_i)Miti≡1(modmi)的一个解。对于任意的n个整数a1,a2,...,a...原创 2019-04-28 14:20:19 · 139 阅读 · 0 评论 -
Prime Distance[POJ2689]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门POJ传送门(数据更强)思路L,R\operatorname{L,R}L,R范围很大,所以直接线性筛是过不了的。但是L-R\operatorname{L-R}L-R的范围只有10610^6106,对于[L,R]\operatorname{[L,R]}[L,R]区间内的任意合数nnn,一定存在小于n\sqrt{n}n...原创 2019-04-24 10:05:21 · 114 阅读 · 0 评论 -
阶乘分解[CH3101]
题面描述传送门思路算术基本定理:N=p1c1p2c2......pmcmN=p_1^{c_1}p_2^{c_2}......p_m^{c_m}N=p1c1p2c2......pmcm题意让我们求N!\operatorname{N!}N!的每一个质因子pip_ipi的个数cic_iciN!\operatorname{N!}N!中质因子ppp的个数实际上就是[1,n][...原创 2019-04-24 11:09:32 · 215 阅读 · 0 评论 -
高中数学必备知识
∑i=1n\huge\sum_{i=1}^ni=1∑n即int ans=0;for(int i=1;i<=n;i++)ans+=i;∏i=1n\Huge\prod_{i=1}^ni=1∏n即int ans=0;for(int i=1;i<=n;i++)ans*=i;x∣y x整除y,即...原创 2019-04-19 19:53:02 · 415 阅读 · 4 评论 -
反素数[BZOJ1053]
题面描述传送门思路三条引理:设mmm是1~N中约数个数最多的数中最小的一个。根据m的定义,m显然满足:∀x<m,g(x)<g(m)\forall x<m,g(x)<g(m)∀x<m,g(x)<g(m)∀x>m,g(x)≤g(m)\forall x>m,g(x)\le g(m)∀x>...原创 2019-04-24 14:13:13 · 176 阅读 · 0 评论 -
石头游戏[CH3401]
题面描述传送门思路题目可以理解为是一个n∗mn*mn∗m的矩阵,在里面进行一些玄学的操作。不难发现,操作序列的长度不超过6,那么1~6的最小公倍数是60,即每经过60秒,所有操作序列都会重新处于最开始的字符处。那么第k(1≤k≤60)k(1\le k\le 60)k(1≤k≤60)秒,第k+60k+60k+60秒执行的字符与第kkk秒一定是相同的。得到了这个结论,我们就可以很容易想到...原创 2019-04-25 18:25:51 · 1130 阅读 · 0 评论