C20181503csy的博客

通分：xyx+y=n!\frac{xy}{x+y}=n!x+yxy​=n!移项并打开括号：xy=n!x+n!yxy=n!x+n!yxy=n!x+n!y合并同类项：(y−n!)x=n!y⟹x=n!yy−n!(y-n!)x=n!y \quad \Longrightarrow \quad x=\frac{n!y}{y-n!}(y−n!)x=n!y⟹x=y−n!n!y​令t=y−n!t=y...

【BZOJ 3667】第一行：CAS,代表数据组数（不大于350），以下CAS行，每行一个数字，保证在64位长整形范围内，并且没有负数。你需要对于每个数字：第一，检验是否是质数，是质数就输出Prime第二，如果不是质数，输出它最大的质因子是哪个。检验一个数是否是质数的朴素方法是O(n)O(\sqrt{n})O(n​)的，并不能满足这道题的需求。其实直到现在，人们对一个数的素性判定依旧没...

【题目描述】对于给出的nnn个询问，每次求有多少个数对(x,y)(x,y)(x,y)，满足a≤x≤b，c≤y≤da≤x≤b，c≤y≤da≤x≤b，c≤y≤d，且gcd(x,y)=kgcd(x,y) = kgcd(x,y)=k.【分析】记f(a,b)=∑i=1a∑j=1b[gcd(i,j)=k]f(a, b)=\sum_{i=1}^a\sum_{j=1}^b[gcd(i,j)=k]f(a,b...

【题目描述】根据一些书上的记载，上帝的一次失败的创世经历是这样的：第一天， 上帝创造了一个世界的基本元素，称做“元”。第二天， 上帝创造了一个新的元素，称作“α”。“α”被定义为“元”构成的集合。容易发现，一共有两种不同的“α”。第三天， 上帝又创造了一个新的元素，称作“β”。“β”被定义为“α”构成的集合。容易发现，一共有四种不同的“β”。第四天， 上帝创造了新的元素“γ”，“γ”被定...

在平常的练习中，我们经常需要求a的b次幂模p的结果。这次就专门讲讲这个。之前，b的规模不是很大。所以我们都是这样做的：for(int i = 1; i <= b; i++) ans *= a, ans %= p;然而，当b很大时，如10^9时，这种方法就大大滴超时了。。。为了搞定b较大时的幂取模运算，在这里提一下快速幂取模的方法。以27^11为例:1、将11转化为二

Lucas定理是用来求C(n,m)%p的值，其中n与m均为非负整数，p是素数。先说说这个值应该怎么算。设Lucas(n, m, p)为上式的值，则Lucas(n, m, p) = cn(n % p, m % p) * Lucas(n / p, m / p, p), Lucas(x, 0, p) = 1;其中cm(n, m) = a! * (b! * (a - b)!) ^ (p - 2

题意:求∑i=1n∑j=1nlcm(i,j)\sum_{i=1}^n\sum_{j=1}^nlcm(i, j)∑i=1n​∑j=1n​lcm(i,j)题解：原式可化为：∑i=1n∑j=1nijgcd(i,j)\sum_{i=1}^n\sum_{j=1}^n\frac{ij}{gcd(i, j)}i=1∑n​j=1∑n​gcd(i,j)ij​∑i=1n∑j=1n∑d=1nijd[d=gcd(...

先来康康这个问题。给定三个正整数A,B,PA,B,PA,B,P，求关于xxx的方程的最小正整数解：Ax≡B(mod P)A^x\equiv B(mod\ P)Ax≡B(mod P)我们令x=a⌈P⌉−b,a∈N,b∈Nx=a\lceil \sqrt{P}\rceil-b,a \in N,b \in Nx=a⌈P​⌉−b,a∈N,b∈N由欧拉定理：Aϕ(P)≡1(mod...