- 博客(25)
- 问答 (2)
- 收藏
- 关注
RSS订阅原创 Mathematica绘制图形
绘制平面图形的命令Plot 显函数绘图ParametricPlot 参数方程绘图PolarPlot 极坐标方程绘图ContourPlot 隐函数绘图ListPlot 绘制点列图例Plot[x Sin[1/x],{x,-1,1}]例绘制三维图形Plot3D 绘制二元函数对应的空间曲面ParametricPlot3D 绘制空间曲线、曲面参数方程对应空间曲线与曲面例常...
2020-04-14 21:43:22
6247
1
原创 poj2104 K-th number(模板)
/************************************************┆ ┏┓ ┏┓ ┆┆┏┛┻━━━┛┻┓ ┆┆┃ ┃ ┆┆┃ ━ ┃ ┆┆┃ ┳┛ ┗┳ ┃ ┆┆┃ ┃ ┆┆┃ ┻ ┃ ┆┆┗━┓ ┏━┛ ┆┆ ┃ ┃ ┆ ┆ ┃ ┗━━━┓ ┆┆ ┃ AC代马...
2019-07-30 14:17:58
111
原创 An Easy Physics Problem
首先分成和圆碰撞和不与圆碰撞两种情况：1.点A和圆不碰撞，这个很好判断上面是射线的轨迹方程，下面是圆的轨迹方程，合并，求解t，若t无解或仅有一个解（即bb-4a*c<=0），直接判断B...
2019-07-17 19:47:13
299
1
原创 The Tower HDU - 6559
题目不是很难，把公式推出来以后直接代入，/************************************************┆ ┏┓ ┏┓ ┆┆┏┛┻━━━┛┻┓ ┆┆┃ ┃ ┆┆┃ ━ ┃ ┆┆┃ ┳┛ ┗┳ ┃ ┆┆┃ ┃ ┆┆┃ ┻ ┃ ┆┆┗━┓ ┏━┛ ┆┆ ┃ ┃ ┆ ┆ ┃ ...
2019-07-17 14:54:39
143
原创 归并求逆序对【模板】
const int maxn = 5e5 + 10;int arr[maxn];int tmp[maxn];ll merge(int l,int r){ if(l == r) return 0; ll ans=0; int mid = l + r >> 1; ans+=merge(l,mid),ans+=merge(mid+1,r); int i = l,j = mi...
2019-07-12 15:16:21
215
原创 二叉查找树【模板】
const int maxn = 1e6 + 10;struct Treap{ int l,r; int val,dat; int cnt,size;}a[maxn];int tot,root,INF=0x7fffffff;int New(int val){ a[++tot].val=val; a[tot].dat=rand(); a[tot].cnt=a[tot].size...
2019-07-11 16:26:07
94
原创 二叉查找树【初模板】
//二叉查找树BST const int maxn = 2e5 + 10;int tot,root,INF=1<<30;struct BST{ int l,r; int val;}a[maxn];int New(int val){ a[++tot].val=val; //a[tot].l=a[tot].r=0; return tot;}void build(){...
2019-07-10 20:57:25
79
原创 树状数组【模板】
const int maxn = 2e5 + 10;int c[maxn];int ask(int x){ int ans=0; for(;x;x-=x&-x) ans+=c[x]; return ans;}void add(int n,int x,int y){ for(;x<=n;x+=x&-x){ c[x]+=y; }}
2019-07-10 19:14:42
72
原创 线段树【模板】//待检验
const int MAXN=50010;int a[MAXN],ans[MAXN<<2],lazy[MAXN<<2];//a[]为原序列信息，ans[]模拟线段树维护区间和，lazy[]为懒惰标记void PushUp(int rt){ ans[rt]=ans[rt<<1]+ans[rt<<1|1];}void Build(in...
2019-07-10 16:32:49
107
原创 可持久化线段树【模板】poj_2104
/************************************************┆ ┏┓ ┏┓ ┆┆┏┛┻━━━┛┻┓ ┆┆┃ ┃ ┆┆┃ ━ ┃ ┆┆┃ ┳┛ ┗┳ ┃ ┆┆┃ ┃ ┆┆┃ ┻ ┃ ┆┆┗━┓ ┏━┛ ┆┆ ┃ ┃ ┆ ┆ ┃ ┗━━━┓ ┆┆ ┃ AC代马...
2019-07-10 16:07:38
95
原创 权值线段树【模板】
const int N = 2e4+100;int T[N<<2];void update(int p,int v,int rt,int l,int r)//值为p的个数加v(-1表示删除，1表示添加) { T[rt]+=v; if(l==r) return ; int m=(l+r>>1); if(p<=m) update(p,v,rt<<...
2019-07-10 16:02:39
213
原创 后缀数组（倍增）【模板】
const int maxn = 2e5+10;char s[maxn];int y[maxn],x[maxn],c[maxn],sa[maxn],rk[maxn],height[maxn];void get_SA(int n,int m)//n表示待求串s的长度，下标由1开始，m表示max(s[i]) { for(int i=1;i<=m;++i) c[i]=0;//桶的初始化 ...
2019-07-10 15:57:21
225
原创 ST表【模板】
const int maxn = 1e5+10;const int maxn2=log2(maxn)+1;//maxn2的值最好直接计算得出 int ST[maxn][maxn2];void Init(int n){ for(int j=1;j<=maxn2;j++){ for(int i=1;i+(1<<j)-1<=n;i++){ ST[i][j]=...
2019-07-10 15:54:09
102
原创 S-Nim POJ - 2960
/************************************************┆ ┏┓ ┏┓ ┆┆┏┛┻━━━┛┻┓ ┆┆┃ ┃ ┆┆┃ ━ ┃ ┆┆┃ ┳┛ ┗┳ ┃ ┆┆┃ ┃ ┆┆┃ ┻ ┃ ┆┆┗━┓ ┏━┛ ┆┆ ┃ ┃ ┆ ┆ ┃ ┗━━━┓ ┆┆ ┃ AC代马...
2019-04-18 23:51:54
110
原创 Magical GCD UVA - 1642
紫书上的数论经典题。题意大概是这样的，给一段正整数序列，问gcd(a i,a i+1,…,a j)*(j-i+1)的最大值。其实不能，就是右端点从前往右扫，遍历整个区间，左端点从右端点往左扫，直到起点。数据维护每个gcd,index.反正核心就是相同gcd的删掉，时间复杂度大概就是nlong(ai)/**********************************************...
2019-04-14 20:32:40
145
转载 高斯消元法解方程组（可以判断有无整数解）
// 高斯消元法解方程组(Gauss-Jordan elimination).(-2表示有浮点数解，但无整数解，//-1表示无解，0表示唯一解，大于0表示无穷解，并返回自由变元的个数)//有equ个方程，var个变元。增广矩阵行数为equ,分别为0到equ-1,列数为var+1,分别为0到var.const int MAXN=50;int a[MAXN][MAXN];//增广矩阵int...
2019-04-03 20:08:33
1788
转载 高斯消元解异或方程
const int maxn = 300;int a[maxn][maxn];int Free_x[maxn]; //自由变元int x[maxn]; //解集 int Free_num = 0; //自由变元数int Guass(int Row, int Column) //系数矩阵的行和列{ int row ...
2019-04-03 14:53:07
403
翻译 高斯消元模板
欢迎使用Markdown编辑器//高斯消元const int N = 20;const double EPS=1e-7;double a[N][N],x[N];double c[N][N];//m表示行数，n表示列数（下标从0开始储存，常数列不算进列数）//返回值：-1表示无解，0表示唯一解，大于0表示无穷解，并返回自由变元的个数int Gauss(int m,int n){ ...
2019-04-01 20:39:45
95
转载 递推算法之平面分割问题总结
这是一类问题，首先由直线划分区域到折线划分区域，再延伸到封闭图形划分区域，最后在推广为平面划分空间的问题。一、n条直线最多分平面问题题目大致如:n条直线，最多可以把平面分为多少个区域。 析:可能你以前就见过这题目，这充其量是一道初中的思考题。当有n-1条直线时，平面最多被分成了f（n-1）个区域。则第n条直线要是切成的区域数最多，就必须与每条直线相交且不能有同一交点。 这样就会得到n-1...
2019-02-24 20:09:42
324
原创 【模板】拓展卢卡斯定理
//拓展卢卡斯定理int pow_mod(int a, int n, int mod)//a^n%mod(a的n次方求余）{ if (n == 0) return 1; int x = pow_mod(a, n / 2, mod); int ans = (long long)x * x % mod; if (n % 2) ans = ((long long)ans * a) %...
2019-02-20 14:26:26
213
原创 逆元【模板】
（https://www.cnblogs.com/zzqc/p/7192436.html）的简写方程ax≡1(mod p),的解称为a关于模p的逆，当gcd(a,p)==1（即a，p互质）时，方程有唯一解，否则无解。费马小定理内容：假如p是质数，且gcd(a,p)=1，那么 a(p-1)≡1（mod p）a关于p的逆元就是a^(p-2) % p,直接快速幂处理拓展欧几里得算法ax+b...
2019-02-19 14:05:56
294
原创 [模板]拓展中国剩余定理
#define maxn 100//模方程数量int a[maxn], r[maxn];void Init(int k) //初始化函数，有k个等式，读入除数a和余数r{for (int i = 0; i &lt; k; i++) scanf("%d", &amp;a[i]);for (int i = 0; i &lt; k; i++) scanf("%d", &
2019-02-18 22:39:40
113
原创 卢卡斯定理（模板）
const int MAXN = 1e5 + 10;//MAXN要大于modint f[MAXN];int pow_mod(int a, int n, int mod)//费马小定理，a^(mod-2)%mod(mod是质数）{if (n == 0)return 1;int x = pow_mod(a, n / 2, mod);int ans = (long long)x * x %...
2019-02-02 17:47:17
94
空空如也
USB红外测温摄像头内核驱动，video_device注册失败
2022-06-17
UVC内核驱动注册失败
2022-06-17
TA创建的收藏夹 TA关注的收藏夹
TA关注的人