hello，good code.

原创 C++实现bam文件读写参考资料

为啥要自己实现这个？因为samtools这种工具读写bam文件太慢了。

原创 完美解决Error: error:0308010C:digital envelope routines::unsupported【windows&linux】

注意，windows在cmd和powershell设置值的命令是不一样的。在运行node.js时出现了如下报错。node.js版本不兼容。下载对应的版本，并安装。

原创 c++ stl 优先队列怎么输出全部数据，同时又获取top数据？

方法：使用优先队列的构造函数构造一个和当前优先队列一样的数据，然后输出。有时，我们在调试时可能会出现需要输出优先队列的全部数据来观看数据情况。例如，我们需要在输入top数据之前输出当前优先队列的数据。总所周知，优先队列是一种常用的数据结构。

原创 Samtools - Duplicate Marking

对于主要读取，此定义与 Picard (v2.10.3) 和 Biobambam2 (bamstreamingmarkduplicates v2.0.57) 中使用的相同。这些工具在确定重复时都不使用补充数据。然而，与 Picard 不同的是，复制模板中的补充读取默认情况下不会在 Samtools 中修改其标志。重复被定义为具有多个对齐的 5' 坐标匹配的模板。对于双端模板，这需要两个主要读取具有匹配的 5' 坐标。当检测到重复时，整体最高质量的模板将被保留，所有其他模板都设置重复标志。

原创 使用podman创建debian/centos/ubuntu系统容器

是一个无守护程序、开源的Linux原生工具，旨在使用Open Containers Initiative轻松查找、运行、构建、共享和部署应用程序 (Podman专注于帮助您维护和修改OCI容器图像的所有命令和功能，例如拉取和标记。它允许您在生产环境中创建、运行和维护这些容器和容器映像。这使得Podman创建的运行容器与任何其他常见容器引擎创建的容器几乎没有区别。Podman控制下的容器既可以由超级用户运行，也可以由非特权用户运行。Podman管理整个容器生态系统，其中包括Pod、容器、容器映像和容器卷。

原创 常见的容器仓库

最近在使用podman，相对于docker优点很明显，不需要root。下面介绍几个常见的容器仓库。

原创 使用picard工具去重自己的sam/bam数据

自 2.0.1 版（2015 年 11 月）起，Picard 需要 Java 1.8 (jdk8u66)。支持 Java 1.7 的最后一个版本是 1.141 版。picard是一组用于操作高通量测序 (HTS) 数据和格式的 Java 命令行工具。Picard 现在在 Java 11 下构建并通过测试。实现的，以支持访问通常用于高通量测序数据的文件格式，例如。Picard 是使用 HTSJDK Java 库。去重自己的sam/bam数据。

原创 207. 课程表 (拓扑排序）

题目链接：https://leetcode-cn.com/problems/course-schedule/解题代码：class Solution {public: bool canFinish(int numCourses, vector<vector<int>>& prerequisites) { vector<vector<int> > edges; edges.resize(numCours.

原创 什么是λ矩阵？什么是多项式矩阵？

为多项式矩阵或λ矩阵。

原创 集合F是啥？

F是一个集合名，如果是在矩阵中一般指代数域F那么数域是啥？设P是由一些复数组成的集合，其中包括0与1，如果P中任意两个数的和、差、积、商（除数不为0）仍是P中的数，则称P为一个数域。常见数域： 复数域C；实数域R；有理数域Q。说明：1）若数集P中任意两个数作某一运算的结果仍在P中，则说数集P对这个运算是封闭的。2）数域的等价定义：如果一个包含0，1在内的数集P对于加法，减法，乘法与除法（除数不为0）是封闭的，则称数集P为一个数域。...

原创 矩阵多项式

原创 矩阵范数

原创 什么是实践意义?和理论意义?

理论意义：比前人同类命题增加了什么内容,有什么新见,与同时代相比,有什么突破之处,在某个理论体系中起了什么作用,对后代理论发展又有什么意义,为一种理论体系的完善增补了什么,对后人有什么启示.所谓理论就是在人类认识论的角度对有关命题的知识范围和适用性进行的研究并得出的成果.比如,数学的算法“＋,－,X,...”是有理论意义的；数与数之间的关系的理论也是有理论意义的.世界是由物质的还是由思想决定的?这类问题的回答,都是由理论意义的.或者说,理论是一种思想的发展,属于行而上的范畴.实践意义就在于可以利用其中

原创 Python的深浅拷贝讲解！

前言在很多语言中都存在深浅拷贝两种拷贝数据的方式，Python中也不例外。本文中详细介绍了Python中的深浅拷贝的相关知识，文章的内容包含： 对象、数据类型、引用 赋值 浅拷贝 深拷贝 一、Python对象我们经常听到：在Python中一切皆对象。其实，说的就是我们在Python中构造的任何数据类型都是一个对象，不管是数字、字符串、字典等常见的数据结构，还是函数，甚至是我们导入的模块等，Python都会把它当做是一个对象来处理。所有的Python对象

原创 Mysql 忘记root密码的完美解决方法

通常在使用Mysql数据库时，如果长时间没有登陆，或者由于工作交接完成度不高，会导致数据库root登陆密码忘记，本文给大家介绍一种当忘记mysql root密码时的解决办法，一起看看吧一、更改my.cnf配置文件1.用命令编辑/etc/my.cnf配置文件，即：vim /etc/my.cnf 或者 vi /etc/my.cnf2.在[mysqld]下添加skip-grant-tables，然后保存并退出3.重启mysql服务：service mysqld restart二、更改r

原创 gcc和g++是什么，有什么区别？

目前（2020-09）GCC 编译器已经更新至 10.2版本，其功能也由最初仅能编译 C 语言，扩增至可以编译多种编程语言，其中就包括 C++ 。除此之外，当下的 GCC 编译器还支持编译 Go、Objective-C，Objective-C ++，Fortran，Ada，D 和 BRIG（HSAIL）等程序，甚至于 GCC 6 以及之前的版本还支持编译 Java 程序。但本教程主要讲解如何使用 GCC 编译器编译运行 C 和 C++ 程序，因此有关其它编程语言如何使用 GCC 编译器编译，将不再做具

原创 Fast Food hdu 1227 dp动态规划

题目大意快餐连锁mcburger拥有公路沿线几家餐馆。最近，他们决定在高速公路上建造几处仓库，每一个位于一个餐厅，并提供一些所需的食材的餐馆。当然，这些仓库应放置，使一个餐厅和它分配的仓库之间的平均距离最小化。你要写一个程序来计算仓库的最佳位置和作业。为了使其更精确，对mcburger管理具有以下规范发布：你将获得位置的N餐厅沿着公路n D1，D2，…，DN（这些都是从公司总部，这恰好

原创 三种编程命名规则（匈牙利法，小驼峰法，大驼峰法）

1. 匈牙利命名：开头字母用变量类型的缩写，其余部分用变量的英文或英文的缩写，要求单词第一个字母大写。For example: long lsum = 0;"l"是类型的缩写； 2. 小驼峰式：(little camel-case)第一个单词首字母小写，后面其他单词首字母大写。For example: string firstName = strin

原创 hdu 1016 素数环 深搜，回溯

#includeusing namespace std;#include#include#includeint n,a[23],vis[23],Case=0;bool isp(int n){for(int i=2;i<=sqrt(n);i++) if(n%i==0) return 0;return 1;}void dfs(int cur){if(cur==n&&isp

原创 全排列 深搜dfs

#include#includeusing namespace std;int a[100],vis[100],n;void dfs(int t){if(n<t){ for(int i=1;i<=n;i++) printf("%d ",a[i]); puts(""); return;}for(int i=1;i<=n;i++)if(vis[i

原创 codevs1048 石子归并 动态规划dp

#include#includeusing namespace std;int n,a[1000]={0},dp[1000][1000]={0},s[1000]={0};int main(){cin>>n;memset(dp,11,sizeof(dp));for(int i=1;i<=n;i++){ cin>>a[i]; s[i]=a[i]+s[i-1]; d

原创 codevs1063 合并果子 贪心

#include#includeusing namespace std;long long a[10005];int main(){int n,i,t;long long ans,p;scanf("%d",&n);for(i=1;i<=n;i++)scanf("%d",&a[i]);sort(a+1,a+n+1);ans=0;for(i=1;i<

原创 codevs1576 最长严格上升子序列 动态规划dp

#include#includeusing namespace std;int main(){ int n,m=-1; cin>>n; long a[1000],dp[1000]={0}; for(int i=1;i<=n;i++) cin>>a[i]; for(int i=1;i<=n;i++) dp[i]=1;

原创 codevs1154 能量项链 入门dp

#include #include #include #include using namespace std;int n,maxx,sz[220],dp[220][220];int main(){scanf("%d",&n);for (int i = 1;i <= n;i++){scanf("%d",&sz[i]);sz[i+n] = sz[i];}for(int j

原创 poj1503 大整数加法

题目：    大整数加法，每个整数至多包含100个数字，最多100个整数相加，计算并输出结果。输入：     最多输入100行数据，每行输入一个正整数，输入0结束输入。输出：    输出所有正整数累加的结果。注意输入01111111的情况。代码如下#includeusing namespace std;#include#includeint main(){

原创 历届试题 九宫重排 （搜索+哈希）

问题描述　　如下面第一个图的九宫格中，放着 1~8 的数字卡片，还有一个格子空着。与空格子相邻的格子中的卡片可以移动到空格中。经过若干次移动，可以形成第二个图所示的局面。　　我们把第一个图的局面记为：12345678.　　把第二个图的局面记为：123.46758　　显然是按从上到下，从左到右的顺序记录数字，空格记为句点。　　本题目的任务是已知九宫的初态和终态，求最少经过

原创 codevs 1215 迷宫 广搜bfs 深搜dfs

题目链接：http://codevs.cn/problem/1215/一个纯广搜的题目.广搜和深搜都可以广搜代码如下#include#include#include#includeusing namespace std;int main(){ int num,n,x0,y0,x1,y1,x2,y2,x3,y3,flag,dir[4][2]={{1,0},{

原创 poj 1230 贪心法

题意 :  有一块棋盘式的场地和一个魔术师，场地中有n道墙，魔术师一次最多能穿越k道墙。求移除最少数量的墙，使得魔术师可以在场地的任意一列表演都可以成功。用二维数组来表示场地思路 : 如果穿不过，贪心找出从该列起向右延伸最长的m道墙，移除这m道墙使得该列能穿过。#include#include#includeusing namespace std;

原创 poj1258 最小生成树，prim

#include#include#includeusing namespace std;int main(){ // freopen("1.txt","r",stdin);int maps[110][110],vis[110]={0},dis[110],n,minn,x,y,ans,k;while(cin>>n&&n){for(int i=1;i<=n;i++) for

原创 poj1488 字符串处理

#include#includeusing namespace std;int main(){string a;unsigned int pos,cnt=1;while(getline(cin,a)){    while((pos=a.find('\"'))!=string::npos){            if(cnt%2==1)        a.rep

原创 poj2408 字符串处理，排序

#include#include#include#include#include#includeusing namespace std;struct node{string a,b;}a[30005],g[30005];struct nod{string a,b;int n;}c[30005];bool Cmp(node a,node b){return a.a<b.

原创 hdu1874 dijkstra

dijkstra算法，呵呵代码#include#include#includeusing namespace std;int main(){ // freopen("1.txt","r",stdin);int t1,t2,s,t,t3,n,m,dis[205],a[205][205];bool vis[205];while(cin>>n>>m){

原创 hdu3549 网络流

#include#include#include#includeusing namespace std;int pre[1005],s=1,t,g[1005][1005],vis[1005],m;bool bfs(){ memset(pre,0,sizeof(pre)); memset(vis,0,sizeof(vis));int k;queue q;q.push

原创 poj 1273 网络流

#include#include#include#includeusing namespace std;int res[202][202],minn,n,m,pre[202],a[202];bool bfs(){queueq;q.push(1);int k;memset(a,0,sizeof(a));a[1]=0x3f3f3f3f;while(!q.empty()){

原创 spfa

粗略讲讲SPFA算法的原理，SPFA算法是1994年西安交通大学段凡丁提出是一种求单源最短路的算法算法中需要用到的主要变量int n; //表示n个点，从1到n标号int s,t; //s为源点，t为终点int d[N]; //d[i]表示源点s到点i的最短路int p[N]; //记录路径（或者说记录前驱）queue q; //

原创 最长公共递增子序列

#include#include#include#includeusingnamespacestd;#define BUF_SIZE 100000#define OUT_SIZE 100000#define ll long long//fread->readboolIOerror=0;inlinecharnc(){staticcharbuf[BUF_SIZE],*p1=buf

原创 codevs 1028 花店橱窗布置

网络流二分图匹配题目描述 Description假设以最美观的方式布置花店的橱窗，有F束花，V个花瓶，我们用美学值（一个整数）表示每束花放入每个花瓶所产生的美学效果。为了取得最佳的美学效果，必须使花的摆放取得最大的美学值。输入描述 Input Description第一行为两个整数F,V(F接下来F行每行V个整数，第i行第j个数表示第i束花放入第j个花瓶的美学值。输出描述 Output Desc

原创 16进制转8进制

#include#include#includeusing namespace std;const int Max=100000+2;int main(){char a[Max];int b[Max*4],n;cin>>n;while(n--){cin>>a;int k=1;int len=strlen(a);for(int i=0;i<len;i++){int t;

原创 基础练习 字母图形

#include#includeusing namespace std;int main(){int n,m;cin>>n>>m;char t;for(int i=1;i<=n;i++){for(int j=i-1,f=1;j>0&&f<=m;j--,f++){t='A'+j;cout<<t;}for(int j=0;j<=m-i;j++){t='A'+j;cout<<

原创 java之大数的处理

Java的提供了操作大数的类,即java.math.BigInter与java.math.BigDecimal类。Big.math.BigInter是针对大整数的处理类，而BigDecimal针对大 小数的处理BigInteger类常用的几种运算方法：public BigInteger （String val) :构造BigIntegerpublic Big