算法模板【最后更新-05】.docVIP

算法模板 网络122 王硕 PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 1 目录 1.数学 1.1 矩阵 ………………………2 1.2 一次方程组的解 ……………3 1.3 矩阵的逆 …………………4 1.4 最大公约数 …………………6 1.5 欧几里得扩展 ……………6 1.6 素数表 ………………………7 1.7 判断素数 …………………8 1.8 分解质因数 …………………9 1.9 欧拉函数 …………………10 1.10 欧拉函数表 …………………11 1.11 mobius函数 …………………12 1.12 数值积分 …………………12 1.13 数值积分（精确） ………13 1.14 快速幂求余 …………………14 1.15 进制转换 …………………15 1.16 格雷码 ………………………16 1.17 高精度类 …………………16 1.18 Fibonacci数列 ……………22 1.19 FFT ………………………23 2.图论 2.1 拓扑排序 …………………24 2.2 dijkstra …………………25 2.3 floyd-warshall ……………26 2.4 kruskal …………………26 3.序列 3.1 快速排序 …………………28 3.2 逆序对 ………………………28 3.3 最长不降子序列长度 ………30 3.4 最长公共子序列长度 ………31 3.5 最长公共上升子序列长度 …31 3.6 最长公共上升子序列 ………32 4.字符串 4.1 Sunday ………………………34 4.2 子串清除 …………………34 4.3 KMP ………………………37 5.数据结构 5.1 并查集 ………………………38 5.2 树状数组 …………………39 5.3 左偏树 ………………………40 5.4 哈希 ………………………41 5.5 RMQ线段树 …………………41 6.其他 6.1 多边形面积 …………………43 6.2 幻方构造 …………………43 6.3 奇数阶次幻方 ……………45 1.1矩阵 #include iostream #include cstring using namespace std; const int MAXN=100; const int MAXM=100; struct Matrix { int n,m; int a[MAXN][MAXM]; void clear() { m=n=0; memset(a,0,sizeof(a)); } Matrix operator +(const Matrix b) const { Matrix tmp; tmp.n=n; tmp.m=m; for(int i=0;in;i++) for(int j=0;jm;j++) tmp.a[i][j]=a[i][j]+b.a[i][j]; return tmp; } Matrix operator *(const Matrix b) const { Matrix tmp; tmp.clear(); tmp.n=n; tmp.m=m; for(int i=0;in;i++) for(int j=0;jb.m;j++) for(int k=0;km;k++) tmp.a[i][j]+=a[i][k]*b.a[k][j]; return tmp; } }; 1.2一次方程组的解 #include iostream #include cstring #include cmath #include algorithm #define MAX 10 #define EPS 1e-8 using namespace std; double a[MAX][MAX+1]; bool l[MAX]; double ans[MAX]; void swap(double x,double y) { double t; t=x; x=y; y=t; } inline int solve(double a[MAX][MAX+1],bool l[],double ans[],const int n) { int res=0,r=0; int i,j,k; memset(ans,0,sizeof(a)); memset(l,false,n); for(i=0;in;i++) { for(j=r;jn;j++) if(fabs(a[i][j])EPS) {