- 博客(25)
- 资源 (1)
- 收藏
- 关注
原创 容斥定理与错排问题
关于错排问题给定 nnn 个数字,求数字 1,2,⋯ ,n1,2,\cdots,n1,2,⋯,n 不排列在第 1,2,⋯ ,n1,2,\cdots,n1,2,⋯,n 上的排列数量称为错排问题。容斥定理X1,⋯ ,XnX_1,\cdots,X_nX1,⋯,Xn 为可数集合,那么 ∣X1∪⋯∪Xn∣|X_1\cup\cdots\cup X_n|∣X1∪⋯∪Xn∣ 为∣X1∪⋯∪Xn∣=∑k=1,1⩽i1<⋯<ik⩽nn(−1)k+1∣Xi1∩⋯∩Xik∣|X_1\cup\cdots
2022-04-21 17:48:55
338
2
原创 旋转公式是怎么推导的
GAMES103 L03 旋转公式是怎么推导的GAMES103 L03考虑代表旋转的四元数 qqq,那么对于参考坐标中的任意一点 v\mathbf vv,其在世界坐标中的朝向为 qvq∗q\mathbf v q^*qvq∗。此时这个物体拥有角速度 w\mathbf ww,角速度的方向就是物体旋转的朝向,角速度的模就是物体的角速度。考虑一小段时间内,假设 w\mathbf ww 在这段时间内不发生变化,那么物体就在这段时间内沿着 w\mathbf ww 确定的方向旋转了 Δθ=∥w∥Δt\Delt
2022-03-15 12:52:41
841
原创 关于转动张量的一部分推导
关于惯性张量的推导对于刚体 AAA,其角动量为E=∫A(r×r˙) dmE=\int_A(r\times\dot r)\,\mathrm dmE=∫A(r×r˙)dm显然该点的线速度与角速度之间有如下关系r˙=ω×r\dot r=\omega\times rr˙=ω×r那么该物体的角动量E=∫Ar×(ω×r) dmE=\int_Ar\times(\omega\times r)\,\mathrm dmE=∫Ar×(ω×r)dm根据向量三重积关系 a×(b×c)=(a⋅c)b−(a
2022-03-13 14:48:35
1715
原创 函数拟合中涉及的矩阵运算
函数拟合中涉及的矩阵运算简便起见,采用幂函数作为基函数。插值插值问题要求对于数据 {(xi,yi):i=0,⋯ ,n−1}\{(x_i,y_i):i=0,\cdots,n-1\}{(xi,yi):i=0,⋯,n−1},找到 nnn 个基函数 bi(x)b_i(x)bi(x) 的线性组合参数 w0,⋯ ,wn−1w_0,\cdots,w_{n-1}w0,⋯,wn−1 使得函数f(x)=w0+w1x+⋯+wn−1xn−1f(x)=w_0+w_1x+\cdots+w_{n-1}x^{n-1}
2022-02-16 22:23:42
751
原创 一个实用的编译技术
一个实用的编译技术使用编译参数结合宏定义实现了在编译时指定不同的输出等级(makefile + c)// defs.h#include <stdio.h>// Output level: ERROR(31) > WARN(93) > INFO(34) > DEBUG(32) > TRACE(90)// `fmt` can only be face-value#ifdef TRACE_OUTPUT #define trace(fmt, arg...)
2021-10-19 20:07:04
114
原创 【图形学】GAMES101 Assignment3 作业框架分析
GAMES101 Assignment3 作业框架分析文章目录GAMES101 Assignment3 作业框架分析rasterizer::drawrasterizer::rasterize_triangle写作业之前看了一眼代码框架,发现疑问挺多的,多方查证之后基本得到了解决。本文主要分析一下 rst::rasterizer::draw 和 rst::rasterizer::rasterize_triangle 方法的内容。这一篇分析参考了网上的资料,代码部分参考了github以为大佬的代码。但是我
2021-09-13 19:37:07
1639
5
原创 【图形学】关于透视校正
关于透视校正使用如下的举证进行透视投影:proj=[n0000n0000n+fnf00−10]\mathrm{proj}=\begin{bmatrix}n&0&0&0\\0&n&0&0\\0&0&n+f&nf\\0&0&-1&0\end{bmatrix}proj=⎣⎢⎢⎡n0000n0000n+f−100nf0⎦⎥⎥⎤假设有两个点 P(x1,y1,z1,w1),Q(x2,y2,z2,w2
2021-09-07 12:07:47
508
原创 【图形学】关于透视投影
文章目录透视投影透视投影我们认定此时已经进行了 module 和 view 变换,剩下的只需要进行投影变换。此外,我们认为摄像机位于原点,面向 −z-z−z 方向,以 yyy 方向为上方向,n,fn,fn,f 为正数,采用齐次坐标。已知投影透视的三条约束:近平面保持不动;远平面的 zzz 轴不动,只有 x,yx,yx,y 的缩放;位于 x=0,y=0x=0,y=0x=0,y=0 的点保持不动。设投影变换的矩阵为 MMM,由于 z<−nz<-nz<−n 的任意一点 (x′
2021-08-17 23:24:41
1416
原创 【计算理论导引】 课后习题1.48
【计算理论导引】 课后习题1.48题干设 Σ={0,1}\Sigma=\{0,1\}Σ={0,1} 以及D={w∣w中子串01与10的出现次数相等},D=\{w|w中子串01与10的出现次数相等\},D={w∣w中子串01与10的出现次数相等},由于 101101101 含有一个 010101 和一个 101010,因此 101∈D101\in D101∈D。而 101010101010 含有两个 101010 和一个 010101,所以 1010∉D1010\not\in D1010∈D
2021-06-03 23:58:25
1293
1
原创 【C++】拷贝构造函数于临时变量小探
文章目录【C++】拷贝构造函数于临时变量小探【C++】拷贝构造函数于临时变量小探今天遇到了一个关于拷贝构造于临时变量的小问题,在这里记录一下。这是代码:#include <iostream>using namespace std;class s1mple {public: int i; s1mple() : i(23) { cout << "constructor" << endl; } s1mple(c
2021-05-13 23:17:49
353
原创 【fft】使用快速傅里叶变换计算大整数乘法
文章目录使用快速傅里叶变换计算大整数乘法1. 从多项式入手2. 换一种方式看待多项式3. 单位复根4. 求点值表达5. 求系数6. 代码使用快速傅里叶变换计算大整数乘法这是一道网络教室的题目,本来是要用分块乘法解决的,但是快速傅里叶也是一种解法。1. 从多项式入手实数可以表示成为多项式的形式,举一个例子143=3×100+4×101+1×102143=3\times10^0+4\times10^1+1\times10^2143=3×100+4×101+1×102这说明 143143143 是多项
2021-04-16 21:43:27
3773
7
原创 【Leetcode】63. 不同路径 II
【Leetcode】63. 不同路径 II题目一个机器人位于一个 m x n 网格的左上角 (起始点在下图中标记为“Start” )。机器人每次只能向下或者向右移动一步。机器人试图达到网格的右下角(在下图中标记为“Finish”)。现在考虑网格中有障碍物。那么从左上角到右下角将会有多少条不同的路径?网格中的障碍物和空位置分别用 1 和 0 来表示。来源:力扣(LeetCode)链接:https://leetcode-cn.com/problems/unique-paths-ii著作权归领扣网
2021-03-28 17:25:07
80
原创 Latex standalone tips
Latex standalone tips在使用 standalone 输出 tikz 包的图形的时候出现了文档错位的问题:检查之后发现是使用了 geometry 包的缘故,将其删去后图形正常。附上代码:\documentclass[UTF8]{standalone}\usepackage{amsmath,amssymb}\usepackage{tikz}%\usepackage{geometry}\usetikzlibrary{graphs}\begin{document}
2021-02-11 18:47:00
1227
原创 伽马gamma函数
伽马函数 gamma function定义如下Γ(x)=∫0+∞ux−1e−udx,x∈R\Gamma(x)=\int_0^{+\infty}{u^{x-1}e^{-u}\mathrm dx},x\in\mathbb RΓ(x)=∫0+∞ux−1e−udx,x∈R根据分部积分有如下性质Γ(x)=(x−1) Γ(x)\Gamma(x)=(x-1)\,\Gamma(x)Γ(x)=(x−1)Γ(x)以下计算 Γ(12)\Gamma(\frac12)Γ(21)Γ(1/2)=∫0+∞1t e−
2020-12-27 11:05:11
1362
原创 植物大战僵尸资源文件提取 总结
PvZ资源文件提取 总结参考植物大战僵尸资源文件(main.pak):格式/提取/工具下载[转贴]植物大战僵尸pak资源文件分析植物大战僵尸资源提取C/C++中判断某一文件或目录是否存在PNG文件格式详解C语言创建文件夹Big Endian 和 Little Endian 详解资源文件格式资源文件名字为main.pak,与主程序PlantsVsZombies.exe位于同一目录下,遵循以下加密方式:资源文件分为两个部分,前一部分记录了每一个资源(包括图片、音乐)的名称长度、文件名称(
2020-12-26 13:53:05
8021
原创 二分查找的平均查找长度
二分查找的平均查找长度对二分查找的平均查找长度进行简单分析。向作出假设:要查找的元素在数组内,数组长度为 nnn. 约定对长度为 nnn 的数组,平均查找长度为随机变量 CnC_nCn,随机变量 InI_nIn 定义如下In={1要查找的关键字是这个长度为n的数组的中间的关键字0otherwise,I_n=\left\{\begin{matrix}1 & 要查找的关键字是这个长度为 n 的数组的中间的关键字\\0 & \mathrm{otherwise}\end{matr
2020-12-19 09:44:02
7800
2
原创 对称差的结合律
对称差的结合律对于集合 A,BA,BA,B,有(A⊕B)⊕C=((A−B)∪(B−A))⊕C=((A∪B)∩(Aˉ∪Bˉ))⊕C=((A∪B)∩(Aˉ∪Bˉ)∪C)∩((Aˉ∩Bˉ)∪(A∩B)∪Cˉ)=((A∪B∪C)∩(Aˉ∪Bˉ∪C))∩(((Aˉ∪B)∩(A∪Bˉ))∪Cˉ)=((A∪B∪C)∩(Aˉ∪Bˉ∪C))∩((Aˉ∪B∪Cˉ)∩(A∪Bˉ∪Cˉ))=(A∪B∪C)∩(Aˉ∪Bˉ∪C)∩(Aˉ∪B∪Cˉ)∩(A∪Bˉ∪Cˉ)\begin{aligned}(A\oplus B)\opl
2020-11-25 20:52:20
8633
1
原创 判断出栈序列
判断出栈序列题目如下请判断:指定的序列能否仅由 入栈 和 出栈 操作得到。输入:有若干组数据输入每组数据中,第一行为两个个整数 n 和 m。n 表示需要依次从 1~n 入栈,m 表示这组数据有 m 个出栈序列需要判断,当 n=0 且 m=0 时停止。接下来有行,每行表示一个出栈序列输出:对每一个出栈序列,如果能正常出栈,则输出 Yes,否则输出 No。sample:input:5 21 2 3 4 55 4 1 2 36 16 5 4 3 2 10 0output:Ye
2020-10-24 22:18:11
431
原创 任意精度的小数运算
任意精度的小数运算问题的大致内容如下给定正整数 n,m,k(n<m)n,m,k(n<m)n,m,k(n<m),求 n/mn/mn/m 的值,保留 kkk 位小数。这道题的我分为两步解决计算 n/mn/mn/m 的任意位小数;通过计算的结果构建链表或类似循环链表的结构;根据链表或类似循环链表的结构输出。计算任意精度的小数有穷小数重温除法的竖式计算过程,本质上是整除与取余数的过程,在除法的每一步中商作为结果对应位的数,余数乘以十后作为下一步的被除数。以 2/72/72
2020-10-11 12:16:35
606
原创 RB-DELETE:红黑树删除操作
RB-DELETE:红黑树删除操作参考文章1、没有子节点1)红色直接删除。2)黑色删除,之后进行自平衡。由于这个节点是黑色的,那么必然存在黑色兄弟节点或者在兄弟节点的位置的子树的黑高为一。2、有一个子节点将这个节点沿着唯一的子节点方向向下交换移动,直到:1)移动到没有子节点参考1进行操作。2)移动到拥有两个子节点参考3进行操作。3、有两个子节点将这个节点的后继节点与这个节点交换值,但颜色不变。然后在原先后继节点的位置删除原本要删除的节点,首先,这个节点必然是其父节点的左子节点。
2020-10-07 23:03:13
297
原创 北京理工大学小学期乐学 t23树上统计
写在前头玩一下csdn的markdown。这道题的大意是给定一棵树,两个相邻点之间的距离为111,对于相邻距离为222的点之间可以新建一条边,要求每个点到其他所有点点的距离之和的和。看似是一道搜索题,但是我多日后再想起来感觉还有动态规划的思路在里头。t23参考帖子:http://lexue.bit.edu.cn/mod/forum/discuss.php?d=126124,阅读前请先看一看这篇文章。如果看不到原帖,这里大致概括本贴中用到的原帖中的结论,对于相邻的节点 i,ji,ji,j,原
2020-10-07 22:16:06
1526
原创 三种交换变量的方法
交换变量文章目录交换变量第三者交换法加法交换法异或交换法总结以下约定要交换的两个变量为 x , y 。第三者交换法定义第三个变量,用于交换两个变量的值。var tmp = x;x = y;y = x;加法交换法这种方法通过加法的运算来交换变量,但是要注意这种方法有局限性不能交换非数字类型的变量对于浮点类型的变量可能会出现问题加法过程中可能会出现超出上限的问题x = x + y;y = x - y;x = x - y;异或交换法通过异或运算的性质得到的一种交换方法,这
2020-10-07 22:04:23
1777
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人