weixin_54010759的博客

原创 CoFSM基于共现尺度空间的多模态遥感图像匹配方法--论文阅读记录

除了论文中提出的经典的尺度空间构造等细节，下图关于base_image的创建也是值得注意的，根据3.3 Frequency of sampling in the spatial domain第二段，有时候可以对输入图片进行长宽扩大一倍达到增大等效第一层的尺度sigma的效果（因为sigma越大卷积耗时越久）Distinctive Image Features from Scale-Invariant Keypoints，作者：David G. Lowe。另外，特征匹配阶段的最近邻次近邻比也值得注意，

原创 sift 解释

转载。

原创 matlab conv2

MATLAB卷积conv、conv2、convn详解-CSDN博客

原创 量化多因子

【代码】量化多因子。

原创 python依次运行多个文件

【代码】python依次运行多个文件。

原创 clickhouse sql语句

对于v1表，选出两种类型，只有CS，只发生CS和NULL转变的两种cik-ticker，得到v2表。select语句，group by，

原创 pytorch GPU版本安装 python windows

查找如下pytorch版本，得到安装命令（测试使用conda命令无法成功安装，使用pip3命令进行安装，命令如下）查找CUDA版本是12.0，查找方式，win + r输入cmd进入命令行模式，输入nvidia-smi，如下，如果在虚拟环境cmd命令界面，输入python进行python编辑模式。在安装上述版本前，已经创建了pytorch19,python版本是3.9，安装命令如下。此外查询torchvision，如果没有安装，安装即可，torch。安装其他库的命令，如下。发现成功安装，如下，

原创 python pandas中时间类型

pre_date = pd.Timestamp(str(t_date)) - pd.Timedelta(days=1)：这行代码计算了t_date的前一天日期。pd.Timestamp(str(t_date))将t_date转换为Timestamp对象，然后使用pd.Timedelta(days=1)创建一个时间增量，表示1天的时间间隔。t_date是一个日期字符串，通过str()函数将其转换为字符串形式，然后使用pd.to_datetime()将其转换为Timestamp对象。#t_date的前一天。

原创 远程VPN登录，IPsec，VPN，win10

windows10 完美解决L2TP无法连接问题。

原创 int类型最大的数测试

原创 数组模拟循环链表

【代码】数组模拟循环链表。

原创 字符串匹配--模拟

【代码】字符串匹配--模拟。

原创 无向图-已知根节点求高度

深搜板子题，无向图，加边加两个，dfs输入两个参数变量，一个是当前深搜节点，另一个是父节点（避免重复搜索父节点），恢复现场。

原创 拓扑排序--代码

在安排每个人的顺序时，有 M个要求，每个要求包含两个整数 a,b，表示小朋友 a 要排在小朋友 bb 的前面。拓扑排序 —— 模板题 AcWing 848. 有向图的拓扑序列。当符合条件的排队顺序不唯一时，编号更小的小朋友尽量更靠前。时间复杂度 O(n+m), n 表示点数，m 表示边数。按排好队列从前到后的顺序在一行内输出每个小朋友的编号。N 个小朋友，编号 1∼N，要排成一队。接下来 M行，每行包含两个整数 a,b。请你找出符合所有要求的排队顺序。保证数对 (a,b)各不相同。第一行包含整数 N,M。

原创 数位dp

数位DP学习整理（数位DP看完这篇你就会了）_unique_pursuit的博客-CSDN博客

原创 生成最小整数

对于数字n，删除k个数，考虑最后一次删除的数的位置，如果已知数字n的位数为m，那么最后一次删除的数的位置可能在第i位，其中，i=1，2，3，4...，m，对于某个情况，设得到的数是nn，问题就转化为了对于数字nn，删除k-1个数，最后的新整数最小值，再回到函数f(n,k)，求min(nni,k-1)，其中nni=删除位置分别是i=1，2，3，。临界条件是k=0，k=0返回对应的数字。现在从m位中删除k位(0

原创 记录字符串压缩

如果限定字符有空格，大小写字母和数字，直接输出压缩结果即可，注意，应该需要将。，如果直接输出出现的次数，输出耗时较长。注：可能上述做法也不是正确解决方案。

原创 动态规划--输出路径06.25

由于需要求解最小字典序，尝试输入数据时逆序输入，其他不改变，状态含义不变，如下，10/11。f[i][j]表示前i个物品，重为j时的最大价值。如何输出路径，并且保证是最小字典序，

原创 双端队列，6.24

求滑动窗口内的最大值，单调队列做法。

原创 数组模拟邻接表

【代码】数组模拟邻接表。

原创 数字组合 自然数的拆分

2.f[i][j]=f[i-1][j]+f[j-wi]，更新第i层的j需要用到上一层的j和j-wi，从尾部遍历，即可，具体，使用f[j]=f[j]+f[j-wi]，j从m开始递减，到wi结束，因为j-wi要大于等于0，取前i个数中的若干个数，和为0的方案有1中，就是一个数都不选这种状态（对应二进制压缩每个二进制位都是0）。1.第i层只和上一层有关系，使用滚动数组求解，少一点判断，j从i开始。

原创 动态规划---三题 石子归并，多边形 金字塔 0607

【代码】动态规划---三题 石子归并，公式求和 金字塔 0607。

原创 子串分值--子串分值和 模拟，找规律

res+=(next[ii]-i)*(i-pre[ii])// next[ii]=I,，此处加的都是0，错误的。因此，需要记录i处字符a[i]的next和pre数组，初始化pre为-1，next为字符串长度。对于得到的s[i]，ii=s[i],next[ii]=i，下一个位置就是当前位置。4c贡献 1=1 c----1*5 next c 没有 5；Next[i]记录s[i]的下一处，pre[i]记录s[i]的下一处，i-p[s[i]] 贡献次数 n[s[i]]-i。

原创 子串--子字符串 0528

设置i，j分别指向头和尾，对于每一个i，遍历j，保证ia[j]，可以，如果a[i]==a[j]，查看下一个，即如果a[i+1]>a[j-1]，可以，若相等，查看下一个，O(n^3)，这里实际在找的是什么，写出代码，看看能不能优化。思考，对于f(i,j)存储区间[i,j]中，最小的ii，使得a[ii]>a[jj]，ii+jj=i+j。a[i]

原创 平方差 找规律-推导-模拟暴力 子串模拟暴力

观察，没有出现只能被2整除却不能被4整除的数，考虑奇偶性。暴力搜索，如下，输出1-100的可以平方差分解的数。

原创 同余方程 乘法逆元 扩展欧几里得算法 5.26

用扩展算法，求ax+my=gcd(a,m)的解把x乘以b/gcd(a,m)得到原方程的特解。ax=b(mod m)得到ax=m(-y)+b，得到ax+my=b;给定整数a,b,m，求解同余方程ax=b(mod m).例： 8x=4(mod 6).整数解x=2；由裴属定理，当gcd(a,m)|b时有解。同余方程 乘法逆元 扩展欧几里得算法。如果x存在整数解，则输出任意一个；扩展欧几里得求解线性同余方程，把同余方程转化为求解不定方程，如果不存在，则输出none。

原创 字符串---第一部分 序列、字串；上升，公共

b[j]，分为三种情况，分别是a[i]不在最长公共子序列中，f(i,j)=f(i-1,j),b[j]不在最长公共子序列中，f(i,j)=f(i,j-1),a[i]和b[j]都不在最长公共子序列中，f(i,j)=f(i-1,j-1),其中第三种可以看成前两种的特例，如果a[i]!遍历s[0]...s[i-1]，如果记录最大的j，使得s[j]

原创 5.25 费解的开关

考虑，第一行开关状态固定，能够影响到第一行的只有第二行，其他行无法影响第一行，所以，对于第一行中的第k列，如果第k列是不亮的，为了使该灯变亮，需要第二行的第k列按钮按下；第三行的开关操作由第二行的灯的状态确定，第四行的开关操作由第三行的灯的状态确定，第五行的操作由第四行的灯的状态确定，第五行是最后一行，到达枚举的终点，如果第五行的灯都是亮的，说明当前枚举的开关操作是合理的，记录开关操作次数，否则说明当前枚举的开关操作不能使灯全部都亮，操作不满足要求， 更新需要的最少的开关的操作次数。

原创 【无标题】

有一根27厘米的细木杆，在第3厘米、7厘米、11厘米、17厘米、23厘米这五个位置上各有一只蚂蚁。木杆很细，只能同时通过一只蚂蚁。开始时，蚂蚁的头朝左还是朝右是任意的，它们只会朝前走或调头，但不会后退。当任意两只蚂蚁碰头时，两只蚂蚁会同时调头朝反方向走。假设蚂蚁们每秒钟可以走1厘米的距离。编写程序，求所有蚂蚁都离开木杆的最小时间和最大时间。思路：当两只蚂蚁相撞时，可以看成两只蚂蚁都没有掉头而直接向前爬行，最小时间即沿蚂蚁朝向到达杆外面的最短时间，最长时间即沿蚂蚁朝向到达杆外面的最长时间。

原创 opencv视频操作

视频读取：1.打开视频文件 case1：通过设备摄像头进行视频获取 case2：打开视频文件 2.获取视频每秒的帧数fps，即frame per second，case1：本地视频，使用opencv中的常数 case2：未知视频，计算，读取时间T内的视频帧数N，如果单线程解决视频处理问题并写视频，得到的处理后的视频不会出什么问题，如果多个线程读取视频帧，视频帧的写入顺序需要注意，怎么才能写好，不发生错位。3.视频逐帧读取，视频可视化显示，视频句柄释放。基于opencv，使用C++/C进行相关操作。

原创 0524 基础题

【代码】0524 基础题。

原创 5.24 基础题目

state_k是state除去j之后的集合，state_k要包含k，k是state_k二进制表示中为1的下标，枚举出k。j++){//n点数。最短Hamilton路径。

原创 二分查找5.18

具有对称性，退出条件L+1==r，并排退出，L和r初始为开区间端点，默认数据已经从小到大存储在a[1]-a[n]中，退出循环返回的是可行区域内的指针值。

原创 5.18 离散线段树-扫描线

5.18 离散线段树-扫描线。

原创 5.17 普通平衡树 splay算法

5.17 普通平衡树。

原创 线段树--数组-可持久化-trie-5.16

可持久化线段树，存储历史版本。

原创 每日一题 最大异或对--数组模拟邻接表

【代码】每日一题。

原创 每日一练2023.4.17-2023.4.18

手写队列，队列有两个指针-front和rear，对于这两个指针，操作都是+1就是上移，如果需要构造循环队列，则需要进行取余对队列长度queueLength进行操作，队列为空等价于front不等于rear，即队首队尾指针不同，稍微整理代码如下，也可能不是，不是的话，利用该队头元素扩展与该节点有边相连且没有被访问的节点，将新的节点入队 ，直到队为空，所有图中相连的节点都被访问过。它表示一个迷宫，其中的1表示墙壁，0表示可以走的路，只能横着走或竖着走，不能斜着走，要求编程序找出从左上角到右下角的最短路线。

原创 链接收藏收录

比如，想要在test1.cpp和test2.cpp文件中使用变量int flag，任选这两个文件的一个，比如选择test1.cpp，在其中定义extern int flag;然后在test2.cpp文件中定义int flag即可，注意，不可以在这两个文件中同时写extern int flag;Linux系统下挂载硬盘不再困扰你！目标是根据画的多边形生成掩膜（即白色区域）图像，生成的图像大小与原图大小保持一致，在对应位置生成掩膜。1、把画的多边形的各顶点坐标读出，要的是在图像中的像素坐标。

原创 opencv部分函数

目标是根据画的多边形生成掩膜（即白色区域）图像，生成的图像大小与原图大小保持一致，在对应位置生成掩膜。3、根据cv2库的fillPoly函数在全0或全1 图的对应多边形位置进行填充。左图，传入大正方形和小正方形的顶点即可；右图，传入小中大三个正方形的顶点即可。1、把画的多边形的各顶点坐标读出，要的是在图像中的像素坐标。2、生成与原图大小一致的全为0或1的图用于后续使用。4、输出填充后得到掩膜的图像cv2.imwrite。进一步，如果需要生成环状的，如下图。opencv画图函数源代码。