轨道力学中二体问题下求解兰伯特方程. 老外写的Matlab程序,我把它转成了Fortran程序. !***************************************************************** subroutine solve_lambert(r1,r2,tt,GM,lw,N,nBranch,v1,v2) implicit real(8)(A-H,O-Z) dimension r1(3),r2(3),v1(3),v2(3),tmp3(3),wih(3),r1…

https://zh.wikipedia.org/zh-hans/李亞普諾夫函數 李雅普诺夫函数(Lyapunov function)是用来证明一动力系统或自治微分方程稳定性的函数.其名称来自俄罗斯数学家亚历山大·李亚普诺夫(Александр Михайлович Ляпунов).李亚普诺夫函数在稳定性理论及控制理论中相当重要. 若一函数可能可以证明系统在某平衡点的稳定性,此函数称为李亚普诺夫候选函数(Lyapunov-candidate-function).不过目前还找不到一般性的方式可建…

coreseek实战(二):windows下mysql数据源部分配置说明 关于coreseek在windows使用mysql数据源的配置,以及中文分词的详细说明,请参考官方文档: mysql数据源配置:http://www.coreseek.cn/products-install/mysql/ 中文分词:http://www.coreseek.cn/products-install/coreseek_mmseg/ 我这里的配置根据dedecms的dede_archives/dede_addona…

pro:给定N*M的矩阵,每次操作一个位置,它会增加2,周围4个位置会增加1.给定初始状态,求一种方案,使得最后的数都为0:(%3意义下. sol:(N*M)^3的复杂度的居然过了.          好像标程是M^3的,因为第一排确定了,后面的都确定了.所以我们只要设关于第一排的方程,那么跑下来,第N+1排的都为0,则合法. (此题由于3的特殊性,x关于3的逆元=x:所以不用求逆元 #include<bits/stdc++.h> #define rep(i,a,b) for(int i=a;…

最近做的项目,由于预算有限,公司决定不采购Windows服务器,而采购基于Linux的服务器. 一般的VPS服务器,如果使用Windows系统,那么Windows Server2012\2016安装好后,就占用1GB的内存了:只要运行起来,2G的内存就用掉了:装上SqlServer.跑个.net,4G内存是标配,用户数量稍微多些4G内存也捉襟见肘了:分配2个VCpu,在桌面状态下不运行任何程序,每个核心的占用率在15%~20%左右. 大家知道64位的CentOS7装好后内存占用是多少吗,用fre…

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5627 题意:Bestcoder的一道题,让你求&意义下的最大生成树. 解法: 贪心,我们从高位到低位贪心,如果含有这一位的边能够构成一棵树的话,我们就可以直接把其他不含有这一位的边全部去掉 然后重复这个行为 这个贪心显然正确啦,至于判断能否构成一颗树,就用并查集就好啦 //HDU 5627 #include <bits/stdc++.h> using namespace std; cons…

题意:给定p=1e9+7,A,B.  求一对X,Y,满足(X+Y)%P=A; 且(X*Y)%P=B: 思路:即,X^2-BX+CΞ0;  那么X=[B+-sqrt(B^2-4C)]/2: 全部部分都要在modP意义下,所以求一个x满足x^2%p=B^2-4C,这个用二次剩余求即可. 套了模板. #include <bits/stdc++.h> using namespace std; ; typedef long long ll; int k; ll a,p,w; struct T{ll x,…

题目链接 首先利用组合数学知识,枚举两人的总胜场数容易得到 这还不是卷积的形式,直接搞的话复杂度大概是O(n^2)的,肯定会TLE.但似乎和卷积有点像?想半天没想出来..多谢Q巨提醒,才知道可以用下面这个公式进行转化 最后,化得的公式为 另外注意,上式右边是一个卷积的形式,但是,所得和的第一项是不需要加上的(不过图中公式没有体现).结合实际意义大概就是,i==0&&j==0时,gcd(i,j)不存在约数d,虽然0可以被任意正整数整除 & 第一项不为0 #include<bit…

//写在前面 单就FFT算法来说的话,下面只给出个人认为比较重要的推导,详细的介绍可参考 FFT算法学习笔记 令v[n]是长度为2N的实序列,V[k]表示该实序列的2N点DFT.定义两个长度为N的实序列g[n]和h[n]为 g[n]=v[2n], h[n]=v[2n+1], 0<=n<N 则可进行如下推导 这里所用的FFT算法能够实现O(nlogn)复杂度的离散傅里叶变换和上面最后所得的关系密切相关. 下面进入正题——模意义下的FFT 还是需要先了解一下关于 复序列的DFT的对称性质及一些补充…

链接:https://www.nowcoder.com/acm/contest/80/B 题目描述 在文某路学车中学高一新生军训中,Jxc正站在太阳下站着军姿,对于这样的酷热的阳光,Jxc 表示非常不爽. Jxc将天空看做一个n*n的矩阵,此时天上有m朵云,这些云会随机分布在m个不同的位置,同时太阳会随机出现在一个位置,Jxc想知道他被太阳晒到的概率是多少,由于他仍在站军姿,所以这个有趣的问题就交给了你.考虑到精度问题,Jxc只需要知道这个概率在对998244353取模意义下的值. Tips:一…