Siyao模板库.docVIP

Siyao@Sphinx 模板库 目录 Number Theroy 1 1、 mul_mod 同余乘：避免乘法溢出 1 2、 GCD算法：求两数最大公约数 已经测试 2 3、 EX_GCD扩展欧几里德：递归求ax+by=gcd(a,b)=d的一组解x,y.效率：o(logn). 2 4、 CF_GCD连分数方法解 ax=gcd(a,b) mod(b) 以及 ax+by=gcd(a,b)=d. 2 5、 msolve(),inv()模线性方程和逆元:解ax=b(mod n). 2 6、 pow_mod()快速幂取模：返回a^n mod p 3 7、 同余一元方程组(2个)：返回是否有解，解为x=m(mod a) 已经测试 FZU 1402 3 8、 同余一元方程组(剩余定理)：x=ai (mod bi),bi两两互素。得到x=x0(mod b1……bn) 3 9、 同余一元方程组(一般情况)：不停合并两方程 3 10、 log()离散对数(模为素数):返回满足a^x=b(mod n)的最小解，无解返回-1 4 11、 mlist()筛法得素数表：int prim[]中保存素数,pnum保存个数，hash可以用bool型 4 12、 pfact()分解素因子：int pf[]中保存素因子int lim[]保存幂次,pfnum保存个数 5 13、 phi(),phis()欧拉函数：phi求1-n中与n互素的数的个数；phis求1-n所有数的phi 5 14、 square_free()筛法生成Square-free数：int sf[]中存结果，shh[]中存因子数 (以便进行容斥) 6 15、 序列中互质的m元组(容斥原理square-free)：m元组最大公约数为1 已经测试 TJU3121 6 16、 区间与数互质(容斥原理dfs)：求[1,r]区间中与x互质的数的个数 已经测试 TJU3317 6 17、 区间与区间互质(容斥原理square-free)：求[1,x],[1,y]互质的对数 7 18、 区间与区间互质(容斥原理dfs)：求[1,x],[1,y]互质的对数 7 19、 枚举约数及对应的欧拉函数：已经测试 HDU2588 7 20、 欧拉函数加速递归函数取模： 已经测试 POJ2720 TJU3313 7 21、 阶乘中素数的阶 已经测试 FZU 1753 8 22、 C(N,M) mod P组合数取模 已经测试 HIT 2813 8 23、 n!(cop) mod P阶乘中与素数互素的部分模p 9 24、 n!n!n!/n!n!n!n! 多阶乘乘除模k(k为square-free数) 9 25、 n!(con) inv mod k 阶乘与k互素部分ff以及逆元fv打表 已经测试 FZU1833 9 26、 sum_mod 等比数列求和取模 10 27、 floor,ceil高斯函数：取上下整 已经测试 SGU 106 10 28、 一次不等式求整数解：L=ax+b=R 已经测试 SGU 106 10 29、 二元一次方程求整数解对：ax+by+c=0 xl=x=xr yl=y=yr 已经测试 SGU 106 10 30、 素数测试之miller-rabin： 返回false则一定是合数 已经测试 PKU1811 10 31、 素数测试函数：isprime(n) 11 32、 Pollard-rho算法 对合数找到一个因数 已经测试 PKU1811 11 33、 分解为素因子的积：利用 Pollard-rho算法 (n为1特殊处理了) 12 34、 Farey逼近分数：得到分母n中与目标分数最接近分数 已测试ECNU 2649 12 35、 Pell方程：x^2-dy^2=1的解 12 36、 带限制可重排列： 13 Number Theroy mul_mod 同余乘：避免乘法溢出 typedef __int64 ll; typedef const __int64 cll; ll mul_mod(ll a,ll b,cll n){ //二进制思想 已经测试 PKU1811 ll ans(0),tmp((a%n+n)%n); b=(b%n+n)%n; //b%=n; while(b){ if(b1) if((ans+=tmp)=n) ans-=n; if((tmp=1)=n) tmp-=n; b=1; } return ans; } GCD算法：求两数最大公约数 已经测试 欧几里得算法.适用非负整数.效率：o(logn). ll gcd(ll a,ll b) { if(!a)return b; if(!b)re